You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
New config setting drift:horizontal_diffusivity, providing time-step independent diffusion, in contrast to drift:current_uncertainty and drift:wind_uncertainty
* New config setting ``drift:horizontal_diffusivity``, providing time-step independent diffusion, in contrast to ``drift:current_uncertainty`` and ``drift:wind_uncertainty``
so I have been planning to use the non-timestep-dependent diffusion figuring that was better. Plus it has the advantage of being able to come up with a semi-physically-based estimate for the horizontal diffusivity of the velocity not being captured in a model times the horizontal grid resolution.
However, now I am getting set up to do some oil simulations and I see it is hard-wired to have:
Yes, current_uncertainty and wind_uncertainty is the legacy from early on (adding random perturbations to winds and currents with given standard deviations). It was first implemented because I thought it had a natural interpretation. In OpenOil the hardwired values remain, but they should be replaced with horizontal_diffusivity, to avoid the time step dependency.
However, the values used are anyway a bit hand-waving (the default uncertainties are randomly chosen and can be modified), so therefore it may not be strictly wrong to use the uncertainties (standard deviation perturbations). Also, one can also interpret the uncertainties differently ("ocean/atmospheric model error and bias") than diffusivity ("unresolved"), and it might make sense to use both mechanisms. Also, there is a plan to implement higher order diffusivity as a Markov chain, where particles have some "memory".
There was a interesting similar discussion in this thread: #1159
I see the configuration setting for horizontal diffusivity
and the note
from
opendrift/history.rst
Line 290 in 34fd858
However, now I am getting set up to do some oil simulations and I see it is hard-wired to have:
in OpenOil initialization.
My questions are:
horizontal_diffusion
(or perhaps they were just set up before the newer capability was available)?horizontal_diffusion
vscurrent_uncertainty
andwind_uncertainty
?Thank you.
The text was updated successfully, but these errors were encountered: