/[MITgcm]/MITgcm/pkg/ecco/ecco_cost.h
ViewVC logotype

Diff of /MITgcm/pkg/ecco/ecco_cost.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.5 by heimbach, Fri May 27 22:10:27 2005 UTC revision 1.17 by heimbach, Mon Mar 6 19:16:53 2006 UTC
# Line 8  c       fields and counters associated w Line 8  c       fields and counters associated w
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     ==================================================================
# Line 274  c     objf_scatxm - time-mean zonal SCAT Line 272  c     objf_scatxm - time-mean zonal SCAT
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
# Line 287  c                  function contribution Line 284  c                  function contribution
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,
# Line 320  c                  function contribution Line 323  c                  function contribution
323       &                objf_scatym,       &                objf_scatym,
324       &                objf_atemp,       &                objf_atemp,
325       &                objf_aqh,       &                objf_aqh,
326         &                objf_precip,
327         &                objf_swflux,
328         &                objf_swdown,
329       &                objf_uwind,       &                objf_uwind,
330       &                objf_vwind,       &                objf_vwind,
331         &                objf_atempm,
332         &                objf_aqhm,
333         &                objf_precipm,
334         &                objf_swfluxm,
335         &                objf_swdownm,
336         &                objf_uwindm,
337         &                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,
# Line 329  c                  function contribution Line 349  c                  function contribution
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,
# Line 338  c                  function contribution Line 358  c                  function contribution
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)
# Line 371  c                  function contribution Line 397  c                  function contribution
397        _RL  objf_scatym(nsx,nsy)        _RL  objf_scatym(nsx,nsy)
398        _RL  objf_atemp(nsx,nsy)        _RL  objf_atemp(nsx,nsy)
399        _RL  objf_aqh  (nsx,nsy)        _RL  objf_aqh  (nsx,nsy)
400          _RL  objf_precip(nsx,nsy)
401          _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)
406          _RL  objf_aqhm  (nsx,nsy)
407          _RL  objf_precipm(nsx,nsy)
408          _RL  objf_swfluxm(nsx,nsy)
409          _RL  objf_swdownm(nsx,nsy)
410          _RL  objf_uwindm(nsx,nsy)
411          _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)
# Line 380  c                  function contribution Line 423  c                  function contribution
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)
# Line 424  c                  function contribution Line 466  c                  function contribution
466       &                num_scatym,       &                num_scatym,
467       &                num_atemp,       &                num_atemp,
468       &                num_aqh,       &                num_aqh,
469         &                num_precip,
470         &                num_swflux,
471         &                num_swdown,
472       &                num_uwind,       &                num_uwind,
473       &                num_vwind,       &                num_vwind,
474         &                num_atempm,
475         &                num_aqhm,
476         &                num_precipm,
477         &                num_swfluxm,
478         &                num_swdownm,
479         &                num_uwindm,
480         &                num_vwindm,
481       &                num_obcsn,       &                num_obcsn,
482       &                num_obcss,       &                num_obcss,
483       &                num_obcsw,       &                num_obcsw,
# Line 433  c                  function contribution Line 485  c                  function contribution
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,
# Line 475  c                  function contribution Line 527  c                  function contribution
527        _RL  num_scatym(nsx,nsy)        _RL  num_scatym(nsx,nsy)
528        _RL  num_atemp(nsx,nsy)        _RL  num_atemp(nsx,nsy)
529        _RL  num_aqh  (nsx,nsy)        _RL  num_aqh  (nsx,nsy)
530          _RL  num_precip(nsx,nsy)
531          _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)
536          _RL  num_aqhm  (nsx,nsy)
537          _RL  num_precipm(nsx,nsy)
538          _RL  num_swfluxm(nsx,nsy)
539          _RL  num_swdownm(nsx,nsy)
540          _RL  num_uwindm(nsx,nsy)
541          _RL  num_vwindm(nsx,nsy)
542        _RL  num_obcsn(nsx,nsy)        _RL  num_obcsn(nsx,nsy)
543        _RL  num_obcss(nsx,nsy)        _RL  num_obcss(nsx,nsy)
544        _RL  num_obcsw(nsx,nsy)        _RL  num_obcsw(nsx,nsy)
# Line 484  c                  function contribution Line 546  c                  function contribution
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)
# Line 494  c                  function contribution Line 555  c                  function contribution
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,
# Line 520  c                  function contribution Line 583  c                  function contribution
583       &                    mult_scaty,       &                    mult_scaty,
584       &                    mult_atemp,       &                    mult_atemp,
585       &                    mult_aqh,       &                    mult_aqh,
586         &                    mult_precip,
587         &                    mult_swflux,
588         &                    mult_swdown,
589       &                    mult_uwind,       &                    mult_uwind,
590       &                    mult_vwind,       &                    mult_vwind,
591       &                    mult_obcsn,       &                    mult_obcsn,
# Line 529  c                  function contribution Line 595  c                  function contribution
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
# Line 563  c                  function contribution Line 632  c                  function contribution
632        _RL  mult_scaty        _RL  mult_scaty
633        _RL  mult_atemp        _RL  mult_atemp
634        _RL  mult_aqh        _RL  mult_aqh
635          _RL  mult_precip
636          _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
# Line 572  c                  function contribution Line 644  c                  function contribution
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     ==========================================================
# Line 588  c                the current model integ Line 661  c                the current model integ
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:
# Line 638  c     velerrfile            - representa Line 715  c     velerrfile            - representa
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,
724         &                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
# Line 664  c     velerrfile            - representa Line 747  c     velerrfile            - representa
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
756          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    
# Line 698  c     wvdrift    - weight for mean merid Line 787  c     wvdrift    - weight for mean merid
787       &                      frame,       &                      frame,
788       &                      cosphi,       &                      cosphi,
789       &                      whflux,wsflux,wtauu,wtauv,       &                      whflux,wsflux,wtauu,wtauv,
790       &                      watemp,waqh,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,
# Line 727  c     wvdrift    - weight for mean merid Line 817  c     wvdrift    - weight for mean merid
817        _RL wscaty  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wscaty  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
818        _RL watemp  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL watemp  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
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)
821          _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)
# Line 739  c     wvdrift    - weight for mean merid Line 832  c     wvdrift    - weight for mean merid
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)
# Line 749  c     wvdrift    - weight for mean merid Line 843  c     wvdrift    - weight for mean merid
843        _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)
844        _RL wbaro        _RL wbaro
845    
846          common /ecco_cost_weights_0_r/
847         &        whflux0, wsflux0, wtau0,
848         &        watemp0, waqh0, wprecip0, wswflux0, wswdown0, wwind0
849          _RL whflux0
850          _RL wsflux0
851          _RL wtau0
852          _RL watemp0
853          _RL waqh0
854          _RL wprecip0
855          _RL wswflux0
856          _RL wswdown0
857          _RL wwind0
858    
859          common /ecco_cost_weights_mean_r/
860         &        wmean_hflux, wmean_sflux, wmean_tau,
861         &        wmean_atemp, wmean_aqh,
862         &        wmean_precip, wmean_swflux, wmean_swdown, wmean_wind
863          _RL wmean_hflux
864          _RL wmean_sflux
865          _RL wmean_tau
866          _RL wmean_atemp
867          _RL wmean_aqh
868          _RL wmean_precip
869          _RL wmean_swflux
870          _RL wmean_swdown
871          _RL wmean_wind
872                                                    
873        common /ecco_cost_weights_2_r/        common /ecco_cost_weights_2_r/
874       &                      whflux2,wsflux2,wtauu2,wtauv2       &                      whflux2,wsflux2,wtauu2,wtauv2
875        _RL whflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL whflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
# Line 840  c     vdriftdat  - drifters meridional v Line 961  c     vdriftdat  - drifters meridional v
961       &                     tpmask,       &                     tpmask,
962       &                     ersobs,       &                     ersobs,
963       &                     ersmask,       &                     ersmask,
964         &                     gfoobs,
965         &                     gfomask,
966       &                     ctdtobs,       &                     ctdtobs,
967       &                     ctdsobs,       &                     ctdsobs,
968       &                     xbtobs,       &                     xbtobs,
# Line 870  c     vdriftdat  - drifters meridional v Line 993  c     vdriftdat  - drifters meridional v
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)
# Line 915  c     driftfile     - reference data fil Line 1040  c     driftfile     - reference data fil
1040       &                     topexmeanfile,       &                     topexmeanfile,
1041       &                     topexfile,       &                     topexfile,
1042       &                     ersfile,       &                     ersfile,
1043         &                     gfofile,
1044       &                     ctdtfile,       &                     ctdtfile,
1045       &                     ctdsfile,       &                     ctdsfile,
1046       &                     ctdtclimfile,       &                     ctdtclimfile,
# Line 937  c     driftfile     - reference data fil Line 1063  c     driftfile     - reference data fil
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
# Line 958  c     using_ers - flag that indicates th Line 1085  c     using_ers - flag that indicates th
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     ==========================================
# Line 980  c     sshperiod      - sampling interval Line 1109  c     sshperiod      - sampling interval
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)
# Line 990  c     sshperiod      - sampling interval Line 1120  c     sshperiod      - sampling interval
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,
# Line 1006  c     sshperiod      - sampling interval Line 1137  c     sshperiod      - sampling interval
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
# Line 1027  c     sshperiod      - sampling interval Line 1156  c     sshperiod      - sampling interval
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
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22