/[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.1 by heimbach, Thu Nov 6 22:10:08 2003 UTC revision 1.12 by heimbach, Fri Oct 14 19:54:49 2005 UTC
# Line 7  c     o Header for averaged temperature, Line 7  c     o Header for averaged temperature,
7  c       fields and counters associated with the averaging.  c       fields and counters associated with the averaging.
8  c  c
9  c     started: Christian Eckert eckert@mit.edu  30-Jun-1999  c     started: Christian Eckert eckert@mit.edu  30-Jun-1999
 c     changed: Christian Eckert eckert@mit.edu  01-Mar-2000  
 c              - Restructured the code in order to create a package  
 c                for the MITgcmUV.  
 c     added sea-ice term: menemenlis@jpl.nasa.gov 26-Feb-2003  
10  c  c
11  c     ==================================================================  c     ==================================================================
12  c     HEADER AVERAGES  c     HEADER AVERAGES
# Line 80  c             intantaneous field. Line 76  c             intantaneous field.
76       &                    wlmean,       &                    wlmean,
77       &                    Sfmean,       &                    Sfmean,
78       &                    Tfmean,                         &                    Tfmean,                  
79       &                    wfmean,       &                    sbar_gen,
80       &                    sbar_gen,       &                    tbar_gen,
81       &                    tbar_gen       &                    wfmean
82    
83  #if (defined (ALLOW_THETA_COST_CONTRIBUTION) || \  #if (defined (ALLOW_THETA_COST_CONTRIBUTION) || \
      defined (ALLOW_SST_COST_CONTRIBUTION) || \  
84       defined (ALLOW_CTDT_COST_CONTRIBUTION) || \       defined (ALLOW_CTDT_COST_CONTRIBUTION) || \
      defined (ALLOW_CTDTCLIM_COST_CONTRIBUTION) || \  
85       defined (ALLOW_XBT_COST_CONTRIBUTION) || \       defined (ALLOW_XBT_COST_CONTRIBUTION) || \
86       defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \       defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \
87       defined (ALLOW_OBCS_COST_CONTRIBUTION))       defined (ALLOW_OBCS_COST_CONTRIBUTION))
88        _RL tbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL tbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
89  #else  #else
90    #ifdef ALLOW_SST_COST_CONTRIBUTION
91          _RL tbar  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
92    #else
93        _RL tbar        _RL tbar
94  #endif  #endif
95    #endif
96    #ifdef GENERIC_BAR_MONTH
97          _RL tbar_gen  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
98          _RL sbar_gen  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
99    #else
100          _RL tbar_gen
101          _RL sbar_gen
102    #endif
103    
104  #if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \  #if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \
      defined (ALLOW_SSS_COST_CONTRIBUTION) || \  
105       defined (ALLOW_CTDS_COST_CONTRIBUTION) || \       defined (ALLOW_CTDS_COST_CONTRIBUTION) || \
      defined (ALLOW_CTDSCLIM_COST_CONTRIBUTION) || \  
106       defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \       defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \
107       defined (ALLOW_OBCS_COST_CONTRIBUTION))       defined (ALLOW_OBCS_COST_CONTRIBUTION))
108        _RL sbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL sbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
109  #else  #else
110    #ifdef ALLOW_SSS_COST_CONTRIBUTION
111          _RL sbar  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
112    #else
113        _RL sbar        _RL sbar
114  #endif  #endif
   
 #ifdef GENERIC_BAR_MONTH  
       _RL tbar_gen  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
       _RL sbar_gen  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
 #else  
       _RL tbar_gen  
       _RL sbar_gen  
115  #endif  #endif
116    
117  #ifdef ALLOW_SSH_COST_CONTRIBUTION  #ifdef ALLOW_SSH_COST_CONTRIBUTION
# Line 122  c             intantaneous field. Line 121  c             intantaneous field.
121  #endif  #endif
122    
123  #if (defined (ALLOW_DRIFTER_COST_CONTRIBUTION) || \  #if (defined (ALLOW_DRIFTER_COST_CONTRIBUTION) || \
124       defined (ALLOW_CURMTR_COST_CONTRIBUTION))       defined (ALLOW_OBCS_COST_CONTRIBUTION))
125        _RL ubar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL ubar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
126        _RL vbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL vbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
127  #else  #else
# Line 259  c     objf_temp  - Temperature contribut Line 258  c     objf_temp  - Temperature contribut
258  c     objf_salt  - Salinity contribution.  c     objf_salt  - Salinity contribution.
259  c     objf_temp0 - Initial conditions Temperature contribution.  c     objf_temp0 - Initial conditions Temperature contribution.
260  c     objf_salt0 - Initial conditions Salinity contribution.  c     objf_salt0 - Initial conditions Salinity contribution.
 c     objf_tmi   - Sea surface temperature contribution.  
261  c     objf_sst   - Sea surface temperature contribution.  c     objf_sst   - Sea surface temperature contribution.
262    c     objf_tmi   - Sea surface temperature contribution.
263  c     objf_sss   - Sea surface salinity contribution.  c     objf_sss   - Sea surface salinity contribution.
 c     objf_atl   - Meridional heat transport in the N-Atlantic  
264  c     objf_ctdt  - Temperature measurements from Woce CTD  c     objf_ctdt  - Temperature measurements from Woce CTD
265  c     objf_ctds  - Salinity measurements from Woce CTD  c     objf_ctds  - Salinity measurements from Woce CTD
266  c     objf_ctdtclim - Temperature from Woce CTD without timetag  c     objf_ctdtclim - Temperature measurements from Woce CTD without timetag
267  c     objf_ctdsclim - Salinity from Woce CTD without timetag  c     objf_ctdsclim - Salinity measurements from Woce CTD without timetag
268  c     objf_xbt   - XBT temperature data  c     objf_xbt   - XBT temperature data
269  c     objf_argot - ARGO temperature profiles  c     objf_argot - ARGO temperature profiles
270  c     objf_argos - ARGO salt profiles  c     objf_argos - ARGO salt profiles
# 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 300  c                  function contribution Line 297  c                  function contribution
297       &                objf_salt,       &                objf_salt,
298       &                objf_temp0,       &                objf_temp0,
299       &                objf_salt0,       &                objf_salt0,
      &                objf_tmi,  
300       &                objf_sst,       &                objf_sst,
301         &                objf_tmi,
302       &                objf_sss,       &                objf_sss,
303       &                objf_ctdt,       &                objf_ctdt,
304       &                objf_ctds,       &                objf_ctds,
# Line 320  c                  function contribution Line 317  c                  function contribution
317       &                objf_scatym,       &                objf_scatym,
318       &                objf_atemp,       &                objf_atemp,
319       &                objf_aqh,       &                objf_aqh,
320         &                objf_precip,
321         &                objf_swflux,
322         &                objf_swdown,
323       &                objf_uwind,       &                objf_uwind,
324       &                objf_vwind,       &                objf_vwind,
325         &                objf_atempm,
326         &                objf_aqhm,
327         &                objf_precipm,
328         &                objf_swfluxm,
329         &                objf_swdownm,
330         &                objf_uwindm,
331         &                objf_vwindm,
332       &                objf_obcsn,       &                objf_obcsn,
333       &                objf_obcss,       &                objf_obcss,
334       &                objf_obcsw,       &                objf_obcsw,
# Line 329  c                  function contribution Line 336  c                  function contribution
336       &                objf_obcsvol,       &                objf_obcsvol,
337       &                objf_curmtr,       &                objf_curmtr,
338       &                objf_ageos,       &                objf_ageos,
339       &                objf_ice,       &                objf_kapgm,
340       &                objf_diffkr,       &                objf_diffkr,
341       &                objf_theta_ini_fin,       &                objf_theta_ini_fin,
342       &                objf_salt_ini_fin       &                objf_salt_ini_fin,
343        _RL  objf_hmean       &                objf_eddytau
344        _RL  objf_obcsvol  
345        _RL  objf_hflux  (nsx,nsy)        _RL  objf_hflux  (nsx,nsy)
346        _RL  objf_hfluxm (nsx,nsy)        _RL  objf_hfluxm (nsx,nsy)
347        _RL  objf_hfluxmm(nsx,nsy)        _RL  objf_hfluxmm(nsx,nsy)
# Line 345  c                  function contribution Line 352  c                  function contribution
352        _RL  objf_tauum  (nsx,nsy)        _RL  objf_tauum  (nsx,nsy)
353        _RL  objf_tauv   (nsx,nsy)        _RL  objf_tauv   (nsx,nsy)
354        _RL  objf_tauvm  (nsx,nsy)        _RL  objf_tauvm  (nsx,nsy)
355          _RL  objf_hmean
356        _RL  objf_h    (nsx,nsy)        _RL  objf_h    (nsx,nsy)
357        _RL  objf_temp (nsx,nsy)        _RL  objf_temp (nsx,nsy)
358        _RL  objf_salt (nsx,nsy)        _RL  objf_salt (nsx,nsy)
359        _RL  objf_temp0(nsx,nsy)        _RL  objf_temp0(nsx,nsy)
360        _RL  objf_salt0(nsx,nsy)        _RL  objf_salt0(nsx,nsy)
       _RL  objf_tmi  (nsx,nsy)  
361        _RL  objf_sst  (nsx,nsy)        _RL  objf_sst  (nsx,nsy)
362          _RL  objf_tmi  (nsx,nsy)
363        _RL  objf_sss  (nsx,nsy)        _RL  objf_sss  (nsx,nsy)
364        _RL  objf_ctdt (nsx,nsy)        _RL  objf_ctdt (nsx,nsy)
365        _RL  objf_ctds (nsx,nsy)        _RL  objf_ctds (nsx,nsy)
# Line 370  c                  function contribution Line 378  c                  function contribution
378        _RL  objf_scatym(nsx,nsy)        _RL  objf_scatym(nsx,nsy)
379        _RL  objf_atemp(nsx,nsy)        _RL  objf_atemp(nsx,nsy)
380        _RL  objf_aqh  (nsx,nsy)        _RL  objf_aqh  (nsx,nsy)
381          _RL  objf_precip(nsx,nsy)
382          _RL  objf_swflux(nsx,nsy)
383          _RL  objf_swdown(nsx,nsy)
384        _RL  objf_uwind(nsx,nsy)        _RL  objf_uwind(nsx,nsy)
385        _RL  objf_vwind(nsx,nsy)        _RL  objf_vwind(nsx,nsy)
386          _RL  objf_atempm(nsx,nsy)
387          _RL  objf_aqhm  (nsx,nsy)
388          _RL  objf_precipm(nsx,nsy)
389          _RL  objf_swfluxm(nsx,nsy)
390          _RL  objf_swdownm(nsx,nsy)
391          _RL  objf_uwindm(nsx,nsy)
392          _RL  objf_vwindm(nsx,nsy)
393        _RL  objf_obcsn(nsx,nsy)        _RL  objf_obcsn(nsx,nsy)
394        _RL  objf_obcss(nsx,nsy)        _RL  objf_obcss(nsx,nsy)
395        _RL  objf_obcsw(nsx,nsy)        _RL  objf_obcsw(nsx,nsy)
396        _RL  objf_obcse(nsx,nsy)        _RL  objf_obcse(nsx,nsy)
397          _RL  objf_obcsvol
398        _RL  objf_curmtr(nsx,nsy)        _RL  objf_curmtr(nsx,nsy)
399        _RL  objf_ageos(nsx,nsy)        _RL  objf_ageos(nsx,nsy)
       _RL  objf_ice  (nsx,nsy)  
400        _RL  objf_kapgm(nsx,nsy)        _RL  objf_kapgm(nsx,nsy)
401        _RL  objf_diffkr(nsx,nsy)        _RL  objf_diffkr(nsx,nsy)
402        _RL  objf_theta_ini_fin(nsx,nsy)        _RL  objf_theta_ini_fin(nsx,nsy)
403        _RL  objf_salt_ini_fin(nsx,nsy)        _RL  objf_salt_ini_fin(nsx,nsy)
404          _RL  objf_eddytau(nsx,nsy)
405    
406          common /ecco_cost_num/
407         &                num_hflux,
408         &                num_hfluxm,
409         &                num_hfluxmm,
410         &                num_sflux,
411         &                num_sfluxm,
412         &                num_sfluxmm,
413         &                num_tauu,
414         &                num_tauum,
415         &                num_tauv,
416         &                num_tauvm,
417         &                num_hmean,
418         &                num_h,
419         &                num_temp,
420         &                num_salt,
421         &                num_temp0,
422         &                num_salt0,
423         &                num_sst,
424         &                num_tmi,
425         &                num_sss,
426         &                num_ctdt,
427         &                num_ctds,
428         &                num_ctdtclim,
429         &                num_ctdsclim,
430         &                num_xbt,
431         &                num_argot,
432         &                num_argos,
433         &                num_drift,
434         &                num_tdrift,
435         &                num_sdrift,
436         &                num_wdrift,
437         &                num_scatx,
438         &                num_scaty,
439         &                num_scatxm,
440         &                num_scatym,
441         &                num_atemp,
442         &                num_aqh,
443         &                num_precip,
444         &                num_swflux,
445         &                num_swdown,
446         &                num_uwind,
447         &                num_vwind,
448         &                num_atempm,
449         &                num_aqhm,
450         &                num_precipm,
451         &                num_swfluxm,
452         &                num_swdownm,
453         &                num_uwindm,
454         &                num_vwindm,
455         &                num_obcsn,
456         &                num_obcss,
457         &                num_obcsw,
458         &                num_obcse,
459         &                num_obcsvol,
460         &                num_curmtr,
461         &                num_ageos,
462         &                num_kapgm,
463         &                num_diffkr,
464         &                num_theta_ini_fin,
465         &                num_salt_ini_fin,
466         &                num_eddytau
467    
468          _RL  num_hflux  (nsx,nsy)
469          _RL  num_hfluxm (nsx,nsy)
470          _RL  num_hfluxmm(nsx,nsy)
471          _RL  num_sflux  (nsx,nsy)
472          _RL  num_sfluxm (nsx,nsy)
473          _RL  num_sfluxmm(nsx,nsy)
474          _RL  num_tauu   (nsx,nsy)
475          _RL  num_tauum  (nsx,nsy)
476          _RL  num_tauv   (nsx,nsy)
477          _RL  num_tauvm  (nsx,nsy)
478          _RL  num_hmean
479          _RL  num_h    (nsx,nsy)
480          _RL  num_temp (nsx,nsy)
481          _RL  num_salt (nsx,nsy)
482          _RL  num_temp0(nsx,nsy)
483          _RL  num_salt0(nsx,nsy)
484          _RL  num_sst  (nsx,nsy)
485          _RL  num_tmi  (nsx,nsy)
486          _RL  num_sss  (nsx,nsy)
487          _RL  num_ctdt (nsx,nsy)
488          _RL  num_ctds (nsx,nsy)
489          _RL  num_ctdtclim (nsx,nsy)
490          _RL  num_ctdsclim (nsx,nsy)
491          _RL  num_xbt  (nsx,nsy)
492          _RL  num_argot(nsx,nsy)
493          _RL  num_argos(nsx,nsy)
494          _RL  num_drift(nsx,nsy)
495          _RL  num_tdrift(nsx,nsy)
496          _RL  num_sdrift(nsx,nsy)
497          _RL  num_wdrift(nsx,nsy)
498          _RL  num_scatx(nsx,nsy)
499          _RL  num_scaty(nsx,nsy)
500          _RL  num_scatxm(nsx,nsy)
501          _RL  num_scatym(nsx,nsy)
502          _RL  num_atemp(nsx,nsy)
503          _RL  num_aqh  (nsx,nsy)
504          _RL  num_precip(nsx,nsy)
505          _RL  num_swflux(nsx,nsy)
506          _RL  num_swdown(nsx,nsy)
507          _RL  num_uwind(nsx,nsy)
508          _RL  num_vwind(nsx,nsy)
509          _RL  num_atempm(nsx,nsy)
510          _RL  num_aqhm  (nsx,nsy)
511          _RL  num_precipm(nsx,nsy)
512          _RL  num_swfluxm(nsx,nsy)
513          _RL  num_swdownm(nsx,nsy)
514          _RL  num_uwindm(nsx,nsy)
515          _RL  num_vwindm(nsx,nsy)
516          _RL  num_obcsn(nsx,nsy)
517          _RL  num_obcss(nsx,nsy)
518          _RL  num_obcsw(nsx,nsy)
519          _RL  num_obcse(nsx,nsy)
520          _RL  num_obcsvol
521          _RL  num_curmtr(nsx,nsy)
522          _RL  num_ageos(nsx,nsy)
523          _RL  num_kapgm(nsx,nsy)
524          _RL  num_diffkr(nsx,nsy)
525          _RL  num_theta_ini_fin(nsx,nsy)
526          _RL  num_salt_ini_fin(nsx,nsy)
527          _RL  num_eddytau(nsx,nsy)
528    
529        common /ecco_cost_aux_r/        common /ecco_cost_aux_r/
530       &                    mult_hflux,       &                    mult_hflux,
# Line 395  c                  function contribution Line 537  c                  function contribution
537       &                    mult_salt,       &                    mult_salt,
538       &                    mult_temp0,       &                    mult_temp0,
539       &                    mult_salt0,       &                    mult_salt0,
      &                    mult_tmi,  
540       &                    mult_sst,       &                    mult_sst,
541         &                    mult_tmi,
542       &                    mult_sss,       &                    mult_sss,
543       &                    mult_ctdt,       &                    mult_ctdt,
544       &                    mult_ctds,       &                    mult_ctds,
# Line 413  c                  function contribution Line 555  c                  function contribution
555       &                    mult_scaty,       &                    mult_scaty,
556       &                    mult_atemp,       &                    mult_atemp,
557       &                    mult_aqh,       &                    mult_aqh,
558         &                    mult_precip,
559         &                    mult_swflux,
560         &                    mult_swdown,
561       &                    mult_uwind,       &                    mult_uwind,
562       &                    mult_vwind,       &                    mult_vwind,
563       &                    mult_obcsn,       &                    mult_obcsn,
# Line 422  c                  function contribution Line 567  c                  function contribution
567       &                    mult_obcsvol,       &                    mult_obcsvol,
568       &                    mult_curmtr,       &                    mult_curmtr,
569       &                    mult_ageos,       &                    mult_ageos,
      &                    mult_ice,  
570       &                    mult_kapgm,       &                    mult_kapgm,
571       &                    mult_diffkr,       &                    mult_diffkr,
572       &                    mult_ini_fin       &                    mult_ini_fin,
573         &                    mult_eddytau
574    
575        _RL  mult_hflux        _RL  mult_hflux
576        _RL  mult_sflux        _RL  mult_sflux
577        _RL  mult_tauu        _RL  mult_tauu
# Line 436  c                  function contribution Line 582  c                  function contribution
582        _RL  mult_salt        _RL  mult_salt
583        _RL  mult_temp0        _RL  mult_temp0
584        _RL  mult_salt0        _RL  mult_salt0
       _RL  mult_tmi  
585        _RL  mult_sst        _RL  mult_sst
586          _RL  mult_tmi
587        _RL  mult_sss        _RL  mult_sss
588        _RL  mult_ctdt        _RL  mult_ctdt
589        _RL  mult_ctds        _RL  mult_ctds
# Line 454  c                  function contribution Line 600  c                  function contribution
600        _RL  mult_scaty        _RL  mult_scaty
601        _RL  mult_atemp        _RL  mult_atemp
602        _RL  mult_aqh        _RL  mult_aqh
603          _RL  mult_precip
604          _RL  mult_swflux
605          _RL  mult_swdown
606        _RL  mult_uwind        _RL  mult_uwind
607        _RL  mult_vwind        _RL  mult_vwind
608        _RL  mult_obcsn        _RL  mult_obcsn
# Line 463  c                  function contribution Line 612  c                  function contribution
612        _RL  mult_obcsvol        _RL  mult_obcsvol
613        _RL  mult_curmtr        _RL  mult_curmtr
614        _RL  mult_ageos        _RL  mult_ageos
       _RL  mult_ice  
615        _RL  mult_kapgm        _RL  mult_kapgm
616        _RL  mult_diffkr        _RL  mult_diffkr
617        _RL  mult_ini_fin        _RL  mult_ini_fin
618          _RL  mult_eddytau
619    
620  c     Record counters relevant for the cost function evaluation.  c     Record counters relevant for the cost function evaluation.
621  c     ==========================================================  c     ==========================================================
# Line 479  c                the current model integ Line 627  c                the current model integ
627    
628        common /ecco_cost_i/        common /ecco_cost_i/
629       &                nmonsrec,       &                nmonsrec,
630       &                ndaysrec       &                ndaysrec,
631         &                nnztbar,
632         &                nnzsbar
633        integer nmonsrec        integer nmonsrec
634        integer ndaysrec        integer ndaysrec
635          integer nnztbar
636          integer nnzsbar
637    
638    
639  c     Data files for the weights used in the cost function:  c     Data files for the weights used in the cost function:
# Line 529  c     velerrfile            - representa Line 681  c     velerrfile            - representa
681       &                salterrfile,       &                salterrfile,
682       &                temperrfile,       &                temperrfile,
683       &                velerrfile,       &                velerrfile,
684         &                salt0errfile,
685         &                temp0errfile,
686         &                vel0errfile,
687       &                atemp_errfile,       &                atemp_errfile,
688       &                aqh_errfile,       &                aqh_errfile,
689         &                precip_errfile,
690         &                swflux_errfile,
691         &                swdown_errfile,
692       &                uwind_errfile,       &                uwind_errfile,
693       &                vwind_errfile       &                vwind_errfile
694        character*(MAX_LEN_FNAM) hflux_errfile        character*(MAX_LEN_FNAM) hflux_errfile
# Line 555  c     velerrfile            - representa Line 713  c     velerrfile            - representa
713        character*(MAX_LEN_FNAM) salterrfile        character*(MAX_LEN_FNAM) salterrfile
714        character*(MAX_LEN_FNAM) temperrfile        character*(MAX_LEN_FNAM) temperrfile
715        character*(MAX_LEN_FNAM) velerrfile        character*(MAX_LEN_FNAM) velerrfile
716          character*(MAX_LEN_FNAM) salt0errfile
717          character*(MAX_LEN_FNAM) temp0errfile
718          character*(MAX_LEN_FNAM) vel0errfile
719        character*(MAX_LEN_FNAM) atemp_errfile        character*(MAX_LEN_FNAM) atemp_errfile
720        character*(MAX_LEN_FNAM) aqh_errfile        character*(MAX_LEN_FNAM) aqh_errfile
721          character*(MAX_LEN_FNAM) precip_errfile
722          character*(MAX_LEN_FNAM) swflux_errfile
723          character*(MAX_LEN_FNAM) swdown_errfile
724        character*(MAX_LEN_FNAM) uwind_errfile        character*(MAX_LEN_FNAM) uwind_errfile
725        character*(MAX_LEN_FNAM) vwind_errfile        character*(MAX_LEN_FNAM) vwind_errfile
726    
# Line 584  c     wctdt      - weight for CTD temper Line 748  c     wctdt      - weight for CTD temper
748  c     wctds      - weight for CTD salinity.  c     wctds      - weight for CTD salinity.
749  c     wudrift    - weight for mean zonal velocity from drifters.  c     wudrift    - weight for mean zonal velocity from drifters.
750  c     wvdrift    - weight for mean meridional velocity from drifters.  c     wvdrift    - weight for mean meridional velocity from drifters.
 c     wcurrent   - weight for  
 c     wbaro      - weight for barotropic velocity adjustments.  
751    
752        common /ecco_cost_weights_r/        common /ecco_cost_weights_r/
753       &                      frame,       &                      frame,
754       &                      cosphi,       &                      cosphi,
755       &                      whflux,wsflux,wtauu,wtauv,       &                      whflux,wsflux,wtauu,wtauv,
756       &                      watemp,waqh,wuwind,wvwind,       &                      watemp,waqh,wprecip,wswflux,wswdown,
757         &                      wuwind,wvwind,
758       &                      wscatx,wscaty,       &                      wscatx,wscaty,
759       &                      wtheta,wtheta2,wthetaLev,       &                      wtheta,wtheta2,wthetaLev,
760       &                      wsalt,wsalt2,wsaltLev,       &                      wsalt,wsalt2,wsaltLev,
# Line 620  c     wbaro      - weight for barotropic Line 783  c     wbaro      - weight for barotropic
783        _RL wscaty  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wscaty  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
784        _RL watemp  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL watemp  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
785        _RL waqh    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL waqh    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
786          _RL wprecip (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
787          _RL wswflux (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
788          _RL wswdown (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
789        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
790        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
791        _RL wtheta  (                            nr,nsx,nsy)        _RL wtheta  (                            nr,nsx,nsy)
# Line 642  c     wbaro      - weight for barotropic Line 808  c     wbaro      - weight for barotropic
808        _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)
809        _RL wbaro        _RL wbaro
810    
811          common /ecco_cost_weights_0_r/
812         &        whflux0, wsflux0, wtau0,
813         &        watemp0, waqh0, wprecip0, wswflux0, wswdown0, wwind0
814          _RL whflux0
815          _RL wsflux0
816          _RL wtau0
817          _RL watemp0
818          _RL waqh0
819          _RL wprecip0
820          _RL wswflux0
821          _RL wswdown0
822          _RL wwind0
823    
824          common /ecco_cost_weights_mean_r/
825         &        wmean_hflux, wmean_sflux, wmean_tau,
826         &        wmean_atemp, wmean_aqh,
827         &        wmean_precip, wmean_swflux, wmean_swdown, wmean_wind
828          _RL wmean_hflux
829          _RL wmean_sflux
830          _RL wmean_tau
831          _RL wmean_atemp
832          _RL wmean_aqh
833          _RL wmean_precip
834          _RL wmean_swflux
835          _RL wmean_swdown
836          _RL wmean_wind
837                                                    
838          common /ecco_cost_weights_2_r/
839         &                      whflux2,wsflux2,wtauu2,wtauv2
840          _RL whflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
841          _RL wsflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
842          _RL wtauu2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
843          _RL wtauv2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
844                                                    
845                          
846    
847  #if (defined (ALLOW_OBCSN_COST_CONTRIBUTION) || \  #if (defined (ALLOW_OBCSN_COST_CONTRIBUTION) || \
848       defined (ALLOW_OBCSN_CONTROL))       defined (ALLOW_OBCSN_CONTROL))
849        common /ecco_cost_weights_obcsn/        common /ecco_cost_weights_obcsn/
# Line 671  c     wbaro      - weight for barotropic Line 873  c     wbaro      - weight for barotropic
873        _RL wobcseLev  (1-oly:sny+oly,nr,nsx,nsy,nobcs)        _RL wobcseLev  (1-oly:sny+oly,nr,nsx,nsy,nobcs)
874  #endif  #endif
875    
876    
877  c     Arrays that contain observations for the model-data comparison:  c     Arrays that contain observations for the model-data comparison:
878  c     ===============================================================  c     ===============================================================
879  c  c
880  c     tdat       - reference temperature data.  c     tdat       - reference temperature data.
881  c     scatxdat   - reference zonal wind stress.  c     scatxdat   - reference zonal wind stress.
882  c     scatydat   - reference meridional wind stress.  c     scatydat   - reference meridional wind stress.
 c     tmidat     - reference TMI sea surface temperature data.  
883  c     sstdat     - reference sea surface temperature data.  c     sstdat     - reference sea surface temperature data.
884    c     tmidat     - reference TMI sea surface temperature data.
885  c     sssdat     - reference sea surface temperature data.  c     sssdat     - reference sea surface temperature data.
886  c     tauxmask   - mask for reference wind stress data.  c     tauxmask   - mask for reference wind stress data.
887  c     tauymask   - mask for reference wind stress data.  c     tauymask   - mask for reference wind stress data.
888  c     scatxmask  - mask for scat wind stress data.  c     scatxmask  - mask for scat wind stress data.
889  c     scatymask  - mask for scat wind stress data.  c     scatymask  - mask for scat wind stress data.
 c     tmimask    - mask for reference sea surface temperature data.  
890  c     sstmask    - mask for reference sea surface temperature data.  c     sstmask    - mask for reference sea surface temperature data.
891    c     tmimask    - mask for reference sea surface temperature data.
892  c     sssmask    - mask for reference sea surface temperature data.  c     sssmask    - mask for reference sea surface temperature data.
893  c     sdat       - reference salinity data.  c     sdat       - reference salinity data.
894  c     tpmean     - reference mean sea surface height data.  c     tpmean     - reference mean sea surface height data.
# Line 706  c     vdriftdat  - drifters meridional v Line 909  c     vdriftdat  - drifters meridional v
909       &                     tdat,       &                     tdat,
910       &                     scatxdat,       &                     scatxdat,
911       &                     scatydat,       &                     scatydat,
      &                     tmidat,  
912       &                     sstdat,       &                     sstdat,
913         &                     tmidat,
914       &                     sssdat,       &                     sssdat,
      &                     tmimask,  
915       &                     sstmask,       &                     sstmask,
916         &                     tmimask,
917       &                     sssmask,       &                     sssmask,
918       &                     tauxmask,       &                     tauxmask,
919       &                     tauymask,       &                     tauymask,
# Line 736  c     vdriftdat  - drifters meridional v Line 939  c     vdriftdat  - drifters meridional v
939        _RL tdat      (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL tdat      (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
940        _RL scatxdat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatxdat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
941        _RL scatydat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatydat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
       _RL tmidat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)  
942        _RL sstdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sstdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
943          _RL tmidat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
944        _RL sssdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sssdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
945        _RL tauxmask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tauxmask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
946        _RL tauymask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tauymask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
947        _RL scatxmask (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatxmask (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
948        _RL scatymask (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatymask (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
       _RL tmimask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)  
949        _RL sstmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sstmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
950          _RL tmimask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
951        _RL sssmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sssmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
952        _RL sdat      (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL sdat      (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
953        _RL tpmean    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpmean    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
# Line 763  c     vdriftdat  - drifters meridional v Line 966  c     vdriftdat  - drifters meridional v
966        _RL curmtruobs(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL curmtruobs(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
967        _RL curmtrvobs(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL curmtrvobs(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
968    
969    
970  c     Files that contain obervations:  c     Files that contain obervations:
971  c     ===============================  c     ===============================
972  c  c
# Line 770  c     tdatfile      - reference data fil Line 974  c     tdatfile      - reference data fil
974  c     sdatfile      - reference data file for salinity.  c     sdatfile      - reference data file for salinity.
975  c     scatxdatfile  - reference data file for zonal wind stress.  c     scatxdatfile  - reference data file for zonal wind stress.
976  c     scatydatfile  - reference data file for meridional wind stress.  c     scatydatfile  - reference data file for meridional wind stress.
 c     tmidatfile    - reference data file for TMI sea surface temp.  
977  c     sstdatfile    - reference data file for sea surface temperature.  c     sstdatfile    - reference data file for sea surface temperature.
978    c     tmidatfile    - reference data file for TMI sea surface temperature.
979  c     topexmeanfile - reference data file for mean sea surface height.  c     topexmeanfile - reference data file for mean sea surface height.
980  c     topexfile     - reference data file for sea surface height data  c     topexfile     - reference data file for sea surface height data
981  c                     (TOPEX/POSEIDON).  c                     (TOPEX/POSEIDON).
# Line 780  c                     (ERS). Line 984  c                     (ERS).
984  c ctdtfile, ctdsfile- reference data file for temperature and salinity  c ctdtfile, ctdsfile- reference data file for temperature and salinity
985  c                     from CTD  c                     from CTD
986  c ctdtclimfile, ctdsclimfile- reference data file for temperature  c ctdtclimfile, ctdsclimfile- reference data file for temperature
987  c                     and salinity from CTD with out timetag  c                      and salinity from CTD with out timetag
988  c     xbtfile       - reference data file for xbt  c     xbtfile       - reference data file for xbt
989  c     ARGOtfile     - reference data file for ARGO  c     ARGOtfile     - reference data file for ARGO
990  c     ARGOsfile     - reference data file for ARGO  c     ARGOsfile     - reference data file for ARGO
# Line 791  c     driftfile     - reference data fil Line 995  c     driftfile     - reference data fil
995       &                     sdatfile,       &                     sdatfile,
996       &                     scatxdatfile,       &                     scatxdatfile,
997       &                     scatydatfile,       &                     scatydatfile,
      &                     tmidatfile,  
998       &                     sstdatfile,       &                     sstdatfile,
999         &                     tmidatfile,
1000       &                     sssdatfile,       &                     sssdatfile,
1001       &                     topexmeanfile,       &                     topexmeanfile,
1002       &                     topexfile,       &                     topexfile,
# Line 813  c     driftfile     - reference data fil Line 1017  c     driftfile     - reference data fil
1017        character*(MAX_LEN_FNAM) sdatfile        character*(MAX_LEN_FNAM) sdatfile
1018        character*(MAX_LEN_FNAM) scatxdatfile        character*(MAX_LEN_FNAM) scatxdatfile
1019        character*(MAX_LEN_FNAM) scatydatfile        character*(MAX_LEN_FNAM) scatydatfile
       character*(MAX_LEN_FNAM) tmidatfile  
1020        character*(MAX_LEN_FNAM) sstdatfile        character*(MAX_LEN_FNAM) sstdatfile
1021          character*(MAX_LEN_FNAM) tmidatfile
1022        character*(MAX_LEN_FNAM) sssdatfile        character*(MAX_LEN_FNAM) sssdatfile
1023        character*(MAX_LEN_FNAM) topexmeanfile        character*(MAX_LEN_FNAM) topexmeanfile
1024        character*(MAX_LEN_FNAM) topexfile        character*(MAX_LEN_FNAM) topexfile
# Line 832  c     driftfile     - reference data fil Line 1036  c     driftfile     - reference data fil
1036        character*(MAX_LEN_FNAM) curmtrufile        character*(MAX_LEN_FNAM) curmtrufile
1037        character*(MAX_LEN_FNAM) curmtrvfile        character*(MAX_LEN_FNAM) curmtrvfile
1038    
1039    
1040  c     Flags used in the model-data comparison:  c     Flags used in the model-data comparison:
1041  c     ========================================  c     ========================================
1042  c  c
# Line 846  c     using_ers - flag that indicates th Line 1051  c     using_ers - flag that indicates th
1051  c     Calendar information for the observations:  c     Calendar information for the observations:
1052  c     ==========================================  c     ==========================================
1053  c  c
 c     tmistartdate   - start date of the sea surface temperature data.  
1054  c     sststartdate   - start date of the sea surface temperature data.  c     sststartdate   - start date of the sea surface temperature data.
1055    c     tmistartdate   - start date of the sea surface temperature data.
1056  c     topexstartdate - start date of the sea surface height data.  c     topexstartdate - start date of the sea surface height data.
1057  c     ersstartdate   - start date of the sea surface height data.  c     ersstartdate   - start date of the sea surface height data.
1058  c     sshperiod      - sampling interval for  sea surface height data.  c     sshperiod      - sampling interval for the sea surface height data.
1059    
1060        common /ecco_cost_data_times_i/        common /ecco_cost_data_times_i/
1061       &                           scatxstartdate,       &                           scatxstartdate,
1062       &                           scatystartdate,       &                           scatystartdate,
      &                           tmistartdate,  
1063       &                           sststartdate,       &                           sststartdate,
      &                           sssstartdate,  
1064       &                           argotstartdate,       &                           argotstartdate,
1065       &                           argosstartdate,       &                           argosstartdate,
1066         &                           tmistartdate,
1067         &                           sssstartdate,
1068       &                           topexstartdate,       &                           topexstartdate,
1069       &                           ersstartdate       &                           ersstartdate
1070        integer scatxstartdate(4)        integer scatxstartdate(4)
1071        integer scatystartdate(4)        integer scatystartdate(4)
       integer tmistartdate(4)  
1072        integer sststartdate(4)        integer sststartdate(4)
       integer sssstartdate(4)  
1073        integer argotstartdate(4)        integer argotstartdate(4)
1074        integer argosstartdate(4)        integer argosstartdate(4)
1075          integer tmistartdate(4)
1076          integer sssstartdate(4)
1077        integer topexstartdate(4)        integer topexstartdate(4)
1078        integer ersstartdate(4)        integer ersstartdate(4)
1079    
# Line 888  c     sshperiod      - sampling interval Line 1093  c     sshperiod      - sampling interval
1093       &                           ersstartdate1,       &                           ersstartdate1,
1094       &                           ersstartdate2,       &                           ersstartdate2,
1095       &                           scatstartdate1,       &                           scatstartdate1,
1096       &                           scatstartdate2,       &                           scatstartdate2
      &                           costIceStart1,  
      &                           costIceStart2,  
      &                           costIceEnd1,  
      &                           costIceEnd2  
1097    
1098        integer tmistartdate1        integer tmistartdate1
1099        integer tmistartdate2        integer tmistartdate2
# Line 910  c     sshperiod      - sampling interval Line 1111  c     sshperiod      - sampling interval
1111        integer ersstartdate2        integer ersstartdate2
1112        integer scatstartdate1        integer scatstartdate1
1113        integer scatstartdate2        integer scatstartdate2
       integer costIceStart1  
       integer costIceStart2  
       integer costIceEnd1  
       integer costIceEnd2  
   
 c     costIceStart   - cost_ice start in seconds wrt model startTime  
 c     costIceEnd     - cost_ice end in seconds wrt model startTime  
1114    
1115        common /ecco_cost_data_times_r/        common /ecco_cost_data_times_r/
1116       &                           topexperiod,       &                           topexperiod,
1117       &                           ersperiod,       &                           ersperiod,
1118       &                           scatperiod,       &                           scatperiod
      &                           costIceStart,  
      &                           costIceEnd  
1119        _RL topexperiod        _RL topexperiod
1120        _RL ersperiod        _RL ersperiod
1121        _RL scatperiod        _RL scatperiod
       _RL costIceStart  
       _RL costIceEnd  
   
 c     cost_ice_flag  - cost_ice flag (see cost_ice.F)  
   
       common /ecco_cost_ice_i/ cost_ice_flag  
       integer cost_ice_flag  
1122    
1123  c     ==================================================================  c     ==================================================================
1124  c     END OF HEADER COST  c     END OF HEADER COST
1125  c     ==================================================================  c     ==================================================================
1126    
1127    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22