[[PageOutline]] [[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 ​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. [[BR]] vertical slides of temperature (between à and 1000 meters) over all bassin (0° - 180°), average between 1°S and 1°N. [[BR]] 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 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) 1. and run NEMO like usually (./opa in mpi if you want) 1. run ferret scripts to have gif files and compare with the output files here given for each experience. '''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:[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/EquatorialBox/gphi0-current-temp-1S1N-all.pdf 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: [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/EquatorialBox/e1e2cst-current-temp-1S1N-all.pdf 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:[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/EquatorialBox/TScst-current-temp-1S1N-all.pdf 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: [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/EquatorialBox/tprof-current-temp-1S1N-all.pdf 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: [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/EquatorialBox/Kelv-current-temp-1S1N-all.pdf Kelv-current-temp-1S1N-all.pdf]