OpenDrift / opendrift

Open source framework for ocean trajectory modelling
https://opendrift.github.io
GNU General Public License v2.0
249 stars 120 forks source link

ROMS reader #742

Open efredj opened 3 years ago

efredj commented 3 years ago

I recently download the opendrift package. I was able to install it and run all the examples. Then I setup a simulation over the Mid Atlantic Bight region using the DOPPIO thredds server data following example from the tutorial.

https://opendrift.github.io/gallery/example_roms_native.html

Everything seems to run correctly: 12:46:57 INFO opendrift.models.basemodel: Fallback values will be used for the following variables which have no readers: 12:46:57 INFO opendrift.models.basemodel: x_wind: 0.000000 12:46:57 INFO opendrift.models.basemodel: y_wind: 0.000000 12:46:57 INFO opendrift.models.basemodel: sea_surface_wave_significant_height: 0.000000 12:46:57 INFO opendrift.models.basemodel: sea_surface_wave_stokes_drift_x_velocity: 0.000000 12:46:57 INFO opendrift.models.basemodel: sea_surface_wave_stokes_drift_y_velocity: 0.000000 12:46:57 INFO opendrift.models.basemodel: sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000 12:46:57 INFO opendrift.models.basemodel: sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000 12:46:57 INFO opendrift.models.basemodel: turbulent_kinetic_energy: 0.000000 12:46:57 INFO opendrift.models.basemodel: turbulent_generic_length_scale: 0.000000 12:46:57 INFO opendrift.models.basemodel: ocean_mixed_layer_thickness: 50.000000 12:46:57 INFO opendrift.models.basemodel: Adding a dynamical landmask with max. priority based on assumed maximum speed of 1 m/s. Adding a customised landmask may be faster... 12:46:58 INFO opendrift.models.basemodel: Using existing reader for land_binary_mask 12:46:58 INFO opendrift.models.basemodel: All points are in ocean 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 01:00:00 - step 1 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 02:00:00 - step 2 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 03:00:00 - step 3 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 04:00:00 - step 4 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 05:00:00 - step 5 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 06:00:00 - step 6 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 07:00:00 - step 7 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 08:00:00 - step 8 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 09:00:00 - step 9 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 10:00:00 - step 10 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-16 11:00:00 - step 11 of 48 - 100 active elements (0 deactivated) show more (open the raw output data in a text editor) ...

12:46:58 INFO opendrift.models.basemodel: 2018-11-17 20:00:00 - step 44 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-17 21:00:00 - step 45 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-17 22:00:00 - step 46 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-17 23:00:00 - step 47 of 48 - 100 active elements (0 deactivated) 12:46:58 INFO opendrift.models.basemodel: 2018-11-18 00:00:00 - step 48 of 48 - 100 active elements (0 deactivated)

However the final graph doesn't plot the trajectory see the jupyter-notebook attached?

test_demo_doppio.zip

Thank you for your help

It is a great package.

Erick

knutfrode commented 3 years ago

Hi Erick,

I apologize for the late reply to this issue. The problem here is that you have seeded elements in november 2018, whereas this reader (ROMS dataset) only covers 12 hours on a specific day (23rd October 2021).

The example scripts uses loglevel=20 (only >=INFO), but I recommend changing this to 0, to show also debug output. You will then see lines like this:

16:08:26 DEBUG   opendrift.models.basemodel: Calling reader roms native
16:08:26 DEBUG   opendrift.models.basemodel: ----------------------------------------
16:08:26 DEBUG   opendrift.models.basemodel:    Outside time coverage of reader.
16:08:26 DEBUG   opendrift.models.basemodel: ----------------------------------------

The OceanDrift-model has a default fallback value of 0 m/s current, and the simulation will thus not crash, but will use this value instead.

erickfredj commented 3 years ago

Thank you for your answer