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/EquatorialBox – NEMO
wiki:WorkingGroups/TestCases/ImplementationStrategyAndTemplate/EquatorialBox

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

--

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, ...

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 files:

coordinates.nc

bathy_r025.nc TScst_init_r025.nc flx_r025.nc

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

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

This test case wants to show:

Outils:

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

There are these 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.

TEST_EQ_BOX01:

TEST_EQ_BOX01_cst:

TEST_EQ_BOX01_noCoriolis:

TEST_EQ_BOX01_tprof:

TEST_EQ_BOX01_kelvin:

there 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.

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)

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

jpidta, jpjdta and jpkdta jpiglo and jpjglo '

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

Attachments (5)

Download all attachments as: .zip