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
- input files:
- coordinates.nc
- bathy_r025.nc
- initial state:
- TScst_init_r025.nc
- TS_init_tprof_r025.nc
- forcing file:
- flx_r025.nc
a detailed TEST CASE: EQUATORIAL BOX
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?
- go into (branches/2015/dev_r5102_CNRS11_TestCase/)NEMOGCM/TOOLS/TEST_CASES/CREATE_FILES
- read README file
- run python make_coordinates.py to generate coordinates_r025.nc
- link coordinates.nc to coordinates_r025.nc
- run python make_bathy.py to generate bathy_meter.nc -> bathy_r025.nc
- link bathy_meter.nc to bathy_r025.nc
- run python make_TScst_init.py to generate TScst_init_r025.nc (you do not need link, filename is written in namelist_cfg file)
- 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)
- current-temp-1S1N-all.pdf (290.6 KB) - added by flavoni 9 years ago.
Download all attachments as: .zip