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) renderer = gs.Renderer( output_timestep=dT_OUT, map_filename=mapfile, output_dir=out_dir, formats=['gif'], # ['gif', 'png'] image_size=(1280, 1024), viewport=((4.5, 53.0), (5.5, 53.5))) model.outputters += renderer gs.set_verbose() # get some more info as the model runs print "Doing a full run:" model.full_run()
# from gnome.basic_types import datetime_value_2d from gnome.utilities.projections import GeoProjection # from gnome.utilities.remote_data import get_datafile # from gnome.environment import Wind, Tide # from gnome.maps import MapFromBNA # from gnome.model import Model # from gnome.spill import point_line_release_spill # from gnome.movers import RandomMover, WindMover, CatsMover, ComponentMover # from gnome.outputters import Renderer, NetCDFOutput, KMZOutput # let's get the console log working: gs.set_verbose() # define base directory base_dir = os.path.dirname(__file__) 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
from gnome.utilities.projections import GeoProjection # from gnome.utilities.remote_data import get_datafile # from gnome.environment import Wind, Tide # from gnome.map import MapFromBNA # from gnome.model import Model # from gnome.spill import point_line_release_spill # from gnome.movers import RandomMover, WindMover, CatsMover, ComponentMover # from gnome.outputters import Renderer, NetCDFOutput, KMZOutput # let's get the console log working: gs.set_verbose() # define base directory base_dir = os.path.dirname(__file__) 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 )