1 |
C $Header: /u/gcmpack/MITgcm/pkg/ecco/ecco_cost_legacy.h,v 1.67 2014/10/09 00:50:16 gforget Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
|
5 |
c ================================================================== |
6 |
c HEADER ECCO_legacy |
7 |
c ================================================================== |
8 |
|
9 |
INTEGER NSSHV4COST |
10 |
PARAMETER ( NSSHV4COST=5 ) |
11 |
|
12 |
common /averages_c/ |
13 |
& tbarfile, |
14 |
& sbarfile, |
15 |
#ifdef ALLOW_SIGMAR_COST_CONTRIBUTION |
16 |
& sigmaRbarfile, |
17 |
#endif |
18 |
& sstbarfile, |
19 |
& psbarfile, |
20 |
& bpbarfile, |
21 |
& iestaubarfile, |
22 |
& ubarfile, |
23 |
& vbarfile, |
24 |
& wbarfile, |
25 |
& tauxbarfile, |
26 |
& tauybarfile, |
27 |
& hfluxmeanbarfile, |
28 |
& sfluxmeanbarfile, |
29 |
& costTranspDataFile |
30 |
character*(MAX_LEN_FNAM) tbarfile |
31 |
character*(MAX_LEN_FNAM) sbarfile |
32 |
#ifdef ALLOW_SIGMAR_COST_CONTRIBUTION |
33 |
character*(MAX_LEN_FNAM) sigmaRbarfile |
34 |
#endif |
35 |
character*(MAX_LEN_FNAM) sstbarfile |
36 |
character*(MAX_LEN_FNAM) psbarfile |
37 |
character*(MAX_LEN_FNAM) bpbarfile |
38 |
character*(MAX_LEN_FNAM) iestaubarfile |
39 |
character*(MAX_LEN_FNAM) ubarfile |
40 |
character*(MAX_LEN_FNAM) vbarfile |
41 |
character*(MAX_LEN_FNAM) wbarfile |
42 |
character*(MAX_LEN_FNAM) tauxbarfile |
43 |
character*(MAX_LEN_FNAM) tauybarfile |
44 |
character*(MAX_LEN_FNAM) hfluxmeanbarfile |
45 |
character*(MAX_LEN_FNAM) sfluxmeanbarfile |
46 |
character*(MAX_LEN_FNAM) costTranspDataFile |
47 |
|
48 |
#ifdef ALLOW_TRANSPORT_COST_CONTRIBUTION |
49 |
common /averages_transp_r/ |
50 |
& transpbar |
51 |
& , transpobs |
52 |
& , wtransp |
53 |
_RL transpbar(maxNumDays,nsx,nsy) |
54 |
_RL transpobs(maxNumDays) |
55 |
_RL wtransp(maxNumDays) |
56 |
#endif |
57 |
|
58 |
common /ecco_cost_aux_r/ |
59 |
& mult_hflux, |
60 |
& mult_sflux, |
61 |
& mult_hfluxmm, |
62 |
& mult_sfluxmm, |
63 |
& mult_tauu, |
64 |
& mult_tauv, |
65 |
& mult_hmean, |
66 |
& mult_h, |
67 |
& mult_tp, |
68 |
& mult_ers, |
69 |
& mult_gfo, |
70 |
& mult_sshv4cost, |
71 |
#ifdef ALLOW_SIGMAR_COST_CONTRIBUTION |
72 |
& mult_sigmaR, |
73 |
#endif |
74 |
& mult_temp, |
75 |
& mult_salt, |
76 |
& mult_temp0, |
77 |
& mult_salt0, |
78 |
& mult_etan0, |
79 |
& mult_uvel0, |
80 |
& mult_vvel0, |
81 |
& mult_sst, |
82 |
& mult_tmi, |
83 |
& mult_sss, |
84 |
& mult_bp, |
85 |
& mult_ies, |
86 |
& mult_ctdt, |
87 |
& mult_ctds, |
88 |
& mult_ctdtclim, |
89 |
& mult_ctdsclim, |
90 |
& mult_xbt, |
91 |
& mult_argot, |
92 |
& mult_argos, |
93 |
& mult_usercost, |
94 |
& mult_drift, |
95 |
& mult_tdrift, |
96 |
& mult_sdrift, |
97 |
& mult_wdrift, |
98 |
& mult_scatx, |
99 |
& mult_scaty, |
100 |
& mult_atemp, |
101 |
& mult_aqh, |
102 |
& mult_precip, |
103 |
& mult_swflux, |
104 |
& mult_swdown, |
105 |
& mult_snowprecip, |
106 |
& mult_lwflux, |
107 |
& mult_lwdown, |
108 |
& mult_evap, |
109 |
& mult_apressure, |
110 |
& mult_runoff, |
111 |
& mult_uwind, |
112 |
& mult_vwind, |
113 |
& mult_curmtr, |
114 |
& mult_kapgm, |
115 |
& mult_kapredi, |
116 |
& mult_diffkr, |
117 |
& mult_ini_fin, |
118 |
& mult_edtau, |
119 |
& mult_bottomdrag, |
120 |
& mult_smooth_ic, |
121 |
& mult_smooth_bc, |
122 |
& mult_transp |
123 |
_RL mult_hflux |
124 |
_RL mult_sflux |
125 |
_RL mult_hfluxmm |
126 |
_RL mult_sfluxmm |
127 |
_RL mult_tauu |
128 |
_RL mult_tauv |
129 |
_RL mult_hmean |
130 |
_RL mult_h |
131 |
_RL mult_tp |
132 |
_RL mult_ers |
133 |
_RL mult_gfo |
134 |
_RL mult_sshv4cost(NSSHV4COST) |
135 |
#ifdef ALLOW_SIGMAR_COST_CONTRIBUTION |
136 |
_RL mult_sigmaR |
137 |
#endif |
138 |
_RL mult_temp |
139 |
_RL mult_salt |
140 |
_RL mult_temp0 |
141 |
_RL mult_salt0 |
142 |
_RL mult_etan0 |
143 |
_RL mult_uvel0 |
144 |
_RL mult_vvel0 |
145 |
_RL mult_sst |
146 |
_RL mult_tmi |
147 |
_RL mult_sss |
148 |
_RL mult_bp |
149 |
_RL mult_ies |
150 |
_RL mult_ctdt |
151 |
_RL mult_ctds |
152 |
_RL mult_ctdtclim |
153 |
_RL mult_ctdsclim |
154 |
_RL mult_xbt |
155 |
_RL mult_argot |
156 |
_RL mult_argos |
157 |
_RL mult_usercost(NUSERCOST) |
158 |
_RL mult_drift |
159 |
_RL mult_tdrift |
160 |
_RL mult_sdrift |
161 |
_RL mult_wdrift |
162 |
_RL mult_scatx |
163 |
_RL mult_scaty |
164 |
_RL mult_atemp |
165 |
_RL mult_aqh |
166 |
_RL mult_precip |
167 |
_RL mult_swflux |
168 |
_RL mult_swdown |
169 |
_RL mult_snowprecip |
170 |
_RL mult_lwflux |
171 |
_RL mult_lwdown |
172 |
_RL mult_evap |
173 |
_RL mult_apressure |
174 |
_RL mult_runoff |
175 |
_RL mult_uwind |
176 |
_RL mult_vwind |
177 |
_RL mult_curmtr |
178 |
_RL mult_kapgm |
179 |
_RL mult_kapredi |
180 |
_RL mult_diffkr |
181 |
_RL mult_ini_fin |
182 |
_RL mult_edtau |
183 |
_RL mult_bottomdrag |
184 |
_RL mult_smooth_ic |
185 |
_RL mult_smooth_bc |
186 |
_RL mult_transp |
187 |
|
188 |
|
189 |
common /ecco_cost_c/ |
190 |
& hflux_errfile, |
191 |
& hfluxm_errfile, |
192 |
& sflux_errfile, |
193 |
& sfluxm_errfile, |
194 |
& tauu_errfile, |
195 |
& tauum_errfile, |
196 |
& tauv_errfile, |
197 |
& tauvm_errfile, |
198 |
& scatx_errfile, |
199 |
& scaty_errfile, |
200 |
& data_errfile, |
201 |
& geoid_errfile, |
202 |
& geoid_covariancefile, |
203 |
& ssh_errfile, |
204 |
& tp_errfile, |
205 |
& ers_errfile, |
206 |
& gfo_errfile, |
207 |
& sshv4cost_scalefile, |
208 |
& sshv4cost_errfile, |
209 |
& ctdt_errfile, |
210 |
& ctds_errfile, |
211 |
& drift_errfile, |
212 |
& udrifterrfile, |
213 |
& vdrifterrfile, |
214 |
#ifdef ALLOW_SIGMAR_COST_CONTRIBUTION |
215 |
& sigmaRerrfile, |
216 |
#endif |
217 |
& salterrfile, |
218 |
& temperrfile, |
219 |
& velerrfile, |
220 |
& salt0errfile, |
221 |
& temp0errfile, |
222 |
& etan0errfile, |
223 |
& uvel0errfile, |
224 |
& vvel0errfile, |
225 |
& vel0errfile, |
226 |
& ssterrfile, |
227 |
& ssserrfile, |
228 |
& bperrfile, |
229 |
& ieserrfile, |
230 |
& atemp_errfile, |
231 |
& aqh_errfile, |
232 |
& precip_errfile, |
233 |
& swflux_errfile, |
234 |
& swdown_errfile, |
235 |
& snowprecip_errfile, |
236 |
& lwflux_errfile, |
237 |
& lwdown_errfile, |
238 |
& evap_errfile, |
239 |
& apressure_errfile, |
240 |
& runoff_errfile, |
241 |
& edtau_errfile, |
242 |
& kapgm_errfile, |
243 |
& kapredi_errfile, |
244 |
& diffkr_errfile, |
245 |
& bottomdrag_errfile, |
246 |
& usercost_errfile, |
247 |
& uwind_errfile, |
248 |
& vwind_errfile |
249 |
character*(MAX_LEN_FNAM) hflux_errfile |
250 |
character*(MAX_LEN_FNAM) sflux_errfile |
251 |
character*(MAX_LEN_FNAM) tauu_errfile |
252 |
character*(MAX_LEN_FNAM) tauv_errfile |
253 |
character*(MAX_LEN_FNAM) hfluxm_errfile |
254 |
character*(MAX_LEN_FNAM) sfluxm_errfile |
255 |
character*(MAX_LEN_FNAM) tauum_errfile |
256 |
character*(MAX_LEN_FNAM) tauvm_errfile |
257 |
character*(MAX_LEN_FNAM) scatx_errfile |
258 |
character*(MAX_LEN_FNAM) scaty_errfile |
259 |
character*(MAX_LEN_FNAM) data_errfile |
260 |
character*(MAX_LEN_FNAM) geoid_errfile |
261 |
character*(MAX_LEN_FNAM) geoid_covariancefile |
262 |
character*(MAX_LEN_FNAM) ssh_errfile |
263 |
character*(MAX_LEN_FNAM) tp_errfile |
264 |
character*(MAX_LEN_FNAM) ers_errfile |
265 |
character*(MAX_LEN_FNAM) gfo_errfile |
266 |
character*(MAX_LEN_FNAM) sshv4cost_scalefile(NSSHV4COST) |
267 |
character*(MAX_LEN_FNAM) sshv4cost_errfile(NSSHV4COST) |
268 |
character*(MAX_LEN_FNAM) ctdt_errfile |
269 |
character*(MAX_LEN_FNAM) ctds_errfile |
270 |
character*(MAX_LEN_FNAM) drift_errfile |
271 |
character*(MAX_LEN_FNAM) udrifterrfile |
272 |
character*(MAX_LEN_FNAM) vdrifterrfile |
273 |
#ifdef ALLOW_SIGMAR_COST_CONTRIBUTION |
274 |
character*(MAX_LEN_FNAM) sigmaRerrfile |
275 |
#endif |
276 |
character*(MAX_LEN_FNAM) salterrfile |
277 |
character*(MAX_LEN_FNAM) temperrfile |
278 |
character*(MAX_LEN_FNAM) velerrfile |
279 |
character*(MAX_LEN_FNAM) salt0errfile |
280 |
character*(MAX_LEN_FNAM) temp0errfile |
281 |
character*(MAX_LEN_FNAM) etan0errfile |
282 |
character*(MAX_LEN_FNAM) uvel0errfile |
283 |
character*(MAX_LEN_FNAM) vvel0errfile |
284 |
character*(MAX_LEN_FNAM) vel0errfile |
285 |
character*(MAX_LEN_FNAM) ssterrfile |
286 |
character*(MAX_LEN_FNAM) ssserrfile |
287 |
character*(MAX_LEN_FNAM) bperrfile |
288 |
character*(MAX_LEN_FNAM) ieserrfile |
289 |
character*(MAX_LEN_FNAM) atemp_errfile |
290 |
character*(MAX_LEN_FNAM) aqh_errfile |
291 |
character*(MAX_LEN_FNAM) precip_errfile |
292 |
character*(MAX_LEN_FNAM) swflux_errfile |
293 |
character*(MAX_LEN_FNAM) swdown_errfile |
294 |
character*(MAX_LEN_FNAM) snowprecip_errfile |
295 |
character*(MAX_LEN_FNAM) lwflux_errfile |
296 |
character*(MAX_LEN_FNAM) lwdown_errfile |
297 |
character*(MAX_LEN_FNAM) evap_errfile |
298 |
character*(MAX_LEN_FNAM) apressure_errfile |
299 |
character*(MAX_LEN_FNAM) runoff_errfile |
300 |
character*(MAX_LEN_FNAM) edtau_errfile |
301 |
character*(MAX_LEN_FNAM) kapgm_errfile |
302 |
character*(MAX_LEN_FNAM) kapredi_errfile |
303 |
character*(MAX_LEN_FNAM) diffkr_errfile |
304 |
character*(MAX_LEN_FNAM) bottomdrag_errfile |
305 |
character*(MAX_LEN_FNAM) usercost_errfile(NUSERCOST) |
306 |
character*(MAX_LEN_FNAM) uwind_errfile |
307 |
character*(MAX_LEN_FNAM) vwind_errfile |
308 |
|
309 |
common /ecco_cost_weights_0_r/ |
310 |
& whflux0, wsflux0, wtau0, |
311 |
& watemp0, waqh0, wprecip0, wsnowprecip0, wwind0, |
312 |
& wswflux0, wswdown0, wlwflux0, wlwdown0, |
313 |
& wevap0, wapressure0, wrunoff0, wkapredi0, |
314 |
& wbottomdrag0,wdiffkr0, wkapgm0, wedtau0 |
315 |
_RL whflux0 |
316 |
_RL wsflux0 |
317 |
_RL wtau0 |
318 |
_RL watemp0 |
319 |
_RL waqh0 |
320 |
_RL wprecip0 |
321 |
_RL wswflux0 |
322 |
_RL wswdown0 |
323 |
_RL wsnowprecip0 |
324 |
_RL wlwflux0 |
325 |
_RL wlwdown0 |
326 |
_RL wevap0 |
327 |
_RL wapressure0 |
328 |
_RL wrunoff0 |
329 |
_RL wbottomdrag0 |
330 |
_RL wwind0 |
331 |
_RL wdiffkr0 |
332 |
_RL wkapgm0 |
333 |
_RL wkapredi0 |
334 |
_RL wedtau0 |
335 |
|
336 |
common /ecco_cost_weights_mean_r/ |
337 |
& wmean_hflux, wmean_sflux, wmean_tau, |
338 |
& wmean_atemp, wmean_aqh, |
339 |
& wmean_precip, wmean_snowprecip, wmean_wind, |
340 |
& wmean_swflux, wmean_swdown, wmean_lwflux, wmean_lwdown, |
341 |
& wmean_evap, wmean_apressure, wmean_runoff |
342 |
_RL wmean_hflux |
343 |
_RL wmean_sflux |
344 |
_RL wmean_tau |
345 |
_RL wmean_atemp |
346 |
_RL wmean_aqh |
347 |
_RL wmean_precip |
348 |
_RL wmean_swflux |
349 |
_RL wmean_swdown |
350 |
_RL wmean_snowprecip |
351 |
_RL wmean_lwflux |
352 |
_RL wmean_lwdown |
353 |
_RL wmean_evap |
354 |
_RL wmean_apressure |
355 |
_RL wmean_runoff |
356 |
_RL wmean_wind |
357 |
|
358 |
common /ecco_cost_data_c/ |
359 |
#ifdef ALLOW_SIGMAR_COST_CONTRIBUTION |
360 |
& sigmaRdatfile, |
361 |
#endif |
362 |
& tdatfile, |
363 |
& sdatfile, |
364 |
& scatxdatfile, |
365 |
& scatydatfile, |
366 |
& sstdatfile, |
367 |
& tmidatfile, |
368 |
& sssdatfile, |
369 |
& bpdatfile, |
370 |
& iesdatfile, |
371 |
& mdtdatfile, |
372 |
& topexfile, |
373 |
& ersfile, |
374 |
& gfofile, |
375 |
& ctdtfile, |
376 |
& ctdsfile, |
377 |
& ctdtclimfile, |
378 |
& ctdsclimfile, |
379 |
& xbtfile, |
380 |
& argotfile, |
381 |
& argosfile, |
382 |
& udriftfile, |
383 |
& vdriftfile, |
384 |
& usercost_datafile, |
385 |
& curmtrufile, |
386 |
& curmtrvfile |
387 |
|
388 |
#ifdef ALLOW_SIGMAR_COST_CONTRIBUTION |
389 |
character*(MAX_LEN_FNAM) sigmaRdatfile |
390 |
#endif |
391 |
character*(MAX_LEN_FNAM) tdatfile |
392 |
character*(MAX_LEN_FNAM) sdatfile |
393 |
character*(MAX_LEN_FNAM) scatxdatfile |
394 |
character*(MAX_LEN_FNAM) scatydatfile |
395 |
character*(MAX_LEN_FNAM) sstdatfile |
396 |
character*(MAX_LEN_FNAM) tmidatfile |
397 |
character*(MAX_LEN_FNAM) sssdatfile |
398 |
character*(MAX_LEN_FNAM) bpdatfile |
399 |
character*(MAX_LEN_FNAM) iesdatfile |
400 |
character*(MAX_LEN_FNAM) mdtdatfile |
401 |
character*(MAX_LEN_FNAM) topexfile |
402 |
character*(MAX_LEN_FNAM) ersfile |
403 |
character*(MAX_LEN_FNAM) gfofile |
404 |
character*(MAX_LEN_FNAM) ctdtfile |
405 |
character*(MAX_LEN_FNAM) ctdsfile |
406 |
character*(MAX_LEN_FNAM) ctdtclimfile |
407 |
character*(MAX_LEN_FNAM) ctdsclimfile |
408 |
character*(MAX_LEN_FNAM) xbtfile |
409 |
character*(MAX_LEN_FNAM) argotfile |
410 |
character*(MAX_LEN_FNAM) argosfile |
411 |
character*(MAX_LEN_FNAM) argofile |
412 |
character*(MAX_LEN_FNAM) usercost_datafile(NUSERCOST) |
413 |
character*(MAX_LEN_FNAM) udriftfile |
414 |
character*(MAX_LEN_FNAM) vdriftfile |
415 |
character*(MAX_LEN_FNAM) curmtrufile |
416 |
character*(MAX_LEN_FNAM) curmtrvfile |
417 |
|
418 |
common /ecco_cost_data_times_i/ |
419 |
& scatxstartdate, |
420 |
& scatystartdate, |
421 |
& sststartdate, |
422 |
& argotstartdate, |
423 |
& argosstartdate, |
424 |
& tmistartdate, |
425 |
& sssstartdate, |
426 |
& bpstartdate, |
427 |
& iesstartdate, |
428 |
& topexstartdate, |
429 |
& ersstartdate, |
430 |
& gfostartdate, |
431 |
& mdtstartdate, |
432 |
& mdtenddate |
433 |
integer scatxstartdate(4) |
434 |
integer scatystartdate(4) |
435 |
integer sststartdate(4) |
436 |
integer argotstartdate(4) |
437 |
integer argosstartdate(4) |
438 |
integer tmistartdate(4) |
439 |
integer sssstartdate(4) |
440 |
integer bpstartdate(4) |
441 |
integer iesstartdate(4) |
442 |
integer topexstartdate(4) |
443 |
integer ersstartdate(4) |
444 |
integer gfostartdate(4) |
445 |
integer mdtstartdate(4) |
446 |
integer mdtenddate(4) |
447 |
|
448 |
common /ecco_cost_data_aux_i/ |
449 |
& tmistartdate1, |
450 |
& tmistartdate2, |
451 |
& sststartdate1, |
452 |
& sststartdate2, |
453 |
& sssstartdate1, |
454 |
& sssstartdate2, |
455 |
& bpstartdate1, |
456 |
& bpstartdate2, |
457 |
& iesstartdate1, |
458 |
& iesstartdate2, |
459 |
& argotstartdate1, |
460 |
& argotstartdate2, |
461 |
& argosstartdate1, |
462 |
& argosstartdate2, |
463 |
& topexstartdate1, |
464 |
& topexstartdate2, |
465 |
& ersstartdate1, |
466 |
& ersstartdate2, |
467 |
& gfostartdate1, |
468 |
& gfostartdate2, |
469 |
& scatstartdate1, |
470 |
& scatstartdate2, |
471 |
& mdtstartdate1, |
472 |
& mdtstartdate2, |
473 |
& mdtenddate1, |
474 |
& mdtenddate2 |
475 |
|
476 |
integer tmistartdate1 |
477 |
integer tmistartdate2 |
478 |
integer sststartdate1 |
479 |
integer sststartdate2 |
480 |
integer sssstartdate1 |
481 |
integer sssstartdate2 |
482 |
integer bpstartdate1 |
483 |
integer bpstartdate2 |
484 |
integer iesstartdate1 |
485 |
integer iesstartdate2 |
486 |
integer argotstartdate1 |
487 |
integer argotstartdate2 |
488 |
integer argosstartdate1 |
489 |
integer argosstartdate2 |
490 |
integer topexstartdate1 |
491 |
integer topexstartdate2 |
492 |
integer ersstartdate1 |
493 |
integer ersstartdate2 |
494 |
integer gfostartdate1 |
495 |
integer gfostartdate2 |
496 |
integer scatstartdate1 |
497 |
integer scatstartdate2 |
498 |
integer mdtstartdate1 |
499 |
integer mdtstartdate2 |
500 |
integer mdtenddate1 |
501 |
integer mdtenddate2 |
502 |
|
503 |
common /ecco_cost_data_times_r/ |
504 |
& topexperiod, |
505 |
& ersperiod, |
506 |
& gfoperiod, |
507 |
& scatperiod |
508 |
_RL topexperiod |
509 |
_RL ersperiod |
510 |
_RL gfoperiod |
511 |
_RL scatperiod |
512 |
|
513 |
common /ecco_cost_data_detrend/ |
514 |
& topexintercept, |
515 |
& ersintercept, |
516 |
& gfointercept, |
517 |
& topexslope, |
518 |
& ersslope, |
519 |
& gfoslope |
520 |
_RL topexintercept |
521 |
_RL ersintercept |
522 |
_RL gfointercept |
523 |
_RL topexslope |
524 |
_RL ersslope |
525 |
_RL gfoslope |
526 |
|
527 |
common /ecco_cost_errfactor/ |
528 |
& sshv4cost_errfactor |
529 |
_RL sshv4cost_errfactor(NSSHV4COST) |
530 |
|
531 |
common /ecco_ssh_daymask_c/ |
532 |
& tpTimeMaskFile, ersTimeMaskFile, gfoTimeMaskFile |
533 |
character*(MAX_LEN_FNAM) tpTimeMaskFile |
534 |
character*(MAX_LEN_FNAM) ersTimeMaskFile |
535 |
character*(MAX_LEN_FNAM) gfoTimeMaskFile |
536 |
|
537 |
c ================================================================== |
538 |
c END OF HEADER ECCO_legacy |
539 |
c ================================================================== |