Nesting and Periodic boundaries not working correctly

nesting
, ,
1

Hi,

This is an issue that Patrick Marchesiello solved for me ~12 years ago in an old version of ROMS_Agrif.
However, this patch was apparently not applied to the Croco model:

I have downloaded the last version of Croco and the nesting is not coded correctly when using periodic boundary conditions. To illustrate the error, I have a run an experiment of an eddy propagating in a periodic channel (EW_Periodic, North/South walls) using a nested grid (child grid is shown as a magenta frame in the figure below).
As you can see after a couple of days (2nd figure), the EW_PERIODIC is applied to the child grid while it should NOT.

If anyone wants to reproduce this example quickly (runs in a couple of minutes using 4 cpus) a simple cppdefs.h, param.h, AGRIF_FixedGrids.in, croco_grd.nc, croco_grd.nc.1, croco_ini.nc, croco_ini.nc.1, Croco.in, Croco.in.1 can be downloaded from my webpage here:

https://www.vincentcombes.com/TEST_EW_PERIODIC.zip

Initial Condition (SSH):

image description

After a couple of days (Error):

image description

It would be great if it could be fixed in the Croco code

Thanks to the development team

Vincent

2

Hi everyone,

After 1 year of silence, I am back with this work.

Has anyone fixed that issue ?

Thank you

Vincent

3

Hi CROCO developers,

I am still interesting in a solution to this problem.
I have tried many things and thought I was able to fix it, but it seems not …
Especially in a real case experiment.

Thank you for any help on this !

Vincent

4

Hi Vincent, Sorry for not responding earlier to your post. We are not planning to include this feature (periodic nest) in CROCO because it is very intrusive and requests duplicating many routines, which may induce other problems with some modules. For these reasons it is not planned to include that. Sorry.

5

Thank you,

I agree that there will be many files to update if we want this configuration to work for all cases.
Perhaps it could be possible to fix it in a simpler case (i.e., no atmospheric coupling, no sediments, no waves, etc.), focusing only on the ROMS ocean folder.
I have tried quite a lot, but unfortunately, without success.

I am open to any kind of help

Vincent