Exemplo n.º 1
0
o = BivalveLarvae(loglevel=0)  # Set loglevel to 0 for debug information
# o = PelagicEggDrift(loglevel=20)  # Set loglevel to 0 for debug information

####################################################################################################################
# READERS
#
# Atmospheric model for wind
# reader_arome = reader_netCDF_CF_generic.Reader(o.test_data_folder() +
#     '16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')

# Ocean model for current
# reader_norkyst = reader_netCDF_pytCF_generic.Reader(o.test_data_folder() +
#     '16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
# # Creating and adding reader for Nordic 4km current dataset
nordic_native = reader_ROMS_native.Reader(
    o.test_data_folder() +
    '2Feb2016_Nordic_sigma_3d/Nordic-4km_SLEVELS_avg_00_subset2Feb2016.nc')
o.add_reader(nordic_native)

# default values for forcing variables
o.set_config('environment:fallback:ocean_vertical_diffusivity', 0.001)  # m2/s
####################################################################################################################

####################################################################################################################
# SEEDING
# time = datetime.utcnow()
o.seed_elements(lon=12.0,
                lat=68.3,
                z=-3.0,
                radius=100,
                number=100,
Exemplo n.º 2
0
o = BivalveLarvae(loglevel=0)  # Set loglevel to 0 for debug information
# o = PelagicEggDrift(loglevel=20)  # Set loglevel to 0 for debug information

####################################################################################################################
# READERS
# 
# Atmospheric model for wind
# reader_arome = reader_netCDF_CF_generic.Reader(o.test_data_folder() +
#     '16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')

# Ocean model for current
# reader_norkyst = reader_netCDF_pytCF_generic.Reader(o.test_data_folder() +
#     '16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
# # Creating and adding reader for Nordic 4km current dataset
nordic_native = reader_ROMS_native.Reader(o.test_data_folder() +
    '2Feb2016_Nordic_sigma_3d/Nordic-4km_SLEVELS_avg_00_subset2Feb2016.nc')
o.add_reader(nordic_native)

# default values for forcing variables
o.set_config('environment:fallback:ocean_vertical_diffusivity',0.001) # m2/s
####################################################################################################################

####################################################################################################################
# SEEDING
# time = datetime.utcnow()
o.seed_elements(lon=12.0, lat=68.3, z=-3.0,
                radius=100, number=100,
                time=nordic_native.start_time,
                terminal_velocity = 0.0)