8 |
c |
c |
9 |
c started: Christian Eckert eckert@mit.edu 30-Jun-1999 |
c started: Christian Eckert eckert@mit.edu 30-Jun-1999 |
10 |
c |
c |
|
c added sea-ice term: menemenlis@jpl.nasa.gov 26-Feb-2003 |
|
|
c |
|
11 |
c ================================================================== |
c ================================================================== |
12 |
c HEADER AVERAGES |
c HEADER AVERAGES |
13 |
c ================================================================== |
c ================================================================== |
21 |
c monrec - number of averaged theta and salinity records. |
c monrec - number of averaged theta and salinity records. |
22 |
|
|
23 |
common /average_i/ |
common /average_i/ |
24 |
& sum1day,sum1mon, |
& sum1day,sum1mon,sum1year, |
25 |
& dayrec,monrec |
& dayrec,monrec,yearrec |
26 |
integer sum1day |
integer sum1day |
27 |
integer sum1mon |
integer sum1mon |
28 |
|
integer sum1year |
29 |
integer dayrec |
integer dayrec |
30 |
integer monrec |
integer monrec |
31 |
|
integer yearrec |
32 |
|
|
33 |
|
c Number of Generic Cost terms: |
34 |
|
c ============================= |
35 |
|
INTEGER NGENCOST |
36 |
|
PARAMETER ( NGENCOST=10 ) |
37 |
|
|
38 |
c Averaged Fields: |
c Averaged Fields: |
39 |
c ================ |
c ================ |
192 |
& tauxbarfile, |
& tauxbarfile, |
193 |
& tauybarfile, |
& tauybarfile, |
194 |
& hfluxbarfile, |
& hfluxbarfile, |
195 |
& sfluxbarfile |
& sfluxbarfile, |
196 |
|
& costTranspDataFile |
197 |
character*(MAX_LEN_FNAM) tbarfile |
character*(MAX_LEN_FNAM) tbarfile |
198 |
character*(MAX_LEN_FNAM) sbarfile |
character*(MAX_LEN_FNAM) sbarfile |
199 |
character*(MAX_LEN_FNAM) psbarfile |
character*(MAX_LEN_FNAM) psbarfile |
204 |
character*(MAX_LEN_FNAM) tauybarfile |
character*(MAX_LEN_FNAM) tauybarfile |
205 |
character*(MAX_LEN_FNAM) hfluxbarfile |
character*(MAX_LEN_FNAM) hfluxbarfile |
206 |
character*(MAX_LEN_FNAM) sfluxbarfile |
character*(MAX_LEN_FNAM) sfluxbarfile |
207 |
|
character*(MAX_LEN_FNAM) costTranspDataFile |
208 |
|
|
209 |
|
#ifdef ALLOW_TRANSPORT_COST_CONTRIBUTION |
210 |
|
INTEGER maxNumDays |
211 |
|
PARAMETER ( maxNumDays = 5480 ) |
212 |
|
common /averages_transp_r/ |
213 |
|
& transpbar |
214 |
|
& , transpobs |
215 |
|
& , wtransp |
216 |
|
_RL transpbar(maxNumDays,nsx,nsy) |
217 |
|
_RL transpobs(maxNumDays) |
218 |
|
_RL wtransp(maxNumDays) |
219 |
|
#endif |
220 |
|
|
221 |
c file precision and field type |
c file precision and field type |
222 |
|
|
274 |
c objf_sfluxmm - Global time-mean Salt flux contribution. |
c objf_sfluxmm - Global time-mean Salt flux contribution. |
275 |
c objf_hmean - Mean sea surface height contribution. |
c objf_hmean - Mean sea surface height contribution. |
276 |
c objf_h - Residual sea surface height contribution. |
c objf_h - Residual sea surface height contribution. |
277 |
|
c objf_tp - Residual sea surface height contribution from T/P |
278 |
|
c objf_ers - Residual sea surface height contribution from T/P |
279 |
|
c objf_gfo - Residual sea surface height contribution from T/P |
280 |
c objf_temp - Temperature contribution. |
c objf_temp - Temperature contribution. |
281 |
c objf_salt - Salinity contribution. |
c objf_salt - Salinity contribution. |
282 |
c objf_temp0 - Initial conditions Temperature contribution. |
c objf_temp0 - Initial conditions Temperature contribution. |
291 |
c objf_xbt - XBT temperature data |
c objf_xbt - XBT temperature data |
292 |
c objf_argot - ARGO temperature profiles |
c objf_argot - ARGO temperature profiles |
293 |
c objf_argos - ARGO salt profiles |
c objf_argos - ARGO salt profiles |
294 |
|
c objf_gencost - gencost user defined contribution |
295 |
c objf_scatxm - time-mean zonal SCAT contribution |
c objf_scatxm - time-mean zonal SCAT contribution |
296 |
c objf_scatym - time-mean meridional SCAT contribution |
c objf_scatym - time-mean meridional SCAT contribution |
297 |
c objf_scatx - zonal SCAT contribution |
c objf_scatx - zonal SCAT contribution |
298 |
c objf_scaty - meridional SCAT contribution |
c objf_scaty - meridional SCAT contribution |
|
c objf_ice - sea-ice volume |
|
299 |
c objf_kapgm - kappa GM contribution |
c objf_kapgm - kappa GM contribution |
300 |
c objf_diffkr - diffusion contribution |
c objf_diffkr - diffusion contribution |
301 |
c objf_theta_ini_fin - final vs. initial theta misfit |
c objf_theta_ini_fin - final vs. initial theta misfit |
302 |
c objf_salt_ini_fin - final vs. initial salt misfit |
c objf_salt_ini_fin - final vs. initial salt misfit |
303 |
|
c objf_eddytau - eddy stress contribution |
304 |
|
c objf_bottomdrag - bottom drag contribution |
305 |
c |
c |
306 |
c mult_"var" - multipliers for the individual cost |
c mult_"var" - multipliers for the individual cost |
307 |
c function contributions. |
c function contributions. |
310 |
& objf_hflux, |
& objf_hflux, |
311 |
& objf_hfluxm, |
& objf_hfluxm, |
312 |
& objf_hfluxmm, |
& objf_hfluxmm, |
313 |
|
& objf_hfluxsmoo, |
314 |
& objf_sflux, |
& objf_sflux, |
315 |
& objf_sfluxm, |
& objf_sfluxm, |
316 |
& objf_sfluxmm, |
& objf_sfluxmm, |
317 |
|
& objf_sfluxsmoo, |
318 |
& objf_tauu, |
& objf_tauu, |
319 |
& objf_tauum, |
& objf_tauum, |
320 |
|
& objf_tauusmoo, |
321 |
& objf_tauv, |
& objf_tauv, |
322 |
& objf_tauvm, |
& objf_tauvm, |
323 |
|
& objf_tauvsmoo, |
324 |
& objf_hmean, |
& objf_hmean, |
325 |
& objf_h, |
& objf_h, |
326 |
|
& objf_tp, |
327 |
|
& objf_ers, |
328 |
|
& objf_gfo, |
329 |
& objf_temp, |
& objf_temp, |
330 |
& objf_salt, |
& objf_salt, |
331 |
& objf_temp0, |
& objf_temp0, |
332 |
& objf_salt0, |
& objf_salt0, |
333 |
|
& objf_temp0smoo, |
334 |
|
& objf_salt0smoo, |
335 |
& objf_sst, |
& objf_sst, |
336 |
& objf_tmi, |
& objf_tmi, |
337 |
& objf_sss, |
& objf_sss, |
342 |
& objf_xbt, |
& objf_xbt, |
343 |
& objf_argot, |
& objf_argot, |
344 |
& objf_argos, |
& objf_argos, |
345 |
|
& objf_gencost, |
346 |
& objf_drift, |
& objf_drift, |
347 |
& objf_tdrift, |
& objf_tdrift, |
348 |
& objf_sdrift, |
& objf_sdrift, |
355 |
& objf_aqh, |
& objf_aqh, |
356 |
& objf_precip, |
& objf_precip, |
357 |
& objf_swflux, |
& objf_swflux, |
358 |
|
& objf_swdown, |
359 |
|
& objf_snowprecip, |
360 |
|
& objf_lwflux, |
361 |
|
& objf_lwdown, |
362 |
|
& objf_evap, |
363 |
|
& objf_apressure, |
364 |
|
& objf_runoff, |
365 |
& objf_uwind, |
& objf_uwind, |
366 |
& objf_vwind, |
& objf_vwind, |
367 |
& objf_atempm, |
& objf_atempm, |
368 |
& objf_aqhm, |
& objf_aqhm, |
369 |
& objf_precipm, |
& objf_precipm, |
370 |
& objf_swfluxm, |
& objf_swfluxm, |
371 |
|
& objf_swdownm, |
372 |
|
& objf_snowprecipm, |
373 |
|
& objf_lwfluxm, |
374 |
|
& objf_lwdownm, |
375 |
|
& objf_evapm, |
376 |
|
& objf_apressurem, |
377 |
|
& objf_runoffm, |
378 |
& objf_uwindm, |
& objf_uwindm, |
379 |
& objf_vwindm, |
& objf_vwindm, |
380 |
|
& objf_atempsmoo, |
381 |
|
& objf_aqhsmoo, |
382 |
|
& objf_precipsmoo, |
383 |
|
& objf_swfluxsmoo, |
384 |
|
& objf_swdownsmoo, |
385 |
|
& objf_snowprecipsmoo, |
386 |
|
& objf_lwfluxsmoo, |
387 |
|
& objf_lwdownsmoo, |
388 |
|
& objf_evapsmoo, |
389 |
|
& objf_apressuresmoo, |
390 |
|
& objf_runoffsmoo, |
391 |
|
& objf_uwindsmoo, |
392 |
|
& objf_vwindsmoo, |
393 |
& objf_obcsn, |
& objf_obcsn, |
394 |
& objf_obcss, |
& objf_obcss, |
395 |
& objf_obcsw, |
& objf_obcsw, |
397 |
& objf_obcsvol, |
& objf_obcsvol, |
398 |
& objf_curmtr, |
& objf_curmtr, |
399 |
& objf_ageos, |
& objf_ageos, |
400 |
& objf_ice, |
& objf_kapgm, |
401 |
& objf_diffkr, |
& objf_diffkr, |
402 |
& objf_theta_ini_fin, |
& objf_theta_ini_fin, |
403 |
& objf_salt_ini_fin, |
& objf_salt_ini_fin, |
404 |
& objf_eddytau |
& objf_eddytau, |
405 |
|
& objf_bottomdrag, |
406 |
|
& objf_transp |
407 |
_RL objf_hflux (nsx,nsy) |
_RL objf_hflux (nsx,nsy) |
408 |
_RL objf_hfluxm (nsx,nsy) |
_RL objf_hfluxm (nsx,nsy) |
409 |
_RL objf_hfluxmm(nsx,nsy) |
_RL objf_hfluxmm(nsx,nsy) |
410 |
|
_RL objf_hfluxsmoo (nsx,nsy) |
411 |
_RL objf_sflux (nsx,nsy) |
_RL objf_sflux (nsx,nsy) |
412 |
_RL objf_sfluxm (nsx,nsy) |
_RL objf_sfluxm (nsx,nsy) |
413 |
_RL objf_sfluxmm(nsx,nsy) |
_RL objf_sfluxmm |
414 |
|
_RL objf_sfluxsmoo (nsx,nsy) |
415 |
_RL objf_tauu (nsx,nsy) |
_RL objf_tauu (nsx,nsy) |
416 |
_RL objf_tauum (nsx,nsy) |
_RL objf_tauum (nsx,nsy) |
417 |
|
_RL objf_tauusmoo (nsx,nsy) |
418 |
_RL objf_tauv (nsx,nsy) |
_RL objf_tauv (nsx,nsy) |
419 |
_RL objf_tauvm (nsx,nsy) |
_RL objf_tauvm (nsx,nsy) |
420 |
|
_RL objf_tauvsmoo (nsx,nsy) |
421 |
_RL objf_hmean |
_RL objf_hmean |
422 |
_RL objf_h (nsx,nsy) |
_RL objf_h (nsx,nsy) |
423 |
|
_RL objf_tp (nsx,nsy) |
424 |
|
_RL objf_ers (nsx,nsy) |
425 |
|
_RL objf_gfo (nsx,nsy) |
426 |
_RL objf_temp (nsx,nsy) |
_RL objf_temp (nsx,nsy) |
427 |
_RL objf_salt (nsx,nsy) |
_RL objf_salt (nsx,nsy) |
428 |
_RL objf_temp0(nsx,nsy) |
_RL objf_temp0(nsx,nsy) |
429 |
_RL objf_salt0(nsx,nsy) |
_RL objf_salt0(nsx,nsy) |
430 |
|
_RL objf_temp0smoo(nsx,nsy) |
431 |
|
_RL objf_salt0smoo(nsx,nsy) |
432 |
_RL objf_sst (nsx,nsy) |
_RL objf_sst (nsx,nsy) |
433 |
_RL objf_tmi (nsx,nsy) |
_RL objf_tmi (nsx,nsy) |
434 |
_RL objf_sss (nsx,nsy) |
_RL objf_sss (nsx,nsy) |
439 |
_RL objf_xbt (nsx,nsy) |
_RL objf_xbt (nsx,nsy) |
440 |
_RL objf_argot(nsx,nsy) |
_RL objf_argot(nsx,nsy) |
441 |
_RL objf_argos(nsx,nsy) |
_RL objf_argos(nsx,nsy) |
442 |
|
_RL objf_gencost(NGENCOST,nsx,nsy) |
443 |
_RL objf_drift(nsx,nsy) |
_RL objf_drift(nsx,nsy) |
444 |
_RL objf_tdrift(nsx,nsy) |
_RL objf_tdrift(nsx,nsy) |
445 |
_RL objf_sdrift(nsx,nsy) |
_RL objf_sdrift(nsx,nsy) |
452 |
_RL objf_aqh (nsx,nsy) |
_RL objf_aqh (nsx,nsy) |
453 |
_RL objf_precip(nsx,nsy) |
_RL objf_precip(nsx,nsy) |
454 |
_RL objf_swflux(nsx,nsy) |
_RL objf_swflux(nsx,nsy) |
455 |
|
_RL objf_swdown(nsx,nsy) |
456 |
|
_RL objf_snowprecip(nsx,nsy) |
457 |
|
_RL objf_lwflux(nsx,nsy) |
458 |
|
_RL objf_lwdown(nsx,nsy) |
459 |
|
_RL objf_evap(nsx,nsy) |
460 |
|
_RL objf_apressure(nsx,nsy) |
461 |
|
_RL objf_runoff(nsx,nsy) |
462 |
_RL objf_uwind(nsx,nsy) |
_RL objf_uwind(nsx,nsy) |
463 |
_RL objf_vwind(nsx,nsy) |
_RL objf_vwind(nsx,nsy) |
464 |
_RL objf_atempm(nsx,nsy) |
_RL objf_atempm(nsx,nsy) |
465 |
_RL objf_aqhm (nsx,nsy) |
_RL objf_aqhm (nsx,nsy) |
466 |
_RL objf_precipm(nsx,nsy) |
_RL objf_precipm(nsx,nsy) |
467 |
_RL objf_swfluxm(nsx,nsy) |
_RL objf_swfluxm(nsx,nsy) |
468 |
|
_RL objf_swdownm(nsx,nsy) |
469 |
|
_RL objf_snowprecipm(nsx,nsy) |
470 |
|
_RL objf_lwfluxm(nsx,nsy) |
471 |
|
_RL objf_lwdownm(nsx,nsy) |
472 |
|
_RL objf_evapm(nsx,nsy) |
473 |
|
_RL objf_apressurem(nsx,nsy) |
474 |
|
_RL objf_runoffm(nsx,nsy) |
475 |
_RL objf_uwindm(nsx,nsy) |
_RL objf_uwindm(nsx,nsy) |
476 |
_RL objf_vwindm(nsx,nsy) |
_RL objf_vwindm(nsx,nsy) |
477 |
|
_RL objf_atempsmoo(nsx,nsy) |
478 |
|
_RL objf_aqhsmoo (nsx,nsy) |
479 |
|
_RL objf_precipsmoo(nsx,nsy) |
480 |
|
_RL objf_swfluxsmoo(nsx,nsy) |
481 |
|
_RL objf_swdownsmoo(nsx,nsy) |
482 |
|
_RL objf_snowprecipsmoo(nsx,nsy) |
483 |
|
_RL objf_lwfluxsmoo(nsx,nsy) |
484 |
|
_RL objf_lwdownsmoo(nsx,nsy) |
485 |
|
_RL objf_evapsmoo(nsx,nsy) |
486 |
|
_RL objf_apressuresmoo(nsx,nsy) |
487 |
|
_RL objf_runoffsmoo(nsx,nsy) |
488 |
|
_RL objf_uwindsmoo(nsx,nsy) |
489 |
|
_RL objf_vwindsmoo(nsx,nsy) |
490 |
_RL objf_obcsn(nsx,nsy) |
_RL objf_obcsn(nsx,nsy) |
491 |
_RL objf_obcss(nsx,nsy) |
_RL objf_obcss(nsx,nsy) |
492 |
_RL objf_obcsw(nsx,nsy) |
_RL objf_obcsw(nsx,nsy) |
494 |
_RL objf_obcsvol |
_RL objf_obcsvol |
495 |
_RL objf_curmtr(nsx,nsy) |
_RL objf_curmtr(nsx,nsy) |
496 |
_RL objf_ageos(nsx,nsy) |
_RL objf_ageos(nsx,nsy) |
|
_RL objf_ice (nsx,nsy) |
|
497 |
_RL objf_kapgm(nsx,nsy) |
_RL objf_kapgm(nsx,nsy) |
498 |
_RL objf_diffkr(nsx,nsy) |
_RL objf_diffkr(nsx,nsy) |
499 |
_RL objf_theta_ini_fin(nsx,nsy) |
_RL objf_theta_ini_fin(nsx,nsy) |
500 |
_RL objf_salt_ini_fin(nsx,nsy) |
_RL objf_salt_ini_fin(nsx,nsy) |
501 |
_RL objf_eddytau(nsx,nsy) |
_RL objf_eddytau(nsx,nsy) |
502 |
|
_RL objf_bottomdrag(nsx,nsy) |
503 |
|
_RL objf_transp |
504 |
|
|
505 |
common /ecco_cost_num/ |
common /ecco_cost_num/ |
506 |
& num_hflux, |
& num_hflux, |
515 |
& num_tauvm, |
& num_tauvm, |
516 |
& num_hmean, |
& num_hmean, |
517 |
& num_h, |
& num_h, |
518 |
|
& num_tp, |
519 |
|
& num_ers, |
520 |
|
& num_gfo, |
521 |
& num_temp, |
& num_temp, |
522 |
& num_salt, |
& num_salt, |
523 |
& num_temp0, |
& num_temp0, |
532 |
& num_xbt, |
& num_xbt, |
533 |
& num_argot, |
& num_argot, |
534 |
& num_argos, |
& num_argos, |
535 |
|
& num_gencost, |
536 |
& num_drift, |
& num_drift, |
537 |
& num_tdrift, |
& num_tdrift, |
538 |
& num_sdrift, |
& num_sdrift, |
545 |
& num_aqh, |
& num_aqh, |
546 |
& num_precip, |
& num_precip, |
547 |
& num_swflux, |
& num_swflux, |
548 |
|
& num_swdown, |
549 |
|
& num_snowprecip, |
550 |
|
& num_lwflux, |
551 |
|
& num_lwdown, |
552 |
|
& num_evap, |
553 |
|
& num_apressure, |
554 |
|
& num_runoff, |
555 |
& num_uwind, |
& num_uwind, |
556 |
& num_vwind, |
& num_vwind, |
557 |
& num_atempm, |
& num_atempm, |
558 |
& num_aqhm, |
& num_aqhm, |
559 |
& num_precipm, |
& num_precipm, |
560 |
& num_swfluxm, |
& num_swfluxm, |
561 |
|
& num_swdownm, |
562 |
|
& num_snowprecipm, |
563 |
|
& num_lwfluxm, |
564 |
|
& num_lwdownm, |
565 |
|
& num_evapm, |
566 |
|
& num_apressurem, |
567 |
|
& num_runoffm, |
568 |
& num_uwindm, |
& num_uwindm, |
569 |
& num_vwindm, |
& num_vwindm, |
570 |
& num_obcsn, |
& num_obcsn, |
574 |
& num_obcsvol, |
& num_obcsvol, |
575 |
& num_curmtr, |
& num_curmtr, |
576 |
& num_ageos, |
& num_ageos, |
577 |
& num_ice, |
& num_kapgm, |
578 |
& num_diffkr, |
& num_diffkr, |
579 |
& num_theta_ini_fin, |
& num_theta_ini_fin, |
580 |
& num_salt_ini_fin, |
& num_salt_ini_fin, |
581 |
& num_eddytau |
& num_eddytau, |
582 |
|
& num_bottomdrag, |
583 |
|
& num_transp |
584 |
|
|
585 |
_RL num_hflux (nsx,nsy) |
_RL num_hflux (nsx,nsy) |
586 |
_RL num_hfluxm (nsx,nsy) |
_RL num_hfluxm (nsx,nsy) |
594 |
_RL num_tauvm (nsx,nsy) |
_RL num_tauvm (nsx,nsy) |
595 |
_RL num_hmean |
_RL num_hmean |
596 |
_RL num_h (nsx,nsy) |
_RL num_h (nsx,nsy) |
597 |
|
_RL num_tp (nsx,nsy) |
598 |
|
_RL num_ers (nsx,nsy) |
599 |
|
_RL num_gfo (nsx,nsy) |
600 |
_RL num_temp (nsx,nsy) |
_RL num_temp (nsx,nsy) |
601 |
_RL num_salt (nsx,nsy) |
_RL num_salt (nsx,nsy) |
602 |
_RL num_temp0(nsx,nsy) |
_RL num_temp0(nsx,nsy) |
611 |
_RL num_xbt (nsx,nsy) |
_RL num_xbt (nsx,nsy) |
612 |
_RL num_argot(nsx,nsy) |
_RL num_argot(nsx,nsy) |
613 |
_RL num_argos(nsx,nsy) |
_RL num_argos(nsx,nsy) |
614 |
|
_RL num_gencost(NGENCOST,nsx,nsy) |
615 |
_RL num_drift(nsx,nsy) |
_RL num_drift(nsx,nsy) |
616 |
_RL num_tdrift(nsx,nsy) |
_RL num_tdrift(nsx,nsy) |
617 |
_RL num_sdrift(nsx,nsy) |
_RL num_sdrift(nsx,nsy) |
624 |
_RL num_aqh (nsx,nsy) |
_RL num_aqh (nsx,nsy) |
625 |
_RL num_precip(nsx,nsy) |
_RL num_precip(nsx,nsy) |
626 |
_RL num_swflux(nsx,nsy) |
_RL num_swflux(nsx,nsy) |
627 |
|
_RL num_swdown(nsx,nsy) |
628 |
|
_RL num_snowprecip(nsx,nsy) |
629 |
|
_RL num_lwflux(nsx,nsy) |
630 |
|
_RL num_lwdown(nsx,nsy) |
631 |
|
_RL num_evap(nsx,nsy) |
632 |
|
_RL num_apressure(nsx,nsy) |
633 |
|
_RL num_runoff(nsx,nsy) |
634 |
_RL num_uwind(nsx,nsy) |
_RL num_uwind(nsx,nsy) |
635 |
_RL num_vwind(nsx,nsy) |
_RL num_vwind(nsx,nsy) |
636 |
_RL num_atempm(nsx,nsy) |
_RL num_atempm(nsx,nsy) |
637 |
_RL num_aqhm (nsx,nsy) |
_RL num_aqhm (nsx,nsy) |
638 |
_RL num_precipm(nsx,nsy) |
_RL num_precipm(nsx,nsy) |
639 |
_RL num_swfluxm(nsx,nsy) |
_RL num_swfluxm(nsx,nsy) |
640 |
|
_RL num_swdownm(nsx,nsy) |
641 |
|
_RL num_snowprecipm(nsx,nsy) |
642 |
|
_RL num_lwfluxm(nsx,nsy) |
643 |
|
_RL num_lwdownm(nsx,nsy) |
644 |
|
_RL num_evapm(nsx,nsy) |
645 |
|
_RL num_apressurem(nsx,nsy) |
646 |
|
_RL num_runoffm(nsx,nsy) |
647 |
_RL num_uwindm(nsx,nsy) |
_RL num_uwindm(nsx,nsy) |
648 |
_RL num_vwindm(nsx,nsy) |
_RL num_vwindm(nsx,nsy) |
649 |
_RL num_obcsn(nsx,nsy) |
_RL num_obcsn(nsx,nsy) |
653 |
_RL num_obcsvol |
_RL num_obcsvol |
654 |
_RL num_curmtr(nsx,nsy) |
_RL num_curmtr(nsx,nsy) |
655 |
_RL num_ageos(nsx,nsy) |
_RL num_ageos(nsx,nsy) |
|
_RL num_ice (nsx,nsy) |
|
656 |
_RL num_kapgm(nsx,nsy) |
_RL num_kapgm(nsx,nsy) |
657 |
_RL num_diffkr(nsx,nsy) |
_RL num_diffkr(nsx,nsy) |
658 |
_RL num_theta_ini_fin(nsx,nsy) |
_RL num_theta_ini_fin(nsx,nsy) |
659 |
_RL num_salt_ini_fin(nsx,nsy) |
_RL num_salt_ini_fin(nsx,nsy) |
660 |
_RL num_eddytau(nsx,nsy) |
_RL num_eddytau(nsx,nsy) |
661 |
|
_RL num_bottomdrag(nsx,nsy) |
662 |
|
_RL num_transp |
663 |
|
|
664 |
common /ecco_cost_aux_r/ |
common /ecco_cost_aux_r/ |
665 |
& mult_hflux, |
& mult_hflux, |
666 |
& mult_sflux, |
& mult_sflux, |
667 |
|
& mult_hfluxmm, |
668 |
|
& mult_sfluxmm, |
669 |
& mult_tauu, |
& mult_tauu, |
670 |
& mult_tauv, |
& mult_tauv, |
671 |
& mult_hmean, |
& mult_hmean, |
672 |
& mult_h, |
& mult_h, |
673 |
|
& mult_tp, |
674 |
|
& mult_ers, |
675 |
|
& mult_gfo, |
676 |
& mult_temp, |
& mult_temp, |
677 |
& mult_salt, |
& mult_salt, |
678 |
& mult_temp0, |
& mult_temp0, |
687 |
& mult_xbt, |
& mult_xbt, |
688 |
& mult_argot, |
& mult_argot, |
689 |
& mult_argos, |
& mult_argos, |
690 |
|
& mult_gencost, |
691 |
& mult_drift, |
& mult_drift, |
692 |
& mult_tdrift, |
& mult_tdrift, |
693 |
& mult_sdrift, |
& mult_sdrift, |
698 |
& mult_aqh, |
& mult_aqh, |
699 |
& mult_precip, |
& mult_precip, |
700 |
& mult_swflux, |
& mult_swflux, |
701 |
|
& mult_swdown, |
702 |
|
& mult_snowprecip, |
703 |
|
& mult_lwflux, |
704 |
|
& mult_lwdown, |
705 |
|
& mult_evap, |
706 |
|
& mult_apressure, |
707 |
|
& mult_runoff, |
708 |
& mult_uwind, |
& mult_uwind, |
709 |
& mult_vwind, |
& mult_vwind, |
710 |
& mult_obcsn, |
& mult_obcsn, |
714 |
& mult_obcsvol, |
& mult_obcsvol, |
715 |
& mult_curmtr, |
& mult_curmtr, |
716 |
& mult_ageos, |
& mult_ageos, |
|
& mult_ice, |
|
717 |
& mult_kapgm, |
& mult_kapgm, |
718 |
& mult_diffkr, |
& mult_diffkr, |
719 |
& mult_ini_fin, |
& mult_ini_fin, |
720 |
& mult_eddytau |
& mult_edtau, |
721 |
|
& mult_bottomdrag, |
722 |
|
& mult_smooth_ic, |
723 |
|
& mult_smooth_bc, |
724 |
|
& mult_transp |
725 |
_RL mult_hflux |
_RL mult_hflux |
726 |
_RL mult_sflux |
_RL mult_sflux |
727 |
|
_RL mult_hfluxmm |
728 |
|
_RL mult_sfluxmm |
729 |
_RL mult_tauu |
_RL mult_tauu |
730 |
_RL mult_tauv |
_RL mult_tauv |
731 |
_RL mult_hmean |
_RL mult_hmean |
732 |
_RL mult_h |
_RL mult_h |
733 |
|
_RL mult_tp |
734 |
|
_RL mult_ers |
735 |
|
_RL mult_gfo |
736 |
_RL mult_temp |
_RL mult_temp |
737 |
_RL mult_salt |
_RL mult_salt |
738 |
_RL mult_temp0 |
_RL mult_temp0 |
747 |
_RL mult_xbt |
_RL mult_xbt |
748 |
_RL mult_argot |
_RL mult_argot |
749 |
_RL mult_argos |
_RL mult_argos |
750 |
|
_RL mult_gencost(NGENCOST) |
751 |
_RL mult_drift |
_RL mult_drift |
752 |
_RL mult_tdrift |
_RL mult_tdrift |
753 |
_RL mult_sdrift |
_RL mult_sdrift |
758 |
_RL mult_aqh |
_RL mult_aqh |
759 |
_RL mult_precip |
_RL mult_precip |
760 |
_RL mult_swflux |
_RL mult_swflux |
761 |
|
_RL mult_swdown |
762 |
|
_RL mult_snowprecip |
763 |
|
_RL mult_lwflux |
764 |
|
_RL mult_lwdown |
765 |
|
_RL mult_evap |
766 |
|
_RL mult_apressure |
767 |
|
_RL mult_runoff |
768 |
_RL mult_uwind |
_RL mult_uwind |
769 |
_RL mult_vwind |
_RL mult_vwind |
770 |
_RL mult_obcsn |
_RL mult_obcsn |
774 |
_RL mult_obcsvol |
_RL mult_obcsvol |
775 |
_RL mult_curmtr |
_RL mult_curmtr |
776 |
_RL mult_ageos |
_RL mult_ageos |
|
_RL mult_ice |
|
777 |
_RL mult_kapgm |
_RL mult_kapgm |
778 |
_RL mult_diffkr |
_RL mult_diffkr |
779 |
_RL mult_ini_fin |
_RL mult_ini_fin |
780 |
_RL mult_eddytau |
_RL mult_edtau |
781 |
|
_RL mult_bottomdrag |
782 |
|
_RL mult_smooth_ic |
783 |
|
_RL mult_smooth_bc |
784 |
|
_RL mult_transp |
785 |
|
|
786 |
c Record counters relevant for the cost function evaluation. |
c Record counters relevant for the cost function evaluation. |
787 |
c ========================================================== |
c ========================================================== |
788 |
c |
c |
789 |
c nmonsrec - number of monthly records that will be generated by |
c nyearsrec - number of yearly records that will be generated by |
790 |
c the current model integration. |
c the current model integration. |
791 |
c ndaysrec - number of daily records that will be generated by |
c nmonsrec - number of monthly records that will be generated by |
792 |
c the current model integration. |
c the current model integration. |
793 |
|
c ndaysrec - number of daily records that will be generated by |
794 |
|
c the current model integration. |
795 |
|
|
796 |
common /ecco_cost_i/ |
common /ecco_cost_i/ |
797 |
|
& nyearsrec, |
798 |
& nmonsrec, |
& nmonsrec, |
799 |
& ndaysrec |
& ndaysrec, |
800 |
|
& nnztbar, |
801 |
|
& nnzsbar |
802 |
|
integer nyearsrec |
803 |
integer nmonsrec |
integer nmonsrec |
804 |
integer ndaysrec |
integer ndaysrec |
805 |
|
integer nnztbar |
806 |
|
integer nnzsbar |
807 |
|
|
808 |
|
|
809 |
c Data files for the weights used in the cost function: |
c Data files for the weights used in the cost function: |
843 |
& geoid_errfile, |
& geoid_errfile, |
844 |
& geoid_covariancefile, |
& geoid_covariancefile, |
845 |
& ssh_errfile, |
& ssh_errfile, |
846 |
|
& tp_errfile, |
847 |
|
& ers_errfile, |
848 |
|
& gfo_errfile, |
849 |
& ctdt_errfile, |
& ctdt_errfile, |
850 |
& ctds_errfile, |
& ctds_errfile, |
851 |
& drift_errfile, |
& drift_errfile, |
854 |
& salterrfile, |
& salterrfile, |
855 |
& temperrfile, |
& temperrfile, |
856 |
& velerrfile, |
& velerrfile, |
857 |
|
& salt0errfile, |
858 |
|
& temp0errfile, |
859 |
|
& vel0errfile, |
860 |
& atemp_errfile, |
& atemp_errfile, |
861 |
& aqh_errfile, |
& aqh_errfile, |
862 |
& precip_errfile, |
& precip_errfile, |
863 |
& swflux_errfile, |
& swflux_errfile, |
864 |
|
& swdown_errfile, |
865 |
|
& snowprecip_errfile, |
866 |
|
& lwflux_errfile, |
867 |
|
& lwdown_errfile, |
868 |
|
& evap_errfile, |
869 |
|
& apressure_errfile, |
870 |
|
& runoff_errfile, |
871 |
|
& edtau_errfile, |
872 |
|
& kapgm_errfile, |
873 |
|
& diffkr_errfile, |
874 |
|
& bottomdrag_errfile, |
875 |
|
& gencost_errfile, |
876 |
& uwind_errfile, |
& uwind_errfile, |
877 |
& vwind_errfile |
& vwind_errfile |
878 |
character*(MAX_LEN_FNAM) hflux_errfile |
character*(MAX_LEN_FNAM) hflux_errfile |
889 |
character*(MAX_LEN_FNAM) geoid_errfile |
character*(MAX_LEN_FNAM) geoid_errfile |
890 |
character*(MAX_LEN_FNAM) geoid_covariancefile |
character*(MAX_LEN_FNAM) geoid_covariancefile |
891 |
character*(MAX_LEN_FNAM) ssh_errfile |
character*(MAX_LEN_FNAM) ssh_errfile |
892 |
|
character*(MAX_LEN_FNAM) tp_errfile |
893 |
|
character*(MAX_LEN_FNAM) ers_errfile |
894 |
|
character*(MAX_LEN_FNAM) gfo_errfile |
895 |
character*(MAX_LEN_FNAM) ctdt_errfile |
character*(MAX_LEN_FNAM) ctdt_errfile |
896 |
character*(MAX_LEN_FNAM) ctds_errfile |
character*(MAX_LEN_FNAM) ctds_errfile |
897 |
character*(MAX_LEN_FNAM) drift_errfile |
character*(MAX_LEN_FNAM) drift_errfile |
900 |
character*(MAX_LEN_FNAM) salterrfile |
character*(MAX_LEN_FNAM) salterrfile |
901 |
character*(MAX_LEN_FNAM) temperrfile |
character*(MAX_LEN_FNAM) temperrfile |
902 |
character*(MAX_LEN_FNAM) velerrfile |
character*(MAX_LEN_FNAM) velerrfile |
903 |
|
character*(MAX_LEN_FNAM) salt0errfile |
904 |
|
character*(MAX_LEN_FNAM) temp0errfile |
905 |
|
character*(MAX_LEN_FNAM) vel0errfile |
906 |
character*(MAX_LEN_FNAM) atemp_errfile |
character*(MAX_LEN_FNAM) atemp_errfile |
907 |
character*(MAX_LEN_FNAM) aqh_errfile |
character*(MAX_LEN_FNAM) aqh_errfile |
908 |
character*(MAX_LEN_FNAM) precip_errfile |
character*(MAX_LEN_FNAM) precip_errfile |
909 |
character*(MAX_LEN_FNAM) swflux_errfile |
character*(MAX_LEN_FNAM) swflux_errfile |
910 |
|
character*(MAX_LEN_FNAM) swdown_errfile |
911 |
|
character*(MAX_LEN_FNAM) snowprecip_errfile |
912 |
|
character*(MAX_LEN_FNAM) lwflux_errfile |
913 |
|
character*(MAX_LEN_FNAM) lwdown_errfile |
914 |
|
character*(MAX_LEN_FNAM) evap_errfile |
915 |
|
character*(MAX_LEN_FNAM) apressure_errfile |
916 |
|
character*(MAX_LEN_FNAM) runoff_errfile |
917 |
|
character*(MAX_LEN_FNAM) edtau_errfile |
918 |
|
character*(MAX_LEN_FNAM) kapgm_errfile |
919 |
|
character*(MAX_LEN_FNAM) diffkr_errfile |
920 |
|
character*(MAX_LEN_FNAM) bottomdrag_errfile |
921 |
|
character*(MAX_LEN_FNAM) gencost_errfile(NGENCOST) |
922 |
character*(MAX_LEN_FNAM) uwind_errfile |
character*(MAX_LEN_FNAM) uwind_errfile |
923 |
character*(MAX_LEN_FNAM) vwind_errfile |
character*(MAX_LEN_FNAM) vwind_errfile |
924 |
|
|
951 |
& frame, |
& frame, |
952 |
& cosphi, |
& cosphi, |
953 |
& whflux,wsflux,wtauu,wtauv, |
& whflux,wsflux,wtauu,wtauv, |
954 |
& watemp,waqh,wprecip,wswflux,wuwind,wvwind, |
& watemp,waqh,wprecip,wsnowprecip, |
955 |
|
& wswflux,wswdown,wlwflux,wlwdown, |
956 |
|
& wevap,wapressure,wrunoff, |
957 |
|
& wbottomdrag, |
958 |
|
& wuwind,wvwind, |
959 |
& wscatx,wscaty, |
& wscatx,wscaty, |
960 |
& wtheta,wtheta2,wthetaLev, |
& wtheta,wtheta2,wthetaLev, |
961 |
& wsalt,wsalt2,wsaltLev, |
& wsalt,wsalt2,wsaltLev, |
962 |
|
& wdiffkr,wdiffkr2,wdiffkrFld, |
963 |
|
& wkapgm,wkapgm2,wkapgmFld, |
964 |
|
& wedtaux,wedtaux2,wedtauxFld, |
965 |
|
& wedtauy,wedtauy2,wedtauyFld, |
966 |
& wsst,wsss, |
& wsst,wsss, |
967 |
& wtp,wers, |
& wtp,wers,wgfo, |
968 |
& wp, |
& wp, |
969 |
& wctdt,wctds, |
& wctdt,wctds, |
970 |
& wudrift,wvdrift, |
& wudrift,wvdrift, |
990 |
_RL waqh (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL waqh (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
991 |
_RL wprecip (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wprecip (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
992 |
_RL wswflux (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wswflux (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
993 |
|
_RL wswdown (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
994 |
|
_RL wsnowprecip (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
995 |
|
_RL wlwflux (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
996 |
|
_RL wlwdown (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
997 |
|
_RL wevap (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
998 |
|
_RL wapressure(1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
999 |
|
_RL wrunoff (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1000 |
|
_RL wbottomdrag (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy) |
1001 |
_RL wuwind (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wuwind (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1002 |
_RL wvwind (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wvwind (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1003 |
_RL wtheta ( nr,nsx,nsy) |
_RL wtheta ( nr,nsx,nsy) |
1010 |
_RL wsss (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wsss (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1011 |
_RL wtp (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wtp (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1012 |
_RL wers (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wers (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1013 |
|
_RL wgfo (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1014 |
_RL wp (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wp (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1015 |
_RL wctdt ( nr,nsx,nsy) |
_RL wctdt ( nr,nsx,nsy) |
1016 |
_RL wctds ( nr,nsx,nsy) |
_RL wctds ( nr,nsx,nsy) |
1020 |
_RL wcurrent2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
_RL wcurrent2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1021 |
_RL wcurrentLev (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
_RL wcurrentLev (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1022 |
_RL wbaro |
_RL wbaro |
1023 |
|
_RL wdiffkr ( nr,nsx,nsy) |
1024 |
|
_RL wdiffkr2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1025 |
|
_RL wdiffkrFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1026 |
|
_RL wkapgm ( nr,nsx,nsy) |
1027 |
|
_RL wkapgm2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1028 |
|
_RL wkapgmFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1029 |
|
_RL wedtaux ( nr,nsx,nsy) |
1030 |
|
_RL wedtaux2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1031 |
|
_RL wedtauxFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1032 |
|
_RL wedtauy ( nr,nsx,nsy) |
1033 |
|
_RL wedtauy2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1034 |
|
_RL wedtauyFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1035 |
|
|
1036 |
common /ecco_cost_weights_0_r/ |
common /ecco_cost_weights_0_r/ |
1037 |
& whflux0, wsflux0, wtau0, |
& whflux0, wsflux0, wtau0, |
1038 |
& watemp0, waqh0, wprecip0, wswflux0, wwind0 |
& watemp0, waqh0, wprecip0, wsnowprecip0, wwind0, |
1039 |
|
& wswflux0, wswdown0, wlwflux0, wlwdown0, |
1040 |
|
& wevap0, wapressure0, wrunoff0, |
1041 |
|
& wbottomdrag0,wdiffkr0, wkapgm0, wedtau0 |
1042 |
_RL whflux0 |
_RL whflux0 |
1043 |
_RL wsflux0 |
_RL wsflux0 |
1044 |
_RL wtau0 |
_RL wtau0 |
1046 |
_RL waqh0 |
_RL waqh0 |
1047 |
_RL wprecip0 |
_RL wprecip0 |
1048 |
_RL wswflux0 |
_RL wswflux0 |
1049 |
|
_RL wswdown0 |
1050 |
|
_RL wsnowprecip0 |
1051 |
|
_RL wlwflux0 |
1052 |
|
_RL wlwdown0 |
1053 |
|
_RL wevap0 |
1054 |
|
_RL wapressure0 |
1055 |
|
_RL wrunoff0 |
1056 |
|
_RL wbottomdrag0 |
1057 |
_RL wwind0 |
_RL wwind0 |
1058 |
|
_RL wdiffkr0 |
1059 |
|
_RL wkapgm0 |
1060 |
|
_RL wedtau0 |
1061 |
|
|
1062 |
common /ecco_cost_weights_mean_r/ |
common /ecco_cost_weights_mean_r/ |
1063 |
& wmean_hflux, wmean_sflux, wmean_tau, |
& wmean_hflux, wmean_sflux, wmean_tau, |
1064 |
& wmean_atemp, wmean_aqh, |
& wmean_atemp, wmean_aqh, |
1065 |
& wmean_precip, wmean_swflux, wmean_wind |
& wmean_precip, wmean_snowprecip, wmean_wind, |
1066 |
|
& wmean_swflux, wmean_swdown, wmean_lwflux, wmean_lwdown, |
1067 |
|
& wmean_evap, wmean_apressure, wmean_runoff |
1068 |
_RL wmean_hflux |
_RL wmean_hflux |
1069 |
_RL wmean_sflux |
_RL wmean_sflux |
1070 |
_RL wmean_tau |
_RL wmean_tau |
1072 |
_RL wmean_aqh |
_RL wmean_aqh |
1073 |
_RL wmean_precip |
_RL wmean_precip |
1074 |
_RL wmean_swflux |
_RL wmean_swflux |
1075 |
|
_RL wmean_swdown |
1076 |
|
_RL wmean_snowprecip |
1077 |
|
_RL wmean_lwflux |
1078 |
|
_RL wmean_lwdown |
1079 |
|
_RL wmean_evap |
1080 |
|
_RL wmean_apressure |
1081 |
|
_RL wmean_runoff |
1082 |
_RL wmean_wind |
_RL wmean_wind |
1083 |
|
|
1084 |
common /ecco_cost_weights_2_r/ |
common /ecco_cost_weights_2_r/ |
1172 |
& tpmask, |
& tpmask, |
1173 |
& ersobs, |
& ersobs, |
1174 |
& ersmask, |
& ersmask, |
1175 |
|
& gfoobs, |
1176 |
|
& gfomask, |
1177 |
& ctdtobs, |
& ctdtobs, |
1178 |
& ctdsobs, |
& ctdsobs, |
1179 |
& xbtobs, |
& xbtobs, |
1204 |
_RL tpmask (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL tpmask (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1205 |
_RL ersobs (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL ersobs (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1206 |
_RL ersmask (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL ersmask (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1207 |
|
_RL gfoobs (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1208 |
|
_RL gfomask (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
1209 |
_RL ctdtobs (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
_RL ctdtobs (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1210 |
_RL ctdsobs (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
_RL ctdsobs (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1211 |
_RL xbtobs (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
_RL xbtobs (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy) |
1251 |
& topexmeanfile, |
& topexmeanfile, |
1252 |
& topexfile, |
& topexfile, |
1253 |
& ersfile, |
& ersfile, |
1254 |
|
& gfofile, |
1255 |
& ctdtfile, |
& ctdtfile, |
1256 |
& ctdsfile, |
& ctdsfile, |
1257 |
& ctdtclimfile, |
& ctdtclimfile, |
1261 |
& argosfile, |
& argosfile, |
1262 |
& udriftfile, |
& udriftfile, |
1263 |
& vdriftfile, |
& vdriftfile, |
1264 |
|
& gencost_datafile, |
1265 |
& curmtrufile, |
& curmtrufile, |
1266 |
& curmtrvfile |
& curmtrvfile |
1267 |
|
|
1275 |
character*(MAX_LEN_FNAM) topexmeanfile |
character*(MAX_LEN_FNAM) topexmeanfile |
1276 |
character*(MAX_LEN_FNAM) topexfile |
character*(MAX_LEN_FNAM) topexfile |
1277 |
character*(MAX_LEN_FNAM) ersfile |
character*(MAX_LEN_FNAM) ersfile |
1278 |
|
character*(MAX_LEN_FNAM) gfofile |
1279 |
character*(MAX_LEN_FNAM) ctdtfile |
character*(MAX_LEN_FNAM) ctdtfile |
1280 |
character*(MAX_LEN_FNAM) ctdsfile |
character*(MAX_LEN_FNAM) ctdsfile |
1281 |
character*(MAX_LEN_FNAM) ctdtclimfile |
character*(MAX_LEN_FNAM) ctdtclimfile |
1284 |
character*(MAX_LEN_FNAM) argotfile |
character*(MAX_LEN_FNAM) argotfile |
1285 |
character*(MAX_LEN_FNAM) argosfile |
character*(MAX_LEN_FNAM) argosfile |
1286 |
character*(MAX_LEN_FNAM) argofile |
character*(MAX_LEN_FNAM) argofile |
1287 |
|
character*(MAX_LEN_FNAM) gencost_datafile(NGENCOST) |
1288 |
character*(MAX_LEN_FNAM) udriftfile |
character*(MAX_LEN_FNAM) udriftfile |
1289 |
character*(MAX_LEN_FNAM) vdriftfile |
character*(MAX_LEN_FNAM) vdriftfile |
1290 |
character*(MAX_LEN_FNAM) curmtrufile |
character*(MAX_LEN_FNAM) curmtrufile |
1298 |
|
|
1299 |
common /ecco_cost_data_flags/ |
common /ecco_cost_data_flags/ |
1300 |
& using_topex, |
& using_topex, |
1301 |
& using_ers |
& using_ers, |
1302 |
|
& using_gfo |
1303 |
logical using_topex |
logical using_topex |
1304 |
logical using_ers |
logical using_ers |
1305 |
|
logical using_gfo |
1306 |
|
|
1307 |
c Calendar information for the observations: |
c Calendar information for the observations: |
1308 |
c ========================================== |
c ========================================== |
1322 |
& tmistartdate, |
& tmistartdate, |
1323 |
& sssstartdate, |
& sssstartdate, |
1324 |
& topexstartdate, |
& topexstartdate, |
1325 |
& ersstartdate |
& ersstartdate, |
1326 |
|
& gfostartdate |
1327 |
integer scatxstartdate(4) |
integer scatxstartdate(4) |
1328 |
integer scatystartdate(4) |
integer scatystartdate(4) |
1329 |
integer sststartdate(4) |
integer sststartdate(4) |
1333 |
integer sssstartdate(4) |
integer sssstartdate(4) |
1334 |
integer topexstartdate(4) |
integer topexstartdate(4) |
1335 |
integer ersstartdate(4) |
integer ersstartdate(4) |
1336 |
|
integer gfostartdate(4) |
1337 |
|
|
1338 |
common /ecco_cost_data_aux_i/ |
common /ecco_cost_data_aux_i/ |
1339 |
& tmistartdate1, |
& tmistartdate1, |
1350 |
& topexstartdate2, |
& topexstartdate2, |
1351 |
& ersstartdate1, |
& ersstartdate1, |
1352 |
& ersstartdate2, |
& ersstartdate2, |
1353 |
|
& gfostartdate1, |
1354 |
|
& gfostartdate2, |
1355 |
& scatstartdate1, |
& scatstartdate1, |
1356 |
& scatstartdate2, |
& scatstartdate2 |
|
& costIceStart1, |
|
|
& costIceStart2, |
|
|
& costIceEnd1, |
|
|
& costIceEnd2 |
|
1357 |
|
|
1358 |
integer tmistartdate1 |
integer tmistartdate1 |
1359 |
integer tmistartdate2 |
integer tmistartdate2 |
1369 |
integer topexstartdate2 |
integer topexstartdate2 |
1370 |
integer ersstartdate1 |
integer ersstartdate1 |
1371 |
integer ersstartdate2 |
integer ersstartdate2 |
1372 |
|
integer gfostartdate1 |
1373 |
|
integer gfostartdate2 |
1374 |
integer scatstartdate1 |
integer scatstartdate1 |
1375 |
integer scatstartdate2 |
integer scatstartdate2 |
|
integer costIceStart1 |
|
|
integer costIceStart2 |
|
|
integer costIceEnd1 |
|
|
integer costIceEnd2 |
|
1376 |
|
|
1377 |
common /ecco_cost_data_times_r/ |
common /ecco_cost_data_times_r/ |
1378 |
& topexperiod, |
& topexperiod, |
1379 |
& ersperiod, |
& ersperiod, |
1380 |
& scatperiod, |
& gfoperiod, |
1381 |
& costIceStart, |
& scatperiod |
|
& costIceEnd |
|
1382 |
_RL topexperiod |
_RL topexperiod |
1383 |
_RL ersperiod |
_RL ersperiod |
1384 |
|
_RL gfoperiod |
1385 |
_RL scatperiod |
_RL scatperiod |
|
_RL costIceStart |
|
|
_RL costIceEnd |
|
|
|
|
|
c cost_ice_flag - cost_ice flag (see cost_ice.F) |
|
1386 |
|
|
1387 |
common /ecco_cost_ice_i/ cost_ice_flag |
common /ecco_cost_data_detrend/ |
1388 |
integer cost_ice_flag |
& topexintercept, |
1389 |
|
& ersintercept, |
1390 |
|
& gfointercept, |
1391 |
|
& topexslope, |
1392 |
|
& ersslope, |
1393 |
|
& gfoslope |
1394 |
|
_RL topexintercept |
1395 |
|
_RL ersintercept |
1396 |
|
_RL gfointercept |
1397 |
|
_RL topexslope |
1398 |
|
_RL ersslope |
1399 |
|
_RL gfoslope |
1400 |
|
|
1401 |
c ================================================================== |
c ================================================================== |
1402 |
c END OF HEADER COST |
c END OF HEADER COST |