/[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.2 by heimbach, Mon Oct 11 16:38:53 2004 UTC revision 1.21 by heimbach, Fri Jul 7 22:53:36 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 23  c     dayrec  - number of averaged surfa Line 21  c     dayrec  - number of averaged surfa
21  c     monrec  - number of averaged theta and salinity records.  c     monrec  - number of averaged theta and salinity records.
22    
23        common /average_i/        common /average_i/
24       &                   sum1day,sum1mon,       &                   sum1day,sum1mon,sum1year,
25       &                   dayrec,monrec       &                   dayrec,monrec,yearrec
26        integer sum1day        integer sum1day
27        integer sum1mon        integer sum1mon
28          integer sum1year
29        integer dayrec        integer dayrec
30        integer monrec        integer monrec
31          integer yearrec
32    
33    
34  c     Averaged Fields:  c     Averaged Fields:
# Line 256  c     objf_hfluxmm    - Global time-mean Line 256  c     objf_hfluxmm    - Global time-mean
256  c     objf_sfluxmm    - Global time-mean Salt flux contribution.  c     objf_sfluxmm    - Global time-mean Salt flux contribution.
257  c     objf_hmean - Mean sea surface height contribution.  c     objf_hmean - Mean sea surface height contribution.
258  c     objf_h     - Residual sea surface height contribution.  c     objf_h     - Residual sea surface height contribution.
259    c     objf_tp    - Residual sea surface height contribution from T/P
260    c     objf_ers   - Residual sea surface height contribution from T/P
261    c     objf_gfo   - Residual sea surface height contribution from T/P
262  c     objf_temp  - Temperature contribution.  c     objf_temp  - Temperature contribution.
263  c     objf_salt  - Salinity contribution.  c     objf_salt  - Salinity contribution.
264  c     objf_temp0 - Initial conditions Temperature contribution.  c     objf_temp0 - Initial conditions Temperature contribution.
# Line 274  c     objf_scatxm - time-mean zonal SCAT Line 277  c     objf_scatxm - time-mean zonal SCAT
277  c     objf_scatym - time-mean meridional SCAT  contribution  c     objf_scatym - time-mean meridional SCAT  contribution
278  c     objf_scatx  - zonal SCAT  contribution  c     objf_scatx  - zonal SCAT  contribution
279  c     objf_scaty  - meridional SCAT  contribution  c     objf_scaty  - meridional SCAT  contribution
 c     objf_ice    - sea-ice volume  
280  c     objf_kapgm  - kappa GM contribution  c     objf_kapgm  - kappa GM contribution
281  c     objf_diffkr - diffusion contribution  c     objf_diffkr - diffusion contribution
282  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 289  c                  function contribution
289       &                objf_hflux,       &                objf_hflux,
290       &                objf_hfluxm,       &                objf_hfluxm,
291       &                objf_hfluxmm,       &                objf_hfluxmm,
292         &                objf_hfluxsmoo,
293       &                objf_sflux,       &                objf_sflux,
294       &                objf_sfluxm,       &                objf_sfluxm,
295       &                objf_sfluxmm,       &                objf_sfluxmm,
296         &                objf_sfluxsmoo,
297       &                objf_tauu,       &                objf_tauu,
298       &                objf_tauum,       &                objf_tauum,
299         &                objf_tauusmoo,
300       &                objf_tauv,       &                objf_tauv,
301       &                objf_tauvm,       &                objf_tauvm,
302         &                objf_tauvsmoo,
303       &                objf_hmean,       &                objf_hmean,
304       &                objf_h,       &                objf_h,
305         &                objf_tp,
306         &                objf_ers,
307         &                objf_gfo,
308       &                objf_temp,       &                objf_temp,
309       &                objf_salt,       &                objf_salt,
310       &                objf_temp0,       &                objf_temp0,
311       &                objf_salt0,       &                objf_salt0,
312         &                objf_temp0smoo,
313         &                objf_salt0smoo,
314       &                objf_sst,       &                objf_sst,
315       &                objf_tmi,       &                objf_tmi,
316       &                objf_sss,       &                objf_sss,
# Line 320  c                  function contribution Line 331  c                  function contribution
331       &                objf_scatym,       &                objf_scatym,
332       &                objf_atemp,       &                objf_atemp,
333       &                objf_aqh,       &                objf_aqh,
334         &                objf_precip,
335         &                objf_swflux,
336         &                objf_swdown,
337       &                objf_uwind,       &                objf_uwind,
338       &                objf_vwind,       &                objf_vwind,
339         &                objf_atempm,
340         &                objf_aqhm,
341         &                objf_precipm,
342         &                objf_swfluxm,
343         &                objf_swdownm,
344         &                objf_uwindm,
345         &                objf_vwindm,
346         &                objf_atempsmoo,
347         &                objf_aqhsmoo,
348         &                objf_precipsmoo,
349         &                objf_swfluxsmoo,
350         &                objf_swdownsmoo,
351         &                objf_uwindsmoo,
352         &                objf_vwindsmoo,
353       &                objf_obcsn,       &                objf_obcsn,
354       &                objf_obcss,       &                objf_obcss,
355       &                objf_obcsw,       &                objf_obcsw,
# Line 329  c                  function contribution Line 357  c                  function contribution
357       &                objf_obcsvol,       &                objf_obcsvol,
358       &                objf_curmtr,       &                objf_curmtr,
359       &                objf_ageos,       &                objf_ageos,
360       &                objf_ice,       &                objf_kapgm,
361       &                objf_diffkr,       &                objf_diffkr,
362       &                objf_theta_ini_fin,       &                objf_theta_ini_fin,
363       &                objf_salt_ini_fin       &                objf_salt_ini_fin,
364         &                objf_eddytau
365        _RL  objf_hflux  (nsx,nsy)        _RL  objf_hflux  (nsx,nsy)
366        _RL  objf_hfluxm (nsx,nsy)        _RL  objf_hfluxm (nsx,nsy)
367        _RL  objf_hfluxmm(nsx,nsy)        _RL  objf_hfluxmm(nsx,nsy)
368          _RL  objf_hfluxsmoo  (nsx,nsy)
369        _RL  objf_sflux  (nsx,nsy)        _RL  objf_sflux  (nsx,nsy)
370        _RL  objf_sfluxm (nsx,nsy)        _RL  objf_sfluxm (nsx,nsy)
371        _RL  objf_sfluxmm(nsx,nsy)        _RL  objf_sfluxmm(nsx,nsy)
372          _RL  objf_sfluxsmoo  (nsx,nsy)
373        _RL  objf_tauu   (nsx,nsy)        _RL  objf_tauu   (nsx,nsy)
374        _RL  objf_tauum  (nsx,nsy)        _RL  objf_tauum  (nsx,nsy)
375          _RL  objf_tauusmoo   (nsx,nsy)
376        _RL  objf_tauv   (nsx,nsy)        _RL  objf_tauv   (nsx,nsy)
377        _RL  objf_tauvm  (nsx,nsy)        _RL  objf_tauvm  (nsx,nsy)
378          _RL  objf_tauvsmoo   (nsx,nsy)
379        _RL  objf_hmean        _RL  objf_hmean
380        _RL  objf_h    (nsx,nsy)        _RL  objf_h    (nsx,nsy)
381          _RL  objf_tp   (nsx,nsy)
382          _RL  objf_ers  (nsx,nsy)
383          _RL  objf_gfo  (nsx,nsy)
384        _RL  objf_temp (nsx,nsy)        _RL  objf_temp (nsx,nsy)
385        _RL  objf_salt (nsx,nsy)        _RL  objf_salt (nsx,nsy)
386        _RL  objf_temp0(nsx,nsy)        _RL  objf_temp0(nsx,nsy)
387        _RL  objf_salt0(nsx,nsy)        _RL  objf_salt0(nsx,nsy)
388          _RL  objf_temp0smoo(nsx,nsy)
389          _RL  objf_salt0smoo(nsx,nsy)
390        _RL  objf_sst  (nsx,nsy)        _RL  objf_sst  (nsx,nsy)
391        _RL  objf_tmi  (nsx,nsy)        _RL  objf_tmi  (nsx,nsy)
392        _RL  objf_sss  (nsx,nsy)        _RL  objf_sss  (nsx,nsy)
# Line 370  c                  function contribution Line 407  c                  function contribution
407        _RL  objf_scatym(nsx,nsy)        _RL  objf_scatym(nsx,nsy)
408        _RL  objf_atemp(nsx,nsy)        _RL  objf_atemp(nsx,nsy)
409        _RL  objf_aqh  (nsx,nsy)        _RL  objf_aqh  (nsx,nsy)
410          _RL  objf_precip(nsx,nsy)
411          _RL  objf_swflux(nsx,nsy)
412          _RL  objf_swdown(nsx,nsy)
413        _RL  objf_uwind(nsx,nsy)        _RL  objf_uwind(nsx,nsy)
414        _RL  objf_vwind(nsx,nsy)        _RL  objf_vwind(nsx,nsy)
415          _RL  objf_atempm(nsx,nsy)
416          _RL  objf_aqhm  (nsx,nsy)
417          _RL  objf_precipm(nsx,nsy)
418          _RL  objf_swfluxm(nsx,nsy)
419          _RL  objf_swdownm(nsx,nsy)
420          _RL  objf_uwindm(nsx,nsy)
421          _RL  objf_vwindm(nsx,nsy)
422          _RL  objf_atempsmoo(nsx,nsy)
423          _RL  objf_aqhsmoo  (nsx,nsy)
424          _RL  objf_precipsmoo(nsx,nsy)
425          _RL  objf_swfluxsmoo(nsx,nsy)
426          _RL  objf_swdownsmoo(nsx,nsy)
427          _RL  objf_uwindsmoo(nsx,nsy)
428          _RL  objf_vwindsmoo(nsx,nsy)
429        _RL  objf_obcsn(nsx,nsy)        _RL  objf_obcsn(nsx,nsy)
430        _RL  objf_obcss(nsx,nsy)        _RL  objf_obcss(nsx,nsy)
431        _RL  objf_obcsw(nsx,nsy)        _RL  objf_obcsw(nsx,nsy)
# Line 379  c                  function contribution Line 433  c                  function contribution
433        _RL  objf_obcsvol        _RL  objf_obcsvol
434        _RL  objf_curmtr(nsx,nsy)        _RL  objf_curmtr(nsx,nsy)
435        _RL  objf_ageos(nsx,nsy)        _RL  objf_ageos(nsx,nsy)
       _RL  objf_ice  (nsx,nsy)  
436        _RL  objf_kapgm(nsx,nsy)        _RL  objf_kapgm(nsx,nsy)
437        _RL  objf_diffkr(nsx,nsy)        _RL  objf_diffkr(nsx,nsy)
438        _RL  objf_theta_ini_fin(nsx,nsy)        _RL  objf_theta_ini_fin(nsx,nsy)
439        _RL  objf_salt_ini_fin(nsx,nsy)        _RL  objf_salt_ini_fin(nsx,nsy)
440          _RL  objf_eddytau(nsx,nsy)
441    
442          common /ecco_cost_num/
443         &                num_hflux,
444         &                num_hfluxm,
445         &                num_hfluxmm,
446         &                num_sflux,
447         &                num_sfluxm,
448         &                num_sfluxmm,
449         &                num_tauu,
450         &                num_tauum,
451         &                num_tauv,
452         &                num_tauvm,
453         &                num_hmean,
454         &                num_h,
455         &                num_tp,
456         &                num_ers,
457         &                num_gfo,
458         &                num_temp,
459         &                num_salt,
460         &                num_temp0,
461         &                num_salt0,
462         &                num_sst,
463         &                num_tmi,
464         &                num_sss,
465         &                num_ctdt,
466         &                num_ctds,
467         &                num_ctdtclim,
468         &                num_ctdsclim,
469         &                num_xbt,
470         &                num_argot,
471         &                num_argos,
472         &                num_drift,
473         &                num_tdrift,
474         &                num_sdrift,
475         &                num_wdrift,
476         &                num_scatx,
477         &                num_scaty,
478         &                num_scatxm,
479         &                num_scatym,
480         &                num_atemp,
481         &                num_aqh,
482         &                num_precip,
483         &                num_swflux,
484         &                num_swdown,
485         &                num_uwind,
486         &                num_vwind,
487         &                num_atempm,
488         &                num_aqhm,
489         &                num_precipm,
490         &                num_swfluxm,
491         &                num_swdownm,
492         &                num_uwindm,
493         &                num_vwindm,
494         &                num_obcsn,
495         &                num_obcss,
496         &                num_obcsw,
497         &                num_obcse,
498         &                num_obcsvol,
499         &                num_curmtr,
500         &                num_ageos,
501         &                num_kapgm,
502         &                num_diffkr,
503         &                num_theta_ini_fin,
504         &                num_salt_ini_fin,
505         &                num_eddytau
506    
507          _RL  num_hflux  (nsx,nsy)
508          _RL  num_hfluxm (nsx,nsy)
509          _RL  num_hfluxmm(nsx,nsy)
510          _RL  num_sflux  (nsx,nsy)
511          _RL  num_sfluxm (nsx,nsy)
512          _RL  num_sfluxmm(nsx,nsy)
513          _RL  num_tauu   (nsx,nsy)
514          _RL  num_tauum  (nsx,nsy)
515          _RL  num_tauv   (nsx,nsy)
516          _RL  num_tauvm  (nsx,nsy)
517          _RL  num_hmean
518          _RL  num_h    (nsx,nsy)
519          _RL  num_tp   (nsx,nsy)
520          _RL  num_ers  (nsx,nsy)
521          _RL  num_gfo  (nsx,nsy)
522          _RL  num_temp (nsx,nsy)
523          _RL  num_salt (nsx,nsy)
524          _RL  num_temp0(nsx,nsy)
525          _RL  num_salt0(nsx,nsy)
526          _RL  num_sst  (nsx,nsy)
527          _RL  num_tmi  (nsx,nsy)
528          _RL  num_sss  (nsx,nsy)
529          _RL  num_ctdt (nsx,nsy)
530          _RL  num_ctds (nsx,nsy)
531          _RL  num_ctdtclim (nsx,nsy)
532          _RL  num_ctdsclim (nsx,nsy)
533          _RL  num_xbt  (nsx,nsy)
534          _RL  num_argot(nsx,nsy)
535          _RL  num_argos(nsx,nsy)
536          _RL  num_drift(nsx,nsy)
537          _RL  num_tdrift(nsx,nsy)
538          _RL  num_sdrift(nsx,nsy)
539          _RL  num_wdrift(nsx,nsy)
540          _RL  num_scatx(nsx,nsy)
541          _RL  num_scaty(nsx,nsy)
542          _RL  num_scatxm(nsx,nsy)
543          _RL  num_scatym(nsx,nsy)
544          _RL  num_atemp(nsx,nsy)
545          _RL  num_aqh  (nsx,nsy)
546          _RL  num_precip(nsx,nsy)
547          _RL  num_swflux(nsx,nsy)
548          _RL  num_swdown(nsx,nsy)
549          _RL  num_uwind(nsx,nsy)
550          _RL  num_vwind(nsx,nsy)
551          _RL  num_atempm(nsx,nsy)
552          _RL  num_aqhm  (nsx,nsy)
553          _RL  num_precipm(nsx,nsy)
554          _RL  num_swfluxm(nsx,nsy)
555          _RL  num_swdownm(nsx,nsy)
556          _RL  num_uwindm(nsx,nsy)
557          _RL  num_vwindm(nsx,nsy)
558          _RL  num_obcsn(nsx,nsy)
559          _RL  num_obcss(nsx,nsy)
560          _RL  num_obcsw(nsx,nsy)
561          _RL  num_obcse(nsx,nsy)
562          _RL  num_obcsvol
563          _RL  num_curmtr(nsx,nsy)
564          _RL  num_ageos(nsx,nsy)
565          _RL  num_kapgm(nsx,nsy)
566          _RL  num_diffkr(nsx,nsy)
567          _RL  num_theta_ini_fin(nsx,nsy)
568          _RL  num_salt_ini_fin(nsx,nsy)
569          _RL  num_eddytau(nsx,nsy)
570    
571        common /ecco_cost_aux_r/        common /ecco_cost_aux_r/
572       &                    mult_hflux,       &                    mult_hflux,
573       &                    mult_sflux,       &                    mult_sflux,
574         &                    mult_hfluxmm,
575         &                    mult_sfluxmm,
576       &                    mult_tauu,       &                    mult_tauu,
577       &                    mult_tauv,       &                    mult_tauv,
578       &                    mult_hmean,       &                    mult_hmean,
579       &                    mult_h,       &                    mult_h,
580         &                    mult_tp,
581         &                    mult_ers,
582         &                    mult_gfo,
583       &                    mult_temp,       &                    mult_temp,
584       &                    mult_salt,       &                    mult_salt,
585       &                    mult_temp0,       &                    mult_temp0,
# Line 414  c                  function contribution Line 602  c                  function contribution
602       &                    mult_scaty,       &                    mult_scaty,
603       &                    mult_atemp,       &                    mult_atemp,
604       &                    mult_aqh,       &                    mult_aqh,
605         &                    mult_precip,
606         &                    mult_swflux,
607         &                    mult_swdown,
608       &                    mult_uwind,       &                    mult_uwind,
609       &                    mult_vwind,       &                    mult_vwind,
610       &                    mult_obcsn,       &                    mult_obcsn,
# Line 423  c                  function contribution Line 614  c                  function contribution
614       &                    mult_obcsvol,       &                    mult_obcsvol,
615       &                    mult_curmtr,       &                    mult_curmtr,
616       &                    mult_ageos,       &                    mult_ageos,
      &                    mult_ice,  
617       &                    mult_kapgm,       &                    mult_kapgm,
618       &                    mult_diffkr,       &                    mult_diffkr,
619       &                    mult_ini_fin       &                    mult_ini_fin,
620         &                    mult_eddytau,
621         &                    mult_smooth_ic,
622         &                    mult_smooth_bc
623        _RL  mult_hflux        _RL  mult_hflux
624        _RL  mult_sflux        _RL  mult_sflux
625          _RL  mult_hfluxmm
626          _RL  mult_sfluxmm
627        _RL  mult_tauu        _RL  mult_tauu
628        _RL  mult_tauv        _RL  mult_tauv
629        _RL  mult_hmean        _RL  mult_hmean
630        _RL  mult_h        _RL  mult_h
631          _RL  mult_tp
632          _RL  mult_ers
633          _RL  mult_gfo
634        _RL  mult_temp        _RL  mult_temp
635        _RL  mult_salt        _RL  mult_salt
636        _RL  mult_temp0        _RL  mult_temp0
# Line 456  c                  function contribution Line 653  c                  function contribution
653        _RL  mult_scaty        _RL  mult_scaty
654        _RL  mult_atemp        _RL  mult_atemp
655        _RL  mult_aqh        _RL  mult_aqh
656          _RL  mult_precip
657          _RL  mult_swflux
658          _RL  mult_swdown
659        _RL  mult_uwind        _RL  mult_uwind
660        _RL  mult_vwind        _RL  mult_vwind
661        _RL  mult_obcsn        _RL  mult_obcsn
# Line 465  c                  function contribution Line 665  c                  function contribution
665        _RL  mult_obcsvol        _RL  mult_obcsvol
666        _RL  mult_curmtr        _RL  mult_curmtr
667        _RL  mult_ageos        _RL  mult_ageos
       _RL  mult_ice  
668        _RL  mult_kapgm        _RL  mult_kapgm
669        _RL  mult_diffkr        _RL  mult_diffkr
670        _RL  mult_ini_fin        _RL  mult_ini_fin
671          _RL  mult_eddytau
672          _RL  mult_smooth_ic
673          _RL  mult_smooth_bc
674    
675  c     Record counters relevant for the cost function evaluation.  c     Record counters relevant for the cost function evaluation.
676  c     ==========================================================  c     ==========================================================
677  c  c
678  c     nmonsrec - number of monthly records that will be generated by  c     nyearsrec - number of yearly records that will be generated by
679  c                the current model integration.  c                 the current model integration.
680  c     ndaysrec - number of  daily  records that will be generated by  c     nmonsrec  - number of monthly records that will be generated by
681  c                the current model integration.  c                 the current model integration.
682    c     ndaysrec  - number of  daily  records that will be generated by
683    c                 the current model integration.
684    
685        common /ecco_cost_i/        common /ecco_cost_i/
686         &                nyearsrec,
687       &                nmonsrec,       &                nmonsrec,
688       &                ndaysrec       &                ndaysrec,
689         &                nnztbar,
690         &                nnzsbar
691          integer nyearsrec
692        integer nmonsrec        integer nmonsrec
693        integer ndaysrec        integer ndaysrec
694          integer nnztbar
695          integer nnzsbar
696    
697    
698  c     Data files for the weights used in the cost function:  c     Data files for the weights used in the cost function:
# Line 531  c     velerrfile            - representa Line 740  c     velerrfile            - representa
740       &                salterrfile,       &                salterrfile,
741       &                temperrfile,       &                temperrfile,
742       &                velerrfile,       &                velerrfile,
743         &                salt0errfile,
744         &                temp0errfile,
745         &                vel0errfile,
746       &                atemp_errfile,       &                atemp_errfile,
747       &                aqh_errfile,       &                aqh_errfile,
748         &                precip_errfile,
749         &                swflux_errfile,
750         &                swdown_errfile,
751       &                uwind_errfile,       &                uwind_errfile,
752       &                vwind_errfile       &                vwind_errfile
753        character*(MAX_LEN_FNAM) hflux_errfile        character*(MAX_LEN_FNAM) hflux_errfile
# Line 557  c     velerrfile            - representa Line 772  c     velerrfile            - representa
772        character*(MAX_LEN_FNAM) salterrfile        character*(MAX_LEN_FNAM) salterrfile
773        character*(MAX_LEN_FNAM) temperrfile        character*(MAX_LEN_FNAM) temperrfile
774        character*(MAX_LEN_FNAM) velerrfile        character*(MAX_LEN_FNAM) velerrfile
775          character*(MAX_LEN_FNAM) salt0errfile
776          character*(MAX_LEN_FNAM) temp0errfile
777          character*(MAX_LEN_FNAM) vel0errfile
778        character*(MAX_LEN_FNAM) atemp_errfile        character*(MAX_LEN_FNAM) atemp_errfile
779        character*(MAX_LEN_FNAM) aqh_errfile        character*(MAX_LEN_FNAM) aqh_errfile
780          character*(MAX_LEN_FNAM) precip_errfile
781          character*(MAX_LEN_FNAM) swflux_errfile
782          character*(MAX_LEN_FNAM) swdown_errfile
783        character*(MAX_LEN_FNAM) uwind_errfile        character*(MAX_LEN_FNAM) uwind_errfile
784        character*(MAX_LEN_FNAM) vwind_errfile        character*(MAX_LEN_FNAM) vwind_errfile
785    
# Line 591  c     wvdrift    - weight for mean merid Line 812  c     wvdrift    - weight for mean merid
812       &                      frame,       &                      frame,
813       &                      cosphi,       &                      cosphi,
814       &                      whflux,wsflux,wtauu,wtauv,       &                      whflux,wsflux,wtauu,wtauv,
815       &                      watemp,waqh,wuwind,wvwind,       &                      watemp,waqh,wprecip,wswflux,wswdown,
816         &                      wuwind,wvwind,
817       &                      wscatx,wscaty,       &                      wscatx,wscaty,
818       &                      wtheta,wtheta2,wthetaLev,       &                      wtheta,wtheta2,wthetaLev,
819       &                      wsalt,wsalt2,wsaltLev,       &                      wsalt,wsalt2,wsaltLev,
820         &                      wdiffkr,wdiffkr2,wdiffkrFld,
821         &                      wkapgm,wkapgm2,wkapgmFld,
822         &                      wedtaux,wedtaux2,wedtauxFld,
823         &                      wedtauy,wedtauy2,wedtauyFld,
824       &                      wsst,wsss,       &                      wsst,wsss,
825       &                      wtp,wers,       &                      wtp,wers,wgfo,
826       &                      wp,       &                      wp,
827       &                      wctdt,wctds,       &                      wctdt,wctds,
828       &                      wudrift,wvdrift,       &                      wudrift,wvdrift,
# Line 620  c     wvdrift    - weight for mean merid Line 846  c     wvdrift    - weight for mean merid
846        _RL wscaty  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wscaty  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
847        _RL watemp  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL watemp  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
848        _RL waqh    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL waqh    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
849          _RL wprecip (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
850          _RL wswflux (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
851          _RL wswdown (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
852        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
853        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
854        _RL wtheta  (                            nr,nsx,nsy)        _RL wtheta  (                            nr,nsx,nsy)
# Line 632  c     wvdrift    - weight for mean merid Line 861  c     wvdrift    - weight for mean merid
861        _RL wsss    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsss    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
862        _RL wtp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
863        _RL wers    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wers    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
864          _RL wgfo    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
865        _RL wp      (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wp      (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
866        _RL wctdt   (                            nr,nsx,nsy)        _RL wctdt   (                            nr,nsx,nsy)
867        _RL wctds   (                            nr,nsx,nsy)        _RL wctds   (                            nr,nsx,nsy)
# Line 641  c     wvdrift    - weight for mean merid Line 871  c     wvdrift    - weight for mean merid
871        _RL wcurrent2   (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL wcurrent2   (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
872        _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)
873        _RL wbaro        _RL wbaro
874          _RL wdiffkr (                            nr,nsx,nsy)
875          _RL wdiffkr2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
876          _RL wdiffkrFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
877          _RL wkapgm (                            nr,nsx,nsy)
878          _RL wkapgm2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
879          _RL wkapgmFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
880          _RL wedtaux (                            nr,nsx,nsy)
881          _RL wedtaux2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
882          _RL wedtauxFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
883          _RL wedtauy (                            nr,nsx,nsy)
884          _RL wedtauy2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
885          _RL wedtauyFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
886    
887          common /ecco_cost_weights_0_r/
888         &        whflux0, wsflux0, wtau0,
889         &        watemp0, waqh0, wprecip0, wswflux0, wswdown0, wwind0,
890         &        wdiffkr0, wkapgm0, wedtau0
891          _RL whflux0
892          _RL wsflux0
893          _RL wtau0
894          _RL watemp0
895          _RL waqh0
896          _RL wprecip0
897          _RL wswflux0
898          _RL wswdown0
899          _RL wwind0
900          _RL wdiffkr0
901          _RL wkapgm0
902          _RL wedtau0
903    
904          common /ecco_cost_weights_mean_r/
905         &        wmean_hflux, wmean_sflux, wmean_tau,
906         &        wmean_atemp, wmean_aqh,
907         &        wmean_precip, wmean_swflux, wmean_swdown, wmean_wind
908          _RL wmean_hflux
909          _RL wmean_sflux
910          _RL wmean_tau
911          _RL wmean_atemp
912          _RL wmean_aqh
913          _RL wmean_precip
914          _RL wmean_swflux
915          _RL wmean_swdown
916          _RL wmean_wind
917                                                    
918          common /ecco_cost_weights_2_r/
919         &                      whflux2,wsflux2,wtauu2,wtauv2
920          _RL whflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
921          _RL wsflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
922          _RL wtauu2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
923          _RL wtauv2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
924                                                    
925                          
926    
927  #if (defined (ALLOW_OBCSN_COST_CONTRIBUTION) || \  #if (defined (ALLOW_OBCSN_COST_CONTRIBUTION) || \
928       defined (ALLOW_OBCSN_CONTROL))       defined (ALLOW_OBCSN_CONTROL))
# Line 724  c     vdriftdat  - drifters meridional v Line 1006  c     vdriftdat  - drifters meridional v
1006       &                     tpmask,       &                     tpmask,
1007       &                     ersobs,       &                     ersobs,
1008       &                     ersmask,       &                     ersmask,
1009         &                     gfoobs,
1010         &                     gfomask,
1011       &                     ctdtobs,       &                     ctdtobs,
1012       &                     ctdsobs,       &                     ctdsobs,
1013       &                     xbtobs,       &                     xbtobs,
# Line 754  c     vdriftdat  - drifters meridional v Line 1038  c     vdriftdat  - drifters meridional v
1038        _RL tpmask    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpmask    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1039        _RL ersobs    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL ersobs    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1040        _RL ersmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL ersmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1041          _RL gfoobs    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1042          _RL gfomask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1043        _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)
1044        _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)
1045        _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 799  c     driftfile     - reference data fil Line 1085  c     driftfile     - reference data fil
1085       &                     topexmeanfile,       &                     topexmeanfile,
1086       &                     topexfile,       &                     topexfile,
1087       &                     ersfile,       &                     ersfile,
1088         &                     gfofile,
1089       &                     ctdtfile,       &                     ctdtfile,
1090       &                     ctdsfile,       &                     ctdsfile,
1091       &                     ctdtclimfile,       &                     ctdtclimfile,
# Line 821  c     driftfile     - reference data fil Line 1108  c     driftfile     - reference data fil
1108        character*(MAX_LEN_FNAM) topexmeanfile        character*(MAX_LEN_FNAM) topexmeanfile
1109        character*(MAX_LEN_FNAM) topexfile        character*(MAX_LEN_FNAM) topexfile
1110        character*(MAX_LEN_FNAM) ersfile        character*(MAX_LEN_FNAM) ersfile
1111          character*(MAX_LEN_FNAM) gfofile
1112        character*(MAX_LEN_FNAM) ctdtfile        character*(MAX_LEN_FNAM) ctdtfile
1113        character*(MAX_LEN_FNAM) ctdsfile        character*(MAX_LEN_FNAM) ctdsfile
1114        character*(MAX_LEN_FNAM) ctdtclimfile        character*(MAX_LEN_FNAM) ctdtclimfile
# Line 842  c     using_ers - flag that indicates th Line 1130  c     using_ers - flag that indicates th
1130    
1131        common /ecco_cost_data_flags/        common /ecco_cost_data_flags/
1132       &                         using_topex,       &                         using_topex,
1133       &                         using_ers       &                         using_ers,
1134         &                         using_gfo
1135        logical using_topex        logical using_topex
1136        logical using_ers        logical using_ers
1137          logical using_gfo
1138    
1139  c     Calendar information for the observations:  c     Calendar information for the observations:
1140  c     ==========================================  c     ==========================================
# Line 864  c     sshperiod      - sampling interval Line 1154  c     sshperiod      - sampling interval
1154       &                           tmistartdate,       &                           tmistartdate,
1155       &                           sssstartdate,       &                           sssstartdate,
1156       &                           topexstartdate,       &                           topexstartdate,
1157       &                           ersstartdate       &                           ersstartdate,
1158         &                           gfostartdate
1159        integer scatxstartdate(4)        integer scatxstartdate(4)
1160        integer scatystartdate(4)        integer scatystartdate(4)
1161        integer sststartdate(4)        integer sststartdate(4)
# Line 874  c     sshperiod      - sampling interval Line 1165  c     sshperiod      - sampling interval
1165        integer sssstartdate(4)        integer sssstartdate(4)
1166        integer topexstartdate(4)        integer topexstartdate(4)
1167        integer ersstartdate(4)        integer ersstartdate(4)
1168          integer gfostartdate(4)
1169    
1170        common /ecco_cost_data_aux_i/        common /ecco_cost_data_aux_i/
1171       &                           tmistartdate1,       &                           tmistartdate1,
# Line 890  c     sshperiod      - sampling interval Line 1182  c     sshperiod      - sampling interval
1182       &                           topexstartdate2,       &                           topexstartdate2,
1183       &                           ersstartdate1,       &                           ersstartdate1,
1184       &                           ersstartdate2,       &                           ersstartdate2,
1185         &                           gfostartdate1,
1186         &                           gfostartdate2,
1187       &                           scatstartdate1,       &                           scatstartdate1,
1188       &                           scatstartdate2,       &                           scatstartdate2
      &                           costIceStart1,  
      &                           costIceStart2,  
      &                           costIceEnd1,  
      &                           costIceEnd2  
1189    
1190        integer tmistartdate1        integer tmistartdate1
1191        integer tmistartdate2        integer tmistartdate2
# Line 911  c     sshperiod      - sampling interval Line 1201  c     sshperiod      - sampling interval
1201        integer topexstartdate2        integer topexstartdate2
1202        integer ersstartdate1        integer ersstartdate1
1203        integer ersstartdate2        integer ersstartdate2
1204          integer gfostartdate1
1205          integer gfostartdate2
1206        integer scatstartdate1        integer scatstartdate1
1207        integer scatstartdate2        integer scatstartdate2
       integer costIceStart1  
       integer costIceStart2  
       integer costIceEnd1  
       integer costIceEnd2  
1208    
1209        common /ecco_cost_data_times_r/        common /ecco_cost_data_times_r/
1210       &                           topexperiod,       &                           topexperiod,
1211       &                           ersperiod,       &                           ersperiod,
1212       &                           scatperiod,       &                           gfoperiod,
1213       &                           costIceStart,       &                           scatperiod
      &                           costIceEnd  
1214        _RL topexperiod        _RL topexperiod
1215        _RL ersperiod        _RL ersperiod
1216          _RL gfoperiod
1217        _RL scatperiod        _RL scatperiod
       _RL costIceStart  
       _RL costIceEnd  
   
 c     cost_ice_flag  - cost_ice flag (see cost_ice.F)  
1218    
1219        common /ecco_cost_ice_i/ cost_ice_flag        common /ecco_cost_data_detrend/
1220        integer cost_ice_flag       &                           topexintercept,
1221         &                           ersintercept,
1222         &                           gfointercept,
1223         &                           topexslope,
1224         &                           ersslope,
1225         &                           gfoslope
1226          _RL topexintercept
1227          _RL ersintercept
1228          _RL gfointercept
1229          _RL topexslope
1230          _RL ersslope
1231          _RL gfoslope
1232    
1233  c     ==================================================================  c     ==================================================================
1234  c     END OF HEADER COST  c     END OF HEADER COST
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22