[[PageOutline]] [[TOC(heading=Robustness-TestCases,Robustness-TestCases/ImplementationStrategyAndTemplate/EquatorialBox/*, depth=2)]] = '''' = = '''TEST CASE: EQUATORIAL BOX''' = = = #adetailedTESTCASE:EQUATORIALBOX in [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] 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 1. read README file 1. run python make_coordinates.py to generate coordinates_r025.nc 1. link coordinates.nc to coordinates_r025.nc 1. run python make_bathy.py to generate bathy_meter.nc -> bathy_r025.nc 1. link bathy_meter.nc to bathy_r025.nc 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) 1. 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. [[BR]]'''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 || = = = 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 || = = = = = 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 || = = = 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 || = = = = = Example 5: Equatorial box stratification and constant wind profile = || '''5. 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 || = = =  Example 6: Equatorial box stratification and variable wind = || '''6. 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:[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/current-temp-1S1N-all.pdf ​current-temp-1S1N-all.pdf]