def make_model(images_dir=os.path.join(base_dir, 'images')):
print 'initializing the model'
# set up the modeling environment
start_time = datetime(2016, 9, 23, 0, 0)
model = Model(start_time=start_time,
duration=timedelta(days=2),
time_step=30 * 60,
uncertain=False)
print 'adding the map'
model.map = GnomeMap() # this is a "water world -- no land anywhere"
# renderere is only top-down view on 2d -- but it's something
renderer = Renderer(output_dir=images_dir,
image_size=(1024, 768),
output_timestep=timedelta(hours=1),
)
renderer.viewport = ((196.14, 71.89), (196.18, 71.93))
print 'adding outputters'
model.outputters += renderer
# Also going to write the results out to a netcdf file
netcdf_file = os.path.join(base_dir, 'script_arctic_plume.nc')
scripting.remove_netcdf(netcdf_file)
model.outputters += NetCDFOutput(netcdf_file,
which_data='most',
# output most of the data associated with the elements
output_timestep=timedelta(hours=2))
print "adding Horizontal and Vertical diffusion"
# Horizontal Diffusion
model.movers += RandomMover(diffusion_coef=500)
# vertical diffusion (different above and below the mixed layer)
model.movers += RandomMover3D(vertical_diffusion_coef_above_ml=5,
vertical_diffusion_coef_below_ml=.11,
mixed_layer_depth=10)
print 'adding Rise Velocity'
# droplets rise as a function of their density and radius
model.movers += TamocRiseVelocityMover()
print 'adding a circular current and eastward current'
fn = 'hycom_glb_regp17_2016092300_subset.nc'
fn_ice = 'hycom-cice_ARCu0.08_046_2016092300_subset.nc'
iconc = IceConcentration.from_netCDF(filename=fn_ice)
ivel = IceVelocity.from_netCDF(filename=fn_ice, grid = iconc.grid)
ic = IceAwareCurrent.from_netCDF(ice_concentration = iconc, ice_velocity= ivel, filename=fn)
model.movers += PyCurrentMover(current = ic)
model.movers += SimpleMover(velocity=(0., 0., 0.))
model.movers += constant_wind_mover(20, 315, units='knots')
# Now to add in the TAMOC "spill"
print "Adding TAMOC spill"
model.spills += tamoc_spill.TamocSpill(release_time=start_time,
start_position=(196.16, 71.91, 40.0),
num_elements=1000,
end_release_time=start_time + timedelta(days=1),
name='TAMOC plume',
TAMOC_interval=None, # how often to re-run TAMOC
)
model.spills[0].data_sources['currents'] = ic
return model
def make_model(images_dir=os.path.join(base_dir, 'images')):
print 'initializing the model'
# set up the modeling environment
start_time = datetime(2016, 9, 18, 1, 0)
model = Model(start_time=start_time,
duration=timedelta(days=3),
time_step=30 * 60,
uncertain=False)
print 'adding the map'
model.map = GnomeMap() # this is a "water world -- no land anywhere"
# renderere is only top-down view on 2d -- but it's something
renderer = Renderer(
output_dir=images_dir,
size=(1024, 768),
output_timestep=timedelta(hours=1),
)
renderer.viewport = ((-87.095, 27.595), (-87.905, 28.405))
print 'adding outputters'
model.outputters += renderer
# Also going to write the results out to a netcdf file
netcdf_file = os.path.join(base_dir, 'gulf_tamoc.nc')
scripting.remove_netcdf(netcdf_file)
model.outputters += NetCDFOutput(
netcdf_file,
which_data='most',
# output most of the data associated with the elements
output_timestep=timedelta(hours=2))
print "adding Horizontal and Vertical diffusion"
# Horizontal Diffusion
model.movers += RandomMover(diffusion_coef=100000)
# vertical diffusion (different above and below the mixed layer)
model.movers += RandomVerticalMover(
vertical_diffusion_coef_above_ml=50,
vertical_diffusion_coef_below_ml=10,
horizontal_diffusion_coef_above_ml=100000,
horizontal_diffusion_coef_below_ml=100,
mixed_layer_depth=10)
print 'adding Rise Velocity'
# droplets rise as a function of their density and radius
model.movers += TamocRiseVelocityMover()
print 'adding the 3D current mover'
gc = GridCurrent.from_netCDF('HYCOM_3d.nc')
model.movers += GridCurrentMover('HYCOM_3d.nc')
# model.movers += SimpleMover(velocity=(0., 0, 0.))
# model.movers += constant_wind_mover(5, 315, units='knots')
# Wind from a buoy
w = Wind(filename='KIKT.osm')
model.movers += WindMover(w)
# Now to add in the TAMOC "spill"
print "Adding TAMOC spill"
model.spills += tamoc_spill.TamocSpill(
release_time=start_time,
start_position=(-87.5, 28.0, 2000),
num_elements=30000,
end_release_time=start_time + timedelta(days=2),
name='TAMOC plume',
TAMOC_interval=None, # how often to re-run TAMOC
)
model.spills[0].data_sources['currents'] = gc
return model
def make_model(images_dir=os.path.join(base_dir, 'images')):
print 'initializing the model'
# set up the modeling environment
start_time = datetime(2004, 12, 31, 13, 0)
model = Model(start_time=start_time,
duration=timedelta(days=3),
time_step=30 * 60,
uncertain=False)
print 'adding the map'
model.map = GnomeMap() # this is a "water world -- no land anywhere"
# renderere is only top-down view on 2d -- but it's something
renderer = Renderer(
output_dir=images_dir,
size=(1024, 768),
output_timestep=timedelta(hours=1),
)
renderer.viewport = ((-.15, -.35), (.15, .35))
print 'adding outputters'
model.outputters += renderer
# Also going to write the results out to a netcdf file
netcdf_file = os.path.join(base_dir, 'script_plume.nc')
scripting.remove_netcdf(netcdf_file)
model.outputters += NetCDFOutput(
netcdf_file,
which_data='most',
# output most of the data associated with the elements
output_timestep=timedelta(hours=2))
print "adding Horizontal and Vertical diffusion"
# Horizontal Diffusion
# model.movers += RandomMover(diffusion_coef=5)
# vertical diffusion (different above and below the mixed layer)
model.movers += RandomVerticalMover(vertical_diffusion_coef_above_ml=5,
vertical_diffusion_coef_below_ml=.11,
mixed_layer_depth=10)
print 'adding Rise Velocity'
# droplets rise as a function of their density and radius
model.movers += RiseVelocityMover()
print 'adding a circular current and eastward current'
# This is .3 m/s south
model.movers += PyCurrentMover(current=vg,
default_num_method='Trapezoid',
extrapolate=True)
model.movers += SimpleMover(velocity=(0., -0.1, 0.))
# Now to add in the TAMOC "spill"
print "Adding TAMOC spill"
model.spills += tamoc_spill.TamocSpill(
release_time=start_time,
start_position=(0, 0, 1000),
num_elements=1000,
end_release_time=start_time + timedelta(days=1),
name='TAMOC plume',
TAMOC_interval=None, # how often to re-run TAMOC
)
return model
def make_model():
"""
Set up a GNOME simulation that uses TAMOC
Set up a spill scenario in GNOME that uses TAMOC to simulate a subsurface
blowout and then pass the TAMOC solution to GNOME for the far-field
particle tracking
"""
# Set up the directory structure for the model
base_dir, images_dir, outfiles_dir = set_directory_structure()
# Set up the modeling environment
print '\n-- Initializing the Model --'
start_time = "2019-06-01T12:00"
model = gs.Model(start_time=start_time,
duration=gs.days(3),
time_step=gs.minutes(30),
uncertain=False)
# Add a map
print '\n-- Adding a Map --'
model.map = gs.GnomeMap()
# Add image output
print '\n-- Adding Image Outputters --'
renderer = gs.Renderer(output_dir=images_dir,
image_size=(1024, 768),
output_timestep=gs.hours(1),
viewport=((-0.15, -0.35), (0.15, 0.35)))
model.outputters += renderer
# Add NetCDF output
print '\n-- Adding NetCDF Outputter --'
if not os.path.exists(outfiles_dir):
os.mkdir(outfiles_dir)
netcdf_file = os.path.join(outfiles_dir, 'script_tamoc.nc')
gs.remove_netcdf(netcdf_file)
file_writer = gs.NetCDFOutput(netcdf_file,
which_data='all',
output_timestep=gs.hours(2))
model.outputters += file_writer
# Add a spill object
print '\n-- Adding a Point Spill --'
end_release_time = model.start_time + gs.hours(12)
point_source = ts.TamocSpill(num_elements=100,
start_position=(0.0, 0.0, 1000.),
release_duration=gs.hours(24),
release_time=start_time,
substance='AD01554',
release_rate=20000.,
units='bbl/day',
gor=500.,
d0=0.5,
phi_0=-np.pi / 2.,
theta_0=0.,
windage_range=(0.01, 0.04),
windage_persist=900,
name='Oil Well Blowout')
model.spills += point_source
# Create an ocean environment
water, wind, waves = base_environment(water_temp=273.15 + 21.,
wind_speed=5.,
wind_dir=225.)
# Add a uniform current in the easterly direction
print '\n-- Adding Currents --'
uniform_current = gs.SimpleMover(velocity=(0.1, 0.0, 0.))
model.movers += uniform_current
# Add a wind mover
wind_mover = gs.WindMover(wind)
model.movers += wind_mover
# Add particle diffusion...note, units are in cm^2/s
print '\n-- Adding Particle Diffusion --'
particle_diffusion = gs.RandomMover3D(
horizontal_diffusion_coef_above_ml=100000.,
horizontal_diffusion_coef_below_ml=10000.,
vertical_diffusion_coef_above_ml=100.,
vertical_diffusion_coef_below_ml=10.,
mixed_layer_depth=15.)
model.movers += particle_diffusion
# Add rise velocity for droplets
print '\n-- Adding Particle Rise Velocity --'
# fixme: we do have code for rise velocity:
# gnome.movers.RiseVelocityMover
# let's test that later
slip_velocity = gs.SimpleMover(velocity=(0.0, 0.0, -0.1))
model.movers += slip_velocity
# Add dissolution
print '\n-- Adding Weathering --'
evaporation = Evaporation(water=water, wind=wind)
model.weatherers += evaporation
dissolution = Dissolution(waves=waves, wind=wind)
model.weatherers += dissolution
return model