/[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.17 by heimbach, Mon Mar 6 19:16:53 2006 UTC revision 1.48 by jmc, Mon Mar 22 02:19:35 2010 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  c     ==================================================================  c     ==================================================================
5  c     HEADER AVERAGES  c     HEADER AVERAGES
# Line 20  c     sum1mon - counter for the monthly Line 22  c     sum1mon - counter for the monthly
22  c     dayrec  - number of averaged surface pressure records.  c     dayrec  - number of averaged surface pressure records.
23  c     monrec  - number of averaged theta and salinity records.  c     monrec  - number of averaged theta and salinity records.
24    
25        common /average_i/        common /average_i/
26       &                   sum1day,sum1mon,       &                   sum1day,sum1mon,sum1year,
27       &                   dayrec,monrec       &                   dayrec,monrec,yearrec
28        integer sum1day        integer sum1day
29        integer sum1mon        integer sum1mon
30          integer sum1year
31        integer dayrec        integer dayrec
32        integer monrec        integer monrec
33          integer yearrec
34    
35    c     Number of sshv4cost Cost terms:
36    c     =============================
37          INTEGER NSSHV4COST
38          PARAMETER ( NSSHV4COST=5 )
39    
40    c     Number of User Cost terms:
41    c     =============================
42          INTEGER NUSERCOST
43          PARAMETER ( NUSERCOST=10 )
44    
45    c     Number of days: (hard-coded to set up some vector dimensions
46    c     =============================
47    c     17 years: 6210
48          INTEGER maxNumDays
49          PARAMETER ( maxNumDays = 6210 )
50    
51    c     Number of Generic Cost terms:
52    c     =============================
53          INTEGER NGENCOST
54          PARAMETER ( NGENCOST=1 )
55    
56    #ifdef ALLOW_GENCOST_CONTRIBUTION
57    c     objf_gencost - gencost user defined contribution
58          common /ecco_gencost_ctrl/
59         &       xx_genbar_dummy
60          _RL  xx_genbar_dummy(NGENCOST)
61    
62          common /ecco_gencost_r_1/
63         &       objf_gencost, num_gencost, mult_gencost,
64         &       gencost_barfld, gencost_modfld, gencost_weight,
65         &       gencost_spmin, gencost_spmax, gencost_spzero
66          _RL  objf_gencost(nsx,nsy,NGENCOST)
67          _RL  num_gencost(nsx,nsy,NGENCOST)
68          _RL  mult_gencost(NGENCOST)
69          _RL  gencost_spmin(NGENCOST)
70          _RL  gencost_spmax(NGENCOST)
71          _RL  gencost_spzero(NGENCOST)
72          _RL  gencost_barfld(1-olx:snx+olx,1-oly:sny+oly,
73         &       nsx,nsy,NGENCOST)
74          _RL  gencost_modfld(1-olx:snx+olx,1-oly:sny+oly,
75         &       nsx,nsy,NGENCOST)
76          _RL  gencost_weight(1-olx:snx+olx,1-oly:sny+oly,
77         &       nsx,nsy,NGENCOST)
78    
79          common /ecco_gencost_r_2/
80         &       gencost_period
81          _RL     gencost_period(NGENCOST)
82    
83          common /ecco_gencost_i_1/
84         &       gencost_nrec
85          integer gencost_nrec(NGENCOST)      
86    
87          common /ecco_gencost_c/
88         &       gencost_errfile,
89         &       gencost_datafile,
90         &       gencost_barfile,
91         &       gencost_avgperiod,
92         &       gencost_mask
93          character*(MAX_LEN_FNAM) gencost_errfile(NGENCOST)
94          character*(MAX_LEN_FNAM) gencost_datafile(NGENCOST)
95          character*(MAX_LEN_FNAM) gencost_barfile(NGENCOST)
96          character*(5)            gencost_avgperiod(NGENCOST)
97          character*(1)            gencost_mask(NGENCOST)
98    
99    #endif /* ALLOW_GENCOST_CONTRIBUTION */
100    
101  c     Averaged Fields:  c     Averaged Fields:
102  c     ================  c     ================
# Line 41  c             intantaneous salinities. Line 110  c             intantaneous salinities.
110  c     psbar - contains the averaged surface pressure field after the call  c     psbar - contains the averaged surface pressure field after the call
111  c             to subroutine POST_DAILY. Before, it accumulates the  c             to subroutine POST_DAILY. Before, it accumulates the
112  c             intantaneous surface pressure field.  c             intantaneous surface pressure field.
113  c     ubar  - contains the averaged zonal velocity component for the  c     ubar  - contains the averaged zonal velocity component for the
114  c             whole integration period. Before, it accumulates the  c             whole integration period. Before, it accumulates the
115  c             intantaneous field.  c             intantaneous field.
116  c     vbar  - contains the averaged zonal velocity component for the  c     vbar  - contains the averaged zonal velocity component for the
# Line 53  c             intantaneous field. Line 122  c             intantaneous field.
122  c     tauybar  - contains the averaged zonal velocity component for the  c     tauybar  - contains the averaged zonal velocity component for the
123  c             whole integration period. Before, it accumulates the  c             whole integration period. Before, it accumulates the
124  c             intantaneous field.  c             intantaneous field.
125  c     hfluxbar  - contains the averaged zonal velocity component for the  c     hfluxmeanbar  - contains the averaged zonal velocity component for the
126  c             whole integration period. Before, it accumulates the  c             whole integration period. Before, it accumulates the
127  c             intantaneous field.  c             intantaneous field.
128  c     sfluxbar  - contains the averaged zonal velocity component for the  c     sfluxmeanbar  - contains the averaged zonal velocity component for the
129  c             whole integration period. Before, it accumulates the  c             whole integration period. Before, it accumulates the
130  c             intantaneous field.  c             intantaneous field.
131    
132        common /averages_r/        common /averages_r/
133       &                    tbar,       &                    tbar,
134       &                    sbar,       &                    sbar,
135         &                    sstbar,
136       &                    psbar,       &                    psbar,
137         &                    bpbar,
138       &                    ubar,       &                    ubar,
139       &                    vbar,       &                    vbar,
140       &                    wbar,       &                    wbar,
141       &                    tauxbar,       &                    tauxbar,
142       &                    tauybar,       &                    tauybar,
143       &                    hfluxbar,       &                    hfluxmeanbar,
144       &                    sfluxbar,       &                    sfluxmeanbar,
145       &                    Slmean,       &                    Slmean,
146       &                    Tlmean,       &                    Tlmean,
147       &                    wlmean,       &                    wlmean,
148       &                    Sfmean,       &                    Sfmean,
149       &                    Tfmean,                         &                    Tfmean,
150       &                    sbar_gen,       &                    sbar_gen,
151       &                    tbar_gen,       &                    tbar_gen,
152       &                    wfmean       &                    wfmean
# Line 101  c             intantaneous field. Line 172  c             intantaneous field.
172        _RL sbar_gen        _RL sbar_gen
173  #endif  #endif
174    
175    #ifdef ALLOW_DAILYSST_COST_CONTRIBUTION
176    cph#ifdef ALLOW_SEAICE_COST_AREASST
177          _RL sstbar  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
178    #else
179          _RL sstbar
180    #endif
181    
182  #if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \  #if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \
183       defined (ALLOW_CTDS_COST_CONTRIBUTION) || \       defined (ALLOW_CTDS_COST_CONTRIBUTION) || \
184       defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \       defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \
# Line 120  c             intantaneous field. Line 198  c             intantaneous field.
198        _RL psbar        _RL psbar
199  #endif  #endif
200    
201    #ifdef ALLOW_BP_COST_CONTRIBUTION
202          _RL bpbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
203    #else
204          _RL bpbar
205    #endif
206    
207  #if (defined (ALLOW_DRIFTER_COST_CONTRIBUTION) || \  #if (defined (ALLOW_DRIFTER_COST_CONTRIBUTION) || \
208       defined (ALLOW_OBCS_COST_CONTRIBUTION))       defined (ALLOW_OBCS_COST_CONTRIBUTION))
209        _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)
# Line 144  c             intantaneous field. Line 228  c             intantaneous field.
228        _RL   Tlmean        _RL   Tlmean
229        _RL   Slmean        _RL   Slmean
230        _RL   Tfmean        _RL   Tfmean
231        _RL   Sfmean                _RL   Sfmean
232  #endif  #endif
233    
234  #ifdef ALLOW_DRIFTW_COST_CONTRIBUTION  #ifdef ALLOW_DRIFTW_COST_CONTRIBUTION
# Line 159  c             intantaneous field. Line 243  c             intantaneous field.
243        _RL tauxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)        _RL tauxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
244        _RL tauybar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)        _RL tauybar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
245  #else  #else
246        _RL tauxbar        _RL tauxbar
247        _RL tauybar        _RL tauybar
248  #endif  #endif
249    
250  #ifdef ALLOW_MEAN_HFLUX_COST_CONTRIBUTION  #ifdef ALLOW_MEAN_HFLUX_COST_CONTRIBUTION
251        _RL hfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)        _RL hfluxmeanbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
252  #else  #else
253        _RL hfluxbar        _RL hfluxmeanbar
254  #endif  #endif
255    
256  #ifdef ALLOW_MEAN_SFLUX_COST_CONTRIBUTION  #ifdef ALLOW_MEAN_SFLUX_COST_CONTRIBUTION
257        _RL sfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)            _RL sfluxmeanbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
258  #else  #else
259        _RL sfluxbar        _RL sfluxmeanbar
260  #endif  #endif
261    
262    
263        common /averages_c/        common /averages_c/
264       &                    tbarfile,       &                    tbarfile,
265       &                    sbarfile,       &                    sbarfile,
266         &                    sstbarfile,
267       &                    psbarfile,       &                    psbarfile,
268         &                    bpbarfile,
269       &                    ubarfile,       &                    ubarfile,
270       &                    vbarfile,       &                    vbarfile,
271       &                    wbarfile,       &                    wbarfile,
272       &                    tauxbarfile,       &                    tauxbarfile,
273       &                    tauybarfile,       &                    tauybarfile,
274       &                    hfluxbarfile,       &                    hfluxmeanbarfile,
275       &                    sfluxbarfile       &                    sfluxmeanbarfile,
276         &                    costTranspDataFile
277        character*(MAX_LEN_FNAM) tbarfile        character*(MAX_LEN_FNAM) tbarfile
278        character*(MAX_LEN_FNAM) sbarfile        character*(MAX_LEN_FNAM) sbarfile
279          character*(MAX_LEN_FNAM) sstbarfile
280        character*(MAX_LEN_FNAM) psbarfile        character*(MAX_LEN_FNAM) psbarfile
281          character*(MAX_LEN_FNAM) bpbarfile
282        character*(MAX_LEN_FNAM) ubarfile        character*(MAX_LEN_FNAM) ubarfile
283        character*(MAX_LEN_FNAM) vbarfile        character*(MAX_LEN_FNAM) vbarfile
284        character*(MAX_LEN_FNAM) wbarfile        character*(MAX_LEN_FNAM) wbarfile
285        character*(MAX_LEN_FNAM) tauxbarfile        character*(MAX_LEN_FNAM) tauxbarfile
286        character*(MAX_LEN_FNAM) tauybarfile        character*(MAX_LEN_FNAM) tauybarfile
287        character*(MAX_LEN_FNAM) hfluxbarfile        character*(MAX_LEN_FNAM) hfluxmeanbarfile
288        character*(MAX_LEN_FNAM) sfluxbarfile        character*(MAX_LEN_FNAM) sfluxmeanbarfile
289          character*(MAX_LEN_FNAM) costTranspDataFile
290    
291    #ifdef ALLOW_TRANSPORT_COST_CONTRIBUTION
292          common /averages_transp_r/
293         &                     transpbar
294         &                   , transpobs
295         &                   , wtransp
296          _RL transpbar(maxNumDays,nsx,nsy)
297          _RL transpobs(maxNumDays)
298          _RL wtransp(maxNumDays)
299    #endif
300    
301  c     file precision and field type  c     file precision and field type
302    
303        common /prec_type_cost/        common /prec_type_cost/
304       &                        cost_iprec,       &                        cost_iprec,
305       &                        cost_yftype       &                        cost_yftype
306    
# Line 254  c     objf_hfluxmm    - Global time-mean Line 354  c     objf_hfluxmm    - Global time-mean
354  c     objf_sfluxmm    - Global time-mean Salt flux contribution.  c     objf_sfluxmm    - Global time-mean Salt flux contribution.
355  c     objf_hmean - Mean sea surface height contribution.  c     objf_hmean - Mean sea surface height contribution.
356  c     objf_h     - Residual sea surface height contribution.  c     objf_h     - Residual sea surface height contribution.
357    c     objf_tp    - Residual sea surface height contribution from T/P
358    c     objf_ers   - Residual sea surface height contribution from T/P
359    c     objf_gfo   - Residual sea surface height contribution from T/P
360  c     objf_temp  - Temperature contribution.  c     objf_temp  - Temperature contribution.
361  c     objf_salt  - Salinity contribution.  c     objf_salt  - Salinity contribution.
362  c     objf_temp0 - Initial conditions Temperature contribution.  c     objf_temp0 - Initial conditions Temperature contribution.
# Line 261  c     objf_salt0 - Initial conditions Sa Line 364  c     objf_salt0 - Initial conditions Sa
364  c     objf_sst   - Sea surface temperature contribution.  c     objf_sst   - Sea surface temperature contribution.
365  c     objf_tmi   - Sea surface temperature contribution.  c     objf_tmi   - Sea surface temperature contribution.
366  c     objf_sss   - Sea surface salinity contribution.  c     objf_sss   - Sea surface salinity contribution.
367  c     objf_ctdt  - Temperature measurements from Woce CTD  c     objf_ctdt  - Temperature measurements from Woce CTD
368  c     objf_ctds  - Salinity measurements from Woce CTD  c     objf_ctds  - Salinity measurements from Woce CTD
369  c     objf_ctdtclim - Temperature measurements from Woce CTD without timetag  c     objf_ctdtclim - Temperature measurements from Woce CTD without timetag
370  c     objf_ctdsclim - Salinity measurements from Woce CTD without timetag  c     objf_ctdsclim - Salinity measurements from Woce CTD without timetag
371  c     objf_xbt   - XBT temperature data  c     objf_xbt   - XBT temperature data
372  c     objf_argot - ARGO temperature profiles  c     objf_argot - ARGO temperature profiles
373  c     objf_argos - ARGO salt profiles  c     objf_argos - ARGO salt profiles
374    c     objf_usercost - user defined cost contribution
375  c     objf_scatxm - time-mean zonal SCAT  contribution  c     objf_scatxm - time-mean zonal SCAT  contribution
376  c     objf_scatym - time-mean meridional SCAT  contribution  c     objf_scatym - time-mean meridional SCAT  contribution
377  c     objf_scatx  - zonal SCAT  contribution  c     objf_scatx  - zonal SCAT  contribution
378  c     objf_scaty  - meridional SCAT  contribution  c     objf_scaty  - meridional SCAT  contribution
379  c     objf_kapgm  - kappa GM contribution  c     objf_kapgm  - kappa GM contribution
380    c     objf_kapredi  - kappa REDI contribution
381  c     objf_diffkr - diffusion contribution  c     objf_diffkr - diffusion contribution
382  c     objf_theta_ini_fin - final vs. initial theta misfit  c     objf_theta_ini_fin - final vs. initial theta misfit
383  c     objf_salt_ini_fin  - final vs. initial salt misfit  c     objf_salt_ini_fin  - final vs. initial salt misfit
384    c     objf_eddytau - eddy streamfunction contribution
385    c     objf_bottomdrag - bottom drag contribution
386  c  c
387  c     mult_"var" - multipliers for the individual cost  c     mult_"var" - multipliers for the individual cost
388  c                  function contributions.  c                  function contributions.
389    
390        common /ecco_cost_objf/        common /ecco_cost_objf/
391       &                objf_hflux,       &     objf_hflux, objf_hfluxm, objf_hfluxmm, objf_hfluxsmoo,
392       &                objf_hfluxm,       &     objf_sflux, objf_sfluxm, objf_sfluxmm, objf_sfluxsmoo,
393       &                objf_hfluxmm,       &     objf_tauu,  objf_tauum,  objf_tauusmoo,
394       &                objf_hfluxsmoo,       &     objf_tauv,  objf_tauvm,  objf_tauvsmoo,
395       &                objf_sflux,       &     objf_hmean,
396       &                objf_sfluxm,       &     objf_h, objf_tp, objf_ers, objf_gfo,
397       &                objf_sfluxmm,       &     objf_sshv4cost,
398       &                objf_sfluxsmoo,       &     objf_temp,      objf_salt,
399       &                objf_tauu,       &     objf_temp0,     objf_salt0,
400       &                objf_tauum,       &     objf_temp0smoo, objf_salt0smoo,
401       &                objf_tauusmoo,       &     objf_sst, objf_tmi, objf_sss,
402       &                objf_tauv,       &     objf_bp,
403       &                objf_tauvm,       &     objf_usercost,
404       &                objf_tauvsmoo,       &     objf_ctdt,      objf_ctds,
405       &                objf_hmean,       &     objf_ctdtclim,  objf_ctdsclim,
406       &                objf_h,       &     objf_xbt, objf_argot,     objf_argos,
407       &                objf_temp,       &     objf_drift, objf_tdrift, objf_sdrift, objf_wdrift,
408       &                objf_salt,       &     objf_scatx, objf_scaty,  objf_scatxm, objf_scatym,
409       &                objf_temp0,       &     objf_atemp,      objf_atempm,      objf_atempsmoo,
410       &                objf_salt0,       &     objf_aqh,        objf_aqhm,        objf_aqhsmoo,
411       &                objf_temp0smoo,       &     objf_precip,     objf_precipm,     objf_precipsmoo,
412       &                objf_salt0smoo,       &     objf_swflux,     objf_swfluxm,     objf_swfluxsmoo,
413       &                objf_sst,       &     objf_swdown,     objf_swdownm,     objf_swdownsmoo,
414       &                objf_tmi,       &     objf_snowprecip, objf_snowprecipm, objf_snowprecipsmoo,
415       &                objf_sss,       &     objf_lwflux,     objf_lwfluxm,     objf_lwfluxsmoo,
416       &                objf_ctdt,       &     objf_lwdown,     objf_lwdownm,     objf_lwdownsmoo,
417       &                objf_ctds,       &     objf_evap,       objf_evapm,       objf_evapsmoo,
418       &                objf_ctdtclim,       &     objf_apressure,  objf_apressurem,  objf_apressuresmoo,
419       &                objf_ctdsclim,       &     objf_runoff,     objf_runoffm,     objf_runoffsmoo,
420       &                objf_xbt,       &     objf_uwind,      objf_uwindm,      objf_uwindsmoo,
421       &                objf_argot,       &     objf_vwind,      objf_vwindm,      objf_vwindsmoo,
422       &                objf_argos,       &     objf_obcsn, objf_obcss, objf_obcsw, objf_obcse,
423       &                objf_drift,       &     objf_obcsvol,
424       &                objf_tdrift,       &     objf_curmtr,
425       &                objf_sdrift,       &     objf_ageos,
426       &                objf_wdrift,       &     objf_kapgm,
427       &                objf_scatx,       &     objf_kapredi,
428       &                objf_scaty,       &     objf_diffkr,
429       &                objf_scatxm,       &     objf_theta_ini_fin, objf_salt_ini_fin,
430       &                objf_scatym,       &     objf_eddytau,
431       &                objf_atemp,       &     objf_bottomdrag,
432       &                objf_aqh,       &     objf_transp
      &                objf_precip,  
      &                objf_swflux,  
      &                objf_swdown,  
      &                objf_uwind,  
      &                objf_vwind,  
      &                objf_atempm,  
      &                objf_aqhm,  
      &                objf_precipm,  
      &                objf_swfluxm,  
      &                objf_swdownm,  
      &                objf_uwindm,  
      &                objf_vwindm,  
      &                objf_atempsmoo,  
      &                objf_aqhsmoo,  
      &                objf_precipsmoo,  
      &                objf_swfluxsmoo,  
      &                objf_swdownsmoo,  
      &                objf_uwindsmoo,  
      &                objf_vwindsmoo,  
      &                objf_obcsn,  
      &                objf_obcss,  
      &                objf_obcsw,  
      &                objf_obcse,  
      &                objf_obcsvol,  
      &                objf_curmtr,  
      &                objf_ageos,  
      &                objf_kapgm,  
      &                objf_diffkr,  
      &                objf_theta_ini_fin,  
      &                objf_salt_ini_fin,  
      &                objf_eddytau  
433    
434        _RL  objf_hflux  (nsx,nsy)        _RL  objf_hflux  (nsx,nsy)
435        _RL  objf_hfluxm (nsx,nsy)        _RL  objf_hfluxm (nsx,nsy)
436        _RL  objf_hfluxmm(nsx,nsy)        _RL  objf_hfluxmm
437        _RL  objf_hfluxsmoo  (nsx,nsy)        _RL  objf_hfluxsmoo  (nsx,nsy)
438        _RL  objf_sflux  (nsx,nsy)        _RL  objf_sflux  (nsx,nsy)
439        _RL  objf_sfluxm (nsx,nsy)        _RL  objf_sfluxm (nsx,nsy)
440        _RL  objf_sfluxmm(nsx,nsy)        _RL  objf_sfluxmm
441        _RL  objf_sfluxsmoo  (nsx,nsy)        _RL  objf_sfluxsmoo  (nsx,nsy)
442        _RL  objf_tauu   (nsx,nsy)        _RL  objf_tauu   (nsx,nsy)
443        _RL  objf_tauum  (nsx,nsy)        _RL  objf_tauum  (nsx,nsy)
# Line 371  c                  function contribution Line 447  c                  function contribution
447        _RL  objf_tauvsmoo   (nsx,nsy)        _RL  objf_tauvsmoo   (nsx,nsy)
448        _RL  objf_hmean        _RL  objf_hmean
449        _RL  objf_h    (nsx,nsy)        _RL  objf_h    (nsx,nsy)
450          _RL  objf_tp   (nsx,nsy)
451          _RL  objf_ers  (nsx,nsy)
452          _RL  objf_gfo  (nsx,nsy)
453          _RL  objf_sshv4cost(NSSHV4COST,nsx,nsy)
454        _RL  objf_temp (nsx,nsy)        _RL  objf_temp (nsx,nsy)
455        _RL  objf_salt (nsx,nsy)        _RL  objf_salt (nsx,nsy)
456        _RL  objf_temp0(nsx,nsy)        _RL  objf_temp0(nsx,nsy)
# Line 379  c                  function contribution Line 459  c                  function contribution
459        _RL  objf_salt0smoo(nsx,nsy)        _RL  objf_salt0smoo(nsx,nsy)
460        _RL  objf_sst  (nsx,nsy)        _RL  objf_sst  (nsx,nsy)
461        _RL  objf_tmi  (nsx,nsy)        _RL  objf_tmi  (nsx,nsy)
462        _RL  objf_sss  (nsx,nsy)        _RL  objf_sss  (nsx,nsy)
463          _RL  objf_bp   (nsx,nsy)
464        _RL  objf_ctdt (nsx,nsy)        _RL  objf_ctdt (nsx,nsy)
465        _RL  objf_ctds (nsx,nsy)        _RL  objf_ctds (nsx,nsy)
466        _RL  objf_ctdtclim (nsx,nsy)        _RL  objf_ctdtclim (nsx,nsy)
467        _RL  objf_ctdsclim (nsx,nsy)        _RL  objf_ctdsclim (nsx,nsy)
468        _RL  objf_xbt  (nsx,nsy)        _RL  objf_xbt  (nsx,nsy)
469        _RL  objf_argot(nsx,nsy)        _RL  objf_argot(nsx,nsy)
470          _RL  objf_usercost(NUSERCOST,nsx,nsy)
471        _RL  objf_argos(nsx,nsy)        _RL  objf_argos(nsx,nsy)
472        _RL  objf_drift(nsx,nsy)        _RL  objf_drift(nsx,nsy)
473        _RL  objf_tdrift(nsx,nsy)        _RL  objf_tdrift(nsx,nsy)
# Line 400  c                  function contribution Line 482  c                  function contribution
482        _RL  objf_precip(nsx,nsy)        _RL  objf_precip(nsx,nsy)
483        _RL  objf_swflux(nsx,nsy)        _RL  objf_swflux(nsx,nsy)
484        _RL  objf_swdown(nsx,nsy)        _RL  objf_swdown(nsx,nsy)
485          _RL  objf_snowprecip(nsx,nsy)
486          _RL  objf_lwflux(nsx,nsy)
487          _RL  objf_lwdown(nsx,nsy)
488          _RL  objf_evap(nsx,nsy)
489          _RL  objf_apressure(nsx,nsy)
490          _RL  objf_runoff(nsx,nsy)
491        _RL  objf_uwind(nsx,nsy)        _RL  objf_uwind(nsx,nsy)
492        _RL  objf_vwind(nsx,nsy)        _RL  objf_vwind(nsx,nsy)
493        _RL  objf_atempm(nsx,nsy)        _RL  objf_atempm(nsx,nsy)
# Line 407  c                  function contribution Line 495  c                  function contribution
495        _RL  objf_precipm(nsx,nsy)        _RL  objf_precipm(nsx,nsy)
496        _RL  objf_swfluxm(nsx,nsy)        _RL  objf_swfluxm(nsx,nsy)
497        _RL  objf_swdownm(nsx,nsy)        _RL  objf_swdownm(nsx,nsy)
498          _RL  objf_snowprecipm(nsx,nsy)
499          _RL  objf_lwfluxm(nsx,nsy)
500          _RL  objf_lwdownm(nsx,nsy)
501          _RL  objf_evapm(nsx,nsy)
502          _RL  objf_apressurem(nsx,nsy)
503          _RL  objf_runoffm(nsx,nsy)
504        _RL  objf_uwindm(nsx,nsy)        _RL  objf_uwindm(nsx,nsy)
505        _RL  objf_vwindm(nsx,nsy)        _RL  objf_vwindm(nsx,nsy)
506        _RL  objf_atempsmoo(nsx,nsy)        _RL  objf_atempsmoo(nsx,nsy)
# Line 414  c                  function contribution Line 508  c                  function contribution
508        _RL  objf_precipsmoo(nsx,nsy)        _RL  objf_precipsmoo(nsx,nsy)
509        _RL  objf_swfluxsmoo(nsx,nsy)        _RL  objf_swfluxsmoo(nsx,nsy)
510        _RL  objf_swdownsmoo(nsx,nsy)        _RL  objf_swdownsmoo(nsx,nsy)
511          _RL  objf_snowprecipsmoo(nsx,nsy)
512          _RL  objf_lwfluxsmoo(nsx,nsy)
513          _RL  objf_lwdownsmoo(nsx,nsy)
514          _RL  objf_evapsmoo(nsx,nsy)
515          _RL  objf_apressuresmoo(nsx,nsy)
516          _RL  objf_runoffsmoo(nsx,nsy)
517        _RL  objf_uwindsmoo(nsx,nsy)        _RL  objf_uwindsmoo(nsx,nsy)
518        _RL  objf_vwindsmoo(nsx,nsy)        _RL  objf_vwindsmoo(nsx,nsy)
519        _RL  objf_obcsn(nsx,nsy)        _RL  objf_obcsn(nsx,nsy)
# Line 424  c                  function contribution Line 524  c                  function contribution
524        _RL  objf_curmtr(nsx,nsy)        _RL  objf_curmtr(nsx,nsy)
525        _RL  objf_ageos(nsx,nsy)        _RL  objf_ageos(nsx,nsy)
526        _RL  objf_kapgm(nsx,nsy)        _RL  objf_kapgm(nsx,nsy)
527          _RL  objf_kapredi(nsx,nsy)
528        _RL  objf_diffkr(nsx,nsy)        _RL  objf_diffkr(nsx,nsy)
529        _RL  objf_theta_ini_fin(nsx,nsy)        _RL  objf_theta_ini_fin(nsx,nsy)
530        _RL  objf_salt_ini_fin(nsx,nsy)        _RL  objf_salt_ini_fin(nsx,nsy)
531        _RL  objf_eddytau(nsx,nsy)        _RL  objf_eddytau(nsx,nsy)
532          _RL  objf_bottomdrag(nsx,nsy)
533          _RL  objf_transp
534    
535        common /ecco_cost_num/        common /ecco_cost_num/
536       &                num_hflux,       &                num_hflux,
# Line 442  c                  function contribution Line 545  c                  function contribution
545       &                num_tauvm,       &                num_tauvm,
546       &                num_hmean,       &                num_hmean,
547       &                num_h,       &                num_h,
548         &                num_tp,
549         &                num_ers,
550         &                num_gfo,
551         &                num_sshv4cost,
552       &                num_temp,       &                num_temp,
553       &                num_salt,       &                num_salt,
554       &                num_temp0,       &                num_temp0,
# Line 449  c                  function contribution Line 556  c                  function contribution
556       &                num_sst,       &                num_sst,
557       &                num_tmi,       &                num_tmi,
558       &                num_sss,       &                num_sss,
559         &                num_bp,
560       &                num_ctdt,       &                num_ctdt,
561       &                num_ctds,       &                num_ctds,
562       &                num_ctdtclim,       &                num_ctdtclim,
# Line 456  c                  function contribution Line 564  c                  function contribution
564       &                num_xbt,       &                num_xbt,
565       &                num_argot,       &                num_argot,
566       &                num_argos,       &                num_argos,
567         &                num_usercost,
568       &                num_drift,       &                num_drift,
569       &                num_tdrift,       &                num_tdrift,
570       &                num_sdrift,       &                num_sdrift,
# Line 469  c                  function contribution Line 578  c                  function contribution
578       &                num_precip,       &                num_precip,
579       &                num_swflux,       &                num_swflux,
580       &                num_swdown,       &                num_swdown,
581         &                num_snowprecip,
582         &                num_lwflux,
583         &                num_lwdown,
584         &                num_evap,
585         &                num_apressure,
586         &                num_runoff,
587       &                num_uwind,       &                num_uwind,
588       &                num_vwind,       &                num_vwind,
589       &                num_atempm,       &                num_atempm,
# Line 476  c                  function contribution Line 591  c                  function contribution
591       &                num_precipm,       &                num_precipm,
592       &                num_swfluxm,       &                num_swfluxm,
593       &                num_swdownm,       &                num_swdownm,
594         &                num_snowprecipm,
595         &                num_lwfluxm,
596         &                num_lwdownm,
597         &                num_evapm,
598         &                num_apressurem,
599         &                num_runoffm,
600       &                num_uwindm,       &                num_uwindm,
601       &                num_vwindm,       &                num_vwindm,
602       &                num_obcsn,       &                num_obcsn,
# Line 486  c                  function contribution Line 607  c                  function contribution
607       &                num_curmtr,       &                num_curmtr,
608       &                num_ageos,       &                num_ageos,
609       &                num_kapgm,       &                num_kapgm,
610         &                num_kapredi,
611       &                num_diffkr,       &                num_diffkr,
612       &                num_theta_ini_fin,       &                num_theta_ini_fin,
613       &                num_salt_ini_fin,       &                num_salt_ini_fin,
614       &                num_eddytau       &                num_eddytau,
615         &                num_bottomdrag,
616         &                num_transp
617    
618        _RL  num_hflux  (nsx,nsy)        _RL  num_hflux  (nsx,nsy)
619        _RL  num_hfluxm (nsx,nsy)        _RL  num_hfluxm (nsx,nsy)
# Line 503  c                  function contribution Line 627  c                  function contribution
627        _RL  num_tauvm  (nsx,nsy)        _RL  num_tauvm  (nsx,nsy)
628        _RL  num_hmean        _RL  num_hmean
629        _RL  num_h    (nsx,nsy)        _RL  num_h    (nsx,nsy)
630          _RL  num_tp   (nsx,nsy)
631          _RL  num_ers  (nsx,nsy)
632          _RL  num_gfo  (nsx,nsy)
633          _RL  num_sshv4cost(NSSHV4COST,nsx,nsy)
634        _RL  num_temp (nsx,nsy)        _RL  num_temp (nsx,nsy)
635        _RL  num_salt (nsx,nsy)        _RL  num_salt (nsx,nsy)
636        _RL  num_temp0(nsx,nsy)        _RL  num_temp0(nsx,nsy)
637        _RL  num_salt0(nsx,nsy)        _RL  num_salt0(nsx,nsy)
638        _RL  num_sst  (nsx,nsy)        _RL  num_sst  (nsx,nsy)
639        _RL  num_tmi  (nsx,nsy)        _RL  num_tmi  (nsx,nsy)
640        _RL  num_sss  (nsx,nsy)        _RL  num_sss  (nsx,nsy)
641          _RL  num_bp   (nsx,nsy)
642        _RL  num_ctdt (nsx,nsy)        _RL  num_ctdt (nsx,nsy)
643        _RL  num_ctds (nsx,nsy)        _RL  num_ctds (nsx,nsy)
644        _RL  num_ctdtclim (nsx,nsy)        _RL  num_ctdtclim (nsx,nsy)
# Line 517  c                  function contribution Line 646  c                  function contribution
646        _RL  num_xbt  (nsx,nsy)        _RL  num_xbt  (nsx,nsy)
647        _RL  num_argot(nsx,nsy)        _RL  num_argot(nsx,nsy)
648        _RL  num_argos(nsx,nsy)        _RL  num_argos(nsx,nsy)
649          _RL  num_usercost(NUSERCOST,nsx,nsy)
650        _RL  num_drift(nsx,nsy)        _RL  num_drift(nsx,nsy)
651        _RL  num_tdrift(nsx,nsy)        _RL  num_tdrift(nsx,nsy)
652        _RL  num_sdrift(nsx,nsy)        _RL  num_sdrift(nsx,nsy)
# Line 530  c                  function contribution Line 660  c                  function contribution
660        _RL  num_precip(nsx,nsy)        _RL  num_precip(nsx,nsy)
661        _RL  num_swflux(nsx,nsy)        _RL  num_swflux(nsx,nsy)
662        _RL  num_swdown(nsx,nsy)        _RL  num_swdown(nsx,nsy)
663          _RL  num_snowprecip(nsx,nsy)
664          _RL  num_lwflux(nsx,nsy)
665          _RL  num_lwdown(nsx,nsy)
666          _RL  num_evap(nsx,nsy)
667          _RL  num_apressure(nsx,nsy)
668          _RL  num_runoff(nsx,nsy)
669        _RL  num_uwind(nsx,nsy)        _RL  num_uwind(nsx,nsy)
670        _RL  num_vwind(nsx,nsy)        _RL  num_vwind(nsx,nsy)
671        _RL  num_atempm(nsx,nsy)        _RL  num_atempm(nsx,nsy)
# Line 537  c                  function contribution Line 673  c                  function contribution
673        _RL  num_precipm(nsx,nsy)        _RL  num_precipm(nsx,nsy)
674        _RL  num_swfluxm(nsx,nsy)        _RL  num_swfluxm(nsx,nsy)
675        _RL  num_swdownm(nsx,nsy)        _RL  num_swdownm(nsx,nsy)
676          _RL  num_snowprecipm(nsx,nsy)
677          _RL  num_lwfluxm(nsx,nsy)
678          _RL  num_lwdownm(nsx,nsy)
679          _RL  num_evapm(nsx,nsy)
680          _RL  num_apressurem(nsx,nsy)
681          _RL  num_runoffm(nsx,nsy)
682        _RL  num_uwindm(nsx,nsy)        _RL  num_uwindm(nsx,nsy)
683        _RL  num_vwindm(nsx,nsy)        _RL  num_vwindm(nsx,nsy)
684        _RL  num_obcsn(nsx,nsy)        _RL  num_obcsn(nsx,nsy)
# Line 547  c                  function contribution Line 689  c                  function contribution
689        _RL  num_curmtr(nsx,nsy)        _RL  num_curmtr(nsx,nsy)
690        _RL  num_ageos(nsx,nsy)        _RL  num_ageos(nsx,nsy)
691        _RL  num_kapgm(nsx,nsy)        _RL  num_kapgm(nsx,nsy)
692          _RL  num_kapredi(nsx,nsy)
693        _RL  num_diffkr(nsx,nsy)        _RL  num_diffkr(nsx,nsy)
694        _RL  num_theta_ini_fin(nsx,nsy)        _RL  num_theta_ini_fin(nsx,nsy)
695        _RL  num_salt_ini_fin(nsx,nsy)        _RL  num_salt_ini_fin(nsx,nsy)
696        _RL  num_eddytau(nsx,nsy)        _RL  num_eddytau(nsx,nsy)
697          _RL  num_bottomdrag(nsx,nsy)
698          _RL  num_transp
699    
700        common /ecco_cost_aux_r/        common /ecco_cost_aux_r/
701       &                    mult_hflux,       &                    mult_hflux,
# Line 561  c                  function contribution Line 706  c                  function contribution
706       &                    mult_tauv,       &                    mult_tauv,
707       &                    mult_hmean,       &                    mult_hmean,
708       &                    mult_h,       &                    mult_h,
709         &                    mult_tp,
710         &                    mult_ers,
711         &                    mult_gfo,
712         &                    mult_sshv4cost,
713       &                    mult_temp,       &                    mult_temp,
714       &                    mult_salt,       &                    mult_salt,
715       &                    mult_temp0,       &                    mult_temp0,
# Line 568  c                  function contribution Line 717  c                  function contribution
717       &                    mult_sst,       &                    mult_sst,
718       &                    mult_tmi,       &                    mult_tmi,
719       &                    mult_sss,       &                    mult_sss,
720         &                    mult_bp,
721       &                    mult_ctdt,       &                    mult_ctdt,
722       &                    mult_ctds,       &                    mult_ctds,
723       &                    mult_ctdtclim,       &                    mult_ctdtclim,
# Line 575  c                  function contribution Line 725  c                  function contribution
725       &                    mult_xbt,       &                    mult_xbt,
726       &                    mult_argot,       &                    mult_argot,
727       &                    mult_argos,       &                    mult_argos,
728         &                    mult_usercost,
729       &                    mult_drift,       &                    mult_drift,
730       &                    mult_tdrift,       &                    mult_tdrift,
731       &                    mult_sdrift,       &                    mult_sdrift,
# Line 586  c                  function contribution Line 737  c                  function contribution
737       &                    mult_precip,       &                    mult_precip,
738       &                    mult_swflux,       &                    mult_swflux,
739       &                    mult_swdown,       &                    mult_swdown,
740         &                    mult_snowprecip,
741         &                    mult_lwflux,
742         &                    mult_lwdown,
743         &                    mult_evap,
744         &                    mult_apressure,
745         &                    mult_runoff,
746       &                    mult_uwind,       &                    mult_uwind,
747       &                    mult_vwind,       &                    mult_vwind,
748       &                    mult_obcsn,       &                    mult_obcsn,
# Line 596  c                  function contribution Line 753  c                  function contribution
753       &                    mult_curmtr,       &                    mult_curmtr,
754       &                    mult_ageos,       &                    mult_ageos,
755       &                    mult_kapgm,       &                    mult_kapgm,
756         &                    mult_kapredi,
757       &                    mult_diffkr,       &                    mult_diffkr,
758       &                    mult_ini_fin,       &                    mult_ini_fin,
759       &                    mult_eddytau,       &                    mult_edtau,
760         &                    mult_bottomdrag,
761       &                    mult_smooth_ic,       &                    mult_smooth_ic,
762       &                    mult_smooth_bc       &                    mult_smooth_bc,
763         &                    mult_transp
764        _RL  mult_hflux        _RL  mult_hflux
765        _RL  mult_sflux        _RL  mult_sflux
766        _RL  mult_hfluxmm        _RL  mult_hfluxmm
# Line 610  c                  function contribution Line 769  c                  function contribution
769        _RL  mult_tauv        _RL  mult_tauv
770        _RL  mult_hmean        _RL  mult_hmean
771        _RL  mult_h        _RL  mult_h
772          _RL  mult_tp
773          _RL  mult_ers
774          _RL  mult_gfo
775          _RL  mult_sshv4cost(NSSHV4COST)
776        _RL  mult_temp        _RL  mult_temp
777        _RL  mult_salt        _RL  mult_salt
778        _RL  mult_temp0        _RL  mult_temp0
# Line 617  c                  function contribution Line 780  c                  function contribution
780        _RL  mult_sst        _RL  mult_sst
781        _RL  mult_tmi        _RL  mult_tmi
782        _RL  mult_sss        _RL  mult_sss
783          _RL  mult_bp
784        _RL  mult_ctdt        _RL  mult_ctdt
785        _RL  mult_ctds        _RL  mult_ctds
786        _RL  mult_ctdtclim        _RL  mult_ctdtclim
# Line 624  c                  function contribution Line 788  c                  function contribution
788        _RL  mult_xbt        _RL  mult_xbt
789        _RL  mult_argot        _RL  mult_argot
790        _RL  mult_argos        _RL  mult_argos
791          _RL  mult_usercost(NUSERCOST)
792        _RL  mult_drift        _RL  mult_drift
793        _RL  mult_tdrift        _RL  mult_tdrift
794        _RL  mult_sdrift        _RL  mult_sdrift
# Line 635  c                  function contribution Line 800  c                  function contribution
800        _RL  mult_precip        _RL  mult_precip
801        _RL  mult_swflux        _RL  mult_swflux
802        _RL  mult_swdown        _RL  mult_swdown
803          _RL  mult_snowprecip
804          _RL  mult_lwflux
805          _RL  mult_lwdown
806          _RL  mult_evap
807          _RL  mult_apressure
808          _RL  mult_runoff
809        _RL  mult_uwind        _RL  mult_uwind
810        _RL  mult_vwind        _RL  mult_vwind
811        _RL  mult_obcsn        _RL  mult_obcsn
# Line 645  c                  function contribution Line 816  c                  function contribution
816        _RL  mult_curmtr        _RL  mult_curmtr
817        _RL  mult_ageos        _RL  mult_ageos
818        _RL  mult_kapgm        _RL  mult_kapgm
819          _RL  mult_kapredi
820        _RL  mult_diffkr        _RL  mult_diffkr
821        _RL  mult_ini_fin        _RL  mult_ini_fin
822        _RL  mult_eddytau        _RL  mult_edtau
823          _RL  mult_bottomdrag
824        _RL  mult_smooth_ic        _RL  mult_smooth_ic
825        _RL  mult_smooth_bc        _RL  mult_smooth_bc
826          _RL  mult_transp
827    
828  c     Record counters relevant for the cost function evaluation.  c     Record counters relevant for the cost function evaluation.
829  c     ==========================================================  c     ==========================================================
830  c  c
831  c     nmonsrec - number of monthly records that will be generated by  c     nyearsrec - number of yearly records that will be generated by
832  c                the current model integration.  c                 the current model integration.
833  c     ndaysrec - number of  daily  records that will be generated by  c     nmonsrec  - number of monthly records that will be generated by
834  c                the current model integration.  c                 the current model integration.
835    c     ndaysrec  - number of  daily  records that will be generated by
836    c                 the current model integration.
837    
838        common /ecco_cost_i/        common /ecco_cost_i/
839         &                nyearsrec,
840       &                nmonsrec,       &                nmonsrec,
841       &                ndaysrec,       &                ndaysrec,
842       &                nnztbar,       &                nnztbar,
843       &                nnzsbar       &                nnzsbar
844          integer nyearsrec
845        integer nmonsrec        integer nmonsrec
846        integer ndaysrec        integer ndaysrec
847        integer nnztbar        integer nnztbar
# Line 692  c     salterrfile           - representa Line 870  c     salterrfile           - representa
870  c     temperrfile           - representation error due unresolved eddies  c     temperrfile           - representation error due unresolved eddies
871  c     velerrfile            - representation error  c     velerrfile            - representation error
872    
873        common /ecco_cost_c/        common /ecco_cost_c/
874       &                hflux_errfile,       &                hflux_errfile,
875       &                hfluxm_errfile,       &                hfluxm_errfile,
876       &                sflux_errfile,       &                sflux_errfile,
# Line 707  c     velerrfile            - representa Line 885  c     velerrfile            - representa
885       &                geoid_errfile,       &                geoid_errfile,
886       &                geoid_covariancefile,       &                geoid_covariancefile,
887       &                ssh_errfile,       &                ssh_errfile,
888         &                tp_errfile,
889         &                ers_errfile,
890         &                gfo_errfile,
891         &                sshv4cost_errfile,
892       &                ctdt_errfile,       &                ctdt_errfile,
893       &                ctds_errfile,       &                ctds_errfile,
894       &                drift_errfile,       &                drift_errfile,
895       &                udrifterrfile,       &                udrifterrfile,
896       &                vdrifterrfile,       &                vdrifterrfile,
897       &                salterrfile,       &                salterrfile,
898       &                temperrfile,       &                temperrfile,
899       &                velerrfile,       &                velerrfile,
900       &                salt0errfile,       &                salt0errfile,
901       &                temp0errfile,       &                temp0errfile,
902       &                vel0errfile,       &                vel0errfile,
903         &                ssterrfile,
904         &                ssserrfile,
905         &                bperrfile,
906       &                atemp_errfile,       &                atemp_errfile,
907       &                aqh_errfile,       &                aqh_errfile,
908       &                precip_errfile,       &                precip_errfile,
909       &                swflux_errfile,       &                swflux_errfile,
910       &                swdown_errfile,       &                swdown_errfile,
911         &                snowprecip_errfile,
912         &                lwflux_errfile,
913         &                lwdown_errfile,
914         &                evap_errfile,
915         &                apressure_errfile,
916         &                runoff_errfile,
917         &                edtau_errfile,
918         &                kapgm_errfile,
919         &                kapredi_errfile,
920         &                diffkr_errfile,
921         &                bottomdrag_errfile,
922         &                usercost_errfile,
923       &                uwind_errfile,       &                uwind_errfile,
924       &                vwind_errfile       &                vwind_errfile
925        character*(MAX_LEN_FNAM) hflux_errfile        character*(MAX_LEN_FNAM) hflux_errfile
# Line 739  c     velerrfile            - representa Line 936  c     velerrfile            - representa
936        character*(MAX_LEN_FNAM) geoid_errfile        character*(MAX_LEN_FNAM) geoid_errfile
937        character*(MAX_LEN_FNAM) geoid_covariancefile        character*(MAX_LEN_FNAM) geoid_covariancefile
938        character*(MAX_LEN_FNAM) ssh_errfile        character*(MAX_LEN_FNAM) ssh_errfile
939        character*(MAX_LEN_FNAM) ctdt_errfile        character*(MAX_LEN_FNAM) tp_errfile
940        character*(MAX_LEN_FNAM) ctds_errfile        character*(MAX_LEN_FNAM) ers_errfile
941          character*(MAX_LEN_FNAM) gfo_errfile
942          character*(MAX_LEN_FNAM) sshv4cost_errfile(NSSHV4COST)
943          character*(MAX_LEN_FNAM) ctdt_errfile
944          character*(MAX_LEN_FNAM) ctds_errfile
945        character*(MAX_LEN_FNAM) drift_errfile        character*(MAX_LEN_FNAM) drift_errfile
946        character*(MAX_LEN_FNAM) udrifterrfile        character*(MAX_LEN_FNAM) udrifterrfile
947        character*(MAX_LEN_FNAM) vdrifterrfile              character*(MAX_LEN_FNAM) vdrifterrfile
948        character*(MAX_LEN_FNAM) salterrfile        character*(MAX_LEN_FNAM) salterrfile
949        character*(MAX_LEN_FNAM) temperrfile        character*(MAX_LEN_FNAM) temperrfile
950        character*(MAX_LEN_FNAM) velerrfile        character*(MAX_LEN_FNAM) velerrfile
951        character*(MAX_LEN_FNAM) salt0errfile        character*(MAX_LEN_FNAM) salt0errfile
952        character*(MAX_LEN_FNAM) temp0errfile        character*(MAX_LEN_FNAM) temp0errfile
953        character*(MAX_LEN_FNAM) vel0errfile        character*(MAX_LEN_FNAM) vel0errfile
954          character*(MAX_LEN_FNAM) ssterrfile
955          character*(MAX_LEN_FNAM) ssserrfile
956          character*(MAX_LEN_FNAM) bperrfile
957        character*(MAX_LEN_FNAM) atemp_errfile        character*(MAX_LEN_FNAM) atemp_errfile
958        character*(MAX_LEN_FNAM) aqh_errfile        character*(MAX_LEN_FNAM) aqh_errfile
959        character*(MAX_LEN_FNAM) precip_errfile        character*(MAX_LEN_FNAM) precip_errfile
960        character*(MAX_LEN_FNAM) swflux_errfile        character*(MAX_LEN_FNAM) swflux_errfile
961        character*(MAX_LEN_FNAM) swdown_errfile        character*(MAX_LEN_FNAM) swdown_errfile
962          character*(MAX_LEN_FNAM) snowprecip_errfile
963          character*(MAX_LEN_FNAM) lwflux_errfile
964          character*(MAX_LEN_FNAM) lwdown_errfile
965          character*(MAX_LEN_FNAM) evap_errfile
966          character*(MAX_LEN_FNAM) apressure_errfile
967          character*(MAX_LEN_FNAM) runoff_errfile
968          character*(MAX_LEN_FNAM) edtau_errfile
969          character*(MAX_LEN_FNAM) kapgm_errfile
970          character*(MAX_LEN_FNAM) kapredi_errfile
971          character*(MAX_LEN_FNAM) diffkr_errfile
972          character*(MAX_LEN_FNAM) bottomdrag_errfile
973          character*(MAX_LEN_FNAM) usercost_errfile(NUSERCOST)
974        character*(MAX_LEN_FNAM) uwind_errfile        character*(MAX_LEN_FNAM) uwind_errfile
975        character*(MAX_LEN_FNAM) vwind_errfile        character*(MAX_LEN_FNAM) vwind_errfile
976    
# Line 783  c     wctds      - weight for CTD salini Line 999  c     wctds      - weight for CTD salini
999  c     wudrift    - weight for mean zonal velocity from drifters.  c     wudrift    - weight for mean zonal velocity from drifters.
1000  c     wvdrift    - weight for mean meridional velocity from drifters.  c     wvdrift    - weight for mean meridional velocity from drifters.
1001    
1002        common /ecco_cost_weights_r/        common /ecco_cost_weights_r/
1003       &                      frame,       &                      frame,
1004       &                      cosphi,       &                      cosphi,
1005       &                      whflux,wsflux,wtauu,wtauv,       &                      whflux,wsflux,wtauu,wtauv,
1006       &                      watemp,waqh,wprecip,wswflux,wswdown,       &                      watemp,waqh,wprecip,wsnowprecip,
1007         &                      wswflux,wswdown,wlwflux,wlwdown,
1008         &                      wevap,wapressure,wrunoff,
1009         &                      wbottomdrag,
1010       &                      wuwind,wvwind,       &                      wuwind,wvwind,
1011       &                      wscatx,wscaty,       &                      wscatx,wscaty,
1012       &                      wtheta,wtheta2,wthetaLev,       &                      wtheta,wtheta2,wthetaLev,
1013       &                      wsalt,wsalt2,wsaltLev,       &                      wsalt,wsalt2,wsaltLev,
1014       &                      wsst,wsss,       &                      wdiffkr,wdiffkr2,wdiffkrFld,
1015         &                      wkapgm,wkapgm2,wkapgmFld,
1016         &                      wkapredi,wkapredi2,wkaprediFld,
1017         &                      wedtaux,wedtaux2,wedtauxFld,
1018         &                      wedtauy,wedtauy2,wedtauyFld,
1019         &                      wsst,wsss,wbp,
1020       &                      wtp,wers,wgfo,       &                      wtp,wers,wgfo,
1021       &                      wp,       &                      wp,wsshv4,
1022       &                      wctdt,wctds,       &                      wctdt,wctds,
1023       &                      wudrift,wvdrift,       &                      wudrift,wvdrift,
1024       &                      whfluxmm,wsfluxmm,       &                      whfluxmm,wsfluxmm,
# Line 820  c     wvdrift    - weight for mean merid Line 1044  c     wvdrift    - weight for mean merid
1044        _RL wprecip (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wprecip (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1045        _RL wswflux (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wswflux (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1046        _RL wswdown (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wswdown (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1047          _RL wsnowprecip (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1048          _RL wlwflux (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1049          _RL wlwdown (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1050          _RL wevap   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1051          _RL wapressure(1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1052          _RL wrunoff (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1053          _RL wbottomdrag (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
1054        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1055        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1056        _RL wtheta  (                            nr,nsx,nsy)        _RL wtheta  (                            nr,nsx,nsy)
# Line 830  c     wvdrift    - weight for mean merid Line 1061  c     wvdrift    - weight for mean merid
1061        _RL wsaltLev  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL wsaltLev  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1062        _RL wsst    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsst    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1063        _RL wsss    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsss    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1064          _RL wbp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1065        _RL wtp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1066          _RL wsshv4  (1-olx:snx+olx,1-oly:sny+oly,NSSHV4COST,nsx,nsy)
1067        _RL wers    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wers    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1068        _RL wgfo    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wgfo    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1069        _RL wp      (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wp      (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
# Line 842  c     wvdrift    - weight for mean merid Line 1075  c     wvdrift    - weight for mean merid
1075        _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)
1076        _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)
1077        _RL wbaro        _RL wbaro
1078          _RL wdiffkr (                            nr,nsx,nsy)
1079          _RL wdiffkr2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1080          _RL wdiffkrFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1081          _RL wkapgm (                            nr,nsx,nsy)
1082          _RL wkapgm2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1083          _RL wkapgmFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1084          _RL wkapredi (                            nr,nsx,nsy)
1085          _RL wkapredi2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1086          _RL wkaprediFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1087          _RL wedtaux (                            nr,nsx,nsy)
1088          _RL wedtaux2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1089          _RL wedtauxFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1090          _RL wedtauy (                            nr,nsx,nsy)
1091          _RL wedtauy2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1092          _RL wedtauyFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1093    
1094        common /ecco_cost_weights_0_r/        common /ecco_cost_weights_0_r/
1095       &        whflux0, wsflux0, wtau0,       &        whflux0, wsflux0, wtau0,
1096       &        watemp0, waqh0, wprecip0, wswflux0, wswdown0, wwind0       &        watemp0, waqh0, wprecip0, wsnowprecip0, wwind0,
1097         &        wswflux0, wswdown0, wlwflux0, wlwdown0,
1098         &        wevap0, wapressure0, wrunoff0, wkapredi0,
1099         &        wbottomdrag0,wdiffkr0, wkapgm0, wedtau0
1100        _RL whflux0        _RL whflux0
1101        _RL wsflux0        _RL wsflux0
1102        _RL wtau0        _RL wtau0
# Line 854  c     wvdrift    - weight for mean merid Line 1105  c     wvdrift    - weight for mean merid
1105        _RL wprecip0        _RL wprecip0
1106        _RL wswflux0        _RL wswflux0
1107        _RL wswdown0        _RL wswdown0
1108          _RL wsnowprecip0
1109          _RL wlwflux0
1110          _RL wlwdown0
1111          _RL wevap0
1112          _RL wapressure0
1113          _RL wrunoff0
1114          _RL wbottomdrag0
1115        _RL wwind0        _RL wwind0
1116          _RL wdiffkr0
1117          _RL wkapgm0
1118          _RL wkapredi0
1119          _RL wedtau0
1120    
1121        common /ecco_cost_weights_mean_r/        common /ecco_cost_weights_mean_r/
1122       &        wmean_hflux, wmean_sflux, wmean_tau,       &        wmean_hflux, wmean_sflux, wmean_tau,
1123       &        wmean_atemp, wmean_aqh,       &        wmean_atemp, wmean_aqh,
1124       &        wmean_precip, wmean_swflux, wmean_swdown, wmean_wind       &        wmean_precip, wmean_snowprecip, wmean_wind,
1125         &        wmean_swflux, wmean_swdown, wmean_lwflux, wmean_lwdown,
1126         &        wmean_evap, wmean_apressure, wmean_runoff
1127        _RL wmean_hflux        _RL wmean_hflux
1128        _RL wmean_sflux        _RL wmean_sflux
1129        _RL wmean_tau        _RL wmean_tau
# Line 868  c     wvdrift    - weight for mean merid Line 1132  c     wvdrift    - weight for mean merid
1132        _RL wmean_precip        _RL wmean_precip
1133        _RL wmean_swflux        _RL wmean_swflux
1134        _RL wmean_swdown        _RL wmean_swdown
1135          _RL wmean_snowprecip
1136          _RL wmean_lwflux
1137          _RL wmean_lwdown
1138          _RL wmean_evap
1139          _RL wmean_apressure
1140          _RL wmean_runoff
1141        _RL wmean_wind        _RL wmean_wind
1142                                                                                                    
1143        common /ecco_cost_weights_2_r/        common /ecco_cost_weights_2_r/
1144       &                      whflux2,wsflux2,wtauu2,wtauv2       &                      whflux2,wsflux2,wtauu2,wtauv2
1145        _RL whflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL whflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1146        _RL wsflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1147        _RL wtauu2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtauu2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1148        _RL wtauv2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtauv2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1149                                                                                                    
1150                                          
1151    
1152  #if (defined (ALLOW_OBCSN_COST_CONTRIBUTION) || \  #if (defined (ALLOW_OBCSN_COST_CONTRIBUTION) || \
1153       defined (ALLOW_OBCSN_CONTROL))       defined (ALLOW_OBCSN_CONTROL))
# Line 918  c     scatydat   - reference meridional Line 1188  c     scatydat   - reference meridional
1188  c     sstdat     - reference sea surface temperature data.  c     sstdat     - reference sea surface temperature data.
1189  c     tmidat     - reference TMI sea surface temperature data.  c     tmidat     - reference TMI sea surface temperature data.
1190  c     sssdat     - reference sea surface temperature data.  c     sssdat     - reference sea surface temperature data.
1191    c     bpdat      - bottom pressure from time-varying GRACE.
1192  c     tauxmask   - mask for reference wind stress data.  c     tauxmask   - mask for reference wind stress data.
1193  c     tauymask   - mask for reference wind stress data.  c     tauymask   - mask for reference wind stress data.
1194  c     scatxmask  - mask for scat wind stress data.  c     scatxmask  - mask for scat wind stress data.
1195  c     scatymask  - mask for scat wind stress data.  c     scatymask  - mask for scat wind stress data.
1196  c     sstmask    - mask for reference sea surface temperature data.  c     sstmask    - mask for reference sea surface temperature data.
1197  c     tmimask    - mask for reference sea surface temperature data.  c     tmimask    - mask for reference sea surface temperature data.
1198  c     sssmask    - mask for reference sea surface temperature data.  c     sssmask    - mask for reference sea surface temperature data.
# Line 934  c     ersobs     - ERS data. Line 1205  c     ersobs     - ERS data.
1205  c     ersmask    - mask for ERS data.  c     ersmask    - mask for ERS data.
1206  c     ctdtobs    - CTD temperature data  c     ctdtobs    - CTD temperature data
1207  c     ctdsobs    - CTD salinity data  c     ctdsobs    - CTD salinity data
1208  c     xbtobs     - XBT data  c     xbtobs     - XBT data
1209  c     argot      - ARGO  temperature data  c     argot      - ARGO  temperature data
1210  c     argos      - ARGO  salt data  c     argos      - ARGO  salt data
1211  c     udriftdat  - drifters zonal velocities  c     udriftdat  - drifters zonal velocities
1212  c     vdriftdat  - drifters meridional velocities  c     vdriftdat  - drifters meridional velocities
1213    
# Line 944  c     vdriftdat  - drifters meridional v Line 1215  c     vdriftdat  - drifters meridional v
1215       &                     tdat,       &                     tdat,
1216       &                     scatxdat,       &                     scatxdat,
1217       &                     scatydat,       &                     scatydat,
1218       &                     sstdat,       &                     bpdat,
      &                     tmidat,  
      &                     sssdat,  
1219       &                     sstmask,       &                     sstmask,
1220       &                     tmimask,       &                     tmimask,
1221       &                     sssmask,       &                     sssmask,
1222         &                     bpmask,
1223       &                     tauxmask,       &                     tauxmask,
1224       &                     tauymask,       &                     tauymask,
1225       &                     scatxmask,       &                     scatxmask,
# Line 972  c     vdriftdat  - drifters meridional v Line 1242  c     vdriftdat  - drifters meridional v
1242       &                     vdriftdat,       &                     vdriftdat,
1243       &                     curmtruobs,       &                     curmtruobs,
1244       &                     curmtrvobs       &                     curmtrvobs
1245        
1246        _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)
1247        _RL scatxdat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatxdat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1248        _RL scatydat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatydat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1249        _RL sstdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL bpdat     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
       _RL tmidat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)  
       _RL sssdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)  
1250        _RL tauxmask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tauxmask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1251        _RL tauymask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tauymask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1252        _RL scatxmask (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatxmask (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
# Line 986  c     vdriftdat  - drifters meridional v Line 1254  c     vdriftdat  - drifters meridional v
1254        _RL sstmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sstmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1255        _RL tmimask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tmimask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1256        _RL sssmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sssmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1257          _RL bpmask    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1258        _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)
1259        _RL tpmean    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpmean    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1260        _RL tpmeanmask(1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpmeanmask(1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
# Line 1020  c     topexfile     - reference data fil Line 1289  c     topexfile     - reference data fil
1289  c                     (TOPEX/POSEIDON).  c                     (TOPEX/POSEIDON).
1290  c     ersfile       - reference data file for sea surface height data  c     ersfile       - reference data file for sea surface height data
1291  c                     (ERS).  c                     (ERS).
1292  c ctdtfile, ctdsfile- reference data file for temperature and salinity  c ctdtfile, ctdsfile- reference data file for temperature and salinity
1293  c                     from CTD  c                     from CTD
1294  c ctdtclimfile, ctdsclimfile- reference data file for temperature  c ctdtclimfile, ctdsclimfile- reference data file for temperature
1295  c                      and salinity from CTD with out timetag  c                      and salinity from CTD with out timetag
1296  c     xbtfile       - reference data file for xbt  c     xbtfile       - reference data file for xbt
1297  c     ARGOtfile     - reference data file for ARGO  c     ARGOtfile     - reference data file for ARGO
1298  c     ARGOsfile     - reference data file for ARGO  c     ARGOsfile     - reference data file for ARGO
1299  c     driftfile     - reference data file for drifter's mean velocities  c     driftfile     - reference data file for drifter mean velocities
1300    
1301        common /ecco_cost_data_c/        common /ecco_cost_data_c/
1302       &                     tdatfile,       &                     tdatfile,
# Line 1037  c     driftfile     - reference data fil Line 1306  c     driftfile     - reference data fil
1306       &                     sstdatfile,       &                     sstdatfile,
1307       &                     tmidatfile,       &                     tmidatfile,
1308       &                     sssdatfile,       &                     sssdatfile,
1309         &                     bpdatfile,
1310       &                     topexmeanfile,       &                     topexmeanfile,
1311       &                     topexfile,       &                     topexfile,
1312       &                     ersfile,       &                     ersfile,
# Line 1048  c     driftfile     - reference data fil Line 1318  c     driftfile     - reference data fil
1318       &                     xbtfile,       &                     xbtfile,
1319       &                     argotfile,       &                     argotfile,
1320       &                     argosfile,       &                     argosfile,
1321       &                     udriftfile,       &                     udriftfile,
1322       &                     vdriftfile,       &                     vdriftfile,
1323         &                     usercost_datafile,
1324       &                     curmtrufile,       &                     curmtrufile,
1325       &                     curmtrvfile       &                     curmtrvfile
1326    
# Line 1060  c     driftfile     - reference data fil Line 1331  c     driftfile     - reference data fil
1331        character*(MAX_LEN_FNAM) sstdatfile        character*(MAX_LEN_FNAM) sstdatfile
1332        character*(MAX_LEN_FNAM) tmidatfile        character*(MAX_LEN_FNAM) tmidatfile
1333        character*(MAX_LEN_FNAM) sssdatfile        character*(MAX_LEN_FNAM) sssdatfile
1334          character*(MAX_LEN_FNAM) bpdatfile
1335        character*(MAX_LEN_FNAM) topexmeanfile        character*(MAX_LEN_FNAM) topexmeanfile
1336        character*(MAX_LEN_FNAM) topexfile        character*(MAX_LEN_FNAM) topexfile
1337        character*(MAX_LEN_FNAM) ersfile        character*(MAX_LEN_FNAM) ersfile
# Line 1072  c     driftfile     - reference data fil Line 1344  c     driftfile     - reference data fil
1344        character*(MAX_LEN_FNAM) argotfile        character*(MAX_LEN_FNAM) argotfile
1345        character*(MAX_LEN_FNAM) argosfile        character*(MAX_LEN_FNAM) argosfile
1346        character*(MAX_LEN_FNAM) argofile        character*(MAX_LEN_FNAM) argofile
1347          character*(MAX_LEN_FNAM) usercost_datafile(NUSERCOST)
1348        character*(MAX_LEN_FNAM) udriftfile        character*(MAX_LEN_FNAM) udriftfile
1349        character*(MAX_LEN_FNAM) vdriftfile              character*(MAX_LEN_FNAM) vdriftfile
1350        character*(MAX_LEN_FNAM) curmtrufile        character*(MAX_LEN_FNAM) curmtrufile
1351        character*(MAX_LEN_FNAM) curmtrvfile        character*(MAX_LEN_FNAM) curmtrvfile
1352    
# Line 1108  c     sshperiod      - sampling interval Line 1381  c     sshperiod      - sampling interval
1381       &                           argosstartdate,       &                           argosstartdate,
1382       &                           tmistartdate,       &                           tmistartdate,
1383       &                           sssstartdate,       &                           sssstartdate,
1384         &                           bpstartdate,
1385       &                           topexstartdate,       &                           topexstartdate,
1386       &                           ersstartdate,       &                           ersstartdate,
1387       &                           gfostartdate       &                           gfostartdate
# Line 1118  c     sshperiod      - sampling interval Line 1392  c     sshperiod      - sampling interval
1392        integer argosstartdate(4)        integer argosstartdate(4)
1393        integer tmistartdate(4)        integer tmistartdate(4)
1394        integer sssstartdate(4)        integer sssstartdate(4)
1395          integer bpstartdate(4)
1396        integer topexstartdate(4)        integer topexstartdate(4)
1397        integer ersstartdate(4)        integer ersstartdate(4)
1398        integer gfostartdate(4)        integer gfostartdate(4)
# Line 1129  c     sshperiod      - sampling interval Line 1404  c     sshperiod      - sampling interval
1404       &                           sststartdate2,       &                           sststartdate2,
1405       &                           sssstartdate1,       &                           sssstartdate1,
1406       &                           sssstartdate2,       &                           sssstartdate2,
1407         &                           bpstartdate1,
1408         &                           bpstartdate2,
1409       &                           argotstartdate1,       &                           argotstartdate1,
1410       &                           argotstartdate2,       &                           argotstartdate2,
1411       &                           argosstartdate1,       &                           argosstartdate1,
# Line 1148  c     sshperiod      - sampling interval Line 1425  c     sshperiod      - sampling interval
1425        integer sststartdate2        integer sststartdate2
1426        integer sssstartdate1        integer sssstartdate1
1427        integer sssstartdate2        integer sssstartdate2
1428          integer bpstartdate1
1429          integer bpstartdate2
1430        integer argotstartdate1        integer argotstartdate1
1431        integer argotstartdate2        integer argotstartdate2
1432        integer argosstartdate1        integer argosstartdate1
# Line 1185  c     sshperiod      - sampling interval Line 1464  c     sshperiod      - sampling interval
1464        _RL ersslope        _RL ersslope
1465        _RL gfoslope        _RL gfoslope
1466    
1467    cgf factor to convert sshv4cost_errfile in m
1468          common /ecco_cost_errfactor/
1469         &         sshv4cost_errfactor
1470          _RL  sshv4cost_errfactor(NSSHV4COST)
1471    
1472    #ifdef ALLOW_SSH_COST_CONTRIBUTION
1473          common /ecco_ssh_daymask_r/
1474         &       tpTimeMask, ersTimeMask, gfoTimeMask
1475          _RL tpTimeMask(maxNumDays)
1476          _RL ersTimeMask(maxNumDays)
1477          _RL gfoTimeMask(maxNumDays)
1478    #endif
1479    
1480          common /ecco_ssh_daymask_c/
1481         &       tpTimeMaskFile, ersTimeMaskFile, gfoTimeMaskFile
1482          character*(MAX_LEN_FNAM) tpTimeMaskFile
1483          character*(MAX_LEN_FNAM) ersTimeMaskFile
1484          character*(MAX_LEN_FNAM) gfoTimeMaskFile
1485    
1486  c     ==================================================================  c     ==================================================================
1487  c     END OF HEADER COST  c     END OF HEADER COST
1488  c     ==================================================================  c     ==================================================================

Legend:
Removed from v.1.17  
changed lines
  Added in v.1.48

  ViewVC Help
Powered by ViewVC 1.1.22