Changes between Version 1 and Version 2 of Working Groups/Test Cases/Implementation Strategy And Template/Equatorial Box


Ignore:
Timestamp:
2015-06-18T11:50:09+02:00 (5 years ago)
Author:
flavoni
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • Working Groups/Test Cases/Implementation Strategy And Template/Equatorial Box

    v1 v2  
    33[[TOC(heading=Robustness-TestCases,Robustness-TestCases/ImplementationStrategyAndTemplate/EquatorialBox/*, depth=2)]] 
    44 
    5 = Title = 
     5= '''''' = 
     6= '''TEST CASE: EQUATORIAL BOX''' = 
     7=  = #adetailedTESTCASE:EQUATORIALBOX 
     8in [https://forge.ipsl.jussieu.fr/nemo/browser/branches/2015/dev_r5102_CNRS11_TestCase/NEMOGCM/TOOLS/TEST_CASES/TEST_EQ_BOX01 ​https://forge.ipsl.jussieu.fr/nemo/browser/branches/2015/dev_r5102_CNRS11_TestCase/NEMOGCM/TOOLS/TEST_CASES/TEST_EQ_BOX01] 
     9 
     10and some examples to create a test case of this type : 
     11 
     12  Equatorial BOX, grid of '''0.25''' ° resolution, between latitude: -30°,30°, and longitude 0°,180°, 
     13 
     14  with zps coordinates, 
     15 
     16  with ubs schema on dynamics and on tracers, with constant TS initial state, with flux only wind (u composant), with flat bathymetry. 
     17 
     18If you want to create a test case what you have to do? 
     19 
     20 1. go into (branches/2015/dev_r5102_CNRS11_TestCase/)NEMOGCM/TOOLS/TEST_CASES/CREATE_FILES 
     21 1. read README file 
     22 1. run python make_coordinates.py to generate coordinates_r025.nc 
     23 1. link coordinates.nc to coordinates_r025.nc 
     24 1. run python make_bathy.py to generate bathy_meter.nc -> bathy_r025.nc 
     25 1. link bathy_meter.nc to bathy_r025.nc 
     26 1. run python make_TScst_init.py to generate TScst_init_r025.nc (you do not need link, filename is written in namelist_cfg file) 
     27 1. run python make_FLX.py to generate flx_r025.nc (you do not need link, filename is written in namelist_cfg file) 
     28 
     29You choose in namelist_cfg cn_exp name and nitend. 
     30 
     31[[BR]]'''NOTA BENE''': if you change resolution (in python scripts) you need to change also in namelist_cfg: 
     32 
     33  '''jpidta, jpjdta and jpkdta''' '''jpiglo and jpjglo''' ''''''''''' 
     34 
     35= Example 1: Equatorial box without Coriolis = 
     36[[BR]] 
     37|| '''1. TEST_EQ_BOX01_noCoriolis''' ''': Equatorial box (180°,60°) : with NO stratification / constant wind trades / NO Coriolis''' || || 
     38|| || || 
     39|| '''options''' || '''values''' || 
     40|| resolution || 0.25° (1/4) || 
     41|| longitude || lon_T=[ 0,180 ] || 
     42|| latitude || lat_T=[-31,31] (to have one point at 0° Equator) || 
     43|| vertical levels || 31 || 
     44|| batymetry || 5000m || 
     45|| forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 
     46|| NO stratification || so=35 / thetao=10 || 
     47|| total time steps || 3600 (=50days) || 
     48|| output frequency || 5 days || 
     49 
     50= Example 2: Equatorial box e1, e2 constant values = 
     51 
     52 
     53|| '''2. TEST_EQ_BOX01_cst''' ''': Equatorial box (180°,60°) : with NO stratification / constant wind trades / e1=e2=constant''' || || 
     54|| || || 
     55|| '''options''' || '''values''' || 
     56|| resolution || 0.25° (1/4) || 
     57|| longitude || lon_T=[ 0,180 ] || 
     58|| latitude || lat_T=[-31,31] (to have one point at 0° Equator) || 
     59|| vertical levels || 31 || 
     60|| batymetry || 5000m || 
     61|| forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 
     62|| NO stratification || so=35 / thetao=10 || 
     63|| total time steps || 3600(=50days) || 
     64|| output frequency || 5 days || 
     65 
     66=  = 
     67= Example 3: Equatorial box with e1, e2 variables with latitudes = 
     68 
     69 
     70|| '''3. TEST_EQ_BOX01_cst''' ''': Equatorial box (180°,60°) : with NO stratification / constant wind trades / e1, e2 variable with longitude''' || || 
     71|| || || 
     72|| '''options''' || '''values''' || 
     73|| resolution || 0.25° (1/4) || 
     74|| longitude || lon_T=[ 0,180 ] || 
     75|| latitude || lat_T=[-31,31] (to have one point at 0° Equator) || 
     76|| vertical levels || 31 || 
     77|| batymetry || 5000m || 
     78|| forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 
     79|| NO stratification || so=35 / thetao=10 || 
     80|| total time steps || 3600(=50days) || 
     81|| output frequency || 5 days || 
     82 
     83= Example 4: Equatorial box with stratification and constant winds = 
     84 
     85 
     86|| '''4. TEST_EQ_BOX01_cst''' ''': Equatorial box (180°,60°) : with stratification / constant wind trades''' || || 
     87|| || || 
     88|| '''options''' || '''values''' || 
     89|| resolution || 0.25° (1/4) || 
     90|| longitude || lon_T=[ 0,180 ] || 
     91|| latitude || lat_T=[-31,31] (to have one point at 0° Equator) || 
     92|| vertical levels || 31 || 
     93|| batymetry || 5000m || 
     94|| forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 
     95|| NO stratification || so=35 / thetao=10 || 
     96|| total time steps || 3600(=50days) || 
     97|| output frequency || 5 days || 
     98 
     99=  = 
     100= Example 5: Equatorial box stratification and constant wind profile = 
     101 
     102 
     103[[BR]] 
     104|| '''5. TEST_EQ_BOX01_cst''' ''': Equatorial box (180°,60°) : with stratification / constant winds (profile)''' || || 
     105|| || || 
     106|| '''options''' || '''values''' || 
     107|| resolution || 0.25° (1/4) || 
     108|| longitude || lon_T=[ 0,180 ] || 
     109|| latitude || lat_T=[-31,31] (to have one point at 0° Equator) || 
     110|| vertical levels || 31 || 
     111|| batymetry || 5000m || 
     112|| forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 
     113|| NO stratification || so=35 / thetao=10 || 
     114|| total time steps || 3600(=50days) || 
     115|| output frequency || 5 days || 
     116 
     117 
     118 
     119= = Example 6: Equatorial box stratification and variable wind  = 
     120 = 
     121 
     122 
     123|| '''6. TEST_EQ_BOX01_cst''' ''': Equatorial box (180°,60°) : with stratification / variable wind (to see Kelvin waves)''' || || 
     124|| || || 
     125|| '''options''' || '''values''' || 
     126|| resolution || 0.25° (1/4) || 
     127|| longitude || lon_T=[ 0,180 ] || 
     128|| latitude || lat_T=[-31,31] (to have one point at 0° Equator) || 
     129|| vertical levels || 31 || 
     130|| batymetry || 5000m || 
     131|| forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 
     132|| NO stratification || so=35 / thetao=10 || 
     133|| total time steps || 3600 (restart from a 1801 time step, restart file from test 6) with winds = 0 || 
     134|| output frequency || 5 days || 
     135 
     136output plot:[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/current-temp-1S1N-all.pdf ​current-temp-1S1N-all.pdf]