Seafloor oilspill rate and total amount released #1085
Replies: 19 comments 40 replies
-
|
You can seed with The total mass is divided by all particles initially. |
Beta Was this translation helpful? Give feedback.
-
|
To continue on the topic. Do backward trajectories work in the oilspill module? |
Beta Was this translation helpful? Give feedback.
-
|
Yes it works, but weathering (eg evaporstion) is then deactivated
…On Tue, Apr 18, 2023, 22:01 Andres Sepulveda ***@***.***> wrote:
To continue on the topic. Do backward trajectories work in the oilspill
module?
—
Reply to this email directly, view it on GitHub
<#1085 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ABH25I76Y633LSBDAQQM3QLXB3XK3ANCNFSM6AAAAAAW6STDCI>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
Beta Was this translation helpful? Give feedback.
-
|
Hi, |
Beta Was this translation helpful? Give feedback.
-
|
Thanks a lot, I am beginner of OpenDrift |
Beta Was this translation helpful? Give feedback.
-
|
Hi, Knut-Frode 1) SIMULATION 1: HOW ADDING THE EFFECT OF SUCH WIND: 2) SIMULATION 2: HOW TO USE THE WIND ALREADY SERVED IN THE HYDRODYNAMIC MODEL AS FORCING Many thanks in advance |
Beta Was this translation helpful? Give feedback.
-
|
So you want to make and compare two simulations with respectively:
To use the wind used to force CROCO, this wind must either be included in the same file as the currents. This you can check with "print r", where r is the created reader (or "readerinfo " from commandline). Alternatively, if you have this wind in a separate file, you can add this as well with another reader. For the run with measured wind, you can add first the reader_timeseries as described above, with your own numbers, and secondly add the CROCO reader. Then the measured wind will be used instead of the wind eventually inlcuded in the CROCO files. |
Beta Was this translation helpful? Give feedback.
-
|
Hi, the wind speed can be found in the BLUK files that CROCO uses as forcing. The variables are named uwnd and vwnd (in m/s) |
Beta Was this translation helpful? Give feedback.
-
|
Thanks Knut-Frode and Andres, it works. Hereafter the code. ``from datetime import timedelta o = OpenOil(loglevel=20) filename_nc = 'croco_Sim01_1month.nc' t0 = datetime(2000, 1, 1) time_seed = [reader_norkyst.start_time, o.set_config('seed:droplet_diameter_min_subsea', 0.0001) o.run(duration=timedelta(days=3), time_step=60, time_step_output=60,outfile="spill_Tunisia_3days.nc")`` |
Beta Was this translation helpful? Give feedback.
-
|
Hi,
|
Beta Was this translation helpful? Give feedback.
-
|
It should be something like this... No need to extract the variables from the BLK file. The only question is if the reader_ROMS_native accept CROCO's uwnd and vwnd as valid. For what I see it expectz Uwind and Vwind. from opendrift.readers import reader_ROMS_native o = Leeway(loglevel=0) filename_his = 'croco_his.nc'; mosa_his = reader_ROMS_native.Reader(filename_his) o.add_reader([mosa_his, variables=['x_sea_water_velocity', 'y_sea_water_velocity']) o.add_reader([mosa_blk, variables=['x_wind', 'y_wind']) |
Beta Was this translation helpful? Give feedback.
-
|
The variable mapping in ROMS reader is hardcoded, as https://github.com/OpenDrift/opendrift/blob/master/opendrift/readers/reader_ROMS_native.py#L37 |
Beta Was this translation helpful? Give feedback.
-
|
I start by modifying a little the "reader_ROMS_native.py" as mentioned above `from datetime import timedelta o = Leeway(loglevel=0) filename_his = 'croco_Sim01_1month.nc' filename_blk = 'croco_blk_ERA5_Y2010M01.nc' o.add_reader([reader_norkyst_his], variables=['time']) However, the following instruction gives an error: I think I have first to read the gridfile ! I have another issue with time, how define the start time and time duration Thank you very much for your time ! |
Beta Was this translation helpful? Give feedback.
-
|
Hi, I tried to do the following to define the reader_norkyst_blk by: This last instruction gave the following error: I did the modification as: Now the error is : I did try to read time from "_his" file: But the problem persist ... |
Beta Was this translation helpful? Give feedback.
-
|
Please show the complete output of the error message? |
Beta Was this translation helpful? Give feedback.
-
File "/home/bejaoui/miniconda3/envs/opendrift/lib/python3.11/site-packages/xarray/core/utils.py", line 455, in getitem During handling of the above exception, another exception occurred: Traceback (most recent call last): File "/home/bejaoui/miniconda3/envs/opendrift/lib/python3.11/site-packages/xarray/core/utils.py", line 455, in getitem
|
Beta Was this translation helpful? Give feedback.
-
|
Here ncdump for _his file: (base) bejaoui@instm-wsz640:/scratch/CROCO/croco_useful_bb/Lagrangian_Andres_Tun/opendrift/examples$ ncdump -h croco_Sim01_1month.nc // global attributes: |
Beta Was this translation helpful? Give feedback.
-
|
Time variable seems to be named |
Beta Was this translation helpful? Give feedback.
-
|
Hi,
The message error is:
During handling of the above exception, another exception occurred: Traceback (most recent call last): |
Beta Was this translation helpful? Give feedback.
-
|
Ok, I think the error is because in CROCO BULk files the time variable is named bulk_time not time. So use NCO scripts to rename ir. This should work ncrename -h -O -v bulk_time,time croco_blk.nc |
Beta Was this translation helpful? Give feedback.
-
|
Hi, When I execute the following after changing bulk_time by time I get a new error ... This is done, knutfrode already updated the code to be able to read bulk_time |
Beta Was this translation helpful? Give feedback.
-
|
The line numbers in the error message do not seem to match the latest version of ROMS reader: Did you manage to update code to the latest version? Or perhaps you have local modifications that make a conflict? You can check with |
Beta Was this translation helpful? Give feedback.
-
|
Hi, Knut, |
Beta Was this translation helpful? Give feedback.
-
|
If you can share the bulk- and the grid-files, then I could have a look at them |
Beta Was this translation helpful? Give feedback.
-
|
yes of course, it took a lot of time. Here the link Here my email adress, B |
Beta Was this translation helpful? Give feedback.
-
|
The unit of time in this file was At least it works for me now, and I can do e.g.
|
Beta Was this translation helpful? Give feedback.
-
|
I now discovered that the file contains a global attribute This attribute could of course be read dynamically, but it does not seem like a good (standard) practice to have the zero time specified this way. |
Beta Was this translation helpful? Give feedback.
-
|
Thank you for your help, very appreciated
After executing the line : **reader_norkyst_blk = reader_ROMS_native.Reader(filename_blk, gridfile='croco_grd.nc') During handling of the above exception, another exception occurred: Traceback (most recent call last): |
Beta Was this translation helpful? Give feedback.
-
|
The line numbers in this error message do not match the latest version of |
Beta Was this translation helpful? Give feedback.
-
|
I am in pwd
Updating 8b4ee8b..b1cc24e What is the terminal command, sorry for the simple question and thanks a lot |
Beta Was this translation helpful? Give feedback.
-
|
If you are happy to discard your local changes (they should mostly be in my corresponding updates), you can do |
Beta Was this translation helpful? Give feedback.
-
|
Of course, with thanks
then
Updating 8b4ee8b..b1cc24e |
Beta Was this translation helpful? Give feedback.
-
|
This seems ok. Then it should hopefully work to read the Croco BLK file. |
Beta Was this translation helpful? Give feedback.
-
|
gives **>>> reader_norkyst_blk = reader_ROMS_native.Reader(filename_blk, gridfile='croco_grd.nc') During handling of the above exception, another exception occurred: Traceback (most recent call last): |
Beta Was this translation helpful? Give feedback.
-
|
This is quite strange, as the same works for me. |
Beta Was this translation helpful? Give feedback.
-
|
Many thanks, I am doing the update ... I will happy to inform you |
Beta Was this translation helpful? Give feedback.
-
|
Hi the command
gives the following in terminal : (base) bejaoui@instm-wsz640:/scratch/CROCO/croco_useful_bb/Lagrangian_Andres_Tun/opendrift$ pip install -e . |
Beta Was this translation helpful? Give feedback.
-
|
I see that you are in "conda base environment".
You must activate the opendrift environment with
`conda activate opendrift`
…On Fri, Apr 28, 2023, 18:22 Béchir Béjaoui ***@***.***> wrote:
Hi the command
pip install -e .
gives the following in terminal :
*(base)
***@***.***:/scratch/CROCO/croco_useful_bb/Lagrangian_Andres_Tun/opendrift$
pip install -e . ERROR: File "setup.py" or "setup.cfg" not found. Directory
cannot be installed in editable mode:
/scratch/CROCO/croco_useful_bb/Lagrangian_Andres_Tun/opendrift (A
"pyproject.toml" file was found, but editable mode currently requires a
setuptools-based build.)*
—
Reply to this email directly, view it on GitHub
<#1085 (reply in thread)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ABH25IZWVRSB535Y4LAAU4DXDPVFHANCNFSM6AAAAAAW6STDCI>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
Beta Was this translation helpful? Give feedback.
-
|
is this means Ok even there is an error ! Attempting uninstall: dask |
Beta Was this translation helpful? Give feedback.
-
|
Hi, after this >>> filename_blk = 'croco_blk_ERA5_Y2010M01_andres.nc' this gives in terminal : >>> reader_norkyst_blk = reader_ROMS_native.Reader(filename_blk, gridfile='croco_grd.nc') is it ok, can I go forward ? |
Beta Was this translation helpful? Give feedback.
-
|
Yes, then it is ok. Then you can do either Btw, NorKyst is a specific Norwegian ocean model used in the examples, so a better name for your reader would be |
Beta Was this translation helpful? Give feedback.
-
|
Hi,
it gives a TypeError, following |
Beta Was this translation helpful? Give feedback.
-
|
Yes, you should add the reader simply like this: |
Beta Was this translation helpful? Give feedback.
-
|
Hi Knute one more clarification, >>> o.set_config('drift:vertical_mixing', True) >>> o.seed_elements(lon=10.124, lat=33.953, z='seafloor', m3_per_hour=30, radius=0, number=7200, 18:10:54 INFO opendrift.models.leeway:209: Seeding elements of object type 1: PIW-1 (Person-in-water (PIW), unknown state (mean values)) |
Beta Was this translation helpful? Give feedback.
-
|
It seems you have created a Leeway simulation object, but this is for
search and rescue.
You should instead use OpenOil.
You will find several examples in the gallery.
…On Fri, Apr 28, 2023, 19:14 Béchir Béjaoui ***@***.***> wrote:
Hi Knute one more clarification,
1.
*>>> o.set_config('drift:vertical_mixing', True)*
Traceback (most recent call last):
File "", line 1, in
File
"/scratch/CROCO/croco_useful_bb/Lagrangian_Andres_Tun/opendrift/opendrift/models/basemodel.py",
line 660, in set_config
raise ValueError('No config setting named %s' % key)
*ValueError: No config setting named drift:vertical_mixing*
1.
*>>> o.seed_elements(lon=10.124, lat=33.953, z='seafloor', m3_per_hour=30,
radius=0, number=7200, ... time=time_seed, oil_type='GAS OIL, EXXON')*
18:10:54 INFO opendrift.models.leeway:209: Seeding elements of object type
1: PIW-1 (Person-in-water (PIW), unknown state (mean values))
Traceback (most recent call last):
File "", line 1, in
File
"/scratch/CROCO/croco_useful_bb/Lagrangian_Andres_Tun/opendrift/opendrift/models/leeway.py",
line 274, in seed_elements
super(Leeway, self).seed_elements(lon, lat,
File
"/scratch/CROCO/croco_useful_bb/Lagrangian_Andres_Tun/opendrift/opendrift/models/basemodel.py",
line 1854, in seed_elements
seafloor_constant = self.get_config(
^^^^^^^^^^^^^^^^
File
"/scratch/CROCO/croco_useful_bb/Lagrangian_Andres_Tun/opendrift/opendrift/models/basemodel.py",
line 706, in get_config
raise ValueError('No config setting named %s' % key)
*ValueError: No config setting named
environment:constant:sea_floor_depth_below_sea_level*
—
Reply to this email directly, view it on GitHub
<#1085 (reply in thread)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ABH25I4UUGDP7HF23YDUOYDXDP3IVANCNFSM6AAAAAAW6STDCI>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
Beta Was this translation helpful? Give feedback.
-
|
Hi Knut,
Hereafter the whole code: o = OpenOil(loglevel=20) # Set loglevel to 0 for debug information filename_his = 'croco_his.nc' filename_blk = 'croco_blk_ERA5_Y2010M01.nc' #Controle the bulk file #plot test of wind speed (do this or this) o.set_config('drift:vertical_mixing', True) time Seeding particlest0_seed = datetime(2010, 1, 1) If there are 720 Tons, and we use 7200 particles, each particle represnts 0.1 Ton or 100 liters of gas oilGas Oil: 850-870 kg/m3Water: 1000 kg/m3Gasoline: 715-780 kg/m3o.seed_elements(lon=10.124, lat=33.953, z='seafloor', m3_per_hour=30, radius=0, number=7200, Setting the range of droplet sizes for the seafloor releaseo.set_config('seed:droplet_diameter_min_subsea', 0.0001) #Running model with a small timestep to resolve the boyant rising #Plot results Slow plot (to save)#o.animation(filename="spill_Tunisia_3days_horizontal.mp4",fast=True) The message in terminal is AttributeError: 'ReaderBlock' object has no attribute 'interpolator1d' ` |
Beta Was this translation helpful? Give feedback.
-
|
18:08:29 INFO opendrift.models.basemodel:1279: ======================== |
Beta Was this translation helpful? Give feedback.
-
|
I am trying this code with the corrections to the reader mentioned before
I had to add the lat_rho and lon_rho variables to the BLK file using NCO, since the option
gave me problems However the output indicates that wind is not to be found 16:04:04 WARNING opendrift.models.basemodel: Missing variables: ['x_wind', 'y_wind', 'sea_surface_wave_significant_height'] What could be the issue? |
Beta Was this translation helpful? Give feedback.
-
|
Dear ALL, o = OpenOil(loglevel=20) # Set loglevel to 0 for debug information #--------------- #filename_his = 'croco_his.nc' #Controle the bulk file #-------------- filename_blk = 'croco_blk_ERA5_Y2010M04.nc' #Controle the bulk file #plot test of wind speed (do this or this) #----------------- o.add_reader([reader_croco_blk,reader_croco_his]) o.set_config('drift:vertical_mixing', True) #time Seeding particles #time_seed = [reader_croco_blk.start_time, #time_seed = [reader_croco_his.start_time, #If there are 720 Tons, and we use 7200 particles, each particle represnts 0.1 Ton or 100 liters of gas oil o.seed_elements(lon=10.124, lat=33.953, z='seafloor', m3_per_hour=30, radius=0, number=7200, #Setting the range of droplet sizes for the seafloor release #Running model with a small timestep to resolve the boyant rising #------------- #Plot results #Plot results with lines and particles colored by depth #Plot vector fields #animation Now could you please tell how :
Many thanks for your support |
Beta Was this translation helpful? Give feedback.
-
|
You can zoom by adding To get bathymetry as background: |
Beta Was this translation helpful? Give feedback.
-
|
Thank you Knut, |
Beta Was this translation helpful? Give feedback.
-
|
This how to add string on the map : |
Beta Was this translation helpful? Give feedback.
-
|
Hi all,
took a lot of time, could you please explain if the came from my workstation or from the calculation ?
|
Beta Was this translation helpful? Give feedback.
-
Hi,
I want to study a case of a seafloor oilspill, a sunken boat with 10 Tons of Gas Oil
How do I define the release rate from the sunken boat?
How do I set the case so that only 10 Tons are released? With the seeding time? Each particle has a volumen that I can relate to the oil density?
Any help is appreciated.
Andres
Beta Was this translation helpful? Give feedback.
All reactions