INGV-2 Action in the 2013 WP

[poddo]

INGV-2 Action in the 2013 WP.
Perform basic tests.
The first test "conservation Box".
The basic Idea is to initialize the model with constant Temperature and Salinity
to check conservation and spurious accelerations depending on different model settings.
The test is done in a flat box and in the sea-mountain case
Starting from Gyre set-up also:
-no-slip condition;
-vorlat;
-energy & enstrophy scheme;
-centered advection;
-bilaplacian operators
have been checked

Step1 : Retrive latest version of the code

svn --username "XXXX" co -r 3819 ​http://forge.ipsl.jussieu.fr/nemo/svn/branches/2012/dev_v3_4_STABLE_2012

Step2: create my Conservation Box configuration

./makenemo -m MyMAC -r GYRE -n CONS_BOX
!!!! Copy and editing some conf files
cp ~OPA_SRC/par_GYRE.h90 ~/MY_SRC/par_CONS_BOX.h90 (attached)
starting from the GYRE_par file here I use:

• jphgr_msh = 1
• ppglam0 = 14.0_wp
• ppgphi0 = 42.0_wp
• ppe1_deg = 0.0625
• ppe2_deg = 0.0625

cp ~/OPA_SRC/par_oce.F90 ~/MY_SRC/par_oce.F90
Simply include the par_CONS_BOX.h90 (not attached)

cp ~/OPA_SRC/DOM/istate.F90 ~/MY_SRC/istate.F90
use constant T and not the profile (philander) (attached)

cp ~/OPA_SRC/DOM/domzgr.F90 ~/MY_SRC/domzgr.F90
change the sea-mountain parameter to work with deg (attached)

Exp1=namelist identical to gyre (attached)
Exp$...
Exp2=namelist with all the differences (attached)
Exp3=as Exp2 but with sea-mountain

From Exp1 to Exp2 the results are consistent, T and S do not vary with time


SSH is zero and there are no velocity

In Exp3 (sea-mountain) See Fig for the bathy

1) I have to change the timestep from 7200 to 300 (see the namelist)
2) there are spurious velocities and bizarre SSH patterns

[poddo]

namelist exp1

[poddo]

namelist exp2

[poddo]

namelist exp3

[poddo]

In the sea mountain case, the results are sensible to the EOS used.
The problem is still there but.
With nn_neos=2
at kt=1 all momentum acc are zero then at kt=2

hpg - Ua: 0.1540743955509789D-32 Va: 0.9190729483525867D-16

With neos=1
from kt=1 to kt=5 zero acc then at kt=6

hpg - Ua: 0.0000000000000000D+00 Va: -0.1724182806109115D-16

With neos=0
at kt=9

hpg - Ua: 0.0000000000000000D+00 Va: -0.1396655148263087D-15

I added a couple of missing "_wp" in eos but this did not solve the problem.