1 | //#define FILTER |
---|
2 | //#define SOLSORYAO |
---|
3 | //#define OPTIMORDER |
---|
4 | //#============================================================================== |
---|
5 | //# fichier de DESCRIPTION sinobad.d |
---|
6 | //#Mohamed Berrada 26/02/2009 |
---|
7 | //#modification (#define) Julien Brajard 06/2011 |
---|
8 | //#modification (filtre) Mohamed Berrada 19/07/2011 |
---|
9 | //#============================================================================== |
---|
10 | //#DEFVAL ---------------------------------------------DEFVAL |
---|
11 | //#definition de valeur |
---|
12 | defval NPCA 15 |
---|
13 | defval NX 32 |
---|
14 | defval NY 22 |
---|
15 | defval NZ 31 |
---|
16 | defval NMAX 800 // ! maximum of iterations for the SOR solver |
---|
17 | defval TA 50 |__DATE__ | number of time steps |
---|
18 | defval TU 2 | start time step |
---|
19 | #ifdef FILTER |
---|
20 | defval OFTL 40 |
---|
21 | defval OFTZ 40 |
---|
22 | defval OFT 80 |
---|
23 | #endif |
---|
24 | #ifndef FILTER |
---|
25 | defval OFT 0 |
---|
26 | #endif |
---|
27 | //#exec disp_valdef |
---|
28 | #ifdef FILTER |
---|
29 | defval K_FILTER |
---|
30 | #endif |
---|
31 | #ifdef SOLSORYAO |
---|
32 | defval K_SOLSORYAO |
---|
33 | #endif |
---|
34 | #ifdef OPTIMORDER |
---|
35 | defval K_OPTIMORDER |
---|
36 | #endif |
---|
37 | |
---|
38 | //#HAT_NAME -----------------------------------------HAT_NAME |
---|
39 | hat_name sinobad |
---|
40 | |
---|
41 | //#OPTION ----------------------------------------------OPTION |
---|
42 | |
---|
43 | option O_M1QN3 | M2QN1 |
---|
44 | //#option O_VARINCR |
---|
45 | option O_GRADTEST |
---|
46 | option O_DBG_NANF |
---|
47 | //#option O_DBG_TING |
---|
48 | option O_EXTOBJ "../object/ncutil.o" |
---|
49 | option O_REAL double |
---|
50 | option o_parallel |
---|
51 | |
---|
52 | //#TRAJECTOIRES (DIM) -----------------------------------TRAJ |
---|
53 | #ifdef FILTER |
---|
54 | traj Tcstf M 0 1 |
---|
55 | traj Tlfil M 1 0 1 OFTL |
---|
56 | traj Tzfil M 1 OFTL 1 OFTZ |
---|
57 | #endif |
---|
58 | traj Tcst M 0 OFT 1 1 |
---|
59 | traj Tsbd M TU OFT 1 TA |
---|
60 | exec disp_traj |
---|
61 | |
---|
62 | #ifdef FILTER |
---|
63 | space S3df M NX NY NZ Tcstf |
---|
64 | space S3dtlf M NX NY NZ Tlfil |
---|
65 | space S3dtzf M NX NY NZ Tzfil |
---|
66 | #endif |
---|
67 | |
---|
68 | space S0d M 1 Tcst |
---|
69 | //#space S3d M NX NY NZ Tcst |
---|
70 | space S3d M NX NY NZ Tcst |
---|
71 | space S3dt M NX NY NZ Tsbd |
---|
72 | #ifdef SOLSORYAO |
---|
73 | space S3dt_b M NX NY NZ Tsbd |
---|
74 | space Ssolsor M NX NY NMAX Tsbd |
---|
75 | #endif |
---|
76 | |
---|
77 | exec disp_space |
---|
78 | |
---|
79 | //####################################################################################### |
---|
80 | //#---------------------------------------- |
---|
81 | //# choisir le vecteur du contrÃŽle (target) |
---|
82 | //#---------------------------------------- |
---|
83 | //#modul pca_ta space S0d noward output NPCA target |
---|
84 | |
---|
85 | modul sshn_c space S3d inter YA1 YA2 noward output 1 |target |
---|
86 | |
---|
87 | //#modul ta_c space S3d input array NPCA output 1 |
---|
88 | #ifndef FILTER |
---|
89 | modul ta_c space S3d noward output 1 target |
---|
90 | #endif |
---|
91 | |
---|
92 | modul sa_c space S3d noward output 1 |target |
---|
93 | modul ua_c space S3d noward output 1 |target |
---|
94 | modul va_c space S3d noward output 1 |target |
---|
95 | |
---|
96 | #ifdef FILTER |
---|
97 | modul nu space S3df noward output 1 target |
---|
98 | modul dta_lfil space S3dtlf input 6 output 1 tempo |
---|
99 | modul ztuv_lfil space S3dtlf input 3 output 2 tempo |
---|
100 | modul dta_zfil space S3dtzf input 4 output 1 tempo |
---|
101 | modul zwy_zfil space S3dtzf input 2 output 1 tempo |
---|
102 | modul ta_c space S3d input 1 output 1 |
---|
103 | #endif |
---|
104 | |
---|
105 | //#------------- |
---|
106 | //# Les modules dans l'ordre du code NEMO/OPA |
---|
107 | //#------------- |
---|
108 | //#boucle sur kt (NEMO/OPA) |
---|
109 | //#step(kt) ---> ---> |
---|
110 | //##START_sbc (debut routine sbc) |
---|
111 | modul sbcgyre_qsr_qns space S3dt inter YA1 YA2 input 1 output 2 tempo |
---|
112 | modul sbcgyre_emp1 space S3dt inter YA1 YA2 input 0 output 1 tempo |
---|
113 | modul sbcgyre_zsumemp_zsurf space S3dt inter YA1 YA2 input 5 output 2 tempo |
---|
114 | modul sbcgyre_emp space S3dt inter YA1 YA2 input 3 output 1 tempo |
---|
115 | modul sbcgyre_emps space S3dt inter YA1 YA2 input 1 output 1 tempo |
---|
116 | modul sbcgyre_utau_vtau space S3dt inter YA1 YA2 input 0 output 2 tempo |
---|
117 | modul sbcgyre_wndm space S3dt inter YA1 YA2 input 4 output 1 tempo |
---|
118 | //##END_sbc (fin routine sbc) |
---|
119 | |
---|
120 | //##START_bn2 (debut routine bn2) (neos=2) |
---|
121 | modul bn2_rn2 space S3dt input 4 output 1 tempo |
---|
122 | //##END_bn2(fin routine bn2) |
---|
123 | |
---|
124 | //##START_zdf (n_evdm=1,nbotfr=2) |
---|
125 | //#mélange vertical |
---|
126 | modul zdf_avt_avmu_avmv space S3dt input 19 output 3 tempo |
---|
127 | modul zdfmxl_hmld space S3dt inter YA1 YA2 input array NZ output 1 tempo |
---|
128 | modul zdfmxl_nmln_hmlp_hmlpt space S3dt inter YA1 YA2 input array NZ output 3 tempo |
---|
129 | //##END_zdf |
---|
130 | |
---|
131 | //##START_trasbc |
---|
132 | //#température et salinté de surface |
---|
133 | modul ta_sa_trasbc space S3dt inter YA1 YA2 input 3 output 2 tempo |
---|
134 | //#//#END_trasbc |
---|
135 | |
---|
136 | //##START_traqsr |
---|
137 | modul ta_traqsr space S3dt input 2 output 1 tempo |
---|
138 | //##END_traqsr |
---|
139 | |
---|
140 | //##START_traadv_cen2 |
---|
141 | //#advection (cen2 second order center scheme) |
---|
142 | modul zcent_hor_traadv_cen2 space S3dt input 5 output 2 tempo |
---|
143 | modul zcens_hor_traadv_cen2 space S3dt clonol zcent_hor_traadv_cen2 |
---|
144 | modul ta_hor_traadv_cen2 space S3dt input 5 output 1 tempo |
---|
145 | modul sa_hor_traadv_cen2 space S3dt input 5 output 1 tempo |
---|
146 | modul zcent_ver_traadv_cen2 space S3dt input 3 output 1 tempo |
---|
147 | modul zcens_ver_traadv_cen2 space S3dt clonol zcent_ver_traadv_cen2 |
---|
148 | modul ta_ver_traadv_cen2 space S3dt input 3 output 1 tempo |
---|
149 | modul sa_ver_traadv_cen2 space S3dt clonol ta_ver_traadv_cen2 |
---|
150 | |
---|
151 | //##END_traadv_cen2 |
---|
152 | |
---|
153 | //##START_traldf_lap |
---|
154 | //#diffusion latérale (laplacian scheme) |
---|
155 | modul ztuv_traldf_lap space S3dt input 3 output 2 tempo |
---|
156 | modul zsuv_traldf_lap space S3dt clonol ztuv_traldf_lap |
---|
157 | modul ta_traldf_lap space S3dt input 5 output 1 tempo |
---|
158 | modul sa_traldf_lap space S3dt clonol ta_traldf_lap |
---|
159 | //##END_traldf_lap |
---|
160 | |
---|
161 | //##START_trazdf_imp |
---|
162 | //#diffusion verticale (schema implicit) |
---|
163 | modul zwisd_trazdf_imp space S3dt input 2 output 3 tempo |
---|
164 | modul zrhs_trazdf_imp space S3dt input 2 output 1 tempo |
---|
165 | modul zrhs_sa_trazdf_imp space S3dt clonol zrhs_trazdf_imp |
---|
166 | modul zwt_trazdf_imp space S3dt input 4 output 1 tempo |
---|
167 | modul ta_de_trazdf_imp space S3dt clonol zwt_trazdf_imp |
---|
168 | modul ta space S3dt input 5 output 1 tempo |
---|
169 | modul sa_de_trazdf_imp space S3dt clonol ta_de_trazdf_imp |
---|
170 | modul sa space S3dt input 5 output 1 tempo |
---|
171 | //##END_trazdf_imp |
---|
172 | |
---|
173 | //##START_tra_nxt |
---|
174 | //#tracer next : filtre asselin |
---|
175 | modul tb space S3dt input 4 output 1 tempo cout |
---|
176 | modul sb space S3dt input 4 output 1 tempo cout |
---|
177 | //##END_tra_nxt |
---|
178 | |
---|
179 | //##START_eos |
---|
180 | //#equation of state, density |
---|
181 | modul rhd space S3dt input 2 output 1 tempo |
---|
182 | modul rhop space S3dt input 1 output 1 tempo |
---|
183 | //##END_eos |
---|
184 | |
---|
185 | //##START_dyn_adv_cen2 |
---|
186 | //#advention |
---|
187 | modul hmf_TF_dynadv_cen2 space S3dt input 6 output 4 tempo |
---|
188 | modul zua_dhmf_dynadv_cen2 space S3dt input 4 output 1 tempo |
---|
189 | modul zva_dhmf_dynadv_cen2 space S3dt input 4 output 1 tempo |
---|
190 | modul zfuv_UW_dynadv_cen2 space S3dt input 7 output 2 tempo |
---|
191 | modul ua_dynadv_cen2 space S3dt input 3 output 1 tempo |
---|
192 | modul va_dynadv_cen2 space S3dt input 3 output 1 tempo |
---|
193 | //##END_dyn_adv_cen2 |
---|
194 | |
---|
195 | //##START_vor |
---|
196 | //#vorticity enstrophy conservative scheme |
---|
197 | modul zwxyz_vor_ene space S3dt input 4 output 3 tempo |
---|
198 | modul ua_vor_ene space S3dt input 7 output 1 tempo |
---|
199 | modul va_vor_ene space S3dt input 7 output 1 tempo |
---|
200 | //##END_vor |
---|
201 | |
---|
202 | //##START_dyn_ldf |
---|
203 | //#lateral diffusion |
---|
204 | modul ua_dyn_ldf_lap space S3dt input 5 output 1 tempo |
---|
205 | modul va_dyn_ldf_lap space S3dt input 5 output 1 tempo |
---|
206 | //##END_dyn_ldf |
---|
207 | |
---|
208 | //##START_dyn_hpg |
---|
209 | //#hydrostratic pressure gradient |
---|
210 | modul zhpij_dyn_hpg_zco space S3dt input 8 output 2 tempo |
---|
211 | modul ua_dyn_hpg_zco space S3dt input 2 output 1 tempo |
---|
212 | modul va_dyn_hpg_zco space S3dt clonol ua_dyn_hpg_zco |
---|
213 | //##END_dyn_hpg |
---|
214 | |
---|
215 | //##START_dyn_zdf |
---|
216 | //#vertical diffusion |
---|
217 | modul zwisd_ua_dynzdf_imp space S3dt input 2 output 3 tempo |
---|
218 | modul zwisd_va_dynzdf_imp space S3dt input 2 output 3 tempo |
---|
219 | modul zrhs_ua_dynzdf_imp space S3dt input 3 output 1 tempo |
---|
220 | modul zrhs_va_dynzdf_imp space S3dt input 3 output 1 tempo |
---|
221 | modul zwt_ua_dynzdf_imp space S3dt clonol zwt_trazdf_imp |
---|
222 | modul zwt_va_dynzdf_imp space S3dt clonol zwt_trazdf_imp |
---|
223 | modul ua_de_dynzdf_imp space S3dt clonol ta_de_trazdf_imp |
---|
224 | modul ua_dynzdf_imp space S3dt input 4 output 1 tempo |
---|
225 | modul va_de_dynzdf_imp space S3dt clonol ta_de_trazdf_imp |
---|
226 | modul va_dynzdf_imp space S3dt clonol ua_dynzdf_imp |
---|
227 | modul ua_norm_dynzdf space S3dt input 2 output 1 tempo |
---|
228 | modul va_norm_dynzdf space S3dt clonol ua_norm_dynzdf |
---|
229 | //##END_dyn_zdf |
---|
230 | |
---|
231 | //##START_dyn_spg |
---|
232 | //#surface presure gradient |
---|
233 | modul spguv_now_dynspg space S3dt inter YA1 YA2 input 3 output 2 tempo |
---|
234 | modul ua_1_dynspg_flt space S3dt input 3 output 1 tempo |
---|
235 | modul va_1_dynspg_flt space S3dt input 3 output 1 tempo |
---|
236 | modul sum_spgu_dynspg_flt space S3dt inter YA1 YA2 input array NZ output 1 tempo |
---|
237 | modul sum_spgv_dynspg_flt space S3dt inter YA1 YA2 input array NZ output 1 tempo |
---|
238 | modul gcb_dynspg_flt space S3dt inter YA1 YA2 input 4 output 1 tempo |
---|
239 | modul gcx_dynspg_flt space S3dt inter YA1 YA2 input 2 output 1 tempo |
---|
240 | #ifndef SOLSORYAO |
---|
241 | modul solsor_dynspg_flt space S3dt inter YA1 YA2 input array 19 output 1 tempo hidjac |
---|
242 | modul gcx2 space S3dt inter YA1 YA2 input 1 output 1 tempo hidjac |
---|
243 | #endif |
---|
244 | |
---|
245 | #ifdef SOLSORYAO |
---|
246 | modul solsor_gx1_x_r space Ssolsor input 8 output 2 tempo |
---|
247 | modul solsor_gx2_x_r space Ssolsor input 9 output 2 tempo |
---|
248 | modul solsor_flag space Ssolsor input 4 output 2 tempo |
---|
249 | #define S3dt S3dt_b |
---|
250 | modul gcx2 space S3dt inter YA1 YA2 input 1 output 1 tempo |
---|
251 | #endif |
---|
252 | |
---|
253 | |
---|
254 | |
---|
255 | modul spguv_dynspg_flt space S3dt inter YA1 YA2 input 3 output 2 tempo |
---|
256 | modul ua space S3dt input 3 output 1 tempo |
---|
257 | modul va space S3dt input 3 output 1 tempo |
---|
258 | modul sshn space S3dt inter YA1 YA2 input 4 output 1 tempo |
---|
259 | modul sshb space S3dt inter YA1 YA2 input 4 output 1 tempo cout |
---|
260 | //##END_dyn_spg |
---|
261 | |
---|
262 | //##START_dyn_nxt |
---|
263 | //#nxt : filtre d'asselin |
---|
264 | modul ub space S3dt input 4 output 1 tempo cout |
---|
265 | modul vb space S3dt input 4 output 1 tempo cout |
---|
266 | //##END_dyn_nxt |
---|
267 | |
---|
268 | //##START_divcur |
---|
269 | //#calcul de div et de rotationnel |
---|
270 | modul hdivn space S3dt input 4 output 1 tempo |
---|
271 | modul rotn space S3dt input 4 output 1 tempo |
---|
272 | //##END_divcur |
---|
273 | |
---|
274 | //##START_wzv |
---|
275 | //#vitesse vertical |
---|
276 | modul wa space S3dt input 2 output 1 tempo |
---|
277 | //##END_wzv |
---|
278 | |
---|
279 | #ifdef SOLSORYAO |
---|
280 | #undef S3dt |
---|
281 | #endif |
---|
282 | |
---|
283 | exec disp_modul |
---|
284 | //############################################################################################# |
---|
285 | |
---|
286 | #ifdef FILTER |
---|
287 | //#---------->dta_lfil |
---|
288 | ctin dta_lfil 1 from nu 1 i j k |
---|
289 | ctin dta_lfil 2 from dta_lfil 1 i j k-1 t-1 |
---|
290 | ctin dta_lfil 3 from ztuv_lfil 1 i-1 j k t |
---|
291 | ctin dta_lfil 4 from ztuv_lfil 1 i j k t |
---|
292 | ctin dta_lfil 5 from ztuv_lfil 2 i j-1 k t |
---|
293 | ctin dta_lfil 6 from ztuv_lfil 2 i j k t |
---|
294 | |
---|
295 | //#---------->ztuv_lfil |
---|
296 | ctin ztuv_lfil 1 from dta_lfil 1 i+1 j k t-1 |
---|
297 | ctin ztuv_lfil 2 from dta_lfil 1 i j k t-1 |
---|
298 | ctin ztuv_lfil 3 from dta_lfil 1 i j+1 k t-1 |
---|
299 | |
---|
300 | //#---------->dta_zfil |
---|
301 | ctin dta_zfil 1 from dta_lfil 1 i j k t |
---|
302 | ctin dta_zfil 2 from dta_zfil 1 i j k t-1 |
---|
303 | ctin dta_zfil 3 from zwy_zfil 1 i j k+1 t |
---|
304 | ctin dta_zfil 4 from zwy_zfil 1 i j k t |
---|
305 | |
---|
306 | //#---------->zwy_zfil |
---|
307 | ctin zwy_zfil 1 from dta_zfil 1 i j k-1 t-1 |
---|
308 | ctin zwy_zfil 2 from dta_zfil 1 i j k t-1 |
---|
309 | |
---|
310 | //#---------->ta_c |
---|
311 | ctin ta_c 1 from dta_zfil 1 i j k-1 t |
---|
312 | #endif |
---|
313 | //#ctin ta_c 1..NPCA from pca_ta 1..NPCA 1 |
---|
314 | |
---|
315 | //exec disp_ct_in |
---|
316 | exec disp_modul |
---|
317 | |
---|
318 | |
---|
319 | //#---------->sbcgyre |
---|
320 | ctin sbcgyre_qsr_qns 1 from tb 1 i j 1 t-1 |
---|
321 | |
---|
322 | ctin sbcgyre_zsumemp_zsurf 1 from sbcgyre_zsumemp_zsurf 1 i-1 j t |
---|
323 | ctin sbcgyre_zsumemp_zsurf 2 from sbcgyre_zsumemp_zsurf 1 NX j-1 t |
---|
324 | ctin sbcgyre_zsumemp_zsurf 3 from sbcgyre_zsumemp_zsurf 2 i-1 j t |
---|
325 | ctin sbcgyre_zsumemp_zsurf 4 from sbcgyre_zsumemp_zsurf 2 NX j-1 t |
---|
326 | ctin sbcgyre_zsumemp_zsurf 5 from sbcgyre_emp1 1 i j t |
---|
327 | |
---|
328 | ctin sbcgyre_emp 1 from sbcgyre_emp1 1 i j t |
---|
329 | ctin sbcgyre_emp 2 from sbcgyre_zsumemp_zsurf 1 NX NY t |
---|
330 | ctin sbcgyre_emp 3 from sbcgyre_zsumemp_zsurf 2 NX NY t |
---|
331 | |
---|
332 | ctin sbcgyre_emps 1 from sbcgyre_emp 1 i j t |
---|
333 | |
---|
334 | ctin sbcgyre_wndm 1 from sbcgyre_utau_vtau 1 i-1 j t |
---|
335 | ctin sbcgyre_wndm 2 from sbcgyre_utau_vtau 2 i j-1 t |
---|
336 | ctin sbcgyre_wndm 3 from sbcgyre_utau_vtau 1 i j t |
---|
337 | ctin sbcgyre_wndm 4 from sbcgyre_utau_vtau 2 i j t |
---|
338 | |
---|
339 | //#---------->bn2 |
---|
340 | ctin bn2_rn2 1 from tb 1 i j k t-1 |
---|
341 | ctin bn2_rn2 2 from tb 1 i j k-1 t-1 |
---|
342 | ctin bn2_rn2 3 from sb 1 i j k t-1 |
---|
343 | ctin bn2_rn2 4 from sb 1 i j k-1 t-1 |
---|
344 | |
---|
345 | //#---------->zdf |
---|
346 | ctin zdf_avt_avmu_avmv 1 from bn2_rn2 1 i j k t |
---|
347 | ctin zdf_avt_avmu_avmv 2 from bn2_rn2 1 i+1 j k t |
---|
348 | ctin zdf_avt_avmu_avmv 3 from bn2_rn2 1 i j+1 k t |
---|
349 | |
---|
350 | ctin zdf_avt_avmu_avmv 4 from ua 1 i-1 j 30 t-1 |
---|
351 | ctin zdf_avt_avmu_avmv 5 from ua 1 i-1 j+1 30 t-1 |
---|
352 | ctin zdf_avt_avmu_avmv 6 from ua 1 i j 30 t-1 |
---|
353 | ctin zdf_avt_avmu_avmv 7 from ua 1 i j+1 30 t-1 |
---|
354 | |
---|
355 | ctin zdf_avt_avmu_avmv 8 from ua 1 i-1 j 2 t-1 |
---|
356 | ctin zdf_avt_avmu_avmv 9 from ua 1 i-1 j+1 2 t-1 |
---|
357 | ctin zdf_avt_avmu_avmv 10 from ua 1 i j 2 t-1 |
---|
358 | ctin zdf_avt_avmu_avmv 11 from ua 1 i j+1 2 t-1 |
---|
359 | |
---|
360 | ctin zdf_avt_avmu_avmv 12 from va 1 i j-1 30 t-1 |
---|
361 | ctin zdf_avt_avmu_avmv 13 from va 1 i j 30 t-1 |
---|
362 | ctin zdf_avt_avmu_avmv 14 from va 1 i+1 j-1 30 t-1 |
---|
363 | ctin zdf_avt_avmu_avmv 15 from va 1 i+1 j 30 t-1 |
---|
364 | |
---|
365 | ctin zdf_avt_avmu_avmv 16 from va 1 i j-1 2 t-1 |
---|
366 | ctin zdf_avt_avmu_avmv 17 from va 1 i j 2 t-1 |
---|
367 | ctin zdf_avt_avmu_avmv 18 from va 1 i+1 j-1 2 t-1 |
---|
368 | ctin zdf_avt_avmu_avmv 19 from va 1 i+1 j 2 t-1 |
---|
369 | |
---|
370 | ctin zdfmxl_hmld 1..NZ from zdf_avt_avmu_avmv 1 i j 1..NZ t |
---|
371 | ctin zdfmxl_nmln_hmlp_hmlpt 1..NZ from rhop 1 i j 1..NZ t-1 |
---|
372 | |
---|
373 | //#---------->trasbc |
---|
374 | ctin ta_sa_trasbc 1 from sbcgyre_qsr_qns 2 i j t |
---|
375 | ctin ta_sa_trasbc 2 from sbcgyre_emps 1 i j t |
---|
376 | ctin ta_sa_trasbc 3 from sa 1 i j 1 t-1 |
---|
377 | |
---|
378 | //#---------->traqsr |
---|
379 | ctin ta_traqsr 1 from sbcgyre_qsr_qns 1 i j t |
---|
380 | ctin ta_traqsr 2 from ta_sa_trasbc 1 i j t |
---|
381 | |
---|
382 | //#---------->traadv_cen2 |
---|
383 | ctin zcent_hor_traadv_cen2 1 from ua 1 i j k t-1 |
---|
384 | ctin zcent_hor_traadv_cen2 2 from va 1 i j k t-1 |
---|
385 | ctin zcent_hor_traadv_cen2 3 from ta 1 i j k t-1 |
---|
386 | ctin zcent_hor_traadv_cen2 4 from ta 1 i j+1 k t-1 |
---|
387 | ctin zcent_hor_traadv_cen2 5 from ta 1 i+1 j k t-1 |
---|
388 | |
---|
389 | ctin zcens_hor_traadv_cen2 1 from ua 1 i j k t-1 |
---|
390 | ctin zcens_hor_traadv_cen2 2 from va 1 i j k t-1 |
---|
391 | ctin zcens_hor_traadv_cen2 3 from sa 1 i j k t-1 |
---|
392 | ctin zcens_hor_traadv_cen2 4 from sa 1 i j+1 k t-1 |
---|
393 | ctin zcens_hor_traadv_cen2 5 from sa 1 i+1 j k t-1 |
---|
394 | |
---|
395 | ctin ta_hor_traadv_cen2 1 from ta_traqsr 1 i j k t |
---|
396 | ctin ta_hor_traadv_cen2 2 from zcent_hor_traadv_cen2 1 i-1 j k t |
---|
397 | ctin ta_hor_traadv_cen2 3 from zcent_hor_traadv_cen2 1 i j k t |
---|
398 | ctin ta_hor_traadv_cen2 4 from zcent_hor_traadv_cen2 2 i j-1 k t |
---|
399 | ctin ta_hor_traadv_cen2 5 from zcent_hor_traadv_cen2 2 i j k t |
---|
400 | |
---|
401 | ctin sa_hor_traadv_cen2 1 from ta_sa_trasbc 2 i j t |
---|
402 | ctin sa_hor_traadv_cen2 2 from zcens_hor_traadv_cen2 1 i-1 j k t |
---|
403 | ctin sa_hor_traadv_cen2 3 from zcens_hor_traadv_cen2 1 i j k t |
---|
404 | ctin sa_hor_traadv_cen2 4 from zcens_hor_traadv_cen2 2 i j-1 k t |
---|
405 | ctin sa_hor_traadv_cen2 5 from zcens_hor_traadv_cen2 2 i j k t |
---|
406 | |
---|
407 | ctin zcent_ver_traadv_cen2 1 from wa 1 i j k t-1 |
---|
408 | ctin zcent_ver_traadv_cen2 2 from ta 1 i j k t-1 |
---|
409 | ctin zcent_ver_traadv_cen2 3 from ta 1 i j k-1 t-1 |
---|
410 | |
---|
411 | ctin zcens_ver_traadv_cen2 1 from wa 1 i j k t-1 |
---|
412 | ctin zcens_ver_traadv_cen2 2 from sa 1 i j k t-1 |
---|
413 | ctin zcens_ver_traadv_cen2 3 from sa 1 i j k-1 t-1 |
---|
414 | |
---|
415 | ctin ta_ver_traadv_cen2 1 from ta_hor_traadv_cen2 1 i j k t |
---|
416 | ctin ta_ver_traadv_cen2 2 from zcent_ver_traadv_cen2 1 i j k t |
---|
417 | ctin ta_ver_traadv_cen2 3 from zcent_ver_traadv_cen2 1 i j k+1 t |
---|
418 | |
---|
419 | ctin sa_ver_traadv_cen2 1 from sa_hor_traadv_cen2 1 i j k t |
---|
420 | ctin sa_ver_traadv_cen2 2 from zcens_ver_traadv_cen2 1 i j k t |
---|
421 | ctin sa_ver_traadv_cen2 3 from zcens_ver_traadv_cen2 1 i j k+1 t |
---|
422 | |
---|
423 | //#---------->traldf_lap |
---|
424 | ctin ztuv_traldf_lap 1 from tb 1 i j k t-1 |
---|
425 | ctin ztuv_traldf_lap 2 from tb 1 i+1 j k t-1 |
---|
426 | ctin ztuv_traldf_lap 3 from tb 1 i j+1 k t-1 |
---|
427 | |
---|
428 | ctin zsuv_traldf_lap 1 from sb 1 i j k t-1 |
---|
429 | ctin zsuv_traldf_lap 2 from sb 1 i+1 j k t-1 |
---|
430 | ctin zsuv_traldf_lap 3 from sb 1 i j+1 k t-1 |
---|
431 | |
---|
432 | ctin ta_traldf_lap 1 from ta_ver_traadv_cen2 1 i j k t |
---|
433 | ctin ta_traldf_lap 2 from ztuv_traldf_lap 1 i j k t |
---|
434 | ctin ta_traldf_lap 3 from ztuv_traldf_lap 1 i-1 j k t |
---|
435 | ctin ta_traldf_lap 4 from ztuv_traldf_lap 2 i j k t |
---|
436 | ctin ta_traldf_lap 5 from ztuv_traldf_lap 2 i j-1 k t |
---|
437 | |
---|
438 | ctin sa_traldf_lap 1 from sa_ver_traadv_cen2 1 i j k t |
---|
439 | ctin sa_traldf_lap 2 from zsuv_traldf_lap 1 i j k t |
---|
440 | ctin sa_traldf_lap 3 from zsuv_traldf_lap 1 i-1 j k t |
---|
441 | ctin sa_traldf_lap 4 from zsuv_traldf_lap 2 i j k t |
---|
442 | ctin sa_traldf_lap 5 from zsuv_traldf_lap 2 i j-1 k t |
---|
443 | |
---|
444 | //#---------->trazdf_imp |
---|
445 | ctin zwisd_trazdf_imp 1 from zdf_avt_avmu_avmv 1 i j k t |
---|
446 | ctin zwisd_trazdf_imp 2 from zdf_avt_avmu_avmv 1 i j k+1 t |
---|
447 | |
---|
448 | ctin zrhs_trazdf_imp 1 from ta_traldf_lap 1 i j k t |
---|
449 | ctin zrhs_trazdf_imp 2 from tb 1 i j k t-1 |
---|
450 | |
---|
451 | ctin zrhs_sa_trazdf_imp 1 from sa_traldf_lap 1 i j k t |
---|
452 | ctin zrhs_sa_trazdf_imp 2 from sb 1 i j k t-1 |
---|
453 | |
---|
454 | ctin zwt_trazdf_imp 1 from zwisd_trazdf_imp 1 i j k t |
---|
455 | ctin zwt_trazdf_imp 2 from zwisd_trazdf_imp 2 i j k-1 t |
---|
456 | ctin zwt_trazdf_imp 3 from zwisd_trazdf_imp 3 i j k t |
---|
457 | ctin zwt_trazdf_imp 4 from zwt_trazdf_imp 1 i j k-1 t |
---|
458 | |
---|
459 | ctin ta_de_trazdf_imp 1 from zwisd_trazdf_imp 1 i j k t |
---|
460 | ctin ta_de_trazdf_imp 2 from ta_de_trazdf_imp 1 i j k-1 t |
---|
461 | ctin ta_de_trazdf_imp 3 from zrhs_trazdf_imp 1 i j k t |
---|
462 | ctin ta_de_trazdf_imp 4 from zwt_trazdf_imp 1 i j k-1 t |
---|
463 | |
---|
464 | ctin ta 1 from zwisd_trazdf_imp 2 i j k t |
---|
465 | ctin ta 2 from ta 1 i j k+1 t |
---|
466 | ctin ta 3 from ta_de_trazdf_imp 1 i j k t |
---|
467 | ctin ta 4 from zwt_trazdf_imp 1 i j k t |
---|
468 | ctin ta 5 from ta_c 1 i j k |
---|
469 | |
---|
470 | ctin sa_de_trazdf_imp 1 from zwisd_trazdf_imp 1 i j k t |
---|
471 | ctin sa_de_trazdf_imp 2 from sa_de_trazdf_imp 1 i j k-1 t |
---|
472 | ctin sa_de_trazdf_imp 3 from zrhs_sa_trazdf_imp 1 i j k t |
---|
473 | ctin sa_de_trazdf_imp 4 from zwt_trazdf_imp 1 i j k-1 t |
---|
474 | |
---|
475 | ctin sa 1 from zwisd_trazdf_imp 2 i j k t |
---|
476 | ctin sa 2 from sa 1 i j k+1 t |
---|
477 | ctin sa 3 from sa_de_trazdf_imp 1 i j k t |
---|
478 | ctin sa 4 from zwt_trazdf_imp 1 i j k t |
---|
479 | ctin sa 5 from sa_c 1 i j k |
---|
480 | |
---|
481 | //#---------->tra_nxt |
---|
482 | ctin tb 1 from tb 1 i j k t-1 |
---|
483 | ctin tb 2 from ta 1 i j k t-1 |
---|
484 | ctin tb 3 from ta 1 i j k t |
---|
485 | ctin tb 4 from ta_c 1 i j k |
---|
486 | |
---|
487 | ctin sb 1 from sb 1 i j k t-1 |
---|
488 | ctin sb 2 from sa 1 i j k t-1 |
---|
489 | ctin sb 3 from sa 1 i j k t |
---|
490 | ctin sb 4 from sa_c 1 i j k |
---|
491 | |
---|
492 | //#---------->eos |
---|
493 | ctin rhd 1 from tb 1 i j k t |
---|
494 | ctin rhd 2 from sb 1 i j k t |
---|
495 | |
---|
496 | ctin rhop 1 from rhd 1 i j k t |
---|
497 | |
---|
498 | |
---|
499 | //#---------->dynadv_cen2 |
---|
500 | ctin hmf_TF_dynadv_cen2 1 from ua 1 i j k t-1 |
---|
501 | ctin hmf_TF_dynadv_cen2 2 from ua 1 i+1 j k t-1 |
---|
502 | ctin hmf_TF_dynadv_cen2 3 from ua 1 i j+1 k t-1 |
---|
503 | ctin hmf_TF_dynadv_cen2 4 from va 1 i j k t-1 |
---|
504 | ctin hmf_TF_dynadv_cen2 5 from va 1 i+1 j k t-1 |
---|
505 | ctin hmf_TF_dynadv_cen2 6 from va 1 i j+1 k t-1 |
---|
506 | |
---|
507 | ctin zua_dhmf_dynadv_cen2 1 from hmf_TF_dynadv_cen2 1 i j k t |
---|
508 | ctin zua_dhmf_dynadv_cen2 2 from hmf_TF_dynadv_cen2 1 i-1 j k t |
---|
509 | ctin zua_dhmf_dynadv_cen2 3 from hmf_TF_dynadv_cen2 2 i j k t |
---|
510 | ctin zua_dhmf_dynadv_cen2 4 from hmf_TF_dynadv_cen2 2 i j-1 k t |
---|
511 | |
---|
512 | ctin zva_dhmf_dynadv_cen2 1 from hmf_TF_dynadv_cen2 3 i j k t |
---|
513 | ctin zva_dhmf_dynadv_cen2 2 from hmf_TF_dynadv_cen2 3 i-1 j k t |
---|
514 | ctin zva_dhmf_dynadv_cen2 3 from hmf_TF_dynadv_cen2 4 i j k t |
---|
515 | ctin zva_dhmf_dynadv_cen2 4 from hmf_TF_dynadv_cen2 4 i j-1 k t |
---|
516 | |
---|
517 | ctin zfuv_UW_dynadv_cen2 1 from wa 1 i j k t-1 |
---|
518 | ctin zfuv_UW_dynadv_cen2 2 from wa 1 i+1 j k t-1 |
---|
519 | ctin zfuv_UW_dynadv_cen2 3 from wa 1 i j+1 k t-1 |
---|
520 | ctin zfuv_UW_dynadv_cen2 4 from ua 1 i j k t-1 |
---|
521 | ctin zfuv_UW_dynadv_cen2 5 from ua 1 i j k-1 t-1 |
---|
522 | ctin zfuv_UW_dynadv_cen2 6 from va 1 i j k t-1 |
---|
523 | ctin zfuv_UW_dynadv_cen2 7 from va 1 i j k-1 t-1 |
---|
524 | |
---|
525 | ctin ua_dynadv_cen2 1 from zfuv_UW_dynadv_cen2 1 i j k t |
---|
526 | ctin ua_dynadv_cen2 2 from zfuv_UW_dynadv_cen2 1 i j k+1 t |
---|
527 | ctin ua_dynadv_cen2 3 from zua_dhmf_dynadv_cen2 1 i j k t |
---|
528 | |
---|
529 | ctin va_dynadv_cen2 1 from zfuv_UW_dynadv_cen2 2 i j k t |
---|
530 | ctin va_dynadv_cen2 2 from zfuv_UW_dynadv_cen2 2 i j k+1 t |
---|
531 | ctin va_dynadv_cen2 3 from zva_dhmf_dynadv_cen2 1 i j k t |
---|
532 | |
---|
533 | //#---------->vor_ene |
---|
534 | ctin zwxyz_vor_ene 1 from ua 1 i j k t-1 |
---|
535 | ctin zwxyz_vor_ene 2 from ua 1 i j+1 k t-1 |
---|
536 | ctin zwxyz_vor_ene 3 from va 1 i j k t-1 |
---|
537 | ctin zwxyz_vor_ene 4 from va 1 i+1 j k t-1 |
---|
538 | |
---|
539 | ctin ua_vor_ene 1 from ua_dynadv_cen2 1 i j k t |
---|
540 | ctin ua_vor_ene 2 from zwxyz_vor_ene 3 i j-1 k t |
---|
541 | ctin ua_vor_ene 3 from zwxyz_vor_ene 3 i j k t |
---|
542 | ctin ua_vor_ene 4 from zwxyz_vor_ene 2 i j-1 k t |
---|
543 | ctin ua_vor_ene 5 from zwxyz_vor_ene 2 i j k t |
---|
544 | ctin ua_vor_ene 6 from zwxyz_vor_ene 2 i+1 j-1 k t |
---|
545 | ctin ua_vor_ene 7 from zwxyz_vor_ene 2 i+1 j k t |
---|
546 | |
---|
547 | ctin va_vor_ene 1 from va_dynadv_cen2 1 i j k t |
---|
548 | ctin va_vor_ene 2 from zwxyz_vor_ene 3 i-1 j k t |
---|
549 | ctin va_vor_ene 3 from zwxyz_vor_ene 3 i j k t |
---|
550 | ctin va_vor_ene 4 from zwxyz_vor_ene 1 i-1 j k t |
---|
551 | ctin va_vor_ene 5 from zwxyz_vor_ene 1 i-1 j+1 k t |
---|
552 | ctin va_vor_ene 6 from zwxyz_vor_ene 1 i j k t |
---|
553 | ctin va_vor_ene 7 from zwxyz_vor_ene 1 i j+1 k t |
---|
554 | |
---|
555 | //#---------->dyn_ldf_lap |
---|
556 | ctin ua_dyn_ldf_lap 1 from ua_vor_ene 1 i j k t |
---|
557 | ctin ua_dyn_ldf_lap 2 from rotn 1 i j-1 k t-2 |
---|
558 | ctin ua_dyn_ldf_lap 3 from rotn 1 i j k t-2 |
---|
559 | ctin ua_dyn_ldf_lap 4 from hdivn 1 i j k t-2 |
---|
560 | ctin ua_dyn_ldf_lap 5 from hdivn 1 i+1 j k t-2 |
---|
561 | |
---|
562 | ctin va_dyn_ldf_lap 1 from va_vor_ene 1 i j k t |
---|
563 | ctin va_dyn_ldf_lap 2 from rotn 1 i-1 j k t-2 |
---|
564 | ctin va_dyn_ldf_lap 3 from rotn 1 i j k t-2 |
---|
565 | ctin va_dyn_ldf_lap 4 from hdivn 1 i j k t-2 |
---|
566 | ctin va_dyn_ldf_lap 5 from hdivn 1 i j+1 k t-2 |
---|
567 | |
---|
568 | //#---------->dyn__hpg_zco |
---|
569 | ctin zhpij_dyn_hpg_zco 1 from rhd 1 i+1 j k t |
---|
570 | ctin zhpij_dyn_hpg_zco 2 from rhd 1 i j+1 k t |
---|
571 | ctin zhpij_dyn_hpg_zco 3 from rhd 1 i j k t |
---|
572 | ctin zhpij_dyn_hpg_zco 4 from rhd 1 i+1 j k-1 t |
---|
573 | ctin zhpij_dyn_hpg_zco 5 from rhd 1 i j+1 k-1 t |
---|
574 | ctin zhpij_dyn_hpg_zco 6 from rhd 1 i j k-1 t |
---|
575 | ctin zhpij_dyn_hpg_zco 7 from zhpij_dyn_hpg_zco 1 i j k-1 t |
---|
576 | ctin zhpij_dyn_hpg_zco 8 from zhpij_dyn_hpg_zco 2 i j k-1 t |
---|
577 | |
---|
578 | ctin ua_dyn_hpg_zco 1 from ua_dyn_ldf_lap 1 i j k t |
---|
579 | ctin ua_dyn_hpg_zco 2 from zhpij_dyn_hpg_zco 1 i j k t |
---|
580 | |
---|
581 | ctin va_dyn_hpg_zco 1 from va_dyn_ldf_lap 1 i j k t |
---|
582 | ctin va_dyn_hpg_zco 2 from zhpij_dyn_hpg_zco 2 i j k t |
---|
583 | |
---|
584 | //#---------->dynzdf_imp |
---|
585 | ctin zwisd_ua_dynzdf_imp 1 from zdf_avt_avmu_avmv 2 i j k t |
---|
586 | ctin zwisd_ua_dynzdf_imp 2 from zdf_avt_avmu_avmv 2 i j k+1 t |
---|
587 | |
---|
588 | ctin zwisd_va_dynzdf_imp 1 from zdf_avt_avmu_avmv 3 i j k t |
---|
589 | ctin zwisd_va_dynzdf_imp 2 from zdf_avt_avmu_avmv 3 i j k+1 t |
---|
590 | |
---|
591 | ctin zrhs_ua_dynzdf_imp 1 from ua_dyn_hpg_zco 1 i j k t |
---|
592 | ctin zrhs_ua_dynzdf_imp 2 from ub 1 i j k t-1 |
---|
593 | ctin zrhs_ua_dynzdf_imp 3 from sbcgyre_utau_vtau 1 i j t |
---|
594 | |
---|
595 | ctin zrhs_va_dynzdf_imp 1 from va_dyn_hpg_zco 1 i j k t |
---|
596 | ctin zrhs_va_dynzdf_imp 2 from vb 1 i j k t-1 |
---|
597 | ctin zrhs_va_dynzdf_imp 3 from sbcgyre_utau_vtau 2 i j t |
---|
598 | |
---|
599 | ctin zwt_ua_dynzdf_imp 1 from zwisd_ua_dynzdf_imp 1 i j k t |
---|
600 | ctin zwt_ua_dynzdf_imp 2 from zwisd_ua_dynzdf_imp 2 i j k-1 t |
---|
601 | ctin zwt_ua_dynzdf_imp 3 from zwisd_ua_dynzdf_imp 3 i j k t |
---|
602 | ctin zwt_ua_dynzdf_imp 4 from zwt_ua_dynzdf_imp 1 i j k-1 t |
---|
603 | |
---|
604 | ctin zwt_va_dynzdf_imp 1 from zwisd_va_dynzdf_imp 1 i j k t |
---|
605 | ctin zwt_va_dynzdf_imp 2 from zwisd_va_dynzdf_imp 2 i j k-1 t |
---|
606 | ctin zwt_va_dynzdf_imp 3 from zwisd_va_dynzdf_imp 3 i j k t |
---|
607 | ctin zwt_va_dynzdf_imp 4 from zwt_va_dynzdf_imp 1 i j k-1 t |
---|
608 | |
---|
609 | ctin ua_de_dynzdf_imp 1 from zwisd_ua_dynzdf_imp 1 i j k t |
---|
610 | ctin ua_de_dynzdf_imp 2 from ua_de_dynzdf_imp 1 i j k-1 t |
---|
611 | ctin ua_de_dynzdf_imp 3 from zrhs_ua_dynzdf_imp 1 i j k t |
---|
612 | ctin ua_de_dynzdf_imp 4 from zwt_ua_dynzdf_imp 1 i j k-1 t |
---|
613 | |
---|
614 | ctin ua_dynzdf_imp 1 from zwisd_ua_dynzdf_imp 2 i j k t |
---|
615 | ctin ua_dynzdf_imp 2 from ua_dynzdf_imp 1 i j k+1 t |
---|
616 | ctin ua_dynzdf_imp 3 from ua_de_dynzdf_imp 1 i j k t |
---|
617 | ctin ua_dynzdf_imp 4 from zwt_ua_dynzdf_imp 1 i j k t |
---|
618 | |
---|
619 | ctin va_de_dynzdf_imp 1 from zwisd_va_dynzdf_imp 1 i j k t |
---|
620 | ctin va_de_dynzdf_imp 2 from va_de_dynzdf_imp 1 i j k-1 t |
---|
621 | ctin va_de_dynzdf_imp 3 from zrhs_va_dynzdf_imp 1 i j k t |
---|
622 | ctin va_de_dynzdf_imp 4 from zwt_va_dynzdf_imp 1 i j k-1 t |
---|
623 | |
---|
624 | ctin va_dynzdf_imp 1 from zwisd_va_dynzdf_imp 2 i j k t |
---|
625 | ctin va_dynzdf_imp 2 from va_dynzdf_imp 1 i j k+1 t |
---|
626 | ctin va_dynzdf_imp 3 from va_de_dynzdf_imp 1 i j k t |
---|
627 | ctin va_dynzdf_imp 4 from zwt_va_dynzdf_imp 1 i j k t |
---|
628 | |
---|
629 | ctin ua_norm_dynzdf 1 from ua_dynzdf_imp 1 i j k t |
---|
630 | ctin ua_norm_dynzdf 2 from ub 1 i j k t-1 |
---|
631 | |
---|
632 | ctin va_norm_dynzdf 1 from va_dynzdf_imp 1 i j k t |
---|
633 | ctin va_norm_dynzdf 2 from vb 1 i j k t-1 |
---|
634 | |
---|
635 | //#---------->ddynspg_flt |
---|
636 | ctin spguv_now_dynspg 1 from sshn 1 i j t-1 |
---|
637 | ctin spguv_now_dynspg 2 from sshn 1 i j+1 t-1 |
---|
638 | ctin spguv_now_dynspg 3 from sshn 1 i+1 j t-1 |
---|
639 | |
---|
640 | ctin ua_1_dynspg_flt 1 from ua_norm_dynzdf 1 i j k t |
---|
641 | ctin ua_1_dynspg_flt 2 from ub 1 i j k t-1 |
---|
642 | ctin ua_1_dynspg_flt 3 from spguv_now_dynspg 1 i j t |
---|
643 | |
---|
644 | ctin va_1_dynspg_flt 1 from va_norm_dynzdf 1 i j k t |
---|
645 | ctin va_1_dynspg_flt 2 from vb 1 i j k t-1 |
---|
646 | ctin va_1_dynspg_flt 3 from spguv_now_dynspg 2 i j t |
---|
647 | |
---|
648 | ctin sum_spgu_dynspg_flt 1..NZ from ua_1_dynspg_flt 1 i j 1..NZ t |
---|
649 | ctin sum_spgv_dynspg_flt 1..NZ from va_1_dynspg_flt 1 i j 1..NZ t |
---|
650 | |
---|
651 | ctin gcb_dynspg_flt 1 from sum_spgu_dynspg_flt 1 i j t |
---|
652 | ctin gcb_dynspg_flt 2 from sum_spgu_dynspg_flt 1 i-1 j t |
---|
653 | ctin gcb_dynspg_flt 3 from sum_spgv_dynspg_flt 1 i j t |
---|
654 | ctin gcb_dynspg_flt 4 from sum_spgv_dynspg_flt 1 i j-1 t |
---|
655 | |
---|
656 | ctin gcx_dynspg_flt 1 from gcx2 1 i j t-1 |
---|
657 | ctin gcx_dynspg_flt 2 from gcx2 1 i j t-2 |
---|
658 | |
---|
659 | ctin spguv_dynspg_flt 1 from gcx2 1 i+1 j t |
---|
660 | ctin spguv_dynspg_flt 2 from gcx2 1 i j t |
---|
661 | ctin spguv_dynspg_flt 3 from gcx2 1 i j+1 t |
---|
662 | |
---|
663 | ctin ua 1 from ua_1_dynspg_flt 1 i j k t |
---|
664 | ctin ua 2 from spguv_dynspg_flt 1 i j t |
---|
665 | ctin ua 3 from ua_c 1 i j k |
---|
666 | |
---|
667 | ctin va 1 from va_1_dynspg_flt 1 i j k t |
---|
668 | ctin va 2 from spguv_dynspg_flt 2 i j t |
---|
669 | ctin va 3 from va_c 1 i j k |
---|
670 | |
---|
671 | ctin sshn 1 from sshb 1 i j t-1 |
---|
672 | ctin sshn 2 from wa 1 i j 1 t-1 |
---|
673 | ctin sshn 3 from sbcgyre_emp 1 i j t |
---|
674 | ctin sshn 4 from sshn_c 1 i j |
---|
675 | |
---|
676 | ctin sshb 1 from sshb 1 i j t-1 |
---|
677 | ctin sshb 2 from sshn 1 i j t-1 |
---|
678 | ctin sshb 3 from sshn 1 i j t |
---|
679 | ctin sshb 4 from sshn_c 1 i j |
---|
680 | |
---|
681 | #ifndef SOLSORYAO |
---|
682 | //############## Ajout de ctin "virtuels" ############## |
---|
683 | ctin solsor_dynspg_flt 1..9 from gcx_dynspg_flt 1 i%~1 j%~1 t |
---|
684 | ctin solsor_dynspg_flt 10..18 from gcb_dynspg_flt 1 i%~1 j%~1 t |
---|
685 | ctin solsor_dynspg_flt 19 from solsor_dynspg_flt 1 i-1 j-1 t |
---|
686 | ctin gcx2 1 from solsor_dynspg_flt 1 i j t |
---|
687 | #endif |
---|
688 | |
---|
689 | #ifdef SOLSORYAO |
---|
690 | ctin solsor_flag 1 from solsor_gx2_x_r 2 i j k t |
---|
691 | ctin solsor_flag 2 from solsor_flag 2 NX NY k-1 t |
---|
692 | ctin solsor_flag 3 from solsor_flag 1 i-1 NY k t |
---|
693 | ctin solsor_flag 4 from solsor_flag 1 i j-1 k t |
---|
694 | |
---|
695 | ctin solsor_gx1_x_r 1 from gcx_dynspg_flt 1 i j t |
---|
696 | ctin solsor_gx1_x_r 2 from gcb_dynspg_flt 1 i j t |
---|
697 | ctin solsor_gx1_x_r 3 from solsor_gx2_x_r 1 i j-1 k-1 t |
---|
698 | ctin solsor_gx1_x_r 4 from solsor_gx2_x_r 1 i-1 j k-1 t |
---|
699 | ctin solsor_gx1_x_r 5 from solsor_gx2_x_r 1 i+1 j k-1 t |
---|
700 | ctin solsor_gx1_x_r 6 from solsor_gx2_x_r 1 i j+1 k-1 t |
---|
701 | ctin solsor_gx1_x_r 7 from solsor_gx2_x_r 1 i j k-1 t |
---|
702 | ctin solsor_gx1_x_r 8 from solsor_flag 2 NX NY k-1 t |
---|
703 | |
---|
704 | ctin solsor_gx2_x_r 1 from gcx_dynspg_flt 1 i j t |
---|
705 | ctin solsor_gx2_x_r 2 from gcb_dynspg_flt 1 i j t |
---|
706 | ctin solsor_gx2_x_r 3 from solsor_gx1_x_r 1 i j-1 k t |
---|
707 | ctin solsor_gx2_x_r 4 from solsor_gx1_x_r 1 i-1 j k t |
---|
708 | ctin solsor_gx2_x_r 5 from solsor_gx1_x_r 1 i+1 j k t |
---|
709 | ctin solsor_gx2_x_r 6 from solsor_gx1_x_r 1 i j+1 k t |
---|
710 | ctin solsor_gx2_x_r 7 from solsor_gx1_x_r 1 i j k t |
---|
711 | ctin solsor_gx2_x_r 8 from solsor_flag 2 NX NY k-1 t |
---|
712 | ctin solsor_gx2_x_r 9 from solsor_gx1_x_r 2 i j k t |
---|
713 | |
---|
714 | ctin gcx2 1 from solsor_gx2_x_r 1 i j NMAX t |
---|
715 | #endif |
---|
716 | |
---|
717 | //#---------->dyn_nxt |
---|
718 | ctin ub 1 from ub 1 i j k t-1 |
---|
719 | ctin ub 2 from ua 1 i j k t-1 |
---|
720 | ctin ub 3 from ua 1 i j k t |
---|
721 | ctin ub 4 from ua_c 1 i j k |
---|
722 | |
---|
723 | ctin vb 1 from vb 1 i j k t-1 |
---|
724 | ctin vb 2 from va 1 i j k t-1 |
---|
725 | ctin vb 3 from va 1 i j k t |
---|
726 | ctin vb 4 from va_c 1 i j k |
---|
727 | |
---|
728 | //#---------->divcur |
---|
729 | ctin hdivn 1 from ua 1 i j k t |
---|
730 | ctin hdivn 2 from ua 1 i-1 j k t |
---|
731 | ctin hdivn 3 from va 1 i j k t |
---|
732 | ctin hdivn 4 from va 1 i j-1 k t |
---|
733 | |
---|
734 | ctin rotn 1 from va 1 i+1 j k t |
---|
735 | ctin rotn 2 from va 1 i j k t |
---|
736 | ctin rotn 3 from ua 1 i j+1 k t |
---|
737 | ctin rotn 4 from ua 1 i j k t |
---|
738 | |
---|
739 | //#---------->wzv |
---|
740 | ctin wa 1 from wa 1 i j k+1 t |
---|
741 | ctin wa 2 from hdivn 1 i j k t |
---|
742 | //#----------> |
---|
743 | |
---|
744 | #ifdef OPTIMORDER |
---|
745 | |
---|
746 | //#---------->sbcgyre + bn2_rn2 + zdf + trasbc + traqsr + traadv + |
---|
747 | //# traldf + trazdf + tra_nxt + eos + dynadv_cen2 + vor_ene |
---|
748 | //# +dyn_ldf + dyn__hpg_zco +dynzdf_imp(-> gcx_dynspg_flt) // |
---|
749 | |
---|
750 | |
---|
751 | order modinspace S3dt |
---|
752 | order YA1 YA2 |
---|
753 | sbcgyre_qsr_qns sbcgyre_emp1 |
---|
754 | spguv_now_dynspg |
---|
755 | gcx_dynspg_flt |
---|
756 | order YA3 |
---|
757 | bn2_rn2 zcent_hor_traadv_cen2 zcens_hor_traadv_cen2 |
---|
758 | zcent_ver_traadv_cen2 zcens_ver_traadv_cen2 |
---|
759 | ztuv_traldf_lap zsuv_traldf_lap |
---|
760 | hmf_TF_dynadv_cen2 zfuv_UW_dynadv_cen2 |
---|
761 | forder |
---|
762 | forder |
---|
763 | |
---|
764 | order YA3 |
---|
765 | order YA1 YA2 |
---|
766 | zua_dhmf_dynadv_cen2 zva_dhmf_dynadv_cen2 zwxyz_vor_ene |
---|
767 | forder |
---|
768 | forder |
---|
769 | |
---|
770 | //#NO PARALLEL |
---|
771 | order YA2 YA1 |
---|
772 | sbcgyre_zsumemp_zsurf |
---|
773 | forder |
---|
774 | order YA1 YA2 |
---|
775 | order YB3 |
---|
776 | zdf_avt_avmu_avmv |
---|
777 | zwisd_trazdf_imp |
---|
778 | zwisd_ua_dynzdf_imp zwisd_va_dynzdf_imp |
---|
779 | forder |
---|
780 | sbcgyre_emp sbcgyre_emps sbcgyre_utau_vtau |
---|
781 | zdfmxl_hmld zdfmxl_nmln_hmlp_hmlpt ta_sa_trasbc |
---|
782 | order YB3 |
---|
783 | ta_traqsr |
---|
784 | ua_dynadv_cen2 va_dynadv_cen2 |
---|
785 | ua_vor_ene va_vor_ene |
---|
786 | ua_dyn_ldf_lap va_dyn_ldf_lap |
---|
787 | ta_hor_traadv_cen2 sa_hor_traadv_cen2 |
---|
788 | ta_ver_traadv_cen2 sa_ver_traadv_cen2 |
---|
789 | ta_traldf_lap sa_traldf_lap |
---|
790 | zrhs_trazdf_imp zrhs_sa_trazdf_imp |
---|
791 | forder |
---|
792 | order YA3 |
---|
793 | zwt_trazdf_imp ta_de_trazdf_imp sa_de_trazdf_imp |
---|
794 | zwt_ua_dynzdf_imp zwt_va_dynzdf_imp |
---|
795 | forder |
---|
796 | order YB3 |
---|
797 | ta sa tb sb rhd rhop |
---|
798 | forder |
---|
799 | forder |
---|
800 | order YA1 YA2 |
---|
801 | order YA3 |
---|
802 | zhpij_dyn_hpg_zco |
---|
803 | ua_dyn_hpg_zco va_dyn_hpg_zco |
---|
804 | zrhs_ua_dynzdf_imp zrhs_va_dynzdf_imp |
---|
805 | ua_de_dynzdf_imp va_de_dynzdf_imp |
---|
806 | forder |
---|
807 | order YB3 |
---|
808 | ua_dynzdf_imp va_dynzdf_imp |
---|
809 | ua_norm_dynzdf va_norm_dynzdf |
---|
810 | ua_1_dynspg_flt va_1_dynspg_flt |
---|
811 | forder |
---|
812 | sum_spgu_dynspg_flt sum_spgv_dynspg_flt |
---|
813 | forder |
---|
814 | |
---|
815 | order YA1 YA2 |
---|
816 | sbcgyre_wndm |
---|
817 | gcb_dynspg_flt |
---|
818 | |
---|
819 | forder |
---|
820 | #ifdef SOLSORYAO |
---|
821 | forder |
---|
822 | #endif |
---|
823 | |
---|
824 | #ifndef SOLSORYAO |
---|
825 | order YA1 YA2 |
---|
826 | solsor_dynspg_flt |
---|
827 | forder |
---|
828 | #endif |
---|
829 | |
---|
830 | #ifdef SOLSORYAO |
---|
831 | order modinspace Ssolsor |
---|
832 | order YA3 |
---|
833 | order YA2 YA1 |
---|
834 | solsor_gx1_x_r |
---|
835 | forder |
---|
836 | |
---|
837 | order YA2 YA1 |
---|
838 | solsor_gx2_x_r |
---|
839 | forder |
---|
840 | |
---|
841 | order YA1 YA2 |
---|
842 | solsor_flag |
---|
843 | forder |
---|
844 | forder |
---|
845 | forder |
---|
846 | #endif |
---|
847 | |
---|
848 | #ifdef SOLSORYAO |
---|
849 | order modinspace S3dt_b |
---|
850 | #endif |
---|
851 | order YA2 YA1 |
---|
852 | gcx2 |
---|
853 | sshn sshb |
---|
854 | forder |
---|
855 | order YA2 YA1 |
---|
856 | spguv_dynspg_flt |
---|
857 | order YA3 |
---|
858 | ua va |
---|
859 | ub vb |
---|
860 | forder |
---|
861 | forder |
---|
862 | |
---|
863 | // #---------->divcur |
---|
864 | order YA3 YA2 YA1 |
---|
865 | hdivn rotn |
---|
866 | forder |
---|
867 | // #---------->wzv |
---|
868 | order YA2 YA1 YB3 |
---|
869 | wa |
---|
870 | |
---|
871 | forder |
---|
872 | |
---|
873 | |
---|
874 | forder |
---|
875 | |
---|
876 | #endif |
---|
877 | |
---|
878 | #ifdef FILTER |
---|
879 | // calcul du filtre |
---|
880 | |
---|
881 | order modinspace S3dtlf |
---|
882 | order YA3 |
---|
883 | order YA2 YA1 |
---|
884 | ztuv_lfil |
---|
885 | forder |
---|
886 | order YA2 YA1 |
---|
887 | dta_lfil |
---|
888 | forder |
---|
889 | forder |
---|
890 | forder |
---|
891 | |
---|
892 | order modinspace S3dtzf |
---|
893 | order YA3 YA2 YA1 |
---|
894 | zwy_zfil |
---|
895 | forder |
---|
896 | order YA3 YA2 YA1 |
---|
897 | dta_zfil |
---|
898 | forder |
---|
899 | forder |
---|
900 | order modinspace S3d |
---|
901 | order YA3 YA2 YA1 |
---|
902 | ta_c |
---|
903 | forder |
---|
904 | forder |
---|
905 | |
---|
906 | #endif |
---|
907 | |
---|
908 | #ifndef OPTIMORDER |
---|
909 | |
---|
910 | order modinspace S3dt |
---|
911 | //#---------->sbcgyre |
---|
912 | order YA2 YA1 |
---|
913 | sbcgyre_qsr_qns sbcgyre_emp1 sbcgyre_zsumemp_zsurf |
---|
914 | forder |
---|
915 | order YA2 YA1 |
---|
916 | sbcgyre_emp sbcgyre_emps sbcgyre_utau_vtau sbcgyre_wndm |
---|
917 | forder |
---|
918 | |
---|
919 | |
---|
920 | //#---------->bn2 |
---|
921 | order YA3 YA2 YA1 |
---|
922 | bn2_rn2 |
---|
923 | forder |
---|
924 | |
---|
925 | //#---------->zdf |
---|
926 | order YA3 YA2 YA1 |
---|
927 | zdf_avt_avmu_avmv |
---|
928 | forder |
---|
929 | order YA2 YA1 |
---|
930 | zdfmxl_hmld zdfmxl_nmln_hmlp_hmlpt |
---|
931 | forder |
---|
932 | //#---------->trasbc |
---|
933 | order YA2 YA1 |
---|
934 | ta_sa_trasbc |
---|
935 | forder |
---|
936 | //#---------->traqsr |
---|
937 | order YA3 YA2 YA1 |
---|
938 | ta_traqsr |
---|
939 | forder |
---|
940 | //#------> traadv_cen2 |
---|
941 | order YA3 |
---|
942 | order YA2 YA1 |
---|
943 | zcent_hor_traadv_cen2 zcens_hor_traadv_cen2 |
---|
944 | forder |
---|
945 | order YA2 YA1 |
---|
946 | ta_hor_traadv_cen2 sa_hor_traadv_cen2 |
---|
947 | forder |
---|
948 | forder |
---|
949 | |
---|
950 | order YA3 YA2 YA1 |
---|
951 | zcent_ver_traadv_cen2 zcens_ver_traadv_cen2 |
---|
952 | forder |
---|
953 | order YA3 YA2 YA1 |
---|
954 | ta_ver_traadv_cen2 sa_ver_traadv_cen2 |
---|
955 | forder |
---|
956 | |
---|
957 | //#------> traldf_lap |
---|
958 | order YA3 YA2 YA1 |
---|
959 | ztuv_traldf_lap zsuv_traldf_lap |
---|
960 | ta_traldf_lap sa_traldf_lap |
---|
961 | forder |
---|
962 | //#------> trazdf_imp |
---|
963 | order YA3 YA2 YA1 |
---|
964 | zwisd_trazdf_imp zrhs_trazdf_imp zrhs_sa_trazdf_imp |
---|
965 | forder |
---|
966 | order YA3 YA2 YA1 |
---|
967 | zwt_trazdf_imp ta_de_trazdf_imp sa_de_trazdf_imp |
---|
968 | forder |
---|
969 | order YB3 YA2 YA1 |
---|
970 | ta sa |
---|
971 | forder |
---|
972 | //#---------->tra_nxt |
---|
973 | order YA3 YA2 YA1 |
---|
974 | tb sb |
---|
975 | forder |
---|
976 | //#---------->eos |
---|
977 | order YA3 YA2 YA1 |
---|
978 | rhd rhop |
---|
979 | forder |
---|
980 | //#------> dynadv_cen2 |
---|
981 | order YA3 |
---|
982 | order YA2 YA1 |
---|
983 | hmf_TF_dynadv_cen2 |
---|
984 | forder |
---|
985 | order YA2 YA1 |
---|
986 | zua_dhmf_dynadv_cen2 zva_dhmf_dynadv_cen2 |
---|
987 | forder |
---|
988 | forder |
---|
989 | |
---|
990 | order YA3 YA2 YA1 |
---|
991 | zfuv_UW_dynadv_cen2 |
---|
992 | forder |
---|
993 | order YA3 YA2 YA1 |
---|
994 | ua_dynadv_cen2 va_dynadv_cen2 |
---|
995 | forder |
---|
996 | //#------> vor_ene && dyn_ldf_lap |
---|
997 | order YA3 |
---|
998 | order YA2 YA1 |
---|
999 | zwxyz_vor_ene |
---|
1000 | forder |
---|
1001 | order YA2 YA1 |
---|
1002 | ua_vor_ene va_vor_ene |
---|
1003 | ua_dyn_ldf_lap va_dyn_ldf_lap |
---|
1004 | forder |
---|
1005 | forder |
---|
1006 | //#---------->dyn__hpg_zco |
---|
1007 | order YA3 YA2 YA1 |
---|
1008 | zhpij_dyn_hpg_zco ua_dyn_hpg_zco va_dyn_hpg_zco |
---|
1009 | forder |
---|
1010 | //#---------->dynzdf_imp |
---|
1011 | order YA3 YA2 YA1 |
---|
1012 | zwisd_ua_dynzdf_imp zwisd_va_dynzdf_imp |
---|
1013 | zrhs_ua_dynzdf_imp zrhs_va_dynzdf_imp |
---|
1014 | forder |
---|
1015 | order YA3 YA2 YA1 |
---|
1016 | zwt_ua_dynzdf_imp ua_de_dynzdf_imp |
---|
1017 | zwt_va_dynzdf_imp va_de_dynzdf_imp |
---|
1018 | forder |
---|
1019 | order YB3 YA2 YA1 |
---|
1020 | ua_dynzdf_imp va_dynzdf_imp |
---|
1021 | forder |
---|
1022 | order YA3 YA2 YA1 |
---|
1023 | ua_norm_dynzdf va_norm_dynzdf |
---|
1024 | forder |
---|
1025 | //#---------->dynspg_flt |
---|
1026 | order YA2 YA1 |
---|
1027 | spguv_now_dynspg |
---|
1028 | forder |
---|
1029 | order YA3 YA2 YA1 |
---|
1030 | ua_1_dynspg_flt va_1_dynspg_flt |
---|
1031 | forder |
---|
1032 | order YA2 YA1 |
---|
1033 | sum_spgu_dynspg_flt sum_spgv_dynspg_flt |
---|
1034 | forder |
---|
1035 | order YA2 YA1 |
---|
1036 | gcb_dynspg_flt gcx_dynspg_flt |
---|
1037 | forder |
---|
1038 | |
---|
1039 | #ifdef SOLSORYAO |
---|
1040 | forder |
---|
1041 | #endif |
---|
1042 | |
---|
1043 | #ifndef SOLSORYAO |
---|
1044 | order YA1 YA2 |
---|
1045 | solsor_dynspg_flt |
---|
1046 | forder |
---|
1047 | #endif |
---|
1048 | |
---|
1049 | #ifdef SOLSORYAO |
---|
1050 | order modinspace Ssolsor |
---|
1051 | order YA3 |
---|
1052 | order YA2 YA1 |
---|
1053 | solsor_gx1_x_r |
---|
1054 | forder |
---|
1055 | |
---|
1056 | order YA2 YA1 |
---|
1057 | solsor_gx2_x_r |
---|
1058 | forder |
---|
1059 | |
---|
1060 | order YA1 YA2 |
---|
1061 | solsor_flag |
---|
1062 | forder |
---|
1063 | forder |
---|
1064 | forder |
---|
1065 | #endif |
---|
1066 | |
---|
1067 | #ifdef SOLSORYAO |
---|
1068 | order modinspace S3dt_b |
---|
1069 | #endif SOLSORYAO |
---|
1070 | |
---|
1071 | order YA2 YA1 |
---|
1072 | gcx2 |
---|
1073 | forder |
---|
1074 | order YA2 YA1 |
---|
1075 | spguv_dynspg_flt |
---|
1076 | sshn sshb |
---|
1077 | order YA3 |
---|
1078 | ua va |
---|
1079 | forder |
---|
1080 | forder |
---|
1081 | // #---------->dyn_nxt |
---|
1082 | order YA3 YA2 YA1 |
---|
1083 | ub vb |
---|
1084 | forder |
---|
1085 | // #---------->divcur |
---|
1086 | order YA3 YA2 YA1 |
---|
1087 | hdivn rotn |
---|
1088 | forder |
---|
1089 | // #---------->wzv |
---|
1090 | order YB3 YA2 YA1 |
---|
1091 | wa |
---|
1092 | forder |
---|
1093 | forder |
---|
1094 | #endif |
---|
1095 | |
---|
1096 | |
---|
1097 | |
---|
1098 | #ifdef SOLSORYAO |
---|
1099 | order spaceintraj Tsbd |
---|
1100 | S3dt Ssolsor S3dt_b |
---|
1101 | forder |
---|
1102 | #endif |
---|
1103 | |
---|
1104 | #ifdef FILTER |
---|
1105 | |
---|
1106 | order spaceintraj Tlfil |
---|
1107 | S3dtlf |
---|
1108 | forder |
---|
1109 | order spaceintraj Tzfil |
---|
1110 | S3dtzf |
---|
1111 | forder |
---|
1112 | order spaceintraj Tcst |
---|
1113 | S3d |
---|
1114 | forder |
---|
1115 | |
---|
1116 | #endif |
---|
1117 | |
---|
1118 | #ifndef SOLSORYAO |
---|
1119 | |
---|
1120 | order spaceintraj Tsbd |
---|
1121 | S3dt |
---|
1122 | forder |
---|
1123 | #endif |
---|
1124 | |
---|
1125 | //#USER FUNCTIONS ---------------------------------INSERT_FCT |
---|
1126 | insert_fct xdisplay |
---|
1127 | insert_fct xcomparYF |
---|
1128 | insert_fct arg xistate_init |
---|
1129 | insert_fct arg true_target_in_tab |
---|
1130 | insert_fct arg xrst_save |
---|
1131 | insert_fct arg xchangesavemode |
---|
1132 | |
---|
1133 | #ifdef FILTER |
---|
1134 | insert_fct arg load_eb |
---|
1135 | #endif |
---|
1136 | |
---|
1137 | insert_fct arg xwriteout |
---|
1138 | insert_fct arg xwritegrad |
---|
1139 | insert_fct arg xinitnc |
---|
1140 | insert_fct arg init_euler |
---|
1141 | //#insert_fct arg load_shape_func |
---|
1142 | //#insert_fct arg load_stdev_pca |
---|
1143 | //#insert_fct arg load_mean |
---|