Opened 11 years ago
Closed 9 years ago
#1127 closed Task (fixed)
INGV-2 Action in the 2013 WP
Reported by: | poddo | Owned by: | nemo |
---|---|---|---|
Priority: | low | Milestone: | |
Component: | OCE | Version: | v3.4 |
Severity: | Keywords: | ||
Cc: |
Description
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
Commit History (0)
(No commits)
Attachments (13)
Change History (15)
Changed 11 years ago by poddo
Changed 11 years ago by poddo
Changed 11 years ago by poddo
Changed 11 years ago by poddo
Changed 11 years ago by poddo
Changed 11 years ago by poddo
Changed 11 years ago by poddo
Changed 11 years ago by poddo
Changed 11 years ago by poddo
Changed 11 years ago by poddo
Changed 11 years ago by poddo
comment:1 Changed 11 years ago by 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.
comment:2 Changed 9 years ago by clevy
- Resolution set to fixed
- Status changed from new to closed
namelist exp1