TOC(heading=Robustness-TestCases,Robustness-TestCases/ImplementationStrategyAndTemplate/EquatorialBox/*, depth=2)?

TEST CASE: EQUATORIAL BOX

Motivations:

This equatorial box is done to :

• show how to create in NEMO a test case without modifying code,
• illustrate equatorial circulation,
• and to illustrate what appens when we remove Coriolis force, stratification, ...

Structure:

The structure of this test case is done like in CONFIGURATION directory:

For this test case we created a directory TEST_EQ_BOX01 in which we give:

• namelist_cfg
• iodef.xml
• cpp keys

We give ( in CREATE_FILES directory ) python scripts to create these files:

• coordinates.nc
• bathy_r025.nc
• TScst_init_r025.nc
• flx_r025.nc

NEMO code read input files, all is setted in namelist and iodef.xml files, and we give ferret scripts to produce plots.

Outils:

In https://forge.ipsl.jussieu.fr/nemo/browser/branches/2015/dev_r5102_CNRS11_TestCase/NEMOGCM/TOOLS/TEST_CASES directory, there are directories:

CREATE_FILES:

In this directory there are PYTHON scripts to create these files :

bathymetry coordinates with e1,e2 constant coordinates with gphi0=constant (to do not have Coriolis) constant fluxes input file for salinity and temperature constant input file for salinity constant and a temperature profile (given by Levitus profile, averaged for pacific ocean)

CREATE_PLOTS:

In this directory there are FERRET scripts to create plots of :

vertical slides of velocities (between à and 1000 meters) over all bassin (0° - 180°), average between 1°S and 1°N.
vertical slides of temperature (between à and 1000 meters) over all bassin (0° - 180°), average between 1°S and 1°N.

And a list of directories with

• namelist_cfg
• iodef.xml
• cpp keys

TEST_EQ_BOX01_noCoriolis

TEST_EQ_BOX01_cst

TEST_EQ_BOX01

TEST_EQ_BOX01_tprof

TEST_EQ_BOX01_kelvin

These directories are 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.

How to run Equatorial Box :

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)

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

• jpidta
• jpjdta
• jpkdta
• jpiglo
• jpjglo

Liste of experiences

Example 1: Equatorial box without Coriolis

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

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

Example 2: Equatorial box e1, e2 constant values

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

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

Example 3: Equatorial box with e1, e2 variables with latitudes

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

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

Example 4: Equatorial box with stratification and constant winds

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

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

Example 5: Equatorial box stratification and variable wind

5. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with stratification / variable wind (to see 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: Kelv-current-temp-1S1N-all.pdf