def simple_model(output_dir):
start_time = "2018-09-20T12:00"
model = gs.Model(
start_time=start_time,
duration=gs.days(1),
time_step=gs.minutes(30),
name="test model for tideflats",
)
model.map = MapFromBNA(bna_file)
model.movers += gs.constant_wind_mover(10, 300, "m/s")
model.spills += gs.surface_point_line_spill(
num_elements=100,
start_position=(5.4, 53.38, 0),
end_position=(5.8, 53.4, 0),
release_time=start_time,
)
model.outputters += gs.Renderer(
output_timestep=gs.hours(1),
map_filename=bna_file,
output_dir=output_dir,
formats=['gif'], # ['gif', 'png']
image_size=(800, 400),
# viewport=((4.5, 53.0),
# (5.5, 53.5)),
)
return model
def simple_model(output_dir):
start_time = "2018-09-20T12:00"
model = gs.Model(start_time=start_time,
duration=gs.days(1),
time_step=gs.minutes(30),
name="test model for tideflats",
)
model.map = MapFromBNA(bna_file)
model.movers += gs.constant_wind_mover(10, 300, "m/s")
model.spills += gs.surface_point_line_spill(num_elements=100,
start_position=(5.4, 53.38, 0),
end_position=(5.8, 53.4, 0),
release_time=start_time,
)
model.outputters += gs.Renderer(output_timestep=gs.hours(1),
map_filename=bna_file,
output_dir=output_dir,
formats=['gif'], # ['gif', 'png']
image_size=(800, 400),
# viewport=((4.5, 53.0),
# (5.5, 53.5)),
)
return model
def make_model():
print 'initializing the model'
model = gs.Model(duration=gs.days(5))
print 'adding outputters'
budget_file = os.path.join(base_dir, 'GNOME_oil_budget.csv')
model.outputters += gs.OilBudgetOutput(budget_file)
print 'adding a spill'
# We need a spill at the very least
spill = gs.point_line_release_spill(num_elements=10, # no need for a lot of elements for a instantaneous release
start_position=(0.0, 0.0, 0.0),
release_time=model.start_time,
amount=1000,
substance='ALASKA NORTH SLOPE (MIDDLE PIPELINE, 1997)',
units='bbl')
model.spills += spill
print 'adding a RandomMover:'
model.movers += gs.RandomMover()
print 'adding a wind mover:'
model.movers += gs.constant_wind_mover(speed=10, direction=0, units="m/s")
model.environment += gs.Water(25, units={"temperature": "C"})
waves = gs.Waves()
model.environment += waves
print 'adding the standard weatherers'
model.add_weathering()
return model
def make_model(images_dir=os.path.join(base_dir, 'images')):
# create the maps:
print 'creating the maps'
mapfile = gs.get_datafile(os.path.join(base_dir, './MassBayMap.bna'))
gnome_map = gs.MapFromBNA(
mapfile,
refloat_halflife=1, # hours
raster_size=2048 * 2048 # about 4 MB
)
renderer = gs.Renderer(mapfile,
images_dir,
image_size=(800, 800),
projection_class=GeoProjection)
print 'initializing the model'
# start_time = datetime(2013, 3, 12, 10, 0)
start_time = "2013-03-12T10:00"
# 15 minutes in seconds
# Default to now, rounded to the nearest hour
model = gs.Model(time_step=gs.minutes(15),
start_time=start_time,
duration=gs.days(1),
map=gnome_map,
uncertain=True)
print 'adding outputters'
model.outputters += renderer
netcdf_file = os.path.join(base_dir, 'script_boston.nc')
gs.remove_netcdf(netcdf_file)
model.outputters += gs.NetCDFOutput(netcdf_file, which_data='all')
model.outputters += gs.KMZOutput(
os.path.join(base_dir, 'script_boston.kmz'))
print 'adding a RandomMover:'
model.movers += gs.RandomMover(diffusion_coef=100000)
print 'adding a wind mover:'
# series = np.zeros((2, ), dtype=datetime_value_2d)
# series[0] = (start_time, (5, 180))
# series[1] = (start_time + timedelta(hours=25), (5, 180))
# w_mover = WindMover(Wind(timeseries=series, units='m/s'))
# model.movers += w_mover
# model.environment += w_mover.wind
w_mover = gs.constant_wind_mover(5, 180, units='m/s')
model.movers += w_mover
print 'adding a cats shio mover:'
curr_file = gs.get_datafile(os.path.join(base_dir, r"./EbbTides.cur"))
tide_file = gs.get_datafile(os.path.join(base_dir, r"./EbbTidesShio.txt"))
c_mover = gs.CatsMover(curr_file, tide=gs.Tide(tide_file))
# this is the value in the file (default)
c_mover.scale_refpoint = (-70.8875, 42.321333)
c_mover.scale = True
c_mover.scale_value = -1
model.movers += c_mover
# TODO: cannot add this till environment base class is created
# model.environment += c_mover.tide
print 'adding a cats ossm mover:'
# ossm_file = get_datafile(os.path.join(base_dir,
# r"./MerrimackMassCoastOSSM.txt"))
curr_file = gs.get_datafile(
os.path.join(base_dir, "MerrimackMassCoast.cur"))
tide_file = gs.get_datafile(
os.path.join(base_dir, "MerrimackMassCoastOSSM.txt"))
c_mover = gs.CatsMover(curr_file, tide=gs.Tide(tide_file))
# but do need to scale (based on river stage)
c_mover.scale = True
c_mover.scale_refpoint = (-70.65, 42.58333)
c_mover.scale_value = 1.
model.movers += c_mover
model.environment += c_mover.tide
print 'adding a cats mover:'
curr_file = gs.get_datafile(os.path.join(base_dir, "MassBaySewage.cur"))
c_mover = gs.CatsMover(curr_file)
# but do need to scale (based on river stage)
c_mover.scale = True
c_mover.scale_refpoint = (-70.78333, 42.39333)
# the scale factor is 0 if user inputs no sewage outfall effects
c_mover.scale_value = .04
model.movers += c_mover
# pat1Angle 315;
# pat1Speed 19.44; pat1SpeedUnits knots;
# pat1ScaleToValue 0.138855
#
# pat2Angle 225;
# pat2Speed 19.44; pat2SpeedUnits knots;
# pat2ScaleToValue 0.05121
#
# scaleBy WindStress
print "adding a component mover:"
component_file1 = gs.get_datafile(os.path.join(base_dir, "WAC10msNW.cur"))
component_file2 = gs.get_datafile(os.path.join(base_dir, "WAC10msSW.cur"))
comp_mover = gs.ComponentMover(component_file1, component_file2,
w_mover.wind)
# todo: callback did not work correctly below - fix!
# comp_mover = ComponentMover(component_file1,
# component_file2,
# Wind(timeseries=series, units='m/s'))
comp_mover.scale_refpoint = (-70.855, 42.275)
comp_mover.pat1_angle = 315
comp_mover.pat1_speed = 19.44
comp_mover.pat1_speed_units = 1
comp_mover.pat1ScaleToValue = .138855
comp_mover.pat2_angle = 225
comp_mover.pat2_speed = 19.44
comp_mover.pat2_speed_units = 1
comp_mover.pat2ScaleToValue = .05121
model.movers += comp_mover
print 'adding a spill'
end_time = gs.asdatetime(start_time) + gs.hours(12)
spill = gs.point_line_release_spill(num_elements=100,
start_position=(-70.911432, 42.369142,
0.0),
release_time=start_time,
end_release_time=end_time)
model.spills += spill
return model
grid_topology=topology)
mapfile = 'Waddensea_ijsselmeer_6.bna'
land_map = gs.MapFromBNA(mapfile)
model.map = TideflatMap(land_map, tideflat)
print 'Adding current'
model.movers += current_mover
print 'Adding RandomMover'
model.movers += gs.RandomMover(diffusion_coeff=50000)
print 'adding a wind mover:'
model.movers += gs.constant_wind_mover(**Wind)
print 'adding spill'
spill = gs.surface_point_line_spill(num_elements=100,
start_position=SpillPosition,
end_position=(SpillPosition[0],
SpillPosition[1] - 0.1, 0),
release_time=asdatetime(start_time),
amount=5000,
units='kg',
name='My spill')
model.spills += spill
model.outputters += gs.NetCDFOutput(os.path.join(out_dir, 'test_output.nc'),
which_data='most',
output_timestep=dT_OUT)
def make_model(images_dir=os.path.join(base_dir, 'images')):
# create the maps:
print 'creating the maps'
mapfile = gs.get_datafile(os.path.join(base_dir, './MassBayMap.bna'))
gnome_map = gs.MapFromBNA(mapfile,
refloat_halflife=1, # hours
raster_size=2048 * 2048 # about 4 MB
)
renderer = gs.Renderer(mapfile,
images_dir,
image_size=(800, 800),
projection_class=GeoProjection)
print 'initializing the model'
# start_time = datetime(2013, 3, 12, 10, 0)
start_time = "2013-03-12T10:00"
# 15 minutes in seconds
# Default to now, rounded to the nearest hour
model = gs.Model(time_step=gs.minute(15),
start_time=start_time,
duration=gs.days(1),
map=gnome_map,
uncertain=True)
print 'adding outputters'
model.outputters += renderer
netcdf_file = os.path.join(base_dir, 'script_boston.nc')
gs.remove_netcdf(netcdf_file)
model.outputters += gs.NetCDFOutput(netcdf_file, which_data='all')
model.outputters += gs.KMZOutput(os.path.join(base_dir, 'script_boston.kmz'))
print 'adding a RandomMover:'
model.movers += gs.RandomMover(diffusion_coef=100000)
print 'adding a wind mover:'
# series = np.zeros((2, ), dtype=datetime_value_2d)
# series[0] = (start_time, (5, 180))
# series[1] = (start_time + timedelta(hours=25), (5, 180))
# w_mover = WindMover(Wind(timeseries=series, units='m/s'))
# model.movers += w_mover
# model.environment += w_mover.wind
w_mover = gs.constant_wind_mover(5, 180, units='m/s')
model.movers += w_mover
print 'adding a cats shio mover:'
curr_file = gs.get_datafile(os.path.join(base_dir, r"./EbbTides.cur"))
tide_file = gs.get_datafile(os.path.join(base_dir, r"./EbbTidesShio.txt"))
c_mover = gs.CatsMover(curr_file, tide=gs.Tide(tide_file))
# this is the value in the file (default)
c_mover.scale_refpoint = (-70.8875, 42.321333)
c_mover.scale = True
c_mover.scale_value = -1
model.movers += c_mover
# TODO: cannot add this till environment base class is created
# model.environment += c_mover.tide
print 'adding a cats ossm mover:'
# ossm_file = get_datafile(os.path.join(base_dir,
# r"./MerrimackMassCoastOSSM.txt"))
curr_file = gs.get_datafile(os.path.join(base_dir,
"MerrimackMassCoast.cur"))
tide_file = gs.get_datafile(os.path.join(base_dir,
"MerrimackMassCoastOSSM.txt"))
c_mover = gs.CatsMover(curr_file, tide=gs.Tide(tide_file))
# but do need to scale (based on river stage)
c_mover.scale = True
c_mover.scale_refpoint = (-70.65, 42.58333)
c_mover.scale_value = 1.
model.movers += c_mover
model.environment += c_mover.tide
print 'adding a cats mover:'
curr_file = gs.get_datafile(os.path.join(base_dir, "MassBaySewage.cur"))
c_mover = gs.CatsMover(curr_file)
# but do need to scale (based on river stage)
c_mover.scale = True
c_mover.scale_refpoint = (-70.78333, 42.39333)
# the scale factor is 0 if user inputs no sewage outfall effects
c_mover.scale_value = .04
model.movers += c_mover
# pat1Angle 315;
# pat1Speed 19.44; pat1SpeedUnits knots;
# pat1ScaleToValue 0.138855
#
# pat2Angle 225;
# pat2Speed 19.44; pat2SpeedUnits knots;
# pat2ScaleToValue 0.05121
#
# scaleBy WindStress
print "adding a component mover:"
component_file1 = gs.get_datafile(os.path.join(base_dir, "WAC10msNW.cur"))
component_file2 = gs.get_datafile(os.path.join(base_dir, "WAC10msSW.cur"))
comp_mover = gs.ComponentMover(component_file1, component_file2, w_mover.wind)
# todo: callback did not work correctly below - fix!
# comp_mover = ComponentMover(component_file1,
# component_file2,
# Wind(timeseries=series, units='m/s'))
comp_mover.scale_refpoint = (-70.855, 42.275)
comp_mover.pat1_angle = 315
comp_mover.pat1_speed = 19.44
comp_mover.pat1_speed_units = 1
comp_mover.pat1ScaleToValue = .138855
comp_mover.pat2_angle = 225
comp_mover.pat2_speed = 19.44
comp_mover.pat2_speed_units = 1
comp_mover.pat2ScaleToValue = .05121
model.movers += comp_mover
print 'adding a spill'
end_time = gs.asdatetime(start_time) + hours(12)
spill = point_line_release_spill(num_elements=100,
start_position=(-70.911432,
42.369142, 0.0),
release_time=start_time,
end_release_time=end_time)
model.spills += spill
return model