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

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

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

revision 1.5 by heimbach, Fri May 27 22:10:27 2005 UTC revision 1.33 by heimbach, Sat Oct 18 12:42:30 2008 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     Averaged Fields:  c     Averaged Fields:
41  c     ================  c     ================
# Line 43  c             intantaneous salinities. Line 49  c             intantaneous salinities.
49  c     psbar - contains the averaged surface pressure field after the call  c     psbar - contains the averaged surface pressure field after the call
50  c             to subroutine POST_DAILY. Before, it accumulates the  c             to subroutine POST_DAILY. Before, it accumulates the
51  c             intantaneous surface pressure field.  c             intantaneous surface pressure field.
52  c     ubar  - contains the averaged zonal velocity component for the  c     ubar  - contains the averaged zonal velocity component for the
53  c             whole integration period. Before, it accumulates the  c             whole integration period. Before, it accumulates the
54  c             intantaneous field.  c             intantaneous field.
55  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 68  c     sfluxbar  - contains the averaged
68  c             whole integration period. Before, it accumulates the  c             whole integration period. Before, it accumulates the
69  c             intantaneous field.  c             intantaneous field.
70    
71        common /averages_r/        common /averages_r/
72       &                    tbar,       &                    tbar,
73       &                    sbar,       &                    sbar,
74         &                    sstbar,
75       &                    psbar,       &                    psbar,
76         &                    bpbar,
77       &                    ubar,       &                    ubar,
78       &                    vbar,       &                    vbar,
79       &                    wbar,       &                    wbar,
# Line 73  c             intantaneous field. Line 81  c             intantaneous field.
81       &                    tauybar,       &                    tauybar,
82       &                    hfluxbar,       &                    hfluxbar,
83       &                    sfluxbar,       &                    sfluxbar,
84       &                    Slmean,       &                    Slmean,
85       &                    Tlmean,       &                    Tlmean,
86       &                    wlmean,       &                    wlmean,
87       &                    Sfmean,       &                    Sfmean,
88       &                    Tfmean,                         &                    Tfmean,
89       &                    sbar_gen,       &                    sbar_gen,
90       &                    tbar_gen,       &                    tbar_gen,
91       &                    wfmean       &                    wfmean
# Line 103  c             intantaneous field. Line 111  c             intantaneous field.
111        _RL sbar_gen        _RL sbar_gen
112  #endif  #endif
113    
114    #ifdef ALLOW_SEAICE_COST_AREASST
115          _RL sstbar  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
116    #else
117          _RL sstbar
118    #endif
119    
120  #if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \  #if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \
121       defined (ALLOW_CTDS_COST_CONTRIBUTION) || \       defined (ALLOW_CTDS_COST_CONTRIBUTION) || \
122       defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \       defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \
# Line 122  c             intantaneous field. Line 136  c             intantaneous field.
136        _RL psbar        _RL psbar
137  #endif  #endif
138    
139    #ifdef ALLOW_BP_COST_CONTRIBUTION
140          _RL bpbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
141    #else
142          _RL bpbar
143    #endif
144    
145  #if (defined (ALLOW_DRIFTER_COST_CONTRIBUTION) || \  #if (defined (ALLOW_DRIFTER_COST_CONTRIBUTION) || \
146       defined (ALLOW_OBCS_COST_CONTRIBUTION))       defined (ALLOW_OBCS_COST_CONTRIBUTION))
147        _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 166  c             intantaneous field.
166        _RL   Tlmean        _RL   Tlmean
167        _RL   Slmean        _RL   Slmean
168        _RL   Tfmean        _RL   Tfmean
169        _RL   Sfmean                _RL   Sfmean
170  #endif  #endif
171    
172  #ifdef ALLOW_DRIFTW_COST_CONTRIBUTION  #ifdef ALLOW_DRIFTW_COST_CONTRIBUTION
# Line 161  c             intantaneous field. Line 181  c             intantaneous field.
181        _RL tauxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)        _RL tauxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
182        _RL tauybar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)        _RL tauybar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
183  #else  #else
184        _RL tauxbar        _RL tauxbar
185        _RL tauybar        _RL tauybar
186  #endif  #endif
187    
188  #ifdef ALLOW_MEAN_HFLUX_COST_CONTRIBUTION  #ifdef ALLOW_MEAN_HFLUX_COST_CONTRIBUTION
189        _RL hfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)        _RL hfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
190  #else  #else
191        _RL hfluxbar        _RL hfluxbar
192  #endif  #endif
193    
194  #ifdef ALLOW_MEAN_SFLUX_COST_CONTRIBUTION  #ifdef ALLOW_MEAN_SFLUX_COST_CONTRIBUTION
195        _RL sfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)            _RL sfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
196  #else  #else
197        _RL sfluxbar        _RL sfluxbar
198  #endif  #endif
199    
200    
201        common /averages_c/        common /averages_c/
202       &                    tbarfile,       &                    tbarfile,
203       &                    sbarfile,       &                    sbarfile,
204         &                    sstbarfile,
205       &                    psbarfile,       &                    psbarfile,
206         &                    bpbarfile,
207       &                    ubarfile,       &                    ubarfile,
208       &                    vbarfile,       &                    vbarfile,
209       &                    wbarfile,       &                    wbarfile,
210       &                    tauxbarfile,       &                    tauxbarfile,
211       &                    tauybarfile,       &                    tauybarfile,
212       &                    hfluxbarfile,       &                    hfluxbarfile,
213       &                    sfluxbarfile       &                    sfluxbarfile,
214         &                    costTranspDataFile
215        character*(MAX_LEN_FNAM) tbarfile        character*(MAX_LEN_FNAM) tbarfile
216        character*(MAX_LEN_FNAM) sbarfile        character*(MAX_LEN_FNAM) sbarfile
217          character*(MAX_LEN_FNAM) sstbarfile
218        character*(MAX_LEN_FNAM) psbarfile        character*(MAX_LEN_FNAM) psbarfile
219          character*(MAX_LEN_FNAM) bpbarfile
220        character*(MAX_LEN_FNAM) ubarfile        character*(MAX_LEN_FNAM) ubarfile
221        character*(MAX_LEN_FNAM) vbarfile        character*(MAX_LEN_FNAM) vbarfile
222        character*(MAX_LEN_FNAM) wbarfile        character*(MAX_LEN_FNAM) wbarfile
# Line 199  c             intantaneous field. Line 224  c             intantaneous field.
224        character*(MAX_LEN_FNAM) tauybarfile        character*(MAX_LEN_FNAM) tauybarfile
225        character*(MAX_LEN_FNAM) hfluxbarfile        character*(MAX_LEN_FNAM) hfluxbarfile
226        character*(MAX_LEN_FNAM) sfluxbarfile        character*(MAX_LEN_FNAM) sfluxbarfile
227          character*(MAX_LEN_FNAM) costTranspDataFile
228    
229    #ifdef ALLOW_TRANSPORT_COST_CONTRIBUTION
230          INTEGER maxNumDays
231          PARAMETER ( maxNumDays = 5480 )
232          common /averages_transp_r/
233         &                     transpbar
234         &                   , transpobs
235         &                   , wtransp
236          _RL transpbar(maxNumDays,nsx,nsy)
237          _RL transpobs(maxNumDays)
238          _RL wtransp(maxNumDays)
239    #endif
240    
241  c     file precision and field type  c     file precision and field type
242    
243        common /prec_type_cost/        common /prec_type_cost/
244       &                        cost_iprec,       &                        cost_iprec,
245       &                        cost_yftype       &                        cost_yftype
246    
# Line 256  c     objf_hfluxmm    - Global time-mean Line 294  c     objf_hfluxmm    - Global time-mean
294  c     objf_sfluxmm    - Global time-mean Salt flux contribution.  c     objf_sfluxmm    - Global time-mean Salt flux contribution.
295  c     objf_hmean - Mean sea surface height contribution.  c     objf_hmean - Mean sea surface height contribution.
296  c     objf_h     - Residual sea surface height contribution.  c     objf_h     - Residual sea surface height contribution.
297    c     objf_tp    - Residual sea surface height contribution from T/P
298    c     objf_ers   - Residual sea surface height contribution from T/P
299    c     objf_gfo   - Residual sea surface height contribution from T/P
300  c     objf_temp  - Temperature contribution.  c     objf_temp  - Temperature contribution.
301  c     objf_salt  - Salinity contribution.  c     objf_salt  - Salinity contribution.
302  c     objf_temp0 - Initial conditions Temperature contribution.  c     objf_temp0 - Initial conditions Temperature contribution.
# Line 263  c     objf_salt0 - Initial conditions Sa Line 304  c     objf_salt0 - Initial conditions Sa
304  c     objf_sst   - Sea surface temperature contribution.  c     objf_sst   - Sea surface temperature contribution.
305  c     objf_tmi   - Sea surface temperature contribution.  c     objf_tmi   - Sea surface temperature contribution.
306  c     objf_sss   - Sea surface salinity contribution.  c     objf_sss   - Sea surface salinity contribution.
307  c     objf_ctdt  - Temperature measurements from Woce CTD  c     objf_ctdt  - Temperature measurements from Woce CTD
308  c     objf_ctds  - Salinity measurements from Woce CTD  c     objf_ctds  - Salinity measurements from Woce CTD
309  c     objf_ctdtclim - Temperature measurements from Woce CTD without timetag  c     objf_ctdtclim - Temperature measurements from Woce CTD without timetag
310  c     objf_ctdsclim - Salinity measurements from Woce CTD without timetag  c     objf_ctdsclim - Salinity measurements from Woce CTD without timetag
311  c     objf_xbt   - XBT temperature data  c     objf_xbt   - XBT temperature data
312  c     objf_argot - ARGO temperature profiles  c     objf_argot - ARGO temperature profiles
313  c     objf_argos - ARGO salt profiles  c     objf_argos - ARGO salt profiles
314    c     objf_gencost - gencost user defined contribution
315  c     objf_scatxm - time-mean zonal SCAT  contribution  c     objf_scatxm - time-mean zonal SCAT  contribution
316  c     objf_scatym - time-mean meridional SCAT  contribution  c     objf_scatym - time-mean meridional SCAT  contribution
317  c     objf_scatx  - zonal SCAT  contribution  c     objf_scatx  - zonal SCAT  contribution
318  c     objf_scaty  - meridional SCAT  contribution  c     objf_scaty  - meridional SCAT  contribution
 c     objf_ice    - sea-ice volume  
319  c     objf_kapgm  - kappa GM contribution  c     objf_kapgm  - kappa GM contribution
320    c     objf_kapredi  - kappa REDI contribution
321  c     objf_diffkr - diffusion contribution  c     objf_diffkr - diffusion contribution
322  c     objf_theta_ini_fin - final vs. initial theta misfit  c     objf_theta_ini_fin - final vs. initial theta misfit
323  c     objf_salt_ini_fin  - final vs. initial salt misfit  c     objf_salt_ini_fin  - final vs. initial salt misfit
324    c     objf_eddytau - eddy streamfunction contribution
325    c     objf_bottomdrag - bottom drag contribution
326  c  c
327  c     mult_"var" - multipliers for the individual cost  c     mult_"var" - multipliers for the individual cost
328  c                  function contributions.  c                  function contributions.
# Line 287  c                  function contribution Line 331  c                  function contribution
331       &                objf_hflux,       &                objf_hflux,
332       &                objf_hfluxm,       &                objf_hfluxm,
333       &                objf_hfluxmm,       &                objf_hfluxmm,
334         &                objf_hfluxsmoo,
335       &                objf_sflux,       &                objf_sflux,
336       &                objf_sfluxm,       &                objf_sfluxm,
337       &                objf_sfluxmm,       &                objf_sfluxmm,
338         &                objf_sfluxsmoo,
339       &                objf_tauu,       &                objf_tauu,
340       &                objf_tauum,       &                objf_tauum,
341         &                objf_tauusmoo,
342       &                objf_tauv,       &                objf_tauv,
343       &                objf_tauvm,       &                objf_tauvm,
344         &                objf_tauvsmoo,
345       &                objf_hmean,       &                objf_hmean,
346       &                objf_h,       &                objf_h,
347         &                objf_tp,
348         &                objf_ers,
349         &                objf_gfo,
350       &                objf_temp,       &                objf_temp,
351       &                objf_salt,       &                objf_salt,
352       &                objf_temp0,       &                objf_temp0,
353       &                objf_salt0,       &                objf_salt0,
354         &                objf_temp0smoo,
355         &                objf_salt0smoo,
356       &                objf_sst,       &                objf_sst,
357       &                objf_tmi,       &                objf_tmi,
358       &                objf_sss,       &                objf_sss,
359         &                objf_bp,
360       &                objf_ctdt,       &                objf_ctdt,
361       &                objf_ctds,       &                objf_ctds,
362       &                objf_ctdtclim,       &                objf_ctdtclim,
# Line 310  c                  function contribution Line 364  c                  function contribution
364       &                objf_xbt,       &                objf_xbt,
365       &                objf_argot,       &                objf_argot,
366       &                objf_argos,       &                objf_argos,
367         &                objf_gencost,
368       &                objf_drift,       &                objf_drift,
369       &                objf_tdrift,       &                objf_tdrift,
370       &                objf_sdrift,       &                objf_sdrift,
# Line 320  c                  function contribution Line 375  c                  function contribution
375       &                objf_scatym,       &                objf_scatym,
376       &                objf_atemp,       &                objf_atemp,
377       &                objf_aqh,       &                objf_aqh,
378         &                objf_precip,
379         &                objf_swflux,
380         &                objf_swdown,
381         &                objf_snowprecip,
382         &                objf_lwflux,
383         &                objf_lwdown,
384         &                objf_evap,
385         &                objf_apressure,
386         &                objf_runoff,
387       &                objf_uwind,       &                objf_uwind,
388       &                objf_vwind,       &                objf_vwind,
389         &                objf_atempm,
390         &                objf_aqhm,
391         &                objf_precipm,
392         &                objf_swfluxm,
393         &                objf_swdownm,
394         &                objf_snowprecipm,
395         &                objf_lwfluxm,
396         &                objf_lwdownm,
397         &                objf_evapm,
398         &                objf_apressurem,
399         &                objf_runoffm,
400         &                objf_uwindm,
401         &                objf_vwindm,
402         &                objf_atempsmoo,
403         &                objf_aqhsmoo,
404         &                objf_precipsmoo,
405         &                objf_swfluxsmoo,
406         &                objf_swdownsmoo,
407         &                objf_snowprecipsmoo,
408         &                objf_lwfluxsmoo,
409         &                objf_lwdownsmoo,
410         &                objf_evapsmoo,
411         &                objf_apressuresmoo,
412         &                objf_runoffsmoo,
413         &                objf_uwindsmoo,
414         &                objf_vwindsmoo,
415       &                objf_obcsn,       &                objf_obcsn,
416       &                objf_obcss,       &                objf_obcss,
417       &                objf_obcsw,       &                objf_obcsw,
# Line 329  c                  function contribution Line 419  c                  function contribution
419       &                objf_obcsvol,       &                objf_obcsvol,
420       &                objf_curmtr,       &                objf_curmtr,
421       &                objf_ageos,       &                objf_ageos,
422       &                objf_ice,       &                objf_kapgm,
423         &                objf_kapredi,
424       &                objf_diffkr,       &                objf_diffkr,
425       &                objf_theta_ini_fin,       &                objf_theta_ini_fin,
426       &                objf_salt_ini_fin,       &                objf_salt_ini_fin,
427       &                objf_eddytau       &                objf_eddytau,
428         &                objf_bottomdrag,
429         &                objf_transp
430        _RL  objf_hflux  (nsx,nsy)        _RL  objf_hflux  (nsx,nsy)
431        _RL  objf_hfluxm (nsx,nsy)        _RL  objf_hfluxm (nsx,nsy)
432        _RL  objf_hfluxmm(nsx,nsy)        _RL  objf_hfluxmm
433          _RL  objf_hfluxsmoo  (nsx,nsy)
434        _RL  objf_sflux  (nsx,nsy)        _RL  objf_sflux  (nsx,nsy)
435        _RL  objf_sfluxm (nsx,nsy)        _RL  objf_sfluxm (nsx,nsy)
436        _RL  objf_sfluxmm(nsx,nsy)        _RL  objf_sfluxmm
437          _RL  objf_sfluxsmoo  (nsx,nsy)
438        _RL  objf_tauu   (nsx,nsy)        _RL  objf_tauu   (nsx,nsy)
439        _RL  objf_tauum  (nsx,nsy)        _RL  objf_tauum  (nsx,nsy)
440          _RL  objf_tauusmoo   (nsx,nsy)
441        _RL  objf_tauv   (nsx,nsy)        _RL  objf_tauv   (nsx,nsy)
442        _RL  objf_tauvm  (nsx,nsy)        _RL  objf_tauvm  (nsx,nsy)
443          _RL  objf_tauvsmoo   (nsx,nsy)
444        _RL  objf_hmean        _RL  objf_hmean
445        _RL  objf_h    (nsx,nsy)        _RL  objf_h    (nsx,nsy)
446          _RL  objf_tp   (nsx,nsy)
447          _RL  objf_ers  (nsx,nsy)
448          _RL  objf_gfo  (nsx,nsy)
449        _RL  objf_temp (nsx,nsy)        _RL  objf_temp (nsx,nsy)
450        _RL  objf_salt (nsx,nsy)        _RL  objf_salt (nsx,nsy)
451        _RL  objf_temp0(nsx,nsy)        _RL  objf_temp0(nsx,nsy)
452        _RL  objf_salt0(nsx,nsy)        _RL  objf_salt0(nsx,nsy)
453          _RL  objf_temp0smoo(nsx,nsy)
454          _RL  objf_salt0smoo(nsx,nsy)
455        _RL  objf_sst  (nsx,nsy)        _RL  objf_sst  (nsx,nsy)
456        _RL  objf_tmi  (nsx,nsy)        _RL  objf_tmi  (nsx,nsy)
457        _RL  objf_sss  (nsx,nsy)        _RL  objf_sss  (nsx,nsy)
458          _RL  objf_bp   (nsx,nsy)
459        _RL  objf_ctdt (nsx,nsy)        _RL  objf_ctdt (nsx,nsy)
460        _RL  objf_ctds (nsx,nsy)        _RL  objf_ctds (nsx,nsy)
461        _RL  objf_ctdtclim (nsx,nsy)        _RL  objf_ctdtclim (nsx,nsy)
# Line 361  c                  function contribution Line 463  c                  function contribution
463        _RL  objf_xbt  (nsx,nsy)        _RL  objf_xbt  (nsx,nsy)
464        _RL  objf_argot(nsx,nsy)        _RL  objf_argot(nsx,nsy)
465        _RL  objf_argos(nsx,nsy)        _RL  objf_argos(nsx,nsy)
466          _RL  objf_gencost(NGENCOST,nsx,nsy)
467        _RL  objf_drift(nsx,nsy)        _RL  objf_drift(nsx,nsy)
468        _RL  objf_tdrift(nsx,nsy)        _RL  objf_tdrift(nsx,nsy)
469        _RL  objf_sdrift(nsx,nsy)        _RL  objf_sdrift(nsx,nsy)
# Line 371  c                  function contribution Line 474  c                  function contribution
474        _RL  objf_scatym(nsx,nsy)        _RL  objf_scatym(nsx,nsy)
475        _RL  objf_atemp(nsx,nsy)        _RL  objf_atemp(nsx,nsy)
476        _RL  objf_aqh  (nsx,nsy)        _RL  objf_aqh  (nsx,nsy)
477          _RL  objf_precip(nsx,nsy)
478          _RL  objf_swflux(nsx,nsy)
479          _RL  objf_swdown(nsx,nsy)
480          _RL  objf_snowprecip(nsx,nsy)
481          _RL  objf_lwflux(nsx,nsy)
482          _RL  objf_lwdown(nsx,nsy)
483          _RL  objf_evap(nsx,nsy)
484          _RL  objf_apressure(nsx,nsy)
485          _RL  objf_runoff(nsx,nsy)
486        _RL  objf_uwind(nsx,nsy)        _RL  objf_uwind(nsx,nsy)
487        _RL  objf_vwind(nsx,nsy)        _RL  objf_vwind(nsx,nsy)
488          _RL  objf_atempm(nsx,nsy)
489          _RL  objf_aqhm  (nsx,nsy)
490          _RL  objf_precipm(nsx,nsy)
491          _RL  objf_swfluxm(nsx,nsy)
492          _RL  objf_swdownm(nsx,nsy)
493          _RL  objf_snowprecipm(nsx,nsy)
494          _RL  objf_lwfluxm(nsx,nsy)
495          _RL  objf_lwdownm(nsx,nsy)
496          _RL  objf_evapm(nsx,nsy)
497          _RL  objf_apressurem(nsx,nsy)
498          _RL  objf_runoffm(nsx,nsy)
499          _RL  objf_uwindm(nsx,nsy)
500          _RL  objf_vwindm(nsx,nsy)
501          _RL  objf_atempsmoo(nsx,nsy)
502          _RL  objf_aqhsmoo  (nsx,nsy)
503          _RL  objf_precipsmoo(nsx,nsy)
504          _RL  objf_swfluxsmoo(nsx,nsy)
505          _RL  objf_swdownsmoo(nsx,nsy)
506          _RL  objf_snowprecipsmoo(nsx,nsy)
507          _RL  objf_lwfluxsmoo(nsx,nsy)
508          _RL  objf_lwdownsmoo(nsx,nsy)
509          _RL  objf_evapsmoo(nsx,nsy)
510          _RL  objf_apressuresmoo(nsx,nsy)
511          _RL  objf_runoffsmoo(nsx,nsy)
512          _RL  objf_uwindsmoo(nsx,nsy)
513          _RL  objf_vwindsmoo(nsx,nsy)
514        _RL  objf_obcsn(nsx,nsy)        _RL  objf_obcsn(nsx,nsy)
515        _RL  objf_obcss(nsx,nsy)        _RL  objf_obcss(nsx,nsy)
516        _RL  objf_obcsw(nsx,nsy)        _RL  objf_obcsw(nsx,nsy)
# Line 380  c                  function contribution Line 518  c                  function contribution
518        _RL  objf_obcsvol        _RL  objf_obcsvol
519        _RL  objf_curmtr(nsx,nsy)        _RL  objf_curmtr(nsx,nsy)
520        _RL  objf_ageos(nsx,nsy)        _RL  objf_ageos(nsx,nsy)
       _RL  objf_ice  (nsx,nsy)  
521        _RL  objf_kapgm(nsx,nsy)        _RL  objf_kapgm(nsx,nsy)
522          _RL  objf_kapredi(nsx,nsy)
523        _RL  objf_diffkr(nsx,nsy)        _RL  objf_diffkr(nsx,nsy)
524        _RL  objf_theta_ini_fin(nsx,nsy)        _RL  objf_theta_ini_fin(nsx,nsy)
525        _RL  objf_salt_ini_fin(nsx,nsy)        _RL  objf_salt_ini_fin(nsx,nsy)
526        _RL  objf_eddytau(nsx,nsy)        _RL  objf_eddytau(nsx,nsy)
527          _RL  objf_bottomdrag(nsx,nsy)
528          _RL  objf_transp
529    
530        common /ecco_cost_num/        common /ecco_cost_num/
531       &                num_hflux,       &                num_hflux,
# Line 400  c                  function contribution Line 540  c                  function contribution
540       &                num_tauvm,       &                num_tauvm,
541       &                num_hmean,       &                num_hmean,
542       &                num_h,       &                num_h,
543         &                num_tp,
544         &                num_ers,
545         &                num_gfo,
546       &                num_temp,       &                num_temp,
547       &                num_salt,       &                num_salt,
548       &                num_temp0,       &                num_temp0,
# Line 407  c                  function contribution Line 550  c                  function contribution
550       &                num_sst,       &                num_sst,
551       &                num_tmi,       &                num_tmi,
552       &                num_sss,       &                num_sss,
553         &                num_bp,
554       &                num_ctdt,       &                num_ctdt,
555       &                num_ctds,       &                num_ctds,
556       &                num_ctdtclim,       &                num_ctdtclim,
# Line 414  c                  function contribution Line 558  c                  function contribution
558       &                num_xbt,       &                num_xbt,
559       &                num_argot,       &                num_argot,
560       &                num_argos,       &                num_argos,
561         &                num_gencost,
562       &                num_drift,       &                num_drift,
563       &                num_tdrift,       &                num_tdrift,
564       &                num_sdrift,       &                num_sdrift,
# Line 424  c                  function contribution Line 569  c                  function contribution
569       &                num_scatym,       &                num_scatym,
570       &                num_atemp,       &                num_atemp,
571       &                num_aqh,       &                num_aqh,
572         &                num_precip,
573         &                num_swflux,
574         &                num_swdown,
575         &                num_snowprecip,
576         &                num_lwflux,
577         &                num_lwdown,
578         &                num_evap,
579         &                num_apressure,
580         &                num_runoff,
581       &                num_uwind,       &                num_uwind,
582       &                num_vwind,       &                num_vwind,
583         &                num_atempm,
584         &                num_aqhm,
585         &                num_precipm,
586         &                num_swfluxm,
587         &                num_swdownm,
588         &                num_snowprecipm,
589         &                num_lwfluxm,
590         &                num_lwdownm,
591         &                num_evapm,
592         &                num_apressurem,
593         &                num_runoffm,
594         &                num_uwindm,
595         &                num_vwindm,
596       &                num_obcsn,       &                num_obcsn,
597       &                num_obcss,       &                num_obcss,
598       &                num_obcsw,       &                num_obcsw,
# Line 433  c                  function contribution Line 600  c                  function contribution
600       &                num_obcsvol,       &                num_obcsvol,
601       &                num_curmtr,       &                num_curmtr,
602       &                num_ageos,       &                num_ageos,
603       &                num_ice,       &                num_kapgm,
604         &                num_kapredi,
605       &                num_diffkr,       &                num_diffkr,
606       &                num_theta_ini_fin,       &                num_theta_ini_fin,
607       &                num_salt_ini_fin,       &                num_salt_ini_fin,
608       &                num_eddytau       &                num_eddytau,
609         &                num_bottomdrag,
610         &                num_transp
611    
612        _RL  num_hflux  (nsx,nsy)        _RL  num_hflux  (nsx,nsy)
613        _RL  num_hfluxm (nsx,nsy)        _RL  num_hfluxm (nsx,nsy)
# Line 451  c                  function contribution Line 621  c                  function contribution
621        _RL  num_tauvm  (nsx,nsy)        _RL  num_tauvm  (nsx,nsy)
622        _RL  num_hmean        _RL  num_hmean
623        _RL  num_h    (nsx,nsy)        _RL  num_h    (nsx,nsy)
624          _RL  num_tp   (nsx,nsy)
625          _RL  num_ers  (nsx,nsy)
626          _RL  num_gfo  (nsx,nsy)
627        _RL  num_temp (nsx,nsy)        _RL  num_temp (nsx,nsy)
628        _RL  num_salt (nsx,nsy)        _RL  num_salt (nsx,nsy)
629        _RL  num_temp0(nsx,nsy)        _RL  num_temp0(nsx,nsy)
630        _RL  num_salt0(nsx,nsy)        _RL  num_salt0(nsx,nsy)
631        _RL  num_sst  (nsx,nsy)        _RL  num_sst  (nsx,nsy)
632        _RL  num_tmi  (nsx,nsy)        _RL  num_tmi  (nsx,nsy)
633        _RL  num_sss  (nsx,nsy)        _RL  num_sss  (nsx,nsy)
634          _RL  num_bp   (nsx,nsy)
635        _RL  num_ctdt (nsx,nsy)        _RL  num_ctdt (nsx,nsy)
636        _RL  num_ctds (nsx,nsy)        _RL  num_ctds (nsx,nsy)
637        _RL  num_ctdtclim (nsx,nsy)        _RL  num_ctdtclim (nsx,nsy)
# Line 465  c                  function contribution Line 639  c                  function contribution
639        _RL  num_xbt  (nsx,nsy)        _RL  num_xbt  (nsx,nsy)
640        _RL  num_argot(nsx,nsy)        _RL  num_argot(nsx,nsy)
641        _RL  num_argos(nsx,nsy)        _RL  num_argos(nsx,nsy)
642          _RL  num_gencost(NGENCOST,nsx,nsy)
643        _RL  num_drift(nsx,nsy)        _RL  num_drift(nsx,nsy)
644        _RL  num_tdrift(nsx,nsy)        _RL  num_tdrift(nsx,nsy)
645        _RL  num_sdrift(nsx,nsy)        _RL  num_sdrift(nsx,nsy)
# Line 475  c                  function contribution Line 650  c                  function contribution
650        _RL  num_scatym(nsx,nsy)        _RL  num_scatym(nsx,nsy)
651        _RL  num_atemp(nsx,nsy)        _RL  num_atemp(nsx,nsy)
652        _RL  num_aqh  (nsx,nsy)        _RL  num_aqh  (nsx,nsy)
653          _RL  num_precip(nsx,nsy)
654          _RL  num_swflux(nsx,nsy)
655          _RL  num_swdown(nsx,nsy)
656          _RL  num_snowprecip(nsx,nsy)
657          _RL  num_lwflux(nsx,nsy)
658          _RL  num_lwdown(nsx,nsy)
659          _RL  num_evap(nsx,nsy)
660          _RL  num_apressure(nsx,nsy)
661          _RL  num_runoff(nsx,nsy)
662        _RL  num_uwind(nsx,nsy)        _RL  num_uwind(nsx,nsy)
663        _RL  num_vwind(nsx,nsy)        _RL  num_vwind(nsx,nsy)
664          _RL  num_atempm(nsx,nsy)
665          _RL  num_aqhm  (nsx,nsy)
666          _RL  num_precipm(nsx,nsy)
667          _RL  num_swfluxm(nsx,nsy)
668          _RL  num_swdownm(nsx,nsy)
669          _RL  num_snowprecipm(nsx,nsy)
670          _RL  num_lwfluxm(nsx,nsy)
671          _RL  num_lwdownm(nsx,nsy)
672          _RL  num_evapm(nsx,nsy)
673          _RL  num_apressurem(nsx,nsy)
674          _RL  num_runoffm(nsx,nsy)
675          _RL  num_uwindm(nsx,nsy)
676          _RL  num_vwindm(nsx,nsy)
677        _RL  num_obcsn(nsx,nsy)        _RL  num_obcsn(nsx,nsy)
678        _RL  num_obcss(nsx,nsy)        _RL  num_obcss(nsx,nsy)
679        _RL  num_obcsw(nsx,nsy)        _RL  num_obcsw(nsx,nsy)
# Line 484  c                  function contribution Line 681  c                  function contribution
681        _RL  num_obcsvol        _RL  num_obcsvol
682        _RL  num_curmtr(nsx,nsy)        _RL  num_curmtr(nsx,nsy)
683        _RL  num_ageos(nsx,nsy)        _RL  num_ageos(nsx,nsy)
       _RL  num_ice  (nsx,nsy)  
684        _RL  num_kapgm(nsx,nsy)        _RL  num_kapgm(nsx,nsy)
685          _RL  num_kapredi(nsx,nsy)
686        _RL  num_diffkr(nsx,nsy)        _RL  num_diffkr(nsx,nsy)
687        _RL  num_theta_ini_fin(nsx,nsy)        _RL  num_theta_ini_fin(nsx,nsy)
688        _RL  num_salt_ini_fin(nsx,nsy)        _RL  num_salt_ini_fin(nsx,nsy)
689        _RL  num_eddytau(nsx,nsy)        _RL  num_eddytau(nsx,nsy)
690          _RL  num_bottomdrag(nsx,nsy)
691          _RL  num_transp
692    
693        common /ecco_cost_aux_r/        common /ecco_cost_aux_r/
694       &                    mult_hflux,       &                    mult_hflux,
695       &                    mult_sflux,       &                    mult_sflux,
696         &                    mult_hfluxmm,
697         &                    mult_sfluxmm,
698       &                    mult_tauu,       &                    mult_tauu,
699       &                    mult_tauv,       &                    mult_tauv,
700       &                    mult_hmean,       &                    mult_hmean,
701       &                    mult_h,       &                    mult_h,
702         &                    mult_tp,
703         &                    mult_ers,
704         &                    mult_gfo,
705       &                    mult_temp,       &                    mult_temp,
706       &                    mult_salt,       &                    mult_salt,
707       &                    mult_temp0,       &                    mult_temp0,
# Line 505  c                  function contribution Line 709  c                  function contribution
709       &                    mult_sst,       &                    mult_sst,
710       &                    mult_tmi,       &                    mult_tmi,
711       &                    mult_sss,       &                    mult_sss,
712         &                    mult_bp,
713       &                    mult_ctdt,       &                    mult_ctdt,
714       &                    mult_ctds,       &                    mult_ctds,
715       &                    mult_ctdtclim,       &                    mult_ctdtclim,
# Line 512  c                  function contribution Line 717  c                  function contribution
717       &                    mult_xbt,       &                    mult_xbt,
718       &                    mult_argot,       &                    mult_argot,
719       &                    mult_argos,       &                    mult_argos,
720         &                    mult_gencost,
721       &                    mult_drift,       &                    mult_drift,
722       &                    mult_tdrift,       &                    mult_tdrift,
723       &                    mult_sdrift,       &                    mult_sdrift,
# Line 520  c                  function contribution Line 726  c                  function contribution
726       &                    mult_scaty,       &                    mult_scaty,
727       &                    mult_atemp,       &                    mult_atemp,
728       &                    mult_aqh,       &                    mult_aqh,
729         &                    mult_precip,
730         &                    mult_swflux,
731         &                    mult_swdown,
732         &                    mult_snowprecip,
733         &                    mult_lwflux,
734         &                    mult_lwdown,
735         &                    mult_evap,
736         &                    mult_apressure,
737         &                    mult_runoff,
738       &                    mult_uwind,       &                    mult_uwind,
739       &                    mult_vwind,       &                    mult_vwind,
740       &                    mult_obcsn,       &                    mult_obcsn,
# Line 529  c                  function contribution Line 744  c                  function contribution
744       &                    mult_obcsvol,       &                    mult_obcsvol,
745       &                    mult_curmtr,       &                    mult_curmtr,
746       &                    mult_ageos,       &                    mult_ageos,
      &                    mult_ice,  
747       &                    mult_kapgm,       &                    mult_kapgm,
748         &                    mult_kapredi,
749       &                    mult_diffkr,       &                    mult_diffkr,
750       &                    mult_ini_fin,       &                    mult_ini_fin,
751       &                    mult_eddytau       &                    mult_edtau,
752         &                    mult_bottomdrag,
753         &                    mult_smooth_ic,
754         &                    mult_smooth_bc,
755         &                    mult_transp
756        _RL  mult_hflux        _RL  mult_hflux
757        _RL  mult_sflux        _RL  mult_sflux
758          _RL  mult_hfluxmm
759          _RL  mult_sfluxmm
760        _RL  mult_tauu        _RL  mult_tauu
761        _RL  mult_tauv        _RL  mult_tauv
762        _RL  mult_hmean        _RL  mult_hmean
763        _RL  mult_h        _RL  mult_h
764          _RL  mult_tp
765          _RL  mult_ers
766          _RL  mult_gfo
767        _RL  mult_temp        _RL  mult_temp
768        _RL  mult_salt        _RL  mult_salt
769        _RL  mult_temp0        _RL  mult_temp0
# Line 548  c                  function contribution Line 771  c                  function contribution
771        _RL  mult_sst        _RL  mult_sst
772        _RL  mult_tmi        _RL  mult_tmi
773        _RL  mult_sss        _RL  mult_sss
774          _RL  mult_bp
775        _RL  mult_ctdt        _RL  mult_ctdt
776        _RL  mult_ctds        _RL  mult_ctds
777        _RL  mult_ctdtclim        _RL  mult_ctdtclim
# Line 555  c                  function contribution Line 779  c                  function contribution
779        _RL  mult_xbt        _RL  mult_xbt
780        _RL  mult_argot        _RL  mult_argot
781        _RL  mult_argos        _RL  mult_argos
782          _RL  mult_gencost(NGENCOST)
783        _RL  mult_drift        _RL  mult_drift
784        _RL  mult_tdrift        _RL  mult_tdrift
785        _RL  mult_sdrift        _RL  mult_sdrift
# Line 563  c                  function contribution Line 788  c                  function contribution
788        _RL  mult_scaty        _RL  mult_scaty
789        _RL  mult_atemp        _RL  mult_atemp
790        _RL  mult_aqh        _RL  mult_aqh
791          _RL  mult_precip
792          _RL  mult_swflux
793          _RL  mult_swdown
794          _RL  mult_snowprecip
795          _RL  mult_lwflux
796          _RL  mult_lwdown
797          _RL  mult_evap
798          _RL  mult_apressure
799          _RL  mult_runoff
800        _RL  mult_uwind        _RL  mult_uwind
801        _RL  mult_vwind        _RL  mult_vwind
802        _RL  mult_obcsn        _RL  mult_obcsn
# Line 572  c                  function contribution Line 806  c                  function contribution
806        _RL  mult_obcsvol        _RL  mult_obcsvol
807        _RL  mult_curmtr        _RL  mult_curmtr
808        _RL  mult_ageos        _RL  mult_ageos
       _RL  mult_ice  
809        _RL  mult_kapgm        _RL  mult_kapgm
810          _RL  mult_kapredi
811        _RL  mult_diffkr        _RL  mult_diffkr
812        _RL  mult_ini_fin        _RL  mult_ini_fin
813        _RL  mult_eddytau        _RL  mult_edtau
814          _RL  mult_bottomdrag
815          _RL  mult_smooth_ic
816          _RL  mult_smooth_bc
817          _RL  mult_transp
818    
819  c     Record counters relevant for the cost function evaluation.  c     Record counters relevant for the cost function evaluation.
820  c     ==========================================================  c     ==========================================================
821  c  c
822  c     nmonsrec - number of monthly records that will be generated by  c     nyearsrec - number of yearly records that will be generated by
823  c                the current model integration.  c                 the current model integration.
824  c     ndaysrec - number of  daily  records that will be generated by  c     nmonsrec  - number of monthly records that will be generated by
825  c                the current model integration.  c                 the current model integration.
826    c     ndaysrec  - number of  daily  records that will be generated by
827    c                 the current model integration.
828    
829        common /ecco_cost_i/        common /ecco_cost_i/
830         &                nyearsrec,
831       &                nmonsrec,       &                nmonsrec,
832       &                ndaysrec       &                ndaysrec,
833         &                nnztbar,
834         &                nnzsbar
835          integer nyearsrec
836        integer nmonsrec        integer nmonsrec
837        integer ndaysrec        integer ndaysrec
838          integer nnztbar
839          integer nnzsbar
840    
841    
842  c     Data files for the weights used in the cost function:  c     Data files for the weights used in the cost function:
# Line 615  c     salterrfile           - representa Line 861  c     salterrfile           - representa
861  c     temperrfile           - representation error due unresolved eddies  c     temperrfile           - representation error due unresolved eddies
862  c     velerrfile            - representation error  c     velerrfile            - representation error
863    
864        common /ecco_cost_c/        common /ecco_cost_c/
865       &                hflux_errfile,       &                hflux_errfile,
866       &                hfluxm_errfile,       &                hfluxm_errfile,
867       &                sflux_errfile,       &                sflux_errfile,
# Line 630  c     velerrfile            - representa Line 876  c     velerrfile            - representa
876       &                geoid_errfile,       &                geoid_errfile,
877       &                geoid_covariancefile,       &                geoid_covariancefile,
878       &                ssh_errfile,       &                ssh_errfile,
879         &                tp_errfile,
880         &                ers_errfile,
881         &                gfo_errfile,
882       &                ctdt_errfile,       &                ctdt_errfile,
883       &                ctds_errfile,       &                ctds_errfile,
884       &                drift_errfile,       &                drift_errfile,
885       &                udrifterrfile,       &                udrifterrfile,
886       &                vdrifterrfile,       &                vdrifterrfile,
887       &                salterrfile,       &                salterrfile,
888       &                temperrfile,       &                temperrfile,
889       &                velerrfile,       &                velerrfile,
890         &                salt0errfile,
891         &                temp0errfile,
892         &                vel0errfile,
893       &                atemp_errfile,       &                atemp_errfile,
894       &                aqh_errfile,       &                aqh_errfile,
895         &                precip_errfile,
896         &                swflux_errfile,
897         &                swdown_errfile,
898         &                snowprecip_errfile,
899         &                lwflux_errfile,
900         &                lwdown_errfile,
901         &                evap_errfile,
902         &                apressure_errfile,
903         &                runoff_errfile,
904         &                edtau_errfile,
905         &                kapgm_errfile,
906         &                kapredi_errfile,
907         &                diffkr_errfile,
908         &                bottomdrag_errfile,
909         &                gencost_errfile,
910       &                uwind_errfile,       &                uwind_errfile,
911       &                vwind_errfile       &                vwind_errfile
912        character*(MAX_LEN_FNAM) hflux_errfile        character*(MAX_LEN_FNAM) hflux_errfile
# Line 656  c     velerrfile            - representa Line 923  c     velerrfile            - representa
923        character*(MAX_LEN_FNAM) geoid_errfile        character*(MAX_LEN_FNAM) geoid_errfile
924        character*(MAX_LEN_FNAM) geoid_covariancefile        character*(MAX_LEN_FNAM) geoid_covariancefile
925        character*(MAX_LEN_FNAM) ssh_errfile        character*(MAX_LEN_FNAM) ssh_errfile
926        character*(MAX_LEN_FNAM) ctdt_errfile        character*(MAX_LEN_FNAM) tp_errfile
927        character*(MAX_LEN_FNAM) ctds_errfile        character*(MAX_LEN_FNAM) ers_errfile
928          character*(MAX_LEN_FNAM) gfo_errfile
929          character*(MAX_LEN_FNAM) ctdt_errfile
930          character*(MAX_LEN_FNAM) ctds_errfile
931        character*(MAX_LEN_FNAM) drift_errfile        character*(MAX_LEN_FNAM) drift_errfile
932        character*(MAX_LEN_FNAM) udrifterrfile        character*(MAX_LEN_FNAM) udrifterrfile
933        character*(MAX_LEN_FNAM) vdrifterrfile              character*(MAX_LEN_FNAM) vdrifterrfile
934        character*(MAX_LEN_FNAM) salterrfile        character*(MAX_LEN_FNAM) salterrfile
935        character*(MAX_LEN_FNAM) temperrfile        character*(MAX_LEN_FNAM) temperrfile
936        character*(MAX_LEN_FNAM) velerrfile        character*(MAX_LEN_FNAM) velerrfile
937          character*(MAX_LEN_FNAM) salt0errfile
938          character*(MAX_LEN_FNAM) temp0errfile
939          character*(MAX_LEN_FNAM) vel0errfile
940        character*(MAX_LEN_FNAM) atemp_errfile        character*(MAX_LEN_FNAM) atemp_errfile
941        character*(MAX_LEN_FNAM) aqh_errfile        character*(MAX_LEN_FNAM) aqh_errfile
942          character*(MAX_LEN_FNAM) precip_errfile
943          character*(MAX_LEN_FNAM) swflux_errfile
944          character*(MAX_LEN_FNAM) swdown_errfile
945          character*(MAX_LEN_FNAM) snowprecip_errfile
946          character*(MAX_LEN_FNAM) lwflux_errfile
947          character*(MAX_LEN_FNAM) lwdown_errfile
948          character*(MAX_LEN_FNAM) evap_errfile
949          character*(MAX_LEN_FNAM) apressure_errfile
950          character*(MAX_LEN_FNAM) runoff_errfile
951          character*(MAX_LEN_FNAM) edtau_errfile
952          character*(MAX_LEN_FNAM) kapgm_errfile
953          character*(MAX_LEN_FNAM) kapredi_errfile
954          character*(MAX_LEN_FNAM) diffkr_errfile
955          character*(MAX_LEN_FNAM) bottomdrag_errfile
956          character*(MAX_LEN_FNAM) gencost_errfile(NGENCOST)
957        character*(MAX_LEN_FNAM) uwind_errfile        character*(MAX_LEN_FNAM) uwind_errfile
958        character*(MAX_LEN_FNAM) vwind_errfile        character*(MAX_LEN_FNAM) vwind_errfile
959    
# Line 694  c     wctds      - weight for CTD salini Line 982  c     wctds      - weight for CTD salini
982  c     wudrift    - weight for mean zonal velocity from drifters.  c     wudrift    - weight for mean zonal velocity from drifters.
983  c     wvdrift    - weight for mean meridional velocity from drifters.  c     wvdrift    - weight for mean meridional velocity from drifters.
984    
985        common /ecco_cost_weights_r/        common /ecco_cost_weights_r/
986       &                      frame,       &                      frame,
987       &                      cosphi,       &                      cosphi,
988       &                      whflux,wsflux,wtauu,wtauv,       &                      whflux,wsflux,wtauu,wtauv,
989       &                      watemp,waqh,wuwind,wvwind,       &                      watemp,waqh,wprecip,wsnowprecip,
990         &                      wswflux,wswdown,wlwflux,wlwdown,
991         &                      wevap,wapressure,wrunoff,
992         &                      wbottomdrag,
993         &                      wuwind,wvwind,
994       &                      wscatx,wscaty,       &                      wscatx,wscaty,
995       &                      wtheta,wtheta2,wthetaLev,       &                      wtheta,wtheta2,wthetaLev,
996       &                      wsalt,wsalt2,wsaltLev,       &                      wsalt,wsalt2,wsaltLev,
997       &                      wsst,wsss,       &                      wdiffkr,wdiffkr2,wdiffkrFld,
998       &                      wtp,wers,       &                      wkapgm,wkapgm2,wkapgmFld,
999         &                      wkapredi,wkapredi2,wkaprediFld,
1000         &                      wedtaux,wedtaux2,wedtauxFld,
1001         &                      wedtauy,wedtauy2,wedtauyFld,
1002         &                      wsst,wsss,wbp,
1003         &                      wtp,wers,wgfo,
1004       &                      wp,       &                      wp,
1005       &                      wctdt,wctds,       &                      wctdt,wctds,
1006       &                      wudrift,wvdrift,       &                      wudrift,wvdrift,
# Line 727  c     wvdrift    - weight for mean merid Line 1024  c     wvdrift    - weight for mean merid
1024        _RL wscaty  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wscaty  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1025        _RL watemp  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL watemp  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1026        _RL waqh    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL waqh    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1027          _RL wprecip (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1028          _RL wswflux (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1029          _RL wswdown (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1030          _RL wsnowprecip (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1031          _RL wlwflux (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1032          _RL wlwdown (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1033          _RL wevap   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1034          _RL wapressure(1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1035          _RL wrunoff (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1036          _RL wbottomdrag (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
1037        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1038        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1039        _RL wtheta  (                            nr,nsx,nsy)        _RL wtheta  (                            nr,nsx,nsy)
# Line 737  c     wvdrift    - weight for mean merid Line 1044  c     wvdrift    - weight for mean merid
1044        _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)
1045        _RL wsst    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsst    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1046        _RL wsss    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsss    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1047          _RL wbp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1048        _RL wtp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1049        _RL wers    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wers    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1050          _RL wgfo    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1051        _RL wp      (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wp      (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1052        _RL wctdt   (                            nr,nsx,nsy)        _RL wctdt   (                            nr,nsx,nsy)
1053        _RL wctds   (                            nr,nsx,nsy)        _RL wctds   (                            nr,nsx,nsy)
# Line 748  c     wvdrift    - weight for mean merid Line 1057  c     wvdrift    - weight for mean merid
1057        _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)
1058        _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)
1059        _RL wbaro        _RL wbaro
1060          _RL wdiffkr (                            nr,nsx,nsy)
1061          _RL wdiffkr2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1062          _RL wdiffkrFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1063          _RL wkapgm (                            nr,nsx,nsy)
1064          _RL wkapgm2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1065          _RL wkapgmFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1066          _RL wkapredi (                            nr,nsx,nsy)
1067          _RL wkapredi2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1068          _RL wkaprediFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1069          _RL wedtaux (                            nr,nsx,nsy)
1070          _RL wedtaux2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1071          _RL wedtauxFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1072          _RL wedtauy (                            nr,nsx,nsy)
1073          _RL wedtauy2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1074          _RL wedtauyFld (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
1075    
1076          common /ecco_cost_weights_0_r/
1077         &        whflux0, wsflux0, wtau0,
1078         &        watemp0, waqh0, wprecip0, wsnowprecip0, wwind0,
1079         &        wswflux0, wswdown0, wlwflux0, wlwdown0,
1080         &        wevap0, wapressure0, wrunoff0, wkapredi0,
1081         &        wbottomdrag0,wdiffkr0, wkapgm0, wedtau0
1082          _RL whflux0
1083          _RL wsflux0
1084          _RL wtau0
1085          _RL watemp0
1086          _RL waqh0
1087          _RL wprecip0
1088          _RL wswflux0
1089          _RL wswdown0
1090          _RL wsnowprecip0
1091          _RL wlwflux0
1092          _RL wlwdown0
1093          _RL wevap0
1094          _RL wapressure0
1095          _RL wrunoff0
1096          _RL wbottomdrag0
1097          _RL wwind0
1098          _RL wdiffkr0
1099          _RL wkapgm0
1100          _RL wkapredi0
1101          _RL wedtau0
1102    
1103          common /ecco_cost_weights_mean_r/
1104         &        wmean_hflux, wmean_sflux, wmean_tau,
1105         &        wmean_atemp, wmean_aqh,
1106         &        wmean_precip, wmean_snowprecip, wmean_wind,
1107         &        wmean_swflux, wmean_swdown, wmean_lwflux, wmean_lwdown,
1108         &        wmean_evap, wmean_apressure, wmean_runoff
1109          _RL wmean_hflux
1110          _RL wmean_sflux
1111          _RL wmean_tau
1112          _RL wmean_atemp
1113          _RL wmean_aqh
1114          _RL wmean_precip
1115          _RL wmean_swflux
1116          _RL wmean_swdown
1117          _RL wmean_snowprecip
1118          _RL wmean_lwflux
1119          _RL wmean_lwdown
1120          _RL wmean_evap
1121          _RL wmean_apressure
1122          _RL wmean_runoff
1123          _RL wmean_wind
1124                                                    
1125        common /ecco_cost_weights_2_r/        common /ecco_cost_weights_2_r/
1126       &                      whflux2,wsflux2,wtauu2,wtauv2       &                      whflux2,wsflux2,wtauu2,wtauv2
1127        _RL whflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL whflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1128        _RL wsflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsflux2 (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1129        _RL wtauu2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtauu2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1130        _RL wtauv2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtauv2  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1131                                                                                                    
1132                                          
1133    
1134  #if (defined (ALLOW_OBCSN_COST_CONTRIBUTION) || \  #if (defined (ALLOW_OBCSN_COST_CONTRIBUTION) || \
1135       defined (ALLOW_OBCSN_CONTROL))       defined (ALLOW_OBCSN_CONTROL))
# Line 797  c     scatydat   - reference meridional Line 1170  c     scatydat   - reference meridional
1170  c     sstdat     - reference sea surface temperature data.  c     sstdat     - reference sea surface temperature data.
1171  c     tmidat     - reference TMI sea surface temperature data.  c     tmidat     - reference TMI sea surface temperature data.
1172  c     sssdat     - reference sea surface temperature data.  c     sssdat     - reference sea surface temperature data.
1173    c     bpdat      - bottom pressure from time-varying GRACE.
1174  c     tauxmask   - mask for reference wind stress data.  c     tauxmask   - mask for reference wind stress data.
1175  c     tauymask   - mask for reference wind stress data.  c     tauymask   - mask for reference wind stress data.
1176  c     scatxmask  - mask for scat wind stress data.  c     scatxmask  - mask for scat wind stress data.
1177  c     scatymask  - mask for scat wind stress data.  c     scatymask  - mask for scat wind stress data.
1178  c     sstmask    - mask for reference sea surface temperature data.  c     sstmask    - mask for reference sea surface temperature data.
1179  c     tmimask    - mask for reference sea surface temperature data.  c     tmimask    - mask for reference sea surface temperature data.
1180  c     sssmask    - mask for reference sea surface temperature data.  c     sssmask    - mask for reference sea surface temperature data.
# Line 813  c     ersobs     - ERS data. Line 1187  c     ersobs     - ERS data.
1187  c     ersmask    - mask for ERS data.  c     ersmask    - mask for ERS data.
1188  c     ctdtobs    - CTD temperature data  c     ctdtobs    - CTD temperature data
1189  c     ctdsobs    - CTD salinity data  c     ctdsobs    - CTD salinity data
1190  c     xbtobs     - XBT data  c     xbtobs     - XBT data
1191  c     argot      - ARGO  temperature data  c     argot      - ARGO  temperature data
1192  c     argos      - ARGO  salt data  c     argos      - ARGO  salt data
1193  c     udriftdat  - drifters zonal velocities  c     udriftdat  - drifters zonal velocities
1194  c     vdriftdat  - drifters meridional velocities  c     vdriftdat  - drifters meridional velocities
1195    
# Line 826  c     vdriftdat  - drifters meridional v Line 1200  c     vdriftdat  - drifters meridional v
1200       &                     sstdat,       &                     sstdat,
1201       &                     tmidat,       &                     tmidat,
1202       &                     sssdat,       &                     sssdat,
1203         &                     bpdat,
1204       &                     sstmask,       &                     sstmask,
1205       &                     tmimask,       &                     tmimask,
1206       &                     sssmask,       &                     sssmask,
1207         &                     bpmask,
1208       &                     tauxmask,       &                     tauxmask,
1209       &                     tauymask,       &                     tauymask,
1210       &                     scatxmask,       &                     scatxmask,
# Line 840  c     vdriftdat  - drifters meridional v Line 1216  c     vdriftdat  - drifters meridional v
1216       &                     tpmask,       &                     tpmask,
1217       &                     ersobs,       &                     ersobs,
1218       &                     ersmask,       &                     ersmask,
1219         &                     gfoobs,
1220         &                     gfomask,
1221       &                     ctdtobs,       &                     ctdtobs,
1222       &                     ctdsobs,       &                     ctdsobs,
1223       &                     xbtobs,       &                     xbtobs,
# Line 849  c     vdriftdat  - drifters meridional v Line 1227  c     vdriftdat  - drifters meridional v
1227       &                     vdriftdat,       &                     vdriftdat,
1228       &                     curmtruobs,       &                     curmtruobs,
1229       &                     curmtrvobs       &                     curmtrvobs
1230        
1231        _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)
1232        _RL scatxdat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatxdat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1233        _RL scatydat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL scatydat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1234        _RL sstdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sstdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1235        _RL tmidat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tmidat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1236        _RL sssdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sssdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1237          _RL bpdat     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1238        _RL tauxmask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tauxmask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1239        _RL tauymask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tauymask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1240        _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 863  c     vdriftdat  - drifters meridional v Line 1242  c     vdriftdat  - drifters meridional v
1242        _RL sstmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sstmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1243        _RL tmimask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tmimask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1244        _RL sssmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sssmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1245          _RL bpmask    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1246        _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)
1247        _RL tpmean    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpmean    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1248        _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 870  c     vdriftdat  - drifters meridional v Line 1250  c     vdriftdat  - drifters meridional v
1250        _RL tpmask    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpmask    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1251        _RL ersobs    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL ersobs    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1252        _RL ersmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL ersmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1253          _RL gfoobs    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1254          _RL gfomask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
1255        _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)
1256        _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)
1257        _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 895  c     topexfile     - reference data fil Line 1277  c     topexfile     - reference data fil
1277  c                     (TOPEX/POSEIDON).  c                     (TOPEX/POSEIDON).
1278  c     ersfile       - reference data file for sea surface height data  c     ersfile       - reference data file for sea surface height data
1279  c                     (ERS).  c                     (ERS).
1280  c ctdtfile, ctdsfile- reference data file for temperature and salinity  c ctdtfile, ctdsfile- reference data file for temperature and salinity
1281  c                     from CTD  c                     from CTD
1282  c ctdtclimfile, ctdsclimfile- reference data file for temperature  c ctdtclimfile, ctdsclimfile- reference data file for temperature
1283  c                      and salinity from CTD with out timetag  c                      and salinity from CTD with out timetag
# Line 912  c     driftfile     - reference data fil Line 1294  c     driftfile     - reference data fil
1294       &                     sstdatfile,       &                     sstdatfile,
1295       &                     tmidatfile,       &                     tmidatfile,
1296       &                     sssdatfile,       &                     sssdatfile,
1297         &                     bpdatfile,
1298       &                     topexmeanfile,       &                     topexmeanfile,
1299       &                     topexfile,       &                     topexfile,
1300       &                     ersfile,       &                     ersfile,
1301         &                     gfofile,
1302       &                     ctdtfile,       &                     ctdtfile,
1303       &                     ctdsfile,       &                     ctdsfile,
1304       &                     ctdtclimfile,       &                     ctdtclimfile,
# Line 922  c     driftfile     - reference data fil Line 1306  c     driftfile     - reference data fil
1306       &                     xbtfile,       &                     xbtfile,
1307       &                     argotfile,       &                     argotfile,
1308       &                     argosfile,       &                     argosfile,
1309       &                     udriftfile,       &                     udriftfile,
1310       &                     vdriftfile,       &                     vdriftfile,
1311         &                     gencost_datafile,
1312       &                     curmtrufile,       &                     curmtrufile,
1313       &                     curmtrvfile       &                     curmtrvfile
1314    
# Line 934  c     driftfile     - reference data fil Line 1319  c     driftfile     - reference data fil
1319        character*(MAX_LEN_FNAM) sstdatfile        character*(MAX_LEN_FNAM) sstdatfile
1320        character*(MAX_LEN_FNAM) tmidatfile        character*(MAX_LEN_FNAM) tmidatfile
1321        character*(MAX_LEN_FNAM) sssdatfile        character*(MAX_LEN_FNAM) sssdatfile
1322          character*(MAX_LEN_FNAM) bpdatfile
1323        character*(MAX_LEN_FNAM) topexmeanfile        character*(MAX_LEN_FNAM) topexmeanfile
1324        character*(MAX_LEN_FNAM) topexfile        character*(MAX_LEN_FNAM) topexfile
1325        character*(MAX_LEN_FNAM) ersfile        character*(MAX_LEN_FNAM) ersfile
1326          character*(MAX_LEN_FNAM) gfofile
1327        character*(MAX_LEN_FNAM) ctdtfile        character*(MAX_LEN_FNAM) ctdtfile
1328        character*(MAX_LEN_FNAM) ctdsfile        character*(MAX_LEN_FNAM) ctdsfile
1329        character*(MAX_LEN_FNAM) ctdtclimfile        character*(MAX_LEN_FNAM) ctdtclimfile
# Line 945  c     driftfile     - reference data fil Line 1332  c     driftfile     - reference data fil
1332        character*(MAX_LEN_FNAM) argotfile        character*(MAX_LEN_FNAM) argotfile
1333        character*(MAX_LEN_FNAM) argosfile        character*(MAX_LEN_FNAM) argosfile
1334        character*(MAX_LEN_FNAM) argofile        character*(MAX_LEN_FNAM) argofile
1335          character*(MAX_LEN_FNAM) gencost_datafile(NGENCOST)
1336        character*(MAX_LEN_FNAM) udriftfile        character*(MAX_LEN_FNAM) udriftfile
1337        character*(MAX_LEN_FNAM) vdriftfile              character*(MAX_LEN_FNAM) vdriftfile
1338        character*(MAX_LEN_FNAM) curmtrufile        character*(MAX_LEN_FNAM) curmtrufile
1339        character*(MAX_LEN_FNAM) curmtrvfile        character*(MAX_LEN_FNAM) curmtrvfile
1340    
# Line 958  c     using_ers - flag that indicates th Line 1346  c     using_ers - flag that indicates th
1346    
1347        common /ecco_cost_data_flags/        common /ecco_cost_data_flags/
1348       &                         using_topex,       &                         using_topex,
1349       &                         using_ers       &                         using_ers,
1350         &                         using_gfo
1351        logical using_topex        logical using_topex
1352        logical using_ers        logical using_ers
1353          logical using_gfo
1354    
1355  c     Calendar information for the observations:  c     Calendar information for the observations:
1356  c     ==========================================  c     ==========================================
# Line 979  c     sshperiod      - sampling interval Line 1369  c     sshperiod      - sampling interval
1369       &                           argosstartdate,       &                           argosstartdate,
1370       &                           tmistartdate,       &                           tmistartdate,
1371       &                           sssstartdate,       &                           sssstartdate,
1372         &                           bpstartdate,
1373       &                           topexstartdate,       &                           topexstartdate,
1374       &                           ersstartdate       &                           ersstartdate,
1375         &                           gfostartdate
1376        integer scatxstartdate(4)        integer scatxstartdate(4)
1377        integer scatystartdate(4)        integer scatystartdate(4)
1378        integer sststartdate(4)        integer sststartdate(4)
# Line 988  c     sshperiod      - sampling interval Line 1380  c     sshperiod      - sampling interval
1380        integer argosstartdate(4)        integer argosstartdate(4)
1381        integer tmistartdate(4)        integer tmistartdate(4)
1382        integer sssstartdate(4)        integer sssstartdate(4)
1383          integer bpstartdate(4)
1384        integer topexstartdate(4)        integer topexstartdate(4)
1385        integer ersstartdate(4)        integer ersstartdate(4)
1386          integer gfostartdate(4)
1387    
1388        common /ecco_cost_data_aux_i/        common /ecco_cost_data_aux_i/
1389       &                           tmistartdate1,       &                           tmistartdate1,
# Line 998  c     sshperiod      - sampling interval Line 1392  c     sshperiod      - sampling interval
1392       &                           sststartdate2,       &                           sststartdate2,
1393       &                           sssstartdate1,       &                           sssstartdate1,
1394       &                           sssstartdate2,       &                           sssstartdate2,
1395         &                           bpstartdate1,
1396         &                           bpstartdate2,
1397       &                           argotstartdate1,       &                           argotstartdate1,
1398       &                           argotstartdate2,       &                           argotstartdate2,
1399       &                           argosstartdate1,       &                           argosstartdate1,
# Line 1006  c     sshperiod      - sampling interval Line 1402  c     sshperiod      - sampling interval
1402       &                           topexstartdate2,       &                           topexstartdate2,
1403       &                           ersstartdate1,       &                           ersstartdate1,
1404       &                           ersstartdate2,       &                           ersstartdate2,
1405         &                           gfostartdate1,
1406         &                           gfostartdate2,
1407       &                           scatstartdate1,       &                           scatstartdate1,
1408       &                           scatstartdate2,       &                           scatstartdate2
      &                           costIceStart1,  
      &                           costIceStart2,  
      &                           costIceEnd1,  
      &                           costIceEnd2  
1409    
1410        integer tmistartdate1        integer tmistartdate1
1411        integer tmistartdate2        integer tmistartdate2
# Line 1019  c     sshperiod      - sampling interval Line 1413  c     sshperiod      - sampling interval
1413        integer sststartdate2        integer sststartdate2
1414        integer sssstartdate1        integer sssstartdate1
1415        integer sssstartdate2        integer sssstartdate2
1416          integer bpstartdate1
1417          integer bpstartdate2
1418        integer argotstartdate1        integer argotstartdate1
1419        integer argotstartdate2        integer argotstartdate2
1420        integer argosstartdate1        integer argosstartdate1
# Line 1027  c     sshperiod      - sampling interval Line 1423  c     sshperiod      - sampling interval
1423        integer topexstartdate2        integer topexstartdate2
1424        integer ersstartdate1        integer ersstartdate1
1425        integer ersstartdate2        integer ersstartdate2
1426          integer gfostartdate1
1427          integer gfostartdate2
1428        integer scatstartdate1        integer scatstartdate1
1429        integer scatstartdate2        integer scatstartdate2
       integer costIceStart1  
       integer costIceStart2  
       integer costIceEnd1  
       integer costIceEnd2  
1430    
1431        common /ecco_cost_data_times_r/        common /ecco_cost_data_times_r/
1432       &                           topexperiod,       &                           topexperiod,
1433       &                           ersperiod,       &                           ersperiod,
1434       &                           scatperiod,       &                           gfoperiod,
1435       &                           costIceStart,       &                           scatperiod
      &                           costIceEnd  
1436        _RL topexperiod        _RL topexperiod
1437        _RL ersperiod        _RL ersperiod
1438          _RL gfoperiod
1439        _RL scatperiod        _RL scatperiod
       _RL costIceStart  
       _RL costIceEnd  
   
 c     cost_ice_flag  - cost_ice flag (see cost_ice.F)  
1440    
1441        common /ecco_cost_ice_i/ cost_ice_flag        common /ecco_cost_data_detrend/
1442        integer cost_ice_flag       &                           topexintercept,
1443         &                           ersintercept,
1444         &                           gfointercept,
1445         &                           topexslope,
1446         &                           ersslope,
1447         &                           gfoslope
1448          _RL topexintercept
1449          _RL ersintercept
1450          _RL gfointercept
1451          _RL topexslope
1452          _RL ersslope
1453          _RL gfoslope
1454    
1455  c     ==================================================================  c     ==================================================================
1456  c     END OF HEADER COST  c     END OF HEADER COST
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22