Limit maximum deep of a particle

[Elisabet Cruz]

Hi everyone!

Has anyone tried to limit the maximum deep that a particle can reach? For example, my particles can not go deeper than 40 m depth.

I tried setting this through the vertical migration flag, using as my probability matrix:

1

1

40

10454400 %maximum time of my simulation, which corresponds to 121 days in seconds

100

But instead of limiting the maximum depth, its moves all my particles to 40 m in the second time index.

Any thoughts on how to limit the maximum deep?

Thanks in advanced,

Elisabet

[Robin Faillettaz]

Hi Elisabet, 

I don’t know if that’s feasible right now in CMS. A workaround would be to cut your netcdf files down to 40 m and set the flag avoid coast to true.

Regards, 

— —
Robin Faillettaz

IFREMER

Unité HALGO — UMR DECOD
Laboratoire de Technologie et de Biologie Halieutique
8, rue François Toullec
56 100 Lorient 
France

Tel: +33 (0)2 97 87 38 12
e-mail : robin.fa...@ifremer.fr 

[nleco...@gmail.com]

Hi Elisabet,

I think you could use the vertical probability matrix, but instead of using one layer, you could put all the layers you want the particles to be. Something like:

1

5

10 20 30 40 50

10454400 %maximum time of my simulation, which corresponds to 121 days in seconds

25

25

25

25

0

Your particles would be in the first 40m (with a ~25% chance to be at 10, 20, 30 or 40 m) and no particles deeper (as the percentage for 50 is set to 0).

Hope this help

[Robin Faillettaz]

Hi, 

That’s a solution too indeed, but this will overwrite the vertical diffusion module so it may not be appropriate depending on what she is trying to model (it forces X% of the particles to be in the defined depth layers). 

Elisabet, we may be able to help you better if you explain us into more details what to are trying to model.

Best, 

[Elisabet Cruz]

Hi!

Thank you for answering so quickly.

I am modeling HAB species, more specifically Dinophysis spp.

I also thought about the solution with the vertical matrix, but the problem is the one that Robin said and also that I am forcing to have a certain number of particles at those depths during all my simulation, which is not so representative for the HAB species. 

Finally, my supervisor and I come about with a solution setting a condition in the "updateloc.f90" program, generating a random number and forcing the particles to do not go deeper than 40±2 m.

Best,

Elisabet