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 ================================================================== |
272 |
c objf_scatym - time-mean meridional SCAT contribution |
c objf_scatym - time-mean meridional SCAT contribution |
273 |
c objf_scatx - zonal SCAT contribution |
c objf_scatx - zonal SCAT contribution |
274 |
c objf_scaty - meridional SCAT contribution |
c objf_scaty - meridional SCAT contribution |
|
c objf_ice - sea-ice volume |
|
275 |
c objf_kapgm - kappa GM contribution |
c objf_kapgm - kappa GM contribution |
276 |
c objf_diffkr - diffusion contribution |
c objf_diffkr - diffusion contribution |
277 |
c objf_theta_ini_fin - final vs. initial theta misfit |
c objf_theta_ini_fin - final vs. initial theta misfit |
284 |
& objf_hflux, |
& objf_hflux, |
285 |
& objf_hfluxm, |
& objf_hfluxm, |
286 |
& objf_hfluxmm, |
& objf_hfluxmm, |
287 |
|
& objf_hfluxsmoo, |
288 |
& objf_sflux, |
& objf_sflux, |
289 |
& objf_sfluxm, |
& objf_sfluxm, |
290 |
& objf_sfluxmm, |
& objf_sfluxmm, |
291 |
|
& objf_sfluxsmoo, |
292 |
& objf_tauu, |
& objf_tauu, |
293 |
& objf_tauum, |
& objf_tauum, |
294 |
|
& objf_tauusmoo, |
295 |
& objf_tauv, |
& objf_tauv, |
296 |
& objf_tauvm, |
& objf_tauvm, |
297 |
|
& objf_tauvsmoo, |
298 |
& objf_hmean, |
& objf_hmean, |
299 |
& objf_h, |
& objf_h, |
300 |
& objf_temp, |
& objf_temp, |
301 |
& objf_salt, |
& objf_salt, |
302 |
& objf_temp0, |
& objf_temp0, |
303 |
& objf_salt0, |
& objf_salt0, |
304 |
|
& objf_temp0smoo, |
305 |
|
& objf_salt0smoo, |
306 |
& objf_sst, |
& objf_sst, |
307 |
& objf_tmi, |
& objf_tmi, |
308 |
& objf_sss, |
& objf_sss, |
325 |
& objf_aqh, |
& objf_aqh, |
326 |
& objf_precip, |
& objf_precip, |
327 |
& objf_swflux, |
& objf_swflux, |
328 |
|
& objf_swdown, |
329 |
& objf_uwind, |
& objf_uwind, |
330 |
& objf_vwind, |
& objf_vwind, |
331 |
& objf_atempm, |
& objf_atempm, |
332 |
& objf_aqhm, |
& objf_aqhm, |
333 |
& objf_precipm, |
& objf_precipm, |
334 |
& objf_swfluxm, |
& objf_swfluxm, |
335 |
|
& objf_swdownm, |
336 |
& objf_uwindm, |
& objf_uwindm, |
337 |
& objf_vwindm, |
& objf_vwindm, |
338 |
|
& objf_atempsmoo, |
339 |
|
& objf_aqhsmoo, |
340 |
|
& objf_precipsmoo, |
341 |
|
& objf_swfluxsmoo, |
342 |
|
& objf_swdownsmoo, |
343 |
|
& objf_uwindsmoo, |
344 |
|
& objf_vwindsmoo, |
345 |
& objf_obcsn, |
& objf_obcsn, |
346 |
& objf_obcss, |
& objf_obcss, |
347 |
& objf_obcsw, |
& objf_obcsw, |
349 |
& objf_obcsvol, |
& objf_obcsvol, |
350 |
& objf_curmtr, |
& objf_curmtr, |
351 |
& objf_ageos, |
& objf_ageos, |
352 |
& objf_ice, |
& objf_kapgm, |
353 |
& objf_diffkr, |
& objf_diffkr, |
354 |
& objf_theta_ini_fin, |
& objf_theta_ini_fin, |
355 |
& objf_salt_ini_fin, |
& objf_salt_ini_fin, |
358 |
_RL objf_hflux (nsx,nsy) |
_RL objf_hflux (nsx,nsy) |
359 |
_RL objf_hfluxm (nsx,nsy) |
_RL objf_hfluxm (nsx,nsy) |
360 |
_RL objf_hfluxmm(nsx,nsy) |
_RL objf_hfluxmm(nsx,nsy) |
361 |
|
_RL objf_hfluxsmoo (nsx,nsy) |
362 |
_RL objf_sflux (nsx,nsy) |
_RL objf_sflux (nsx,nsy) |
363 |
_RL objf_sfluxm (nsx,nsy) |
_RL objf_sfluxm (nsx,nsy) |
364 |
_RL objf_sfluxmm(nsx,nsy) |
_RL objf_sfluxmm(nsx,nsy) |
365 |
|
_RL objf_sfluxsmoo (nsx,nsy) |
366 |
_RL objf_tauu (nsx,nsy) |
_RL objf_tauu (nsx,nsy) |
367 |
_RL objf_tauum (nsx,nsy) |
_RL objf_tauum (nsx,nsy) |
368 |
|
_RL objf_tauusmoo (nsx,nsy) |
369 |
_RL objf_tauv (nsx,nsy) |
_RL objf_tauv (nsx,nsy) |
370 |
_RL objf_tauvm (nsx,nsy) |
_RL objf_tauvm (nsx,nsy) |
371 |
|
_RL objf_tauvsmoo (nsx,nsy) |
372 |
_RL objf_hmean |
_RL objf_hmean |
373 |
_RL objf_h (nsx,nsy) |
_RL objf_h (nsx,nsy) |
374 |
_RL objf_temp (nsx,nsy) |
_RL objf_temp (nsx,nsy) |
375 |
_RL objf_salt (nsx,nsy) |
_RL objf_salt (nsx,nsy) |
376 |
_RL objf_temp0(nsx,nsy) |
_RL objf_temp0(nsx,nsy) |
377 |
_RL objf_salt0(nsx,nsy) |
_RL objf_salt0(nsx,nsy) |
378 |
|
_RL objf_temp0smoo(nsx,nsy) |
379 |
|
_RL objf_salt0smoo(nsx,nsy) |
380 |
_RL objf_sst (nsx,nsy) |
_RL objf_sst (nsx,nsy) |
381 |
_RL objf_tmi (nsx,nsy) |
_RL objf_tmi (nsx,nsy) |
382 |
_RL objf_sss (nsx,nsy) |
_RL objf_sss (nsx,nsy) |
399 |
_RL objf_aqh (nsx,nsy) |
_RL objf_aqh (nsx,nsy) |
400 |
_RL objf_precip(nsx,nsy) |
_RL objf_precip(nsx,nsy) |
401 |
_RL objf_swflux(nsx,nsy) |
_RL objf_swflux(nsx,nsy) |
402 |
|
_RL objf_swdown(nsx,nsy) |
403 |
_RL objf_uwind(nsx,nsy) |
_RL objf_uwind(nsx,nsy) |
404 |
_RL objf_vwind(nsx,nsy) |
_RL objf_vwind(nsx,nsy) |
405 |
_RL objf_atempm(nsx,nsy) |
_RL objf_atempm(nsx,nsy) |
406 |
_RL objf_aqhm (nsx,nsy) |
_RL objf_aqhm (nsx,nsy) |
407 |
_RL objf_precipm(nsx,nsy) |
_RL objf_precipm(nsx,nsy) |
408 |
_RL objf_swfluxm(nsx,nsy) |
_RL objf_swfluxm(nsx,nsy) |
409 |
|
_RL objf_swdownm(nsx,nsy) |
410 |
_RL objf_uwindm(nsx,nsy) |
_RL objf_uwindm(nsx,nsy) |
411 |
_RL objf_vwindm(nsx,nsy) |
_RL objf_vwindm(nsx,nsy) |
412 |
|
_RL objf_atempsmoo(nsx,nsy) |
413 |
|
_RL objf_aqhsmoo (nsx,nsy) |
414 |
|
_RL objf_precipsmoo(nsx,nsy) |
415 |
|
_RL objf_swfluxsmoo(nsx,nsy) |
416 |
|
_RL objf_swdownsmoo(nsx,nsy) |
417 |
|
_RL objf_uwindsmoo(nsx,nsy) |
418 |
|
_RL objf_vwindsmoo(nsx,nsy) |
419 |
_RL objf_obcsn(nsx,nsy) |
_RL objf_obcsn(nsx,nsy) |
420 |
_RL objf_obcss(nsx,nsy) |
_RL objf_obcss(nsx,nsy) |
421 |
_RL objf_obcsw(nsx,nsy) |
_RL objf_obcsw(nsx,nsy) |
423 |
_RL objf_obcsvol |
_RL objf_obcsvol |
424 |
_RL objf_curmtr(nsx,nsy) |
_RL objf_curmtr(nsx,nsy) |
425 |
_RL objf_ageos(nsx,nsy) |
_RL objf_ageos(nsx,nsy) |
|
_RL objf_ice (nsx,nsy) |
|
426 |
_RL objf_kapgm(nsx,nsy) |
_RL objf_kapgm(nsx,nsy) |
427 |
_RL objf_diffkr(nsx,nsy) |
_RL objf_diffkr(nsx,nsy) |
428 |
_RL objf_theta_ini_fin(nsx,nsy) |
_RL objf_theta_ini_fin(nsx,nsy) |
468 |
& num_aqh, |
& num_aqh, |
469 |
& num_precip, |
& num_precip, |
470 |
& num_swflux, |
& num_swflux, |
471 |
|
& num_swdown, |
472 |
& num_uwind, |
& num_uwind, |
473 |
& num_vwind, |
& num_vwind, |
474 |
& num_atempm, |
& num_atempm, |
475 |
& num_aqhm, |
& num_aqhm, |
476 |
& num_precipm, |
& num_precipm, |
477 |
& num_swfluxm, |
& num_swfluxm, |
478 |
|
& num_swdownm, |
479 |
& num_uwindm, |
& num_uwindm, |
480 |
& num_vwindm, |
& num_vwindm, |
481 |
& num_obcsn, |
& num_obcsn, |
485 |
& num_obcsvol, |
& num_obcsvol, |
486 |
& num_curmtr, |
& num_curmtr, |
487 |
& num_ageos, |
& num_ageos, |
488 |
& num_ice, |
& num_kapgm, |
489 |
& num_diffkr, |
& num_diffkr, |
490 |
& num_theta_ini_fin, |
& num_theta_ini_fin, |
491 |
& num_salt_ini_fin, |
& num_salt_ini_fin, |
529 |
_RL num_aqh (nsx,nsy) |
_RL num_aqh (nsx,nsy) |
530 |
_RL num_precip(nsx,nsy) |
_RL num_precip(nsx,nsy) |
531 |
_RL num_swflux(nsx,nsy) |
_RL num_swflux(nsx,nsy) |
532 |
|
_RL num_swdown(nsx,nsy) |
533 |
_RL num_uwind(nsx,nsy) |
_RL num_uwind(nsx,nsy) |
534 |
_RL num_vwind(nsx,nsy) |
_RL num_vwind(nsx,nsy) |
535 |
_RL num_atempm(nsx,nsy) |
_RL num_atempm(nsx,nsy) |
536 |
_RL num_aqhm (nsx,nsy) |
_RL num_aqhm (nsx,nsy) |
537 |
_RL num_precipm(nsx,nsy) |
_RL num_precipm(nsx,nsy) |
538 |
_RL num_swfluxm(nsx,nsy) |
_RL num_swfluxm(nsx,nsy) |
539 |
|
_RL num_swdownm(nsx,nsy) |
540 |
_RL num_uwindm(nsx,nsy) |
_RL num_uwindm(nsx,nsy) |
541 |
_RL num_vwindm(nsx,nsy) |
_RL num_vwindm(nsx,nsy) |
542 |
_RL num_obcsn(nsx,nsy) |
_RL num_obcsn(nsx,nsy) |
546 |
_RL num_obcsvol |
_RL num_obcsvol |
547 |
_RL num_curmtr(nsx,nsy) |
_RL num_curmtr(nsx,nsy) |
548 |
_RL num_ageos(nsx,nsy) |
_RL num_ageos(nsx,nsy) |
|
_RL num_ice (nsx,nsy) |
|
549 |
_RL num_kapgm(nsx,nsy) |
_RL num_kapgm(nsx,nsy) |
550 |
_RL num_diffkr(nsx,nsy) |
_RL num_diffkr(nsx,nsy) |
551 |
_RL num_theta_ini_fin(nsx,nsy) |
_RL num_theta_ini_fin(nsx,nsy) |
555 |
common /ecco_cost_aux_r/ |
common /ecco_cost_aux_r/ |
556 |
& mult_hflux, |
& mult_hflux, |
557 |
& mult_sflux, |
& mult_sflux, |
558 |
|
& mult_hfluxmm, |
559 |
|
& mult_sfluxmm, |
560 |
& mult_tauu, |
& mult_tauu, |
561 |
& mult_tauv, |
& mult_tauv, |
562 |
& mult_hmean, |
& mult_hmean, |
585 |
& mult_aqh, |
& mult_aqh, |
586 |
& mult_precip, |
& mult_precip, |
587 |
& mult_swflux, |
& mult_swflux, |
588 |
|
& mult_swdown, |
589 |
& mult_uwind, |
& mult_uwind, |
590 |
& mult_vwind, |
& mult_vwind, |
591 |
& mult_obcsn, |
& mult_obcsn, |
595 |
& mult_obcsvol, |
& mult_obcsvol, |
596 |
& mult_curmtr, |
& mult_curmtr, |
597 |
& mult_ageos, |
& mult_ageos, |
|
& mult_ice, |
|
598 |
& mult_kapgm, |
& mult_kapgm, |
599 |
& mult_diffkr, |
& mult_diffkr, |
600 |
& mult_ini_fin, |
& mult_ini_fin, |
601 |
& mult_eddytau |
& mult_eddytau, |
602 |
|
& mult_smooth_ic, |
603 |
|
& mult_smooth_bc |
604 |
|
|
605 |
_RL mult_hflux |
_RL mult_hflux |
606 |
_RL mult_sflux |
_RL mult_sflux |
607 |
|
_RL mult_hfluxmm |
608 |
|
_RL mult_sfluxmm |
609 |
_RL mult_tauu |
_RL mult_tauu |
610 |
_RL mult_tauv |
_RL mult_tauv |
611 |
_RL mult_hmean |
_RL mult_hmean |
634 |
_RL mult_aqh |
_RL mult_aqh |
635 |
_RL mult_precip |
_RL mult_precip |
636 |
_RL mult_swflux |
_RL mult_swflux |
637 |
|
_RL mult_swdown |
638 |
_RL mult_uwind |
_RL mult_uwind |
639 |
_RL mult_vwind |
_RL mult_vwind |
640 |
_RL mult_obcsn |
_RL mult_obcsn |
644 |
_RL mult_obcsvol |
_RL mult_obcsvol |
645 |
_RL mult_curmtr |
_RL mult_curmtr |
646 |
_RL mult_ageos |
_RL mult_ageos |
|
_RL mult_ice |
|
647 |
_RL mult_kapgm |
_RL mult_kapgm |
648 |
_RL mult_diffkr |
_RL mult_diffkr |
649 |
_RL mult_ini_fin |
_RL mult_ini_fin |
650 |
_RL mult_eddytau |
_RL mult_eddytau |
651 |
|
_RL mult_smooth_ic |
652 |
|
_RL mult_smooth_bc |
653 |
|
|
654 |
c Record counters relevant for the cost function evaluation. |
c Record counters relevant for the cost function evaluation. |
655 |
c ========================================================== |
c ========================================================== |
661 |
|
|
662 |
common /ecco_cost_i/ |
common /ecco_cost_i/ |
663 |
& nmonsrec, |
& nmonsrec, |
664 |
& ndaysrec |
& ndaysrec, |
665 |
|
& nnztbar, |
666 |
|
& nnzsbar |
667 |
integer nmonsrec |
integer nmonsrec |
668 |
integer ndaysrec |
integer ndaysrec |
669 |
|
integer nnztbar |
670 |
|
integer nnzsbar |
671 |
|
|
672 |
|
|
673 |
c Data files for the weights used in the cost function: |
c Data files for the weights used in the cost function: |
715 |
& salterrfile, |
& salterrfile, |
716 |
& temperrfile, |
& temperrfile, |
717 |
& velerrfile, |
& velerrfile, |
718 |
|
& salt0errfile, |
719 |
|
& temp0errfile, |
720 |
|
& vel0errfile, |
721 |
& atemp_errfile, |
& atemp_errfile, |
722 |
& aqh_errfile, |
& aqh_errfile, |
723 |
& precip_errfile, |
& precip_errfile, |
724 |
& swflux_errfile, |
& swflux_errfile, |
725 |
|
& swdown_errfile, |
726 |
& uwind_errfile, |
& uwind_errfile, |
727 |
& vwind_errfile |
& vwind_errfile |
728 |
character*(MAX_LEN_FNAM) hflux_errfile |
character*(MAX_LEN_FNAM) hflux_errfile |
747 |
character*(MAX_LEN_FNAM) salterrfile |
character*(MAX_LEN_FNAM) salterrfile |
748 |
character*(MAX_LEN_FNAM) temperrfile |
character*(MAX_LEN_FNAM) temperrfile |
749 |
character*(MAX_LEN_FNAM) velerrfile |
character*(MAX_LEN_FNAM) velerrfile |
750 |
|
character*(MAX_LEN_FNAM) salt0errfile |
751 |
|
character*(MAX_LEN_FNAM) temp0errfile |
752 |
|
character*(MAX_LEN_FNAM) vel0errfile |
753 |
character*(MAX_LEN_FNAM) atemp_errfile |
character*(MAX_LEN_FNAM) atemp_errfile |
754 |
character*(MAX_LEN_FNAM) aqh_errfile |
character*(MAX_LEN_FNAM) aqh_errfile |
755 |
character*(MAX_LEN_FNAM) precip_errfile |
character*(MAX_LEN_FNAM) precip_errfile |
756 |
character*(MAX_LEN_FNAM) swflux_errfile |
character*(MAX_LEN_FNAM) swflux_errfile |
757 |
|
character*(MAX_LEN_FNAM) swdown_errfile |
758 |
character*(MAX_LEN_FNAM) uwind_errfile |
character*(MAX_LEN_FNAM) uwind_errfile |
759 |
character*(MAX_LEN_FNAM) vwind_errfile |
character*(MAX_LEN_FNAM) vwind_errfile |
760 |
|
|
787 |
& frame, |
& frame, |
788 |
& cosphi, |
& cosphi, |
789 |
& whflux,wsflux,wtauu,wtauv, |
& whflux,wsflux,wtauu,wtauv, |
790 |
& watemp,waqh,wprecip,wswflux,wuwind,wvwind, |
& watemp,waqh,wprecip,wswflux,wswdown, |
791 |
|
& wuwind,wvwind, |
792 |
& wscatx,wscaty, |
& wscatx,wscaty, |
793 |
& wtheta,wtheta2,wthetaLev, |
& wtheta,wtheta2,wthetaLev, |
794 |
& wsalt,wsalt2,wsaltLev, |
& wsalt,wsalt2,wsaltLev, |
795 |
& wsst,wsss, |
& wsst,wsss, |
796 |
& wtp,wers, |
& wtp,wers,wgfo, |
797 |
& wp, |
& wp, |
798 |
& wctdt,wctds, |
& wctdt,wctds, |
799 |
& wudrift,wvdrift, |
& wudrift,wvdrift, |
819 |
_RL waqh (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL waqh (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
820 |
_RL wprecip (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wprecip (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
821 |
_RL wswflux (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wswflux (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
822 |
|
_RL wswdown (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
823 |
_RL wuwind (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wuwind (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
824 |
_RL wvwind (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wvwind (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
825 |
_RL wtheta ( nr,nsx,nsy) |
_RL wtheta ( nr,nsx,nsy) |
832 |
_RL wsss (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wsss (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
833 |
_RL wtp (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wtp (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
834 |
_RL wers (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wers (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
835 |
|
_RL wgfo (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
836 |
_RL wp (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL wp (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
837 |
_RL wctdt ( nr,nsx,nsy) |
_RL wctdt ( nr,nsx,nsy) |
838 |
_RL wctds ( nr,nsx,nsy) |
_RL wctds ( nr,nsx,nsy) |
845 |
|
|
846 |
common /ecco_cost_weights_0_r/ |
common /ecco_cost_weights_0_r/ |
847 |
& whflux0, wsflux0, wtau0, |
& whflux0, wsflux0, wtau0, |
848 |
& watemp0, waqh0, wprecip0, wswflux0, wwind0 |
& watemp0, waqh0, wprecip0, wswflux0, wswdown0, wwind0 |
849 |
_RL whflux0 |
_RL whflux0 |
850 |
_RL wsflux0 |
_RL wsflux0 |
851 |
_RL wtau0 |
_RL wtau0 |
853 |
_RL waqh0 |
_RL waqh0 |
854 |
_RL wprecip0 |
_RL wprecip0 |
855 |
_RL wswflux0 |
_RL wswflux0 |
856 |
|
_RL wswdown0 |
857 |
_RL wwind0 |
_RL wwind0 |
858 |
|
|
859 |
common /ecco_cost_weights_mean_r/ |
common /ecco_cost_weights_mean_r/ |
860 |
& wmean_hflux, wmean_sflux, wmean_tau, |
& wmean_hflux, wmean_sflux, wmean_tau, |
861 |
& wmean_atemp, wmean_aqh, |
& wmean_atemp, wmean_aqh, |
862 |
& wmean_precip, wmean_swflux, wmean_wind |
& wmean_precip, wmean_swflux, wmean_swdown, wmean_wind |
863 |
_RL wmean_hflux |
_RL wmean_hflux |
864 |
_RL wmean_sflux |
_RL wmean_sflux |
865 |
_RL wmean_tau |
_RL wmean_tau |
867 |
_RL wmean_aqh |
_RL wmean_aqh |
868 |
_RL wmean_precip |
_RL wmean_precip |
869 |
_RL wmean_swflux |
_RL wmean_swflux |
870 |
|
_RL wmean_swdown |
871 |
_RL wmean_wind |
_RL wmean_wind |
872 |
|
|
873 |
common /ecco_cost_weights_2_r/ |
common /ecco_cost_weights_2_r/ |
961 |
& tpmask, |
& tpmask, |
962 |
& ersobs, |
& ersobs, |
963 |
& ersmask, |
& ersmask, |
964 |
|
& gfoobs, |
965 |
|
& gfomask, |
966 |
& ctdtobs, |
& ctdtobs, |
967 |
& ctdsobs, |
& ctdsobs, |
968 |
& xbtobs, |
& xbtobs, |
993 |
_RL tpmask (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL tpmask (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
994 |
_RL ersobs (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL ersobs (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
995 |
_RL ersmask (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
_RL ersmask (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
996 |
|
_RL gfoobs (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
997 |
|
_RL gfomask (1-olx:snx+olx,1-oly:sny+oly, nsx,nsy) |
998 |
_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) |
999 |
_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) |
1000 |
_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) |
1040 |
& topexmeanfile, |
& topexmeanfile, |
1041 |
& topexfile, |
& topexfile, |
1042 |
& ersfile, |
& ersfile, |
1043 |
|
& gfofile, |
1044 |
& ctdtfile, |
& ctdtfile, |
1045 |
& ctdsfile, |
& ctdsfile, |
1046 |
& ctdtclimfile, |
& ctdtclimfile, |
1063 |
character*(MAX_LEN_FNAM) topexmeanfile |
character*(MAX_LEN_FNAM) topexmeanfile |
1064 |
character*(MAX_LEN_FNAM) topexfile |
character*(MAX_LEN_FNAM) topexfile |
1065 |
character*(MAX_LEN_FNAM) ersfile |
character*(MAX_LEN_FNAM) ersfile |
1066 |
|
character*(MAX_LEN_FNAM) gfofile |
1067 |
character*(MAX_LEN_FNAM) ctdtfile |
character*(MAX_LEN_FNAM) ctdtfile |
1068 |
character*(MAX_LEN_FNAM) ctdsfile |
character*(MAX_LEN_FNAM) ctdsfile |
1069 |
character*(MAX_LEN_FNAM) ctdtclimfile |
character*(MAX_LEN_FNAM) ctdtclimfile |
1085 |
|
|
1086 |
common /ecco_cost_data_flags/ |
common /ecco_cost_data_flags/ |
1087 |
& using_topex, |
& using_topex, |
1088 |
& using_ers |
& using_ers, |
1089 |
|
& using_gfo |
1090 |
logical using_topex |
logical using_topex |
1091 |
logical using_ers |
logical using_ers |
1092 |
|
logical using_gfo |
1093 |
|
|
1094 |
c Calendar information for the observations: |
c Calendar information for the observations: |
1095 |
c ========================================== |
c ========================================== |
1109 |
& tmistartdate, |
& tmistartdate, |
1110 |
& sssstartdate, |
& sssstartdate, |
1111 |
& topexstartdate, |
& topexstartdate, |
1112 |
& ersstartdate |
& ersstartdate, |
1113 |
|
& gfostartdate |
1114 |
integer scatxstartdate(4) |
integer scatxstartdate(4) |
1115 |
integer scatystartdate(4) |
integer scatystartdate(4) |
1116 |
integer sststartdate(4) |
integer sststartdate(4) |
1120 |
integer sssstartdate(4) |
integer sssstartdate(4) |
1121 |
integer topexstartdate(4) |
integer topexstartdate(4) |
1122 |
integer ersstartdate(4) |
integer ersstartdate(4) |
1123 |
|
integer gfostartdate(4) |
1124 |
|
|
1125 |
common /ecco_cost_data_aux_i/ |
common /ecco_cost_data_aux_i/ |
1126 |
& tmistartdate1, |
& tmistartdate1, |
1137 |
& topexstartdate2, |
& topexstartdate2, |
1138 |
& ersstartdate1, |
& ersstartdate1, |
1139 |
& ersstartdate2, |
& ersstartdate2, |
1140 |
|
& gfostartdate1, |
1141 |
|
& gfostartdate2, |
1142 |
& scatstartdate1, |
& scatstartdate1, |
1143 |
& scatstartdate2, |
& scatstartdate2 |
|
& costIceStart1, |
|
|
& costIceStart2, |
|
|
& costIceEnd1, |
|
|
& costIceEnd2 |
|
1144 |
|
|
1145 |
integer tmistartdate1 |
integer tmistartdate1 |
1146 |
integer tmistartdate2 |
integer tmistartdate2 |
1156 |
integer topexstartdate2 |
integer topexstartdate2 |
1157 |
integer ersstartdate1 |
integer ersstartdate1 |
1158 |
integer ersstartdate2 |
integer ersstartdate2 |
1159 |
|
integer gfostartdate1 |
1160 |
|
integer gfostartdate2 |
1161 |
integer scatstartdate1 |
integer scatstartdate1 |
1162 |
integer scatstartdate2 |
integer scatstartdate2 |
|
integer costIceStart1 |
|
|
integer costIceStart2 |
|
|
integer costIceEnd1 |
|
|
integer costIceEnd2 |
|
1163 |
|
|
1164 |
common /ecco_cost_data_times_r/ |
common /ecco_cost_data_times_r/ |
1165 |
& topexperiod, |
& topexperiod, |
1166 |
& ersperiod, |
& ersperiod, |
1167 |
& scatperiod, |
& gfoperiod, |
1168 |
& costIceStart, |
& scatperiod |
|
& costIceEnd |
|
1169 |
_RL topexperiod |
_RL topexperiod |
1170 |
_RL ersperiod |
_RL ersperiod |
1171 |
|
_RL gfoperiod |
1172 |
_RL scatperiod |
_RL scatperiod |
|
_RL costIceStart |
|
|
_RL costIceEnd |
|
|
|
|
|
c cost_ice_flag - cost_ice flag (see cost_ice.F) |
|
1173 |
|
|
1174 |
common /ecco_cost_ice_i/ cost_ice_flag |
common /ecco_cost_data_detrend/ |
1175 |
integer cost_ice_flag |
& topexintercept, |
1176 |
|
& ersintercept, |
1177 |
|
& gfointercept, |
1178 |
|
& topexslope, |
1179 |
|
& ersslope, |
1180 |
|
& gfoslope |
1181 |
|
_RL topexintercept |
1182 |
|
_RL ersintercept |
1183 |
|
_RL gfointercept |
1184 |
|
_RL topexslope |
1185 |
|
_RL ersslope |
1186 |
|
_RL gfoslope |
1187 |
|
|
1188 |
c ================================================================== |
c ================================================================== |
1189 |
c END OF HEADER COST |
c END OF HEADER COST |