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