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