New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
WorkingGroups/TestCases/ImplementationStrategyAndTemplate – NEMO
wiki:WorkingGroups/TestCases/ImplementationStrategyAndTemplate

Version 16 (modified by flavoni, 9 years ago) (diff)

--

How to create a TEST CASE: STRATEGY

in https://forge.ipsl.jussieu.fr/nemo/browser/branches/2015/dev_r5102_CNRS11_TestCase/NEMOGCM/TOOLS/TEST_CASES there are:

CREATE_FILES directory : python scripts to create

  1. input files: 
    1. coordinates.nc 
    2. bathy_r025.nc
  2. initial state:
    1. TScst_init_r025.nc
    2. TS_init_tprof_r025.nc
  3. forcing file:
    1. flx_r025.nc

a detailed TEST CASE: EQUATORIAL BOX

in https://forge.ipsl.jussieu.fr/nemo/browser/branches/2015/dev_r5102_CNRS11_TestCase/NEMOGCM/TOOLS/TEST_CASES/TEST_EQ_BOX01

and some examples to create a test case of this type :

Equatorial BOX, grid of 0.25 ° resolution, between latitude: -30°,30°, and longitude 0°,180°,

with zps coordinates,

with ubs schema on dynamics and on tracers, with constant TS initial state, with flux only wind (u composant), with flat bathymetry.

If you want to create a test case what you have to do?

  1. go into (branches/2015/dev_r5102_CNRS11_TestCase/)NEMOGCM/TOOLS/TEST_CASES/CREATE_FILES
  2. read README file
  3. run python make_coordinates.py to generate coordinates_r025.nc
  4. link coordinates.nc to coordinates_r025.nc
  5. run python make_bathy.py to generate bathy_meter.nc -> bathy_r025.nc
  6. link bathy_meter.nc to bathy_r025.nc
  7. run python make_TScst_init.py to generate TScst_init_r025.nc (you do not need link, filename is written in namelist_cfg file)
  8. run python make_FLX.py to generate flx_r025.nc (you do not need link, filename is written in namelist_cfg file)

You choose in namelist_cfg cn_exp name and nitend.


NOTA BENE: if you change resolution (in python scripts) you need to change also in namelist_cfg:

jpidta, jpjdta and jpkdta jpiglo and jpjglo '

1. TEST_EQ_BOX01_noCoriolis : Equatorial box (180°,60°) : with NO stratification / constant wind trades / NO Coriolis
options values
resolution 0.25° (1/4)
longitude lon_T=[ 0,180 ]
latitude lat_T=[-31,31] (to have one point at 0° Equator)
vertical levels 31
batymetry 5000m
forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
NO stratification so=35 / thetao=10
total time steps 3600 (=50days)
output frequency 5 days

-

2. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with NO stratification / constant wind trades / e1=e2=constant
options values
resolution 0.25° (1/4)
longitude lon_T=[ 0,180 ]
latitude lat_T=[-31,31] (to have one point at 0° Equator)
vertical levels 31
batymetry 5000m
forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
NO stratification so=35 / thetao=10
total time steps 3600(=50days)
output frequency 5 days

-

3. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with NO stratification / constant wind trades / e1, e2 variable with longitude
options values
resolution 0.25° (1/4)
longitude lon_T=[ 0,180 ]
latitude lat_T=[-31,31] (to have one point at 0° Equator)
vertical levels 31
batymetry 5000m
forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
NO stratification so=35 / thetao=10
total time steps 3600(=50days)
output frequency 5 days

-

4. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with stratification / constant wind trades
options values
resolution 0.25° (1/4)
longitude lon_T=[ 0,180 ]
latitude lat_T=[-31,31] (to have one point at 0° Equator)
vertical levels 31
batymetry 5000m
forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
NO stratification so=35 / thetao=10
total time steps 3600(=50days)
output frequency 5 days

-

5. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with stratification / constant wind trades
options values
resolution 0.25° (1/4)
longitude lon_T=[ 0,180 ]
latitude lat_T=[-31,31] (to have one point at 0° Equator)
vertical levels 31
batymetry 5000m
forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
NO stratification so=35 / thetao=10
total time steps 3600(=50days)
output frequency 5 days

-

6. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with stratification / constant winds (profile)
options values
resolution 0.25° (1/4)
longitude lon_T=[ 0,180 ]
latitude lat_T=[-31,31] (to have one point at 0° Equator)
vertical levels 31
batymetry 5000m
forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
NO stratification so=35 / thetao=10
total time steps 3600(=50days)
output frequency 5 days

-

7. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with stratification / variable wind (Kelvin waves)
options values
resolution 0.25° (1/4)
longitude lon_T=[ 0,180 ]
latitude lat_T=[-31,31] (to have one point at 0° Equator)
vertical levels 31
batymetry 5000m
forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
NO stratification so=35 / thetao=10
total time steps 3600 (restart from a 1801 time step, restart file from test 6) with winds = 0
output frequency 5 days

output plot:current-temp-1S1N-all.pdf

Attachments (1)

Download all attachments as: .zip