Changes between Version 29 and Version 30 of WorkingGroups/TestCases/ImplementationStrategyAndTemplate/EquatorialBox
- Timestamp:
- 2017-04-27T21:31:03+02:00 (7 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
WorkingGroups/TestCases/ImplementationStrategyAndTemplate/EquatorialBox
v29 v30 1 1 [[PageOutline]] 2 2 3 [[TOC( heading=Robustness-TestCases,Robustness-TestCases/ImplementationStrategyAndTemplate/EquatorialBox/*, depth=2)]]3 [[TOC(WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/EquatorialBox/*, depth=2)]] 4 4 5 = '''TEST CASE: EQUATORIAL BOX''' = 6 == Motivations: == 5 = TEST CASE: EQUATORIAL BOX 6 7 == Motivations: 7 8 This equatorial box is done to : 8 9 9 * show how to create in NEMO a test case without modifying code, 10 10 * illustrate basic features for equatorial circulation (impact of Coriolis force, stratification, ...) 11 11 12 == Structure: == 13 The structure of this test case is done like in CONFIGURATION directory: 14 12 == Structure: 13 The structure of this test case is done like in CONFIGURATION directory.[[BR]] 15 14 For this test case we created a directory TEST_EQ_BOX01 in which we give: 16 17 * namelist_cfg 18 * iodef.xml 15 * `namelist_cfg` 16 * `iodef.xml` 19 17 * cpp keys 20 18 21 19 We give ( in CREATE_FILES directory ) '''python scripts''' to create these files: 20 * `coordinates.nc` 21 * `bathy_r025.nc` 22 * `TScst_init_r025.nc` 23 * `flx_r025.nc` 22 24 23 * coordinates.nc 24 * bathy_r025.nc 25 * TScst_init_r025.nc 26 * flx_r025.nc 27 28 NEMO code read input files, all is setted in namelist and iodef.xml files, and we give ferret scripts to produce plots. 25 NEMO code read input files, all is setted in `namelist` and `iodef.xml` files, and we give ferret scripts to produce plots. 29 26 30 27 == Outils: == 31 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:28 In [browser:branches/2015/dev_r5102_CNRS11_TestCase/NEMOGCM/TOOLS/TEST_CASES TEST_CASES] directory, there are directories: 32 29 33 30 === CREATE_FILES: === … … 42 39 * input file for salinity constant and a temperature profile (given by Levitus profile, averaged for pacific ocean) 43 40 44 === CREATE_PLOTS: ===41 === CREATE_PLOTS: 45 42 In this directory there are ''' FERRET scripts ''' to create plots of : 46 43 … … 48 45 * vertical slides of temperature (between à and 1000 meters) over all bassin (0° - 180°), average between 1°S and 1°N. 49 46 50 == And a list of directories with ==51 * namelist_cfg52 * iodef.xml47 == And a list of directories with 48 * `namelist_cfg` 49 * `iodef.xml` 53 50 * cpp keys 54 51 55 === TEST_EQ_BOX01_noCoriolis : ===52 === TEST_EQ_BOX01_noCoriolis : 56 53 Equatorial box without Coriolis 57 54 58 === TEST_EQ_BOX01_cst ===55 === TEST_EQ_BOX01_cst 59 56 Equatorial box with e1= e2 = constant 60 57 61 === TEST_EQ_BOX01 ===58 === TEST_EQ_BOX01 62 59 Equatorial box with e2 variable with latitude 63 60 64 === TEST_EQ_BOX01_tprof ===61 === TEST_EQ_BOX01_tprof 65 62 Equatorial box with stratification and constant winds 66 63 67 === TEST_EQ_BOX01_kelvin ===64 === TEST_EQ_BOX01_kelvin 68 65 Equatorial box with stratification, start form ''' restart''' of TEST_EQ_BOW01_tprof''' and winds=0 (to see kelvin waves)''' ''''' ''''' ' 69 66 … … 89 86 1. run ferret scripts to have gif files and compare with the pdf output files here given for each of 5 experiences. 90 87 91 ''' NOTA BENE''': if you change resolution (in python scripts) you need to change also in namelist_cfg:88 '''noTA BENE''': if you change resolution (in python scripts) you need to change also in namelist_cfg: 92 89 93 90 * '''jpidta''' … … 99 96 = Liste of experiences : = 100 97 = Example 1: Equatorial box without Coriolis = 101 || '''1. TEST_EQ_BOX01_noCoriolis''' ''': Equatorial box (180°,60°) : with NO stratification / constant wind trades / NOCoriolis''' || ||98 || '''1. TEST_EQ_BOX01_noCoriolis''' ''': Equatorial box (180°,60°) : with no stratification / constant wind trades / no Coriolis''' || || 102 99 || || || 103 100 || '''options''' || '''values''' || … … 108 105 || batymetry || 5000m || 109 106 || forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 110 || NOstratification || so=35 / thetao=10 ||107 || no stratification || so=35 / thetao=10 || 111 108 || total time steps || 3600 (=50days) || 112 109 || output frequency || 5 days || 113 114 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 ] 110 output plot: [[PdfImg(gphi0-current-temp-1S1N-all.pdf, width="800")]] 115 111 116 112 = Example 2: Equatorial box e1, e2 constant values = 117 || '''2. TEST_EQ_BOX01_cst''' ''': Equatorial box (180°,60°) : with NOstratification / constant wind trades / e1=e2=constant''' || ||113 || '''2. TEST_EQ_BOX01_cst''' ''': Equatorial box (180°,60°) : with no stratification / constant wind trades / e1=e2=constant''' || || 118 114 || || || 119 115 || '''options''' || '''values''' || … … 124 120 || batymetry || 5000m || 125 121 || forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 126 || NOstratification || so=35 / thetao=10 ||122 || no stratification || so=35 / thetao=10 || 127 123 || total time steps || 3600(=50days) || 128 124 || output frequency || 5 days || 129 130 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] 125 output plot: [[PdfImg(e1e2cst-current-temp-1S1N-all.pdf, width="800")]] 131 126 132 127 = Example 3: Equatorial box with e1, e2 variables with latitudes = 133 || '''3. TEST_EQ_BOX01_cst''' ''': Equatorial box (180°,60°) : with NOstratification / constant wind trades / e1, e2 variable with longitude''' || ||128 || '''3. TEST_EQ_BOX01_cst''' ''': Equatorial box (180°,60°) : with no stratification / constant wind trades / e1, e2 variable with longitude''' || || 134 129 || || || 135 130 || '''options''' || '''values''' || … … 140 135 || batymetry || 5000m || 141 136 || forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 142 || NOstratification || so=35 / thetao=10 ||137 || no stratification || so=35 / thetao=10 || 143 138 || total time steps || 3600(=50days) || 144 139 || output frequency || 5 days || 145 146 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] 140 output plot: [[PdfImg(TScst-current-temp-1S1N-all.pdf, width="800")]] 147 141 148 142 = Example 4: Equatorial box with stratification and constant winds = … … 156 150 || batymetry || 5000m || 157 151 || forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 158 || NOstratification || so=35 / thetao=vertical profile ||152 || no stratification || so=35 / thetao=vertical profile || 159 153 || total time steps || 3600(=50days) || 160 154 || output frequency || 5 days || 161 162 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] 155 output plot: [[PdfImg(tprof-current-temp-1S1N-all.pdf, width="800")]] 163 156 164 157 = Example 5: Equatorial box stratification and variable wind = … … 172 165 || batymetry || 5000m || 173 166 || forcing files || utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0 || 174 || NOstratification || so=35 / thetao=vertical profile ||167 || no stratification || so=35 / thetao=vertical profile || 175 168 || total time steps || 3600 (restart from a 1801 time step, restart file from test 6) with winds = 0 || 176 169 || output frequency || 5 days || 177 178 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] 170 output plot: [[PdfImg(Kelv-current-temp-1S1N-all.pdf, width="800")]]