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()
Example #2
0
# 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
Example #3
0
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
                              )