/[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.8 by heimbach, Fri Jul 29 05:19:40 2005 UTC revision 1.38 by heimbach, Wed May 13 13:39:38 2009 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  c     ==================================================================  c     ==================================================================
5  c     HEADER AVERAGES  c     HEADER AVERAGES
# Line 8  c       fields and counters associated w Line 10  c       fields and counters associated w
10  c  c
11  c     started: Christian Eckert eckert@mit.edu  30-Jun-1999  c     started: Christian Eckert eckert@mit.edu  30-Jun-1999
12  c  c
 c     added sea-ice term: menemenlis@jpl.nasa.gov 26-Feb-2003  
 c  
13  c     ==================================================================  c     ==================================================================
14  c     HEADER AVERAGES  c     HEADER AVERAGES
15  c     ==================================================================  c     ==================================================================
# Line 22  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 Generic Cost terms:
36    c     =============================
37          INTEGER NGENCOST
38          PARAMETER ( NGENCOST=10 )
39    
40    c     Number of days: (hard-coded to set up some vector dimensions
41    c     =============================
42    c     17 years: 6210
43          INTEGER maxNumDays
44          PARAMETER ( maxNumDays = 6210 )
45    
46  c     Averaged Fields:  c     Averaged Fields:
47  c     ================  c     ================
# Line 43  c             intantaneous salinities. Line 55  c             intantaneous salinities.
55  c     psbar - contains the averaged surface pressure field after the call  c     psbar - contains the averaged surface pressure field after the call
56  c             to subroutine POST_DAILY. Before, it accumulates the  c             to subroutine POST_DAILY. Before, it accumulates the
57  c             intantaneous surface pressure field.  c             intantaneous surface pressure field.
58  c     ubar  - contains the averaged zonal velocity component for the  c     ubar  - contains the averaged zonal velocity component for the
59  c             whole integration period. Before, it accumulates the  c             whole integration period. Before, it accumulates the
60  c             intantaneous field.  c             intantaneous field.
61  c     vbar  - contains the averaged zonal velocity component for the  c     vbar  - contains the averaged zonal velocity component for the
# Line 62  c     sfluxbar  - contains the averaged Line 74  c     sfluxbar  - contains the averaged
74  c             whole integration period. Before, it accumulates the  c             whole integration period. Before, it accumulates the
75  c             intantaneous field.  c             intantaneous field.
76    
77        common /averages_r/        common /averages_r/
78       &                    tbar,       &                    tbar,
79       &                    sbar,       &                    sbar,
80         &                    sstbar,
81       &                    psbar,       &                    psbar,
82         &                    bpbar,
83       &                    ubar,       &                    ubar,
84       &                    vbar,       &                    vbar,
85       &                    wbar,       &                    wbar,
# Line 73  c             intantaneous field. Line 87  c             intantaneous field.
87       &                    tauybar,       &                    tauybar,
88       &                    hfluxbar,       &                    hfluxbar,
89       &                    sfluxbar,       &                    sfluxbar,
90       &                    Slmean,       &                    Slmean,
91       &                    Tlmean,       &                    Tlmean,
92       &                    wlmean,       &                    wlmean,
93       &                    Sfmean,       &                    Sfmean,
94       &                    Tfmean,                         &                    Tfmean,
95       &                    sbar_gen,       &                    sbar_gen,
96       &                    tbar_gen,       &                    tbar_gen,
97       &                    wfmean       &                    wfmean
# Line 103  c             intantaneous field. Line 117  c             intantaneous field.
117        _RL sbar_gen        _RL sbar_gen
118  #endif  #endif
119    
120    #ifdef ALLOW_DAILYSST_COST_CONTRIBUTION
121    cph#ifdef ALLOW_SEAICE_COST_AREASST
122          _RL sstbar  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
123    #else
124          _RL sstbar
125    #endif
126    
127  #if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \  #if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \
128       defined (ALLOW_CTDS_COST_CONTRIBUTION) || \       defined (ALLOW_CTDS_COST_CONTRIBUTION) || \
129       defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \       defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \
# Line 122  c             intantaneous field. Line 143  c             intantaneous field.
143        _RL psbar        _RL psbar
144  #endif  #endif
145    
146    #ifdef ALLOW_BP_COST_CONTRIBUTION
147          _RL bpbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
148    #else
149          _RL bpbar
150    #endif
151    
152  #if (defined (ALLOW_DRIFTER_COST_CONTRIBUTION) || \  #if (defined (ALLOW_DRIFTER_COST_CONTRIBUTION) || \
153       defined (ALLOW_OBCS_COST_CONTRIBUTION))       defined (ALLOW_OBCS_COST_CONTRIBUTION))
154        _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 146  c             intantaneous field. Line 173  c             intantaneous field.
173        _RL   Tlmean        _RL   Tlmean
174        _RL   Slmean        _RL   Slmean
175        _RL   Tfmean        _RL   Tfmean
176        _RL   Sfmean                _RL   Sfmean
177  #endif  #endif
178    
179  #ifdef ALLOW_DRIFTW_COST_CONTRIBUTION  #ifdef ALLOW_DRIFTW_COST_CONTRIBUTION
# Line 161  c             intantaneous field. Line 188  c             intantaneous field.
188        _RL tauxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)        _RL tauxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
189        _RL tauybar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)        _RL tauybar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
190  #else  #else
191        _RL tauxbar        _RL tauxbar
192        _RL tauybar        _RL tauybar
193  #endif  #endif
194    
195  #ifdef ALLOW_MEAN_HFLUX_COST_CONTRIBUTION  #ifdef ALLOW_MEAN_HFLUX_COST_CONTRIBUTION
196        _RL hfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)        _RL hfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
197  #else  #else
198        _RL hfluxbar        _RL hfluxbar
199  #endif  #endif
200    
201  #ifdef ALLOW_MEAN_SFLUX_COST_CONTRIBUTION  #ifdef ALLOW_MEAN_SFLUX_COST_CONTRIBUTION
202        _RL sfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)            _RL sfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
203  #else  #else
204        _RL sfluxbar        _RL sfluxbar
205  #endif  #endif
206    
207    
208        common /averages_c/        common /averages_c/
209       &                    tbarfile,       &                    tbarfile,
210       &                    sbarfile,       &                    sbarfile,
211         &                    sstbarfile,
212       &                    psbarfile,       &                    psbarfile,
213         &                    bpbarfile,
214       &                    ubarfile,       &                    ubarfile,
215       &                    vbarfile,       &                    vbarfile,
216       &                    wbarfile,       &                    wbarfile,
217       &                    tauxbarfile,       &                    tauxbarfile,
218       &                    tauybarfile,       &                    tauybarfile,
219       &                    hfluxbarfile,       &                    hfluxbarfile,
220       &                    sfluxbarfile       &                    sfluxbarfile,
221         &                    costTranspDataFile
222        character*(MAX_LEN_FNAM) tbarfile        character*(MAX_LEN_FNAM) tbarfile
223        character*(MAX_LEN_FNAM) sbarfile        character*(MAX_LEN_FNAM) sbarfile
224          character*(MAX_LEN_FNAM) sstbarfile
225        character*(MAX_LEN_FNAM) psbarfile        character*(MAX_LEN_FNAM) psbarfile
226          character*(MAX_LEN_FNAM) bpbarfile
227        character*(MAX_LEN_FNAM) ubarfile        character*(MAX_LEN_FNAM) ubarfile
228        character*(MAX_LEN_FNAM) vbarfile        character*(MAX_LEN_FNAM) vbarfile
229        character*(MAX_LEN_FNAM) wbarfile        character*(MAX_LEN_FNAM) wbarfile
# Line 199  c             intantaneous field. Line 231  c             intantaneous field.
231        character*(MAX_LEN_FNAM) tauybarfile        character*(MAX_LEN_FNAM) tauybarfile
232        character*(MAX_LEN_FNAM) hfluxbarfile        character*(MAX_LEN_FNAM) hfluxbarfile
233        character*(MAX_LEN_FNAM) sfluxbarfile        character*(MAX_LEN_FNAM) sfluxbarfile
234          character*(MAX_LEN_FNAM) costTranspDataFile
235    
236    #ifdef ALLOW_TRANSPORT_COST_CONTRIBUTION
237          common /averages_transp_r/
238         &                     transpbar
239         &                   , transpobs
240         &                   , wtransp
241          _RL transpbar(maxNumDays,nsx,nsy)
242          _RL transpobs(maxNumDays)
243          _RL wtransp(maxNumDays)
244    #endif
245    
246  c     file precision and field type  c     file precision and field type
247    
248        common /prec_type_cost/        common /prec_type_cost/
249       &                        cost_iprec,       &                        cost_iprec,
250       &                        cost_yftype       &                        cost_yftype
251    
# Line 256  c     objf_hfluxmm    - Global time-mean Line 299  c     objf_hfluxmm    - Global time-mean
299  c     objf_sfluxmm    - Global time-mean Salt flux contribution.  c     objf_sfluxmm    - Global time-mean Salt flux contribution.
300  c     objf_hmean - Mean sea surface height contribution.  c     objf_hmean - Mean sea surface height contribution.
301  c     objf_h     - Residual sea surface height contribution.  c     objf_h     - Residual sea surface height contribution.
302    c     objf_tp    - Residual sea surface height contribution from T/P
303    c     objf_ers   - Residual sea surface height contribution from T/P
304    c     objf_gfo   - Residual sea surface height contribution from T/P
305  c     objf_temp  - Temperature contribution.  c     objf_temp  - Temperature contribution.
306  c     objf_salt  - Salinity contribution.  c     objf_salt  - Salinity contribution.
307  c     objf_temp0 - Initial conditions Temperature contribution.  c     objf_temp0 - Initial conditions Temperature contribution.
# Line 263  c     objf_salt0 - Initial conditions Sa Line 309  c     objf_salt0 - Initial conditions Sa
309  c     objf_sst   - Sea surface temperature contribution.  c     objf_sst   - Sea surface temperature contribution.
310  c     objf_tmi   - Sea surface temperature contribution.  c     objf_tmi   - Sea surface temperature contribution.
311  c     objf_sss   - Sea surface salinity contribution.  c     objf_sss   - Sea surface salinity contribution.
312  c     objf_ctdt  - Temperature measurements from Woce CTD  c     objf_ctdt  - Temperature measurements from Woce CTD
313  c     objf_ctds  - Salinity measurements from Woce CTD  c     objf_ctds  - Salinity measurements from Woce CTD
314  c     objf_ctdtclim - Temperature measurements from Woce CTD without timetag  c     objf_ctdtclim - Temperature measurements from Woce CTD without timetag
315  c     objf_ctdsclim - Salinity measurements from Woce CTD without timetag  c     objf_ctdsclim - Salinity measurements from Woce CTD without timetag
316  c     objf_xbt   - XBT temperature data  c     objf_xbt   - XBT temperature data
317  c     objf_argot - ARGO temperature profiles  c     objf_argot - ARGO temperature profiles
318  c     objf_argos - ARGO salt profiles  c     objf_argos - ARGO salt profiles
319    c     objf_gencost - gencost user defined contribution
320  c     objf_scatxm - time-mean zonal SCAT  contribution  c     objf_scatxm - time-mean zonal SCAT  contribution
321  c     objf_scatym - time-mean meridional SCAT  contribution  c     objf_scatym - time-mean meridional SCAT  contribution
322  c     objf_scatx  - zonal SCAT  contribution  c     objf_scatx  - zonal SCAT  contribution
323  c     objf_scaty  - meridional SCAT  contribution  c     objf_scaty  - meridional SCAT  contribution
 c     objf_ice    - sea-ice volume  
324  c     objf_kapgm  - kappa GM contribution  c     objf_kapgm  - kappa GM contribution
325    c     objf_kapredi  - kappa REDI contribution
326  c     objf_diffkr - diffusion contribution  c     objf_diffkr - diffusion contribution
327  c     objf_theta_ini_fin - final vs. initial theta misfit  c     objf_theta_ini_fin - final vs. initial theta misfit
328  c     objf_salt_ini_fin  - final vs. initial salt misfit  c     objf_salt_ini_fin  - final vs. initial salt misfit
329    c     objf_eddytau - eddy streamfunction contribution
330    c     objf_bottomdrag - bottom drag contribution
331  c  c
332  c     mult_"var" - multipliers for the individual cost  c     mult_"var" - multipliers for the individual cost
333  c                  function contributions.  c                  function contributions.
# Line 287  c                  function contribution Line 336  c                  function contribution
336       &                objf_hflux,       &                objf_hflux,
337       &                objf_hfluxm,       &                objf_hfluxm,
338       &                objf_hfluxmm,       &                objf_hfluxmm,
339         &                objf_hfluxsmoo,
340       &                objf_sflux,       &                objf_sflux,
341       &                objf_sfluxm,       &                objf_sfluxm,
342       &                objf_sfluxmm,       &                objf_sfluxmm,
343         &                objf_sfluxsmoo,
344       &                objf_tauu,       &                objf_tauu,
345       &                objf_tauum,       &                objf_tauum,
346         &                objf_tauusmoo,
347       &                objf_tauv,       &                objf_tauv,
348       &                objf_tauvm,       &                objf_tauvm,
349         &                objf_tauvsmoo,
350       &                objf_hmean,       &                objf_hmean,
351       &                objf_h,       &                objf_h,
352         &                objf_tp,
353         &                objf_ers,
354         &                objf_gfo,
355       &                objf_temp,       &                objf_temp,
356       &                objf_salt,       &                objf_salt,
357       &                objf_temp0,       &                objf_temp0,
358       &                objf_salt0,       &                objf_salt0,
359         &                objf_temp0smoo,
360         &                objf_salt0smoo,
361       &                objf_sst,       &                objf_sst,
362       &                objf_tmi,       &                objf_tmi,
363       &                objf_sss,       &                objf_sss,
364         &                objf_bp,
365       &                objf_ctdt,       &                objf_ctdt,
366       &                objf_ctds,       &                objf_ctds,
367       &                objf_ctdtclim,       &                objf_ctdtclim,
# Line 310  c                  function contribution Line 369  c                  function contribution
369       &                objf_xbt,       &                objf_xbt,
370       &                objf_argot,       &                objf_argot,
371       &                objf_argos,       &                objf_argos,
372         &                objf_gencost,
373       &                objf_drift,       &                objf_drift,
374       &                objf_tdrift,       &                objf_tdrift,
375       &                objf_sdrift,       &                objf_sdrift,
# Line 322  c                  function contribution Line 382  c                  function contribution
382       &                objf_aqh,       &                objf_aqh,
383       &                objf_precip,       &                objf_precip,
384       &                objf_swflux,       &                objf_swflux,
385         &                objf_swdown,
386         &                objf_snowprecip,
387         &                objf_lwflux,
388         &                objf_lwdown,
389         &                objf_evap,
390         &                objf_apressure,
391         &                objf_runoff,
392       &                objf_uwind,       &                objf_uwind,
393       &                objf_vwind,       &                objf_vwind,
394       &                objf_atempm,       &                objf_atempm,
395       &                objf_aqhm,       &                objf_aqhm,
396       &                objf_precipm,       &                objf_precipm,
397       &                objf_swfluxm,       &                objf_swfluxm,
398         &                objf_swdownm,
399         &                objf_snowprecipm,
400         &                objf_lwfluxm,
401         &                objf_lwdownm,
402         &                objf_evapm,
403         &                objf_apressurem,
404         &                objf_runoffm,
405       &                objf_uwindm,       &                objf_uwindm,
406       &                objf_vwindm,       &                objf_vwindm,
407         &                objf_atempsmoo,
408         &                objf_aqhsmoo,
409         &                objf_precipsmoo,
410         &                objf_swfluxsmoo,
411         &                objf_swdownsmoo,
412         &                objf_snowprecipsmoo,
413         &                objf_lwfluxsmoo,
414         &                objf_lwdownsmoo,
415         &                objf_evapsmoo,
416         &                objf_apressuresmoo,
417         &                objf_runoffsmoo,
418         &                objf_uwindsmoo,
419         &                objf_vwindsmoo,
420       &                objf_obcsn,       &                objf_obcsn,
421       &                objf_obcss,       &                objf_obcss,
422       &                objf_obcsw,       &                objf_obcsw,
# Line 337  c                  function contribution Line 424  c                  function contribution
424       &                objf_obcsvol,       &                objf_obcsvol,
425       &                objf_curmtr,       &                objf_curmtr,
426       &                objf_ageos,       &                objf_ageos,
427       &                objf_ice,       &                objf_kapgm,
428         &                objf_kapredi,
429       &                objf_diffkr,       &                objf_diffkr,
430       &                objf_theta_ini_fin,       &                objf_theta_ini_fin,
431       &                objf_salt_ini_fin,       &                objf_salt_ini_fin,
432       &                objf_eddytau       &                objf_eddytau,
433         &                objf_bottomdrag,
434         &                objf_transp
435        _RL  objf_hflux  (nsx,nsy)        _RL  objf_hflux  (nsx,nsy)
436        _RL  objf_hfluxm (nsx,nsy)        _RL  objf_hfluxm (nsx,nsy)
437        _RL  objf_hfluxmm(nsx,nsy)        _RL  objf_hfluxmm
438          _RL  objf_hfluxsmoo  (nsx,nsy)
439        _RL  objf_sflux  (nsx,nsy)        _RL  objf_sflux  (nsx,nsy)
440        _RL  objf_sfluxm (nsx,nsy)        _RL  objf_sfluxm (nsx,nsy)
441        _RL  objf_sfluxmm(nsx,nsy)        _RL  objf_sfluxmm
442          _RL  objf_sfluxsmoo  (nsx,nsy)
443        _RL  objf_tauu   (nsx,nsy)        _RL  objf_tauu   (nsx,nsy)
444        _RL  objf_tauum  (nsx,nsy)        _RL  objf_tauum  (nsx,nsy)
445          _RL  objf_tauusmoo   (nsx,nsy)
446        _RL  objf_tauv   (nsx,nsy)        _RL  objf_tauv   (nsx,nsy)
447        _RL  objf_tauvm  (nsx,nsy)        _RL  objf_tauvm  (nsx,nsy)
448          _RL  objf_tauvsmoo   (nsx,nsy)
449        _RL  objf_hmean        _RL  objf_hmean
450        _RL  objf_h    (nsx,nsy)        _RL  objf_h    (nsx,nsy)
451          _RL  objf_tp   (nsx,nsy)
452          _RL  objf_ers  (nsx,nsy)
453          _RL  objf_gfo  (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)
457        _RL  objf_salt0(nsx,nsy)        _RL  objf_salt0(nsx,nsy)
458          _RL  objf_temp0smoo(nsx,nsy)
459          _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)
# Line 369  c                  function contribution Line 468  c                  function contribution
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_argos(nsx,nsy)        _RL  objf_argos(nsx,nsy)
471          _RL  objf_gencost(NGENCOST,nsx,nsy)
472        _RL  objf_drift(nsx,nsy)        _RL  objf_drift(nsx,nsy)
473        _RL  objf_tdrift(nsx,nsy)        _RL  objf_tdrift(nsx,nsy)
474        _RL  objf_sdrift(nsx,nsy)        _RL  objf_sdrift(nsx,nsy)
# Line 381  c                  function contribution Line 481  c                  function contribution
481        _RL  objf_aqh  (nsx,nsy)        _RL  objf_aqh  (nsx,nsy)
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)
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)
494        _RL  objf_aqhm  (nsx,nsy)        _RL  objf_aqhm  (nsx,nsy)
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)
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)
507          _RL  objf_aqhsmoo  (nsx,nsy)
508          _RL  objf_precipsmoo(nsx,nsy)
509          _RL  objf_swfluxsmoo(nsx,nsy)
510          _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)
518          _RL  objf_vwindsmoo(nsx,nsy)
519        _RL  objf_obcsn(nsx,nsy)        _RL  objf_obcsn(nsx,nsy)
520        _RL  objf_obcss(nsx,nsy)        _RL  objf_obcss(nsx,nsy)
521        _RL  objf_obcsw(nsx,nsy)        _RL  objf_obcsw(nsx,nsy)
# Line 396  c                  function contribution Line 523  c                  function contribution
523        _RL  objf_obcsvol        _RL  objf_obcsvol
524        _RL  objf_curmtr(nsx,nsy)        _RL  objf_curmtr(nsx,nsy)
525        _RL  objf_ageos(nsx,nsy)        _RL  objf_ageos(nsx,nsy)
       _RL  objf_ice  (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 416  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_temp,       &                num_temp,
552       &                num_salt,       &                num_salt,
553       &                num_temp0,       &                num_temp0,
# Line 423  c                  function contribution Line 555  c                  function contribution
555       &                num_sst,       &                num_sst,
556       &                num_tmi,       &                num_tmi,
557       &                num_sss,       &                num_sss,
558         &                num_bp,
559       &                num_ctdt,       &                num_ctdt,
560       &                num_ctds,       &                num_ctds,
561       &                num_ctdtclim,       &                num_ctdtclim,
# Line 430  c                  function contribution Line 563  c                  function contribution
563       &                num_xbt,       &                num_xbt,
564       &                num_argot,       &                num_argot,
565       &                num_argos,       &                num_argos,
566         &                num_gencost,
567       &                num_drift,       &                num_drift,
568       &                num_tdrift,       &                num_tdrift,
569       &                num_sdrift,       &                num_sdrift,
# Line 442  c                  function contribution Line 576  c                  function contribution
576       &                num_aqh,       &                num_aqh,
577       &                num_precip,       &                num_precip,
578       &                num_swflux,       &                num_swflux,
579         &                num_swdown,
580         &                num_snowprecip,
581         &                num_lwflux,
582         &                num_lwdown,
583         &                num_evap,
584         &                num_apressure,
585         &                num_runoff,
586       &                num_uwind,       &                num_uwind,
587       &                num_vwind,       &                num_vwind,
588       &                num_atempm,       &                num_atempm,
589       &                num_aqhm,       &                num_aqhm,
590       &                num_precipm,       &                num_precipm,
591       &                num_swfluxm,       &                num_swfluxm,
592         &                num_swdownm,
593         &                num_snowprecipm,
594         &                num_lwfluxm,
595         &                num_lwdownm,
596         &                num_evapm,
597         &                num_apressurem,
598         &                num_runoffm,
599       &                num_uwindm,       &                num_uwindm,
600       &                num_vwindm,       &                num_vwindm,
601       &                num_obcsn,       &                num_obcsn,
# Line 457  c                  function contribution Line 605  c                  function contribution
605       &                num_obcsvol,       &                num_obcsvol,
606       &                num_curmtr,       &                num_curmtr,
607       &                num_ageos,       &                num_ageos,
608       &                num_ice,       &                num_kapgm,
609         &                num_kapredi,
610       &                num_diffkr,       &                num_diffkr,
611       &                num_theta_ini_fin,       &                num_theta_ini_fin,
612       &                num_salt_ini_fin,       &                num_salt_ini_fin,
613       &                num_eddytau       &                num_eddytau,
614         &                num_bottomdrag,
615         &                num_transp
616    
617        _RL  num_hflux  (nsx,nsy)        _RL  num_hflux  (nsx,nsy)
618        _RL  num_hfluxm (nsx,nsy)        _RL  num_hfluxm (nsx,nsy)
# Line 475  c                  function contribution Line 626  c                  function contribution
626        _RL  num_tauvm  (nsx,nsy)        _RL  num_tauvm  (nsx,nsy)
627        _RL  num_hmean        _RL  num_hmean
628        _RL  num_h    (nsx,nsy)        _RL  num_h    (nsx,nsy)
629          _RL  num_tp   (nsx,nsy)
630          _RL  num_ers  (nsx,nsy)
631          _RL  num_gfo  (nsx,nsy)
632        _RL  num_temp (nsx,nsy)        _RL  num_temp (nsx,nsy)
633        _RL  num_salt (nsx,nsy)        _RL  num_salt (nsx,nsy)
634        _RL  num_temp0(nsx,nsy)        _RL  num_temp0(nsx,nsy)
635        _RL  num_salt0(nsx,nsy)        _RL  num_salt0(nsx,nsy)
636        _RL  num_sst  (nsx,nsy)        _RL  num_sst  (nsx,nsy)
637        _RL  num_tmi  (nsx,nsy)        _RL  num_tmi  (nsx,nsy)
638        _RL  num_sss  (nsx,nsy)        _RL  num_sss  (nsx,nsy)
639          _RL  num_bp   (nsx,nsy)
640        _RL  num_ctdt (nsx,nsy)        _RL  num_ctdt (nsx,nsy)
641        _RL  num_ctds (nsx,nsy)        _RL  num_ctds (nsx,nsy)
642        _RL  num_ctdtclim (nsx,nsy)        _RL  num_ctdtclim (nsx,nsy)
# Line 489  c                  function contribution Line 644  c                  function contribution
644        _RL  num_xbt  (nsx,nsy)        _RL  num_xbt  (nsx,nsy)
645        _RL  num_argot(nsx,nsy)        _RL  num_argot(nsx,nsy)
646        _RL  num_argos(nsx,nsy)        _RL  num_argos(nsx,nsy)
647          _RL  num_gencost(NGENCOST,nsx,nsy)
648        _RL  num_drift(nsx,nsy)        _RL  num_drift(nsx,nsy)
649        _RL  num_tdrift(nsx,nsy)        _RL  num_tdrift(nsx,nsy)
650        _RL  num_sdrift(nsx,nsy)        _RL  num_sdrift(nsx,nsy)
# Line 501  c                  function contribution Line 657  c                  function contribution
657        _RL  num_aqh  (nsx,nsy)        _RL  num_aqh  (nsx,nsy)
658        _RL  num_precip(nsx,nsy)        _RL  num_precip(nsx,nsy)
659        _RL  num_swflux(nsx,nsy)        _RL  num_swflux(nsx,nsy)
660          _RL  num_swdown(nsx,nsy)
661          _RL  num_snowprecip(nsx,nsy)
662          _RL  num_lwflux(nsx,nsy)
663          _RL  num_lwdown(nsx,nsy)
664          _RL  num_evap(nsx,nsy)
665          _RL  num_apressure(nsx,nsy)
666          _RL  num_runoff(nsx,nsy)
667        _RL  num_uwind(nsx,nsy)        _RL  num_uwind(nsx,nsy)
668        _RL  num_vwind(nsx,nsy)        _RL  num_vwind(nsx,nsy)
669        _RL  num_atempm(nsx,nsy)        _RL  num_atempm(nsx,nsy)
670        _RL  num_aqhm  (nsx,nsy)        _RL  num_aqhm  (nsx,nsy)
671        _RL  num_precipm(nsx,nsy)        _RL  num_precipm(nsx,nsy)
672        _RL  num_swfluxm(nsx,nsy)        _RL  num_swfluxm(nsx,nsy)
673          _RL  num_swdownm(nsx,nsy)
674          _RL  num_snowprecipm(nsx,nsy)
675          _RL  num_lwfluxm(nsx,nsy)
676          _RL  num_lwdownm(nsx,nsy)
677          _RL  num_evapm(nsx,nsy)
678          _RL  num_apressurem(nsx,nsy)
679          _RL  num_runoffm(nsx,nsy)
680        _RL  num_uwindm(nsx,nsy)        _RL  num_uwindm(nsx,nsy)
681        _RL  num_vwindm(nsx,nsy)        _RL  num_vwindm(nsx,nsy)
682        _RL  num_obcsn(nsx,nsy)        _RL  num_obcsn(nsx,nsy)
# Line 516  c                  function contribution Line 686  c                  function contribution
686        _RL  num_obcsvol        _RL  num_obcsvol
687        _RL  num_curmtr(nsx,nsy)        _RL  num_curmtr(nsx,nsy)
688        _RL  num_ageos(nsx,nsy)        _RL  num_ageos(nsx,nsy)
       _RL  num_ice  (nsx,nsy)  
689        _RL  num_kapgm(nsx,nsy)        _RL  num_kapgm(nsx,nsy)
690          _RL  num_kapredi(nsx,nsy)
691        _RL  num_diffkr(nsx,nsy)        _RL  num_diffkr(nsx,nsy)
692        _RL  num_theta_ini_fin(nsx,nsy)        _RL  num_theta_ini_fin(nsx,nsy)
693        _RL  num_salt_ini_fin(nsx,nsy)        _RL  num_salt_ini_fin(nsx,nsy)
694        _RL  num_eddytau(nsx,nsy)        _RL  num_eddytau(nsx,nsy)
695          _RL  num_bottomdrag(nsx,nsy)
696          _RL  num_transp
697    
698        common /ecco_cost_aux_r/        common /ecco_cost_aux_r/
699       &                    mult_hflux,       &                    mult_hflux,
700       &                    mult_sflux,       &                    mult_sflux,
701         &                    mult_hfluxmm,
702         &                    mult_sfluxmm,
703       &                    mult_tauu,       &                    mult_tauu,
704       &                    mult_tauv,       &                    mult_tauv,
705       &                    mult_hmean,       &                    mult_hmean,
706       &                    mult_h,       &                    mult_h,
707         &                    mult_tp,
708         &                    mult_ers,
709         &                    mult_gfo,
710       &                    mult_temp,       &                    mult_temp,
711       &                    mult_salt,       &                    mult_salt,
712       &                    mult_temp0,       &                    mult_temp0,
# Line 537  c                  function contribution Line 714  c                  function contribution
714       &                    mult_sst,       &                    mult_sst,
715       &                    mult_tmi,       &                    mult_tmi,
716       &                    mult_sss,       &                    mult_sss,
717         &                    mult_bp,
718       &                    mult_ctdt,       &                    mult_ctdt,
719       &                    mult_ctds,       &                    mult_ctds,
720       &                    mult_ctdtclim,       &                    mult_ctdtclim,
# Line 544  c                  function contribution Line 722  c                  function contribution
722       &                    mult_xbt,       &                    mult_xbt,
723       &                    mult_argot,       &                    mult_argot,
724       &                    mult_argos,       &                    mult_argos,
725         &                    mult_gencost,
726       &                    mult_drift,       &                    mult_drift,
727       &                    mult_tdrift,       &                    mult_tdrift,
728       &                    mult_sdrift,       &                    mult_sdrift,
# Line 554  c                  function contribution Line 733  c                  function contribution
733       &                    mult_aqh,       &                    mult_aqh,
734       &                    mult_precip,       &                    mult_precip,
735       &                    mult_swflux,       &                    mult_swflux,
736         &                    mult_swdown,
737         &                    mult_snowprecip,
738         &                    mult_lwflux,
739         &                    mult_lwdown,
740         &                    mult_evap,
741         &                    mult_apressure,
742         &                    mult_runoff,
743       &                    mult_uwind,       &                    mult_uwind,
744       &                    mult_vwind,       &                    mult_vwind,
745       &                    mult_obcsn,       &                    mult_obcsn,
# Line 563  c                  function contribution Line 749  c                  function contribution
749       &                    mult_obcsvol,       &                    mult_obcsvol,
750       &                    mult_curmtr,       &                    mult_curmtr,
751       &                    mult_ageos,       &                    mult_ageos,
      &                    mult_ice,  
752       &                    mult_kapgm,       &                    mult_kapgm,
753         &                    mult_kapredi,
754       &                    mult_diffkr,       &                    mult_diffkr,
755       &                    mult_ini_fin,       &                    mult_ini_fin,
756       &                    mult_eddytau       &                    mult_edtau,
757         &                    mult_bottomdrag,
758         &                    mult_smooth_ic,
759         &                    mult_smooth_bc,
760         &                    mult_transp
761        _RL  mult_hflux        _RL  mult_hflux
762        _RL  mult_sflux        _RL  mult_sflux
763          _RL  mult_hfluxmm
764          _RL  mult_sfluxmm
765        _RL  mult_tauu        _RL  mult_tauu
766        _RL  mult_tauv        _RL  mult_tauv
767        _RL  mult_hmean        _RL  mult_hmean
768        _RL  mult_h        _RL  mult_h
769          _RL  mult_tp
770          _RL  mult_ers
771          _RL  mult_gfo
772        _RL  mult_temp        _RL  mult_temp
773        _RL  mult_salt        _RL  mult_salt
774        _RL  mult_temp0        _RL  mult_temp0
# Line 582  c                  function contribution Line 776  c                  function contribution
776        _RL  mult_sst        _RL  mult_sst
777        _RL  mult_tmi        _RL  mult_tmi
778        _RL  mult_sss        _RL  mult_sss
779          _RL  mult_bp
780        _RL  mult_ctdt        _RL  mult_ctdt
781        _RL  mult_ctds        _RL  mult_ctds
782        _RL  mult_ctdtclim        _RL  mult_ctdtclim
# Line 589  c                  function contribution Line 784  c                  function contribution
784        _RL  mult_xbt        _RL  mult_xbt
785        _RL  mult_argot        _RL  mult_argot
786        _RL  mult_argos        _RL  mult_argos
787          _RL  mult_gencost(NGENCOST)
788        _RL  mult_drift        _RL  mult_drift
789        _RL  mult_tdrift        _RL  mult_tdrift
790        _RL  mult_sdrift        _RL  mult_sdrift
# Line 599  c                  function contribution Line 795  c                  function contribution
795        _RL  mult_aqh        _RL  mult_aqh
796        _RL  mult_precip        _RL  mult_precip
797        _RL  mult_swflux        _RL  mult_swflux
798          _RL  mult_swdown
799          _RL  mult_snowprecip
800          _RL  mult_lwflux
801          _RL  mult_lwdown
802          _RL  mult_evap
803          _RL  mult_apressure
804          _RL  mult_runoff
805        _RL  mult_uwind        _RL  mult_uwind
806        _RL  mult_vwind        _RL  mult_vwind
807        _RL  mult_obcsn        _RL  mult_obcsn
# Line 608  c                  function contribution Line 811  c                  function contribution
811        _RL  mult_obcsvol        _RL  mult_obcsvol
812        _RL  mult_curmtr        _RL  mult_curmtr
813        _RL  mult_ageos        _RL  mult_ageos
       _RL  mult_ice  
814        _RL  mult_kapgm        _RL  mult_kapgm
815          _RL  mult_kapredi
816        _RL  mult_diffkr        _RL  mult_diffkr
817        _RL  mult_ini_fin        _RL  mult_ini_fin
818        _RL  mult_eddytau        _RL  mult_edtau
819          _RL  mult_bottomdrag
820          _RL  mult_smooth_ic
821          _RL  mult_smooth_bc
822          _RL  mult_transp
823    
824  c     Record counters relevant for the cost function evaluation.  c     Record counters relevant for the cost function evaluation.
825  c     ==========================================================  c     ==========================================================
826  c  c
827  c     nmonsrec - number of monthly records that will be generated by  c     nyearsrec - number of yearly records that will be generated by
828  c                the current model integration.  c                 the current model integration.
829  c     ndaysrec - number of  daily  records that will be generated by  c     nmonsrec  - number of monthly records that will be generated by
830  c                the current model integration.  c                 the current model integration.
831    c     ndaysrec  - number of  daily  records that will be generated by
832    c                 the current model integration.
833    
834        common /ecco_cost_i/        common /ecco_cost_i/
835         &                nyearsrec,
836       &                nmonsrec,       &                nmonsrec,
837       &                ndaysrec       &                ndaysrec,
838         &                nnztbar,
839         &                nnzsbar
840          integer nyearsrec
841        integer nmonsrec        integer nmonsrec
842        integer ndaysrec        integer ndaysrec
843          integer nnztbar
844          integer nnzsbar
845    
846    
847  c     Data files for the weights used in the cost function:  c     Data files for the weights used in the cost function:
# Line 651  c     salterrfile           - representa Line 866  c     salterrfile           - representa
866  c     temperrfile           - representation error due unresolved eddies  c     temperrfile           - representation error due unresolved eddies
867  c     velerrfile            - representation error  c     velerrfile            - representation error
868    
869        common /ecco_cost_c/        common /ecco_cost_c/
870       &                hflux_errfile,       &                hflux_errfile,
871       &                hfluxm_errfile,       &                hfluxm_errfile,
872       &                sflux_errfile,       &                sflux_errfile,
# Line 666  c     velerrfile            - representa Line 881  c     velerrfile            - representa
881       &                geoid_errfile,       &                geoid_errfile,
882       &                geoid_covariancefile,       &                geoid_covariancefile,
883       &                ssh_errfile,       &                ssh_errfile,
884         &                tp_errfile,
885         &                ers_errfile,
886         &                gfo_errfile,
887       &                ctdt_errfile,       &                ctdt_errfile,
888       &                ctds_errfile,       &                ctds_errfile,
889       &                drift_errfile,       &                drift_errfile,
890       &                udrifterrfile,       &                udrifterrfile,
891       &                vdrifterrfile,       &                vdrifterrfile,
892       &                salterrfile,       &                salterrfile,
893       &                temperrfile,       &                temperrfile,
894       &                velerrfile,       &                velerrfile,
895         &                salt0errfile,
896         &                temp0errfile,
897         &                vel0errfile,
898       &                atemp_errfile,       &                atemp_errfile,
899       &                aqh_errfile,       &                aqh_errfile,
900       &                precip_errfile,       &                precip_errfile,
901       &                swflux_errfile,       &                swflux_errfile,
902         &                swdown_errfile,
903         &                snowprecip_errfile,
904         &                lwflux_errfile,
905         &                lwdown_errfile,
906         &                evap_errfile,
907         &                apressure_errfile,
908         &                runoff_errfile,
909         &                edtau_errfile,
910         &                kapgm_errfile,
911         &                kapredi_errfile,
912         &                diffkr_errfile,
913         &                bottomdrag_errfile,
914         &                gencost_errfile,
915       &                uwind_errfile,       &                uwind_errfile,
916       &                vwind_errfile       &                vwind_errfile
917        character*(MAX_LEN_FNAM) hflux_errfile        character*(MAX_LEN_FNAM) hflux_errfile
# Line 694  c     velerrfile            - representa Line 928  c     velerrfile            - representa
928        character*(MAX_LEN_FNAM) geoid_errfile        character*(MAX_LEN_FNAM) geoid_errfile
929        character*(MAX_LEN_FNAM) geoid_covariancefile        character*(MAX_LEN_FNAM) geoid_covariancefile
930        character*(MAX_LEN_FNAM) ssh_errfile        character*(MAX_LEN_FNAM) ssh_errfile
931        character*(MAX_LEN_FNAM) ctdt_errfile        character*(MAX_LEN_FNAM) tp_errfile
932        character*(MAX_LEN_FNAM) ctds_errfile        character*(MAX_LEN_FNAM) ers_errfile
933          character*(MAX_LEN_FNAM) gfo_errfile
934          character*(MAX_LEN_FNAM) ctdt_errfile
935          character*(MAX_LEN_FNAM) ctds_errfile
936        character*(MAX_LEN_FNAM) drift_errfile        character*(MAX_LEN_FNAM) drift_errfile
937        character*(MAX_LEN_FNAM) udrifterrfile        character*(MAX_LEN_FNAM) udrifterrfile
938        character*(MAX_LEN_FNAM) vdrifterrfile              character*(MAX_LEN_FNAM) vdrifterrfile
939        character*(MAX_LEN_FNAM) salterrfile        character*(MAX_LEN_FNAM) salterrfile
940        character*(MAX_LEN_FNAM) temperrfile        character*(MAX_LEN_FNAM) temperrfile
941        character*(MAX_LEN_FNAM) velerrfile        character*(MAX_LEN_FNAM) velerrfile
942          character*(MAX_LEN_FNAM) salt0errfile
943          character*(MAX_LEN_FNAM) temp0errfile
944          character*(MAX_LEN_FNAM) vel0errfile
945        character*(MAX_LEN_FNAM) atemp_errfile        character*(MAX_LEN_FNAM) atemp_errfile
946        character*(MAX_LEN_FNAM) aqh_errfile        character*(MAX_LEN_FNAM) aqh_errfile
947        character*(MAX_LEN_FNAM) precip_errfile        character*(MAX_LEN_FNAM) precip_errfile
948        character*(MAX_LEN_FNAM) swflux_errfile        character*(MAX_LEN_FNAM) swflux_errfile
949          character*(MAX_LEN_FNAM) swdown_errfile
950          character*(MAX_LEN_FNAM) snowprecip_errfile
951          character*(MAX_LEN_FNAM) lwflux_errfile
952          character*(MAX_LEN_FNAM) lwdown_errfile
953          character*(MAX_LEN_FNAM) evap_errfile
954          character*(MAX_LEN_FNAM) apressure_errfile
955          character*(MAX_LEN_FNAM) runoff_errfile
956          character*(MAX_LEN_FNAM) edtau_errfile
957          character*(MAX_LEN_FNAM) kapgm_errfile
958          character*(MAX_LEN_FNAM) kapredi_errfile
959          character*(MAX_LEN_FNAM) diffkr_errfile
960          character*(MAX_LEN_FNAM) bottomdrag_errfile
961          character*(MAX_LEN_FNAM) gencost_errfile(NGENCOST)
962        character*(MAX_LEN_FNAM) uwind_errfile        character*(MAX_LEN_FNAM) uwind_errfile
963        character*(MAX_LEN_FNAM) vwind_errfile        character*(MAX_LEN_FNAM) vwind_errfile
964    
# Line 734  c     wctds      - weight for CTD salini Line 987  c     wctds      - weight for CTD salini
987  c     wudrift    - weight for mean zonal velocity from drifters.  c     wudrift    - weight for mean zonal velocity from drifters.
988  c     wvdrift    - weight for mean meridional velocity from drifters.  c     wvdrift    - weight for mean meridional velocity from drifters.
989    
990        common /ecco_cost_weights_r/        common /ecco_cost_weights_r/
991       &                      frame,       &                      frame,
992       &                      cosphi,       &                      cosphi,
993       &                      whflux,wsflux,wtauu,wtauv,       &                      whflux,wsflux,wtauu,wtauv,
994       &                      watemp,waqh,wprecip,wswflux,wuwind,wvwind,       &                      watemp,waqh,wprecip,wsnowprecip,
995         &                      wswflux,wswdown,wlwflux,wlwdown,
996         &                      wevap,wapressure,wrunoff,
997         &                      wbottomdrag,
998         &                      wuwind,wvwind,
999       &                      wscatx,wscaty,       &                      wscatx,wscaty,
1000       &                      wtheta,wtheta2,wthetaLev,       &                      wtheta,wtheta2,wthetaLev,
1001       &                      wsalt,wsalt2,wsaltLev,       &                      wsalt,wsalt2,wsaltLev,
1002       &                      wsst,wsss,       &                      wdiffkr,wdiffkr2,wdiffkrFld,
1003       &                      wtp,wers,       &                      wkapgm,wkapgm2,wkapgmFld,
1004         &                      wkapredi,wkapredi2,wkaprediFld,
1005         &                      wedtaux,wedtaux2,wedtauxFld,
1006         &                      wedtauy,wedtauy2,wedtauyFld,
1007         &                      wsst,wsss,wbp,
1008         &                      wtp,wers,wgfo,
1009       &                      wp,       &                      wp,
1010       &                      wctdt,wctds,       &                      wctdt,wctds,
1011       &                      wudrift,wvdrift,       &                      wudrift,wvdrift,
# Line 769  c     wvdrift    - weight for mean merid Line 1031  c     wvdrift    - weight for mean merid
1031        _RL waqh    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL waqh    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1032        _RL wprecip (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wprecip (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1033        _RL wswflux (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wswflux (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1034          _RL wswdown (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1035          _RL wsnowprecip (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1036          _RL wlwflux (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1037          _RL wlwdown (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1038          _RL wevap   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1039          _RL wapressure(1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1040          _RL wrunoff (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1041          _RL wbottomdrag (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
1042        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1043        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1044        _RL wtheta  (                            nr,nsx,nsy)        _RL wtheta  (                            nr,nsx,nsy)
# Line 779  c     wvdrift    - weight for mean merid Line 1049  c     wvdrift    - weight for mean merid
1049        _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)
1050        _RL wsst    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsst    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1051        _RL wsss    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsss    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1052          _RL wbp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1053        _RL wtp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1054        _RL wers    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wers    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1055          _RL wgfo    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1056        _RL wp      (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wp      (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1057        _RL wctdt   (                            nr,nsx,nsy)        _RL wctdt   (                            nr,nsx,nsy)
1058        _RL wctds   (                            nr,nsx,nsy)        _RL wctds   (                            nr,nsx,nsy)
# Line 790  c     wvdrift    - weight for mean merid Line 1062  c     wvdrift    - weight for mean merid
1062        _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)
1063        _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)
1064        _RL wbaro        _RL wbaro
1065          _RL wdiffkr (                            nr,nsx,nsy)
1066          _RL wdiffkr2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1067          _RL wdiffkrFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1068          _RL wkapgm (                            nr,nsx,nsy)
1069          _RL wkapgm2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1070          _RL wkapgmFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1071          _RL wkapredi (                            nr,nsx,nsy)
1072          _RL wkapredi2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1073          _RL wkaprediFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1074          _RL wedtaux (                            nr,nsx,nsy)
1075          _RL wedtaux2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1076          _RL wedtauxFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1077          _RL wedtauy (                            nr,nsx,nsy)
1078          _RL wedtauy2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1079          _RL wedtauyFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1080    
1081        common /ecco_cost_weights_0_r/        common /ecco_cost_weights_0_r/
1082       &        whflux0, wsflux0, wtau0,       &        whflux0, wsflux0, wtau0,
1083       &        watemp0, waqh0, wprecip0, wswflux0, wwind0       &        watemp0, waqh0, wprecip0, wsnowprecip0, wwind0,
1084         &        wswflux0, wswdown0, wlwflux0, wlwdown0,
1085         &        wevap0, wapressure0, wrunoff0, wkapredi0,
1086         &        wbottomdrag0,wdiffkr0, wkapgm0, wedtau0
1087        _RL whflux0        _RL whflux0
1088        _RL wsflux0        _RL wsflux0
1089        _RL wtau0        _RL wtau0
# Line 801  c     wvdrift    - weight for mean merid Line 1091  c     wvdrift    - weight for mean merid
1091        _RL waqh0        _RL waqh0
1092        _RL wprecip0        _RL wprecip0
1093        _RL wswflux0        _RL wswflux0
1094          _RL wswdown0
1095          _RL wsnowprecip0
1096          _RL wlwflux0
1097          _RL wlwdown0
1098          _RL wevap0
1099          _RL wapressure0
1100          _RL wrunoff0
1101          _RL wbottomdrag0
1102        _RL wwind0        _RL wwind0
1103          _RL wdiffkr0
1104          _RL wkapgm0
1105          _RL wkapredi0
1106          _RL wedtau0
1107    
1108        common /ecco_cost_weights_mean_r/        common /ecco_cost_weights_mean_r/
1109       &        wmean_hflux, wmean_sflux, wmean_tau,       &        wmean_hflux, wmean_sflux, wmean_tau,
1110       &        wmean_atemp, wmean_aqh,       &        wmean_atemp, wmean_aqh,
1111       &        wmean_precip, wmean_swflux, wmean_wind       &        wmean_precip, wmean_snowprecip, wmean_wind,
1112         &        wmean_swflux, wmean_swdown, wmean_lwflux, wmean_lwdown,
1113         &        wmean_evap, wmean_apressure, wmean_runoff
1114        _RL wmean_hflux        _RL wmean_hflux
1115        _RL wmean_sflux        _RL wmean_sflux
1116        _RL wmean_tau        _RL wmean_tau
# Line 814  c     wvdrift    - weight for mean merid Line 1118  c     wvdrift    - weight for mean merid
1118        _RL wmean_aqh        _RL wmean_aqh
1119        _RL wmean_precip        _RL wmean_precip
1120        _RL wmean_swflux        _RL wmean_swflux
1121          _RL wmean_swdown
1122          _RL wmean_snowprecip
1123          _RL wmean_lwflux
1124          _RL wmean_lwdown
1125          _RL wmean_evap
1126          _RL wmean_apressure
1127          _RL wmean_runoff
1128        _RL wmean_wind        _RL wmean_wind
1129                                                                                                    
1130        common /ecco_cost_weights_2_r/        common /ecco_cost_weights_2_r/
1131       &                      whflux2,wsflux2,wtauu2,wtauv2       &                      whflux2,wsflux2,wtauu2,wtauv2
1132        _RL whflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL whflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1133        _RL wsflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1134        _RL wtauu2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtauu2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1135        _RL wtauv2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtauv2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1136                                                                                                    
1137                                          
1138    
1139  #if (defined (ALLOW_OBCSN_COST_CONTRIBUTION) || \  #if (defined (ALLOW_OBCSN_COST_CONTRIBUTION) || \
1140       defined (ALLOW_OBCSN_CONTROL))       defined (ALLOW_OBCSN_CONTROL))
# Line 864  c     scatydat   - reference meridional Line 1175  c     scatydat   - reference meridional
1175  c     sstdat     - reference sea surface temperature data.  c     sstdat     - reference sea surface temperature data.
1176  c     tmidat     - reference TMI sea surface temperature data.  c     tmidat     - reference TMI sea surface temperature data.
1177  c     sssdat     - reference sea surface temperature data.  c     sssdat     - reference sea surface temperature data.
1178    c     bpdat      - bottom pressure from time-varying GRACE.
1179  c     tauxmask   - mask for reference wind stress data.  c     tauxmask   - mask for reference wind stress data.
1180  c     tauymask   - mask for reference wind stress data.  c     tauymask   - mask for reference wind stress data.
1181  c     scatxmask  - mask for scat wind stress data.  c     scatxmask  - mask for scat wind stress data.
1182  c     scatymask  - mask for scat wind stress data.  c     scatymask  - mask for scat wind stress data.
1183  c     sstmask    - mask for reference sea surface temperature data.  c     sstmask    - mask for reference sea surface temperature data.
1184  c     tmimask    - mask for reference sea surface temperature data.  c     tmimask    - mask for reference sea surface temperature data.
1185  c     sssmask    - mask for reference sea surface temperature data.  c     sssmask    - mask for reference sea surface temperature data.
# Line 880  c     ersobs     - ERS data. Line 1192  c     ersobs     - ERS data.
1192  c     ersmask    - mask for ERS data.  c     ersmask    - mask for ERS data.
1193  c     ctdtobs    - CTD temperature data  c     ctdtobs    - CTD temperature data
1194  c     ctdsobs    - CTD salinity data  c     ctdsobs    - CTD salinity data
1195  c     xbtobs     - XBT data  c     xbtobs     - XBT data
1196  c     argot      - ARGO  temperature data  c     argot      - ARGO  temperature data
1197  c     argos      - ARGO  salt data  c     argos      - ARGO  salt data
1198  c     udriftdat  - drifters zonal velocities  c     udriftdat  - drifters zonal velocities
1199  c     vdriftdat  - drifters meridional velocities  c     vdriftdat  - drifters meridional velocities
1200    
# Line 893  c     vdriftdat  - drifters meridional v Line 1205  c     vdriftdat  - drifters meridional v
1205       &                     sstdat,       &                     sstdat,
1206       &                     tmidat,       &                     tmidat,
1207       &                     sssdat,       &                     sssdat,
1208         &                     bpdat,
1209       &                     sstmask,       &                     sstmask,
1210       &                     tmimask,       &                     tmimask,
1211       &                     sssmask,       &                     sssmask,
1212         &                     bpmask,
1213       &                     tauxmask,       &                     tauxmask,
1214       &                     tauymask,       &                     tauymask,
1215       &                     scatxmask,       &                     scatxmask,
# Line 907  c     vdriftdat  - drifters meridional v Line 1221  c     vdriftdat  - drifters meridional v
1221       &                     tpmask,       &                     tpmask,
1222       &                     ersobs,       &                     ersobs,
1223       &                     ersmask,       &                     ersmask,
1224         &                     gfoobs,
1225         &                     gfomask,
1226       &                     ctdtobs,       &                     ctdtobs,
1227       &                     ctdsobs,       &                     ctdsobs,
1228       &                     xbtobs,       &                     xbtobs,
# Line 916  c     vdriftdat  - drifters meridional v Line 1232  c     vdriftdat  - drifters meridional v
1232       &                     vdriftdat,       &                     vdriftdat,
1233       &                     curmtruobs,       &                     curmtruobs,
1234       &                     curmtrvobs       &                     curmtrvobs
1235        
1236        _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)
1237        _RL scatxdat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatxdat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1238        _RL scatydat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatydat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1239        _RL sstdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sstdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1240        _RL tmidat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tmidat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1241        _RL sssdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sssdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1242          _RL bpdat     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1243        _RL tauxmask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tauxmask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1244        _RL tauymask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tauymask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1245        _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 930  c     vdriftdat  - drifters meridional v Line 1247  c     vdriftdat  - drifters meridional v
1247        _RL sstmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sstmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1248        _RL tmimask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tmimask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1249        _RL sssmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sssmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1250          _RL bpmask    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1251        _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)
1252        _RL tpmean    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpmean    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1253        _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 937  c     vdriftdat  - drifters meridional v Line 1255  c     vdriftdat  - drifters meridional v
1255        _RL tpmask    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpmask    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1256        _RL ersobs    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL ersobs    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1257        _RL ersmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL ersmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1258          _RL gfoobs    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1259          _RL gfomask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1260        _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)
1261        _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)
1262        _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 962  c     topexfile     - reference data fil Line 1282  c     topexfile     - reference data fil
1282  c                     (TOPEX/POSEIDON).  c                     (TOPEX/POSEIDON).
1283  c     ersfile       - reference data file for sea surface height data  c     ersfile       - reference data file for sea surface height data
1284  c                     (ERS).  c                     (ERS).
1285  c ctdtfile, ctdsfile- reference data file for temperature and salinity  c ctdtfile, ctdsfile- reference data file for temperature and salinity
1286  c                     from CTD  c                     from CTD
1287  c ctdtclimfile, ctdsclimfile- reference data file for temperature  c ctdtclimfile, ctdsclimfile- reference data file for temperature
1288  c                      and salinity from CTD with out timetag  c                      and salinity from CTD with out timetag
# Line 979  c     driftfile     - reference data fil Line 1299  c     driftfile     - reference data fil
1299       &                     sstdatfile,       &                     sstdatfile,
1300       &                     tmidatfile,       &                     tmidatfile,
1301       &                     sssdatfile,       &                     sssdatfile,
1302         &                     bpdatfile,
1303       &                     topexmeanfile,       &                     topexmeanfile,
1304       &                     topexfile,       &                     topexfile,
1305       &                     ersfile,       &                     ersfile,
1306         &                     gfofile,
1307       &                     ctdtfile,       &                     ctdtfile,
1308       &                     ctdsfile,       &                     ctdsfile,
1309       &                     ctdtclimfile,       &                     ctdtclimfile,
# Line 989  c     driftfile     - reference data fil Line 1311  c     driftfile     - reference data fil
1311       &                     xbtfile,       &                     xbtfile,
1312       &                     argotfile,       &                     argotfile,
1313       &                     argosfile,       &                     argosfile,
1314       &                     udriftfile,       &                     udriftfile,
1315       &                     vdriftfile,       &                     vdriftfile,
1316         &                     gencost_datafile,
1317       &                     curmtrufile,       &                     curmtrufile,
1318       &                     curmtrvfile       &                     curmtrvfile
1319    
# Line 1001  c     driftfile     - reference data fil Line 1324  c     driftfile     - reference data fil
1324        character*(MAX_LEN_FNAM) sstdatfile        character*(MAX_LEN_FNAM) sstdatfile
1325        character*(MAX_LEN_FNAM) tmidatfile        character*(MAX_LEN_FNAM) tmidatfile
1326        character*(MAX_LEN_FNAM) sssdatfile        character*(MAX_LEN_FNAM) sssdatfile
1327          character*(MAX_LEN_FNAM) bpdatfile
1328        character*(MAX_LEN_FNAM) topexmeanfile        character*(MAX_LEN_FNAM) topexmeanfile
1329        character*(MAX_LEN_FNAM) topexfile        character*(MAX_LEN_FNAM) topexfile
1330        character*(MAX_LEN_FNAM) ersfile        character*(MAX_LEN_FNAM) ersfile
1331          character*(MAX_LEN_FNAM) gfofile
1332        character*(MAX_LEN_FNAM) ctdtfile        character*(MAX_LEN_FNAM) ctdtfile
1333        character*(MAX_LEN_FNAM) ctdsfile        character*(MAX_LEN_FNAM) ctdsfile
1334        character*(MAX_LEN_FNAM) ctdtclimfile        character*(MAX_LEN_FNAM) ctdtclimfile
# Line 1012  c     driftfile     - reference data fil Line 1337  c     driftfile     - reference data fil
1337        character*(MAX_LEN_FNAM) argotfile        character*(MAX_LEN_FNAM) argotfile
1338        character*(MAX_LEN_FNAM) argosfile        character*(MAX_LEN_FNAM) argosfile
1339        character*(MAX_LEN_FNAM) argofile        character*(MAX_LEN_FNAM) argofile
1340          character*(MAX_LEN_FNAM) gencost_datafile(NGENCOST)
1341        character*(MAX_LEN_FNAM) udriftfile        character*(MAX_LEN_FNAM) udriftfile
1342        character*(MAX_LEN_FNAM) vdriftfile              character*(MAX_LEN_FNAM) vdriftfile
1343        character*(MAX_LEN_FNAM) curmtrufile        character*(MAX_LEN_FNAM) curmtrufile
1344        character*(MAX_LEN_FNAM) curmtrvfile        character*(MAX_LEN_FNAM) curmtrvfile
1345    
# Line 1025  c     using_ers - flag that indicates th Line 1351  c     using_ers - flag that indicates th
1351    
1352        common /ecco_cost_data_flags/        common /ecco_cost_data_flags/
1353       &                         using_topex,       &                         using_topex,
1354       &                         using_ers       &                         using_ers,
1355         &                         using_gfo
1356        logical using_topex        logical using_topex
1357        logical using_ers        logical using_ers
1358          logical using_gfo
1359    
1360  c     Calendar information for the observations:  c     Calendar information for the observations:
1361  c     ==========================================  c     ==========================================
# Line 1046  c     sshperiod      - sampling interval Line 1374  c     sshperiod      - sampling interval
1374       &                           argosstartdate,       &                           argosstartdate,
1375       &                           tmistartdate,       &                           tmistartdate,
1376       &                           sssstartdate,       &                           sssstartdate,
1377         &                           bpstartdate,
1378       &                           topexstartdate,       &                           topexstartdate,
1379       &                           ersstartdate       &                           ersstartdate,
1380         &                           gfostartdate
1381        integer scatxstartdate(4)        integer scatxstartdate(4)
1382        integer scatystartdate(4)        integer scatystartdate(4)
1383        integer sststartdate(4)        integer sststartdate(4)
# Line 1055  c     sshperiod      - sampling interval Line 1385  c     sshperiod      - sampling interval
1385        integer argosstartdate(4)        integer argosstartdate(4)
1386        integer tmistartdate(4)        integer tmistartdate(4)
1387        integer sssstartdate(4)        integer sssstartdate(4)
1388          integer bpstartdate(4)
1389        integer topexstartdate(4)        integer topexstartdate(4)
1390        integer ersstartdate(4)        integer ersstartdate(4)
1391          integer gfostartdate(4)
1392    
1393        common /ecco_cost_data_aux_i/        common /ecco_cost_data_aux_i/
1394       &                           tmistartdate1,       &                           tmistartdate1,
# Line 1065  c     sshperiod      - sampling interval Line 1397  c     sshperiod      - sampling interval
1397       &                           sststartdate2,       &                           sststartdate2,
1398       &                           sssstartdate1,       &                           sssstartdate1,
1399       &                           sssstartdate2,       &                           sssstartdate2,
1400         &                           bpstartdate1,
1401         &                           bpstartdate2,
1402       &                           argotstartdate1,       &                           argotstartdate1,
1403       &                           argotstartdate2,       &                           argotstartdate2,
1404       &                           argosstartdate1,       &                           argosstartdate1,
# Line 1073  c     sshperiod      - sampling interval Line 1407  c     sshperiod      - sampling interval
1407       &                           topexstartdate2,       &                           topexstartdate2,
1408       &                           ersstartdate1,       &                           ersstartdate1,
1409       &                           ersstartdate2,       &                           ersstartdate2,
1410         &                           gfostartdate1,
1411         &                           gfostartdate2,
1412       &                           scatstartdate1,       &                           scatstartdate1,
1413       &                           scatstartdate2,       &                           scatstartdate2
      &                           costIceStart1,  
      &                           costIceStart2,  
      &                           costIceEnd1,  
      &                           costIceEnd2  
1414    
1415        integer tmistartdate1        integer tmistartdate1
1416        integer tmistartdate2        integer tmistartdate2
# Line 1086  c     sshperiod      - sampling interval Line 1418  c     sshperiod      - sampling interval
1418        integer sststartdate2        integer sststartdate2
1419        integer sssstartdate1        integer sssstartdate1
1420        integer sssstartdate2        integer sssstartdate2
1421          integer bpstartdate1
1422          integer bpstartdate2
1423        integer argotstartdate1        integer argotstartdate1
1424        integer argotstartdate2        integer argotstartdate2
1425        integer argosstartdate1        integer argosstartdate1
# Line 1094  c     sshperiod      - sampling interval Line 1428  c     sshperiod      - sampling interval
1428        integer topexstartdate2        integer topexstartdate2
1429        integer ersstartdate1        integer ersstartdate1
1430        integer ersstartdate2        integer ersstartdate2
1431          integer gfostartdate1
1432          integer gfostartdate2
1433        integer scatstartdate1        integer scatstartdate1
1434        integer scatstartdate2        integer scatstartdate2
       integer costIceStart1  
       integer costIceStart2  
       integer costIceEnd1  
       integer costIceEnd2  
1435    
1436        common /ecco_cost_data_times_r/        common /ecco_cost_data_times_r/
1437       &                           topexperiod,       &                           topexperiod,
1438       &                           ersperiod,       &                           ersperiod,
1439       &                           scatperiod,       &                           gfoperiod,
1440       &                           costIceStart,       &                           scatperiod
      &                           costIceEnd  
1441        _RL topexperiod        _RL topexperiod
1442        _RL ersperiod        _RL ersperiod
1443          _RL gfoperiod
1444        _RL scatperiod        _RL scatperiod
       _RL costIceStart  
       _RL costIceEnd  
1445    
1446  c     cost_ice_flag  - cost_ice flag (see cost_ice.F)        common /ecco_cost_data_detrend/
1447         &                           topexintercept,
1448         &                           ersintercept,
1449         &                           gfointercept,
1450         &                           topexslope,
1451         &                           ersslope,
1452         &                           gfoslope
1453          _RL topexintercept
1454          _RL ersintercept
1455          _RL gfointercept
1456          _RL topexslope
1457          _RL ersslope
1458          _RL gfoslope
1459    
1460    #ifdef ALLOW_SSH_COST_CONTRIBUTION
1461          common /ecco_ssh_daymask_r/
1462         &       tpTimeMask, ersTimeMask, gfoTimeMask
1463          _RL tpTimeMask(maxNumDays)
1464          _RL ersTimeMask(maxNumDays)
1465          _RL gfoTimeMask(maxNumDays)
1466    #endif
1467    
1468        common /ecco_cost_ice_i/ cost_ice_flag        common /ecco_ssh_daymask_c/
1469        integer cost_ice_flag       &       tpTimeMaskFile, ersTimeMaskFile, gfoTimeMaskFile
1470          character*(MAX_LEN_FNAM) tpTimeMaskFile
1471          character*(MAX_LEN_FNAM) ersTimeMaskFile
1472          character*(MAX_LEN_FNAM) gfoTimeMaskFile
1473    
1474  c     ==================================================================  c     ==================================================================
1475  c     END OF HEADER COST  c     END OF HEADER COST
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.38
  ViewVC Help
Powered by ViewVC 1.1.22