/[MITgcm]/MITgcm/pkg/cost/cost.h
ViewVC logotype

Diff of /MITgcm/pkg/cost/cost.h

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

revision 1.2.6.2 by heimbach, Thu Mar 21 03:55:16 2002 UTC revision 1.2.6.3 by heimbach, Thu Apr 4 10:58:58 2002 UTC
# Line 49  c             intantaneous surface press Line 49  c             intantaneous surface press
49  c     ubar  - contains the averaged zonal velocity component for the  c     ubar  - contains the averaged zonal velocity component for the
50  c             whole integration period. Before, it accumulates the  c             whole integration period. Before, it accumulates the
51  c             intantaneous field.  c             intantaneous field.
52    c     vbar  - contains the averaged zonal velocity component for the
53    c             whole integration period. Before, it accumulates the
54    c             intantaneous field.
55    c     tauxbar  - contains the averaged zonal velocity component for the
56    c             whole integration period. Before, it accumulates the
57    c             intantaneous field.
58    c     tauybar  - contains the averaged zonal velocity component for the
59    c             whole integration period. Before, it accumulates the
60    c             intantaneous field.
61    c     hfluxbar  - contains the averaged zonal velocity component for the
62    c             whole integration period. Before, it accumulates the
63    c             intantaneous field.
64    c     sfluxbar  - contains the averaged zonal velocity component for the
65    c             whole integration period. Before, it accumulates the
66    c             intantaneous field.
67    
68        common /averages_r/        common /averages_r/
69       &                    tbar,       &                    tbar,
70       &                    sbar,       &                    sbar,
71       &                    psbar,       &                    psbar,
72       &                    ubar,       &                    ubar,
73       &                    vbar       &                    vbar,
74         &                    wbar,
75         &                    tauxbar,
76         &                    tauybar,
77         &                    hfluxbar,
78         &                    sfluxbar,
79         &                    Slmean,
80         &                    Tlmean,
81         &                    wlmean,
82         &                    Sfmean,
83         &                    Tfmean,                  
84         &                    wfmean
85    
86  #if (defined (ALLOW_THETA_COST_CONTRIBUTION) || \  #if (defined (ALLOW_THETA_COST_CONTRIBUTION) || \
      defined (ALLOW_SST_COST_CONTRIBUTION) || \  
87       defined (ALLOW_CTDT_COST_CONTRIBUTION) || \       defined (ALLOW_CTDT_COST_CONTRIBUTION) || \
88       defined (ALLOW_XBT_COST_CONTRIBUTION) || \       defined (ALLOW_XBT_COST_CONTRIBUTION) || \
89         defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \
90       defined (ALLOW_OBCS_COST_CONTRIBUTION))       defined (ALLOW_OBCS_COST_CONTRIBUTION))
91        _RL tbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL tbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
92  #else  #else
93    #ifdef ALLOW_SST_COST_CONTRIBUTION
94          _RL tbar  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
95    #else
96        _RL tbar        _RL tbar
97  #endif  #endif
98    #endif
99    
100  #if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \  #if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \
      defined (ALLOW_SSS_COST_CONTRIBUTION) || \  
101       defined (ALLOW_CTDS_COST_CONTRIBUTION) || \       defined (ALLOW_CTDS_COST_CONTRIBUTION) || \
102         defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \
103       defined (ALLOW_OBCS_COST_CONTRIBUTION))       defined (ALLOW_OBCS_COST_CONTRIBUTION))
104        _RL sbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL sbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
105  #else  #else
106    #ifdef ALLOW_SSS_COST_CONTRIBUTION
107          _RL sbar  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
108    #else
109        _RL sbar        _RL sbar
110  #endif  #endif
111    #endif
112    
113  #ifdef ALLOW_SSH_COST_CONTRIBUTION  #ifdef ALLOW_SSH_COST_CONTRIBUTION
114        _RL psbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)        _RL psbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
# Line 91  c             intantaneous field. Line 125  c             intantaneous field.
125        _RL vbar        _RL vbar
126  #endif  #endif
127    
128    #ifdef ALLOW_DRIFTW_COST_CONTRIBUTION
129          _RL wbar  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
130    #else
131          _RL wbar
132    #endif
133    
134    #ifdef ALLOW_DRIFT_COST_CONTRIBUTION
135          _RL   Tlmean(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
136          _RL   Slmean(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
137          _RL   Tfmean(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
138          _RL   Sfmean(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
139    #else
140          _RL   Tlmean
141          _RL   Slmean
142          _RL   Tfmean
143          _RL   Sfmean        
144    #endif
145    
146    #ifdef ALLOW_DRIFTW_COST_CONTRIBUTION
147          _RL   wlmean(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
148          _RL   wfmean(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
149    #else
150          _RL   wlmean
151          _RL   wfmean
152    #endif
153    
154    #ifdef ALLOW_SCAT_COST_CONTRIBUTION
155          _RL tauxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
156          _RL tauybar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
157    #else
158          _RL tauxbar
159          _RL tauybar
160    #endif
161    
162    #ifdef ALLOW_MEAN_HFLUX_COST_CONTRIBUTION
163          _RL hfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)
164    #else
165          _RL hfluxbar
166    #endif
167    
168    #ifdef ALLOW_MEAN_SFLUX_COST_CONTRIBUTION
169          _RL sfluxbar  (1-olx:snx+olx,1-oly:sny+oly,  nsx,nsy)    
170    #else
171          _RL sfluxbar
172    #endif
173    
174    
175        common /averages_c/        common /averages_c/
# Line 98  c             intantaneous field. Line 177  c             intantaneous field.
177       &                    sbarfile,       &                    sbarfile,
178       &                    psbarfile,       &                    psbarfile,
179       &                    ubarfile,       &                    ubarfile,
180       &                    vbarfile       &                    vbarfile,
181         &                    wbarfile,
182         &                    tauxbarfile,
183         &                    tauybarfile,
184         &                    hfluxbarfile,
185         &                    sfluxbarfile
186        character*(MAX_LEN_FNAM) tbarfile        character*(MAX_LEN_FNAM) tbarfile
187        character*(MAX_LEN_FNAM) sbarfile        character*(MAX_LEN_FNAM) sbarfile
188        character*(MAX_LEN_FNAM) psbarfile        character*(MAX_LEN_FNAM) psbarfile
189        character*(MAX_LEN_FNAM) ubarfile        character*(MAX_LEN_FNAM) ubarfile
190        character*(MAX_LEN_FNAM) vbarfile        character*(MAX_LEN_FNAM) vbarfile
191          character*(MAX_LEN_FNAM) wbarfile
192          character*(MAX_LEN_FNAM) tauxbarfile
193          character*(MAX_LEN_FNAM) tauybarfile
194          character*(MAX_LEN_FNAM) hfluxbarfile
195          character*(MAX_LEN_FNAM) sfluxbarfile
196    
197  c     file precision and field type  c     file precision and field type
198    
# Line 153  c     objf_hflux    - Heat flux contribu Line 242  c     objf_hflux    - Heat flux contribu
242  c     objf_sflux    - Salt flux contribution.  c     objf_sflux    - Salt flux contribution.
243  c     objf_tauu  - Zonal wind stress contribution.  c     objf_tauu  - Zonal wind stress contribution.
244  c     objf_tauv  - Meridional wind stress contribution.  c     objf_tauv  - Meridional wind stress contribution.
245    c     objf_hfluxm    - time-mean Heat flux contribution.
246    c     objf_sfluxm    - time-mean Salt flux contribution.
247    c     objf_tauum  - time-mean Zonal wind stress contribution.
248    c     objf_tauvm  - time-mean Meridional wind stress contribution.
249    c     objf_hfluxmm    - Global time-mean Heat flux contribution.
250    c     objf_sfluxmm    - Global time-mean Salt flux contribution.
251  c     objf_hmean - Mean sea surface height contribution.  c     objf_hmean - Mean sea surface height contribution.
252  c     objf_h     - Residual sea surface height contribution.  c     objf_h     - Residual sea surface height contribution.
253  c     objf_temp  - Temperature contribution.  c     objf_temp  - Temperature contribution.
254  c     objf_salt  - Salinity contribution.  c     objf_salt  - Salinity contribution.
255    c     objf_temp0 - Initial conditions Temperature contribution.
256    c     objf_salt0 - Initial conditions Salinity contribution.
257  c     objf_sst   - Sea surface temperature contribution.  c     objf_sst   - Sea surface temperature contribution.
258  c     objf_sss   - Sea surface salinity contribution.  c     objf_sss   - Sea surface salinity contribution.
259  c     objf_atl   - Meridional heat transport in the N-Atlantic  c     objf_atl   - Meridional heat transport in the N-Atlantic
260  c     objf_ctdt  - Temperature measurements from Woce CTD  c     objf_ctdt  - Temperature measurements from Woce CTD
261  c     objf_ctds  - Salinity measurements from Woce CTD  c     objf_ctds  - Salinity measurements from Woce CTD
262  c     objf_xbt   - XBT temperature data  c     objf_xbt   - XBT temperature data
263    c     objf_argot - ARGO temperature profiles
264    c     objf_argos - ARGO salt profiles
265    c     objf_scatxm - time-mean zonal SCAT  contribution
266    c     objf_scatym - time-mean meridional SCAT  contribution
267    c     objf_scatx  - zonal SCAT  contribution
268    c     objf_scaty  - meridional SCAT  contribution
269  c  c
270  c     mult_"var" - multipliers for the individual cost  c     mult_"var" - multipliers for the individual cost
271  c                  function contributions.  c                  function contributions.
# Line 170  c                  function contribution Line 273  c                  function contribution
273        common /cost_r/        common /cost_r/
274       &                fc,       &                fc,
275       &                objf_hflux,       &                objf_hflux,
276         &                objf_hfluxm,
277         &                objf_hfluxmm,
278       &                objf_sflux,       &                objf_sflux,
279         &                objf_sfluxm,
280         &                objf_sfluxmm,
281       &                objf_tauu,       &                objf_tauu,
282         &                objf_tauum,
283       &                objf_tauv,       &                objf_tauv,
284         &                objf_tauvm,
285       &                objf_hmean,       &                objf_hmean,
286       &                objf_h,       &                objf_h,
287       &                objf_temp,       &                objf_temp,
288       &                objf_salt,       &                objf_salt,
289         &                objf_temp0,
290         &                objf_salt0,
291       &                objf_sst,       &                objf_sst,
292       &                objf_sss,       &                objf_sss,
293       &                objf_atl,       &                objf_atl,
294       &                objf_ctdt,       &                objf_ctdt,
295       &                objf_ctds,       &                objf_ctds,
296       &                objf_xbt,       &                objf_xbt,
297         &                objf_argot,
298         &                objf_argos,
299       &                objf_drift,       &                objf_drift,
300         &                objf_tdrift,
301         &                objf_sdrift,
302         &                objf_wdrift,
303         &                objf_scatx,
304         &                objf_scaty,
305         &                objf_scatxm,
306         &                objf_scatym,
307       &                objf_atemp,       &                objf_atemp,
308       &                objf_aqh,       &                objf_aqh,
309       &                objf_uwind,       &                objf_uwind,
# Line 193  c                  function contribution Line 313  c                  function contribution
313       &                objf_obcsw,       &                objf_obcsw,
314       &                objf_obcse       &                objf_obcse
315        _RL  fc        _RL  fc
316        _RL  objf_hflux(nsx,nsy)        _RL  objf_hflux  (nsx,nsy)
317        _RL  objf_sflux(nsx,nsy)        _RL  objf_hfluxm (nsx,nsy)
318        _RL  objf_tauu (nsx,nsy)        _RL  objf_hfluxmm(nsx,nsy)
319        _RL  objf_tauv (nsx,nsy)        _RL  objf_sflux  (nsx,nsy)
320          _RL  objf_sfluxm (nsx,nsy)
321          _RL  objf_sfluxmm(nsx,nsy)
322          _RL  objf_tauu   (nsx,nsy)
323          _RL  objf_tauum  (nsx,nsy)
324          _RL  objf_tauv   (nsx,nsy)
325          _RL  objf_tauvm  (nsx,nsy)
326        _RL  objf_hmean        _RL  objf_hmean
327        _RL  objf_h    (nsx,nsy)        _RL  objf_h    (nsx,nsy)
328        _RL  objf_temp (nsx,nsy)        _RL  objf_temp (nsx,nsy)
329        _RL  objf_salt (nsx,nsy)        _RL  objf_salt (nsx,nsy)
330          _RL  objf_temp0(nsx,nsy)
331          _RL  objf_salt0(nsx,nsy)
332        _RL  objf_sst  (nsx,nsy)        _RL  objf_sst  (nsx,nsy)
333        _RL  objf_sss  (nsx,nsy)        _RL  objf_sss  (nsx,nsy)
334        _RL  objf_atl  (nsx,nsy)        _RL  objf_atl  (nsx,nsy)
335        _RL  objf_ctdt (nsx,nsy)        _RL  objf_ctdt (nsx,nsy)
336        _RL  objf_ctds (nsx,nsy)        _RL  objf_ctds (nsx,nsy)
337        _RL  objf_xbt  (nsx,nsy)        _RL  objf_xbt  (nsx,nsy)
338          _RL  objf_argot(nsx,nsy)
339          _RL  objf_argos(nsx,nsy)
340        _RL  objf_drift(nsx,nsy)        _RL  objf_drift(nsx,nsy)
341          _RL  objf_tdrift(nsx,nsy)
342          _RL  objf_sdrift(nsx,nsy)
343          _RL  objf_wdrift(nsx,nsy)
344          _RL  objf_scatx(nsx,nsy)
345          _RL  objf_scaty(nsx,nsy)
346          _RL  objf_scatxm(nsx,nsy)
347          _RL  objf_scatym(nsx,nsy)
348        _RL  objf_atemp(nsx,nsy)        _RL  objf_atemp(nsx,nsy)
349        _RL  objf_aqh  (nsx,nsy)        _RL  objf_aqh  (nsx,nsy)
350        _RL  objf_uwind(nsx,nsy)        _RL  objf_uwind(nsx,nsy)
# Line 226  c                  function contribution Line 363  c                  function contribution
363       &                    mult_h,       &                    mult_h,
364       &                    mult_temp,       &                    mult_temp,
365       &                    mult_salt,       &                    mult_salt,
366         &                    mult_temp0,
367         &                    mult_salt0,
368       &                    mult_sst,       &                    mult_sst,
369       &                    mult_sss,       &                    mult_sss,
370       &                    mult_atl,       &                    mult_atl,
371       &                    mult_ctdt,       &                    mult_ctdt,
372       &                    mult_ctds,       &                    mult_ctds,
373       &                    mult_xbt,       &                    mult_xbt,
374         &                    mult_argot,
375         &                    mult_argos,
376       &                    mult_drift,       &                    mult_drift,
377         &                    mult_tdrift,
378         &                    mult_sdrift,
379         &                    mult_wdrift,
380         &                    mult_scatx,
381         &                    mult_scaty,
382       &                    mult_atemp,       &                    mult_atemp,
383       &                    mult_aqh,       &                    mult_aqh,
384       &                    mult_uwind,       &                    mult_uwind,
# Line 249  c                  function contribution Line 395  c                  function contribution
395        _RL  mult_h        _RL  mult_h
396        _RL  mult_temp        _RL  mult_temp
397        _RL  mult_salt        _RL  mult_salt
398          _RL  mult_temp0
399          _RL  mult_salt0
400        _RL  mult_sst        _RL  mult_sst
401        _RL  mult_sss        _RL  mult_sss
402        _RL  mult_atl        _RL  mult_atl
403        _RL  mult_ctdt        _RL  mult_ctdt
404        _RL  mult_ctds        _RL  mult_ctds
405        _RL  mult_xbt        _RL  mult_xbt
406          _RL  mult_argot
407          _RL  mult_argos
408        _RL  mult_drift        _RL  mult_drift
409          _RL  mult_tdrift
410          _RL  mult_sdrift
411          _RL  mult_wdrift
412          _RL  mult_scatx
413          _RL  mult_scaty
414        _RL  mult_atemp        _RL  mult_atemp
415        _RL  mult_aqh        _RL  mult_aqh
416        _RL  mult_uwind        _RL  mult_uwind
# Line 287  c Line 442  c
442  c     hflux_errfile         - heat flux error.  c     hflux_errfile         - heat flux error.
443  c     sflux_errfile         - salt flux error.  c     sflux_errfile         - salt flux error.
444  c     tauu_errfile          - zonal wind stress error.  c     tauu_errfile          - zonal wind stress error.
445    c     tauum_errfile         - zonal wind stress error.
446  c     tauv_errfile          - meridional wind stress error.  c     tauv_errfile          - meridional wind stress error.
447    c     tauvm_errfile         - meridional wind stress error.
448    c     tscatx_errfile        - zonal wind stress error.
449    c     tscaty_errfile        - meridional wind stress error.
450  c     data_errfile          - weights for theta, salt, and SST  c     data_errfile          - weights for theta, salt, and SST
451  c     geoid_errfile         - geoid error.  c     geoid_errfile         - geoid error.
452  c     geoid_covariancefile  - geoid error covariance.  c     geoid_covariancefile  - geoid error covariance.
# Line 295  c     ssh_errfile           - sea surfac Line 454  c     ssh_errfile           - sea surfac
454  c     ctdt_errfile          - CTD temperature error.  c     ctdt_errfile          - CTD temperature error.
455  c     ctds_errfile          - CTD salinity error.  c     ctds_errfile          - CTD salinity error.
456  c     drift_errfile         - drifter error.  c     drift_errfile         - drifter error.
457    c     salterrfile           - representation error due unresolved eddies
458    c     temperrfile           - representation error due unresolved eddies
459    
460        common /cost_c/        common /cost_c/
461       &                hflux_errfile,       &                hflux_errfile,
462         &                hfluxm_errfile,
463       &                sflux_errfile,       &                sflux_errfile,
464         &                sfluxm_errfile,
465       &                tauu_errfile,       &                tauu_errfile,
466         &                tauum_errfile,
467       &                tauv_errfile,       &                tauv_errfile,
468         &                tauvm_errfile,
469         &                scatx_errfile,
470         &                scaty_errfile,
471       &                data_errfile,       &                data_errfile,
472       &                geoid_errfile,       &                geoid_errfile,
473       &                geoid_covariancefile,       &                geoid_covariancefile,
# Line 308  c     drift_errfile         - drifter er Line 475  c     drift_errfile         - drifter er
475       &                ctdt_errfile,       &                ctdt_errfile,
476       &                ctds_errfile,       &                ctds_errfile,
477       &                drift_errfile,       &                drift_errfile,
478         &                udrifterrfile,
479         &                vdrifterrfile,
480         &                salterrfile,
481         &                temperrfile,
482       &                atemp_errfile,       &                atemp_errfile,
483       &                aqh_errfile,       &                aqh_errfile,
484       &                uwind_errfile,       &                uwind_errfile,
# Line 316  c     drift_errfile         - drifter er Line 487  c     drift_errfile         - drifter er
487        character*(MAX_LEN_FNAM) sflux_errfile        character*(MAX_LEN_FNAM) sflux_errfile
488        character*(MAX_LEN_FNAM) tauu_errfile        character*(MAX_LEN_FNAM) tauu_errfile
489        character*(MAX_LEN_FNAM) tauv_errfile        character*(MAX_LEN_FNAM) tauv_errfile
490          character*(MAX_LEN_FNAM) hfluxm_errfile
491          character*(MAX_LEN_FNAM) sfluxm_errfile
492          character*(MAX_LEN_FNAM) tauum_errfile
493          character*(MAX_LEN_FNAM) tauvm_errfile
494          character*(MAX_LEN_FNAM) scatx_errfile
495          character*(MAX_LEN_FNAM) scaty_errfile
496        character*(MAX_LEN_FNAM) data_errfile        character*(MAX_LEN_FNAM) data_errfile
497        character*(MAX_LEN_FNAM) geoid_errfile        character*(MAX_LEN_FNAM) geoid_errfile
498        character*(MAX_LEN_FNAM) geoid_covariancefile        character*(MAX_LEN_FNAM) geoid_covariancefile
# Line 323  c     drift_errfile         - drifter er Line 500  c     drift_errfile         - drifter er
500        character*(MAX_LEN_FNAM) ctdt_errfile        character*(MAX_LEN_FNAM) ctdt_errfile
501        character*(MAX_LEN_FNAM) ctds_errfile        character*(MAX_LEN_FNAM) ctds_errfile
502        character*(MAX_LEN_FNAM) drift_errfile        character*(MAX_LEN_FNAM) drift_errfile
503          character*(MAX_LEN_FNAM) udrifterrfile
504          character*(MAX_LEN_FNAM) vdrifterrfile      
505          character*(MAX_LEN_FNAM) salterrfile
506          character*(MAX_LEN_FNAM) temperrfile
507        character*(MAX_LEN_FNAM) atemp_errfile        character*(MAX_LEN_FNAM) atemp_errfile
508        character*(MAX_LEN_FNAM) aqh_errfile        character*(MAX_LEN_FNAM) aqh_errfile
509        character*(MAX_LEN_FNAM) uwind_errfile        character*(MAX_LEN_FNAM) uwind_errfile
# Line 333  c     Arrays where the weights are store Line 514  c     Arrays where the weights are store
514  c     ====================================  c     ====================================
515  c  c
516  c     cosphi     - cosine of latitude.  c     cosphi     - cosine of latitude.
517  c     whflux        - weight for heat flux.  c     whflux     - weight for heat flux.
518  c     wsflux        - weight for salt flux.  c     wsflux     - weight for salt flux.
519  c     wtauu        - weight for zonal wind stress.  c     wtauu      - weight for zonal wind stress.
520  c     wtauu        - weight for meridional wind stress.  c     wtauu      - weight for meridional wind stress.
521  c     wtheta         - weight for temperature.  c     wscatx     - weight for zonal scat stress.
522    c     wscaty     - weight for meridional scat stress.
523    c     wtheta     - weight for temperature.
524    c     wtheta2    - representation error due to unresolved eddies
525  c     wsst       - weight for sea surface temperature.  c     wsst       - weight for sea surface temperature.
526  c     wsss       - weight for sea surface salinity.  c     wsss       - weight for sea surface salinity.
527  c     wsalt         - weight for salinity.  c     wsalt      - weight for salinity.
528    c     wsalt2     - representation error due to unresolved eddies
529  c     wtp        - weight for TOPEX/POSEIDON data.  c     wtp        - weight for TOPEX/POSEIDON data.
530  c     wers       - weight for ERS data.  c     wers       - weight for ERS data.
531  c     wp         - weight for geoid.  c     wp         - weight for geoid.
# Line 353  c     wvdrift    - weight for mean merid Line 538  c     wvdrift    - weight for mean merid
538       &                      cosphi,       &                      cosphi,
539       &                      whflux,wsflux,wtauu,wtauv,       &                      whflux,wsflux,wtauu,wtauv,
540       &                      watemp,waqh,wuwind,wvwind,       &                      watemp,waqh,wuwind,wvwind,
541       &                      wtheta,wsst,wsss,       &                      wscatx,wscaty,
542       &                      wsalt,       &                      wtheta,wtheta2,wsst,wsss,
543         &                      wsalt,wsalt2,
544       &                      wtp,wers,       &                      wtp,wers,
545       &                      wp,       &                      wp,
546       &                      wctdt,wctds,       &                      wctdt,wctds,
547       &                      wudrift,wvdrift       &                      wudrift,wvdrift,
548         &                      whfluxmm,wsfluxmm
549        _RL frame   (1-olx:snx+olx,1-oly:sny+oly           )        _RL frame   (1-olx:snx+olx,1-oly:sny+oly           )
550        _RL cosphi  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL cosphi  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
551        _RL whflux  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL whflux  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
552          _RL whfluxm (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
553          _RL whfluxmm(1-olx:snx+olx,1-oly:sny+oly)
554        _RL wsflux  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsflux  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
555          _RL wsfluxm (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
556          _RL wsfluxmm(1-olx:snx+olx,1-oly:sny+oly)
557        _RL wtauu   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtauu   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
558        _RL wtauv   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtauv   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
559          _RL wtauum  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
560          _RL wtauvm  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
561          _RL wscatx  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
562          _RL wscaty  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
563        _RL watemp  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL watemp  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
564        _RL waqh    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL waqh    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
565        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wuwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
566        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wvwind  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
567        _RL wtheta  (                            nr,nsx,nsy)        _RL wtheta  (                            nr,nsx,nsy)
568        _RL wsalt   (                            nr,nsx,nsy)        _RL wsalt   (                            nr,nsx,nsy)
569          _RL wtheta2 (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
570          _RL wsalt2  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
571        _RL wsst    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsst    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
572        _RL wsss    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wsss    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
573        _RL wtp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL wtp     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
# Line 407  c     Arrays that contain observations f Line 604  c     Arrays that contain observations f
604  c     ===============================================================  c     ===============================================================
605  c  c
606  c     tdat       - reference temperature data.  c     tdat       - reference temperature data.
607    c     scatxdat   - reference zonal wind stress.
608    c     scatydat   - reference meridional wind stress.
609  c     sstdat     - reference sea surface temperature data.  c     sstdat     - reference sea surface temperature data.
610    c     sssdat     - reference sea surface temperature data.
611    c     tauxmask   - mask for reference wind stress data.
612    c     tauymask   - mask for reference wind stress data.
613    c     scatxmask  - mask for scat wind stress data.
614    c     scatymask  - mask for scat wind stress data.
615  c     sstmask    - mask for reference sea surface temperature data.  c     sstmask    - mask for reference sea surface temperature data.
616    c     sssmask    - mask for reference sea surface temperature data.
617  c     sdat       - reference salinity data.  c     sdat       - reference salinity data.
618  c     tpmean     - reference mean sea surface height data.  c     tpmean     - reference mean sea surface height data.
619  c     tpmeanmask - mask for reference mean sea surface height data.  c     tpmeanmask - mask for reference mean sea surface height data.
# Line 419  c     ersmask    - mask for ERS data. Line 624  c     ersmask    - mask for ERS data.
624  c     ctdtobs    - CTD temperature data  c     ctdtobs    - CTD temperature data
625  c     ctdsobs    - CTD salinity data  c     ctdsobs    - CTD salinity data
626  c     xbtobs     - XBT data  c     xbtobs     - XBT data
627    c     argot      - ARGO  temperature data
628    c     argos      - ARGO  salt data
629  c     udriftdat  - drifters zonal velocities  c     udriftdat  - drifters zonal velocities
630  c     vdriftdat  - drifters meridional velocities  c     vdriftdat  - drifters meridional velocities
631    
632        common /cost_data_r/        common /cost_data_r/
633       &                     tdat,       &                     tdat,
634         &                     scatxdat,
635         &                     scatydat,
636       &                     sstdat,       &                     sstdat,
637         &                     sssdat,
638       &                     sstmask,       &                     sstmask,
639         &                     sssmask,
640         &                     tauxmask,
641         &                     tauymask,
642         &                     scatxmask,
643         &                     scatymask,
644       &                     sdat,       &                     sdat,
      &                     sssdat,  
      &                     sssmask,  
645       &                     tpmean,       &                     tpmean,
646       &                     tpmeanmask,       &                     tpmeanmask,
647       &                     tpobs,       &                     tpobs,
# Line 438  c     vdriftdat  - drifters meridional v Line 651  c     vdriftdat  - drifters meridional v
651       &                     ctdtobs,       &                     ctdtobs,
652       &                     ctdsobs,       &                     ctdsobs,
653       &                     xbtobs,       &                     xbtobs,
654         &                     argotobs,
655         &                     argosobs,
656       &                     udriftdat,       &                     udriftdat,
657       &                     vdriftdat       &                     vdriftdat
658            
659        _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)
660          _RL scatxdat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
661          _RL scatydat  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
662        _RL sstdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sstdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
       _RL sstmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)  
       _RL sdat      (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
663        _RL sssdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sssdat    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
664          _RL tauxmask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
665          _RL tauymask  (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
666          _RL scatxmask (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
667          _RL scatymask (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
668          _RL sstmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
669        _RL sssmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL sssmask   (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
670          _RL sdat      (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
671        _RL tpmean    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpmean    (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
672        _RL tpmeanmask(1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpmeanmask(1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
673        _RL tpobs     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL tpobs     (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
# Line 456  c     vdriftdat  - drifters meridional v Line 677  c     vdriftdat  - drifters meridional v
677        _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)
678        _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)
679        _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)
680          _RL argotobs  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
681          _RL argosobs  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
682        _RL udriftdat (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL udriftdat (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
683        _RL vdriftdat (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)        _RL vdriftdat (1-olx:snx+olx,1-oly:sny+oly,   nsx,nsy)
684    
# Line 465  c     =============================== Line 688  c     ===============================
688  c  c
689  c     tdatfile      - reference data file for temperature.  c     tdatfile      - reference data file for temperature.
690  c     sdatfile      - reference data file for salinity.  c     sdatfile      - reference data file for salinity.
691    c     scatxdatfile  - reference data file for zonal wind stress.
692    c     scatydatfile  - reference data file for meridional wind stress.
693  c     sstdatfile    - reference data file for sea surface temperature.  c     sstdatfile    - reference data file for sea surface temperature.
694  c     topexmeanfile - reference data file for mean sea surface height.  c     topexmeanfile - reference data file for mean sea surface height.
695  c     topexfile     - reference data file for sea surface height data  c     topexfile     - reference data file for sea surface height data
# Line 474  c                     (ERS). Line 699  c                     (ERS).
699  c ctdtfile, ctdsfile- reference data file for temperature and salinity  c ctdtfile, ctdsfile- reference data file for temperature and salinity
700  c                     from CTD  c                     from CTD
701  c     xbtfile       - reference data file for xbt  c     xbtfile       - reference data file for xbt
702    c     ARGOtfile     - reference data file for ARGO
703    c     ARGOsfile     - reference data file for ARGO
704  c     driftfile     - reference data file for drifter's mean velocities  c     driftfile     - reference data file for drifter's mean velocities
705    
706        common /cost_data_c/        common /cost_data_c/
707       &                     tdatfile,       &                     tdatfile,
708       &                     sdatfile,       &                     sdatfile,
709         &                     scatxdatfile,
710         &                     scatydatfile,
711       &                     sstdatfile,       &                     sstdatfile,
712       &                     sssdatfile,       &                     sssdatfile,
713       &                     topexmeanfile,       &                     topexmeanfile,
# Line 487  c     driftfile     - reference data fil Line 716  c     driftfile     - reference data fil
716       &                     ctdtfile,       &                     ctdtfile,
717       &                     ctdsfile,       &                     ctdsfile,
718       &                     xbtfile,       &                     xbtfile,
719         &                     argotfile,
720         &                     argosfile,
721       &                     udriftfile,       &                     udriftfile,
722       &                     vdriftfile       &                     vdriftfile
723        character*(MAX_LEN_FNAM) tdatfile        character*(MAX_LEN_FNAM) tdatfile
724        character*(MAX_LEN_FNAM) sdatfile        character*(MAX_LEN_FNAM) sdatfile
725          character*(MAX_LEN_FNAM) scatxdatfile
726          character*(MAX_LEN_FNAM) scatydatfile
727        character*(MAX_LEN_FNAM) sstdatfile        character*(MAX_LEN_FNAM) sstdatfile
728        character*(MAX_LEN_FNAM) sssdatfile        character*(MAX_LEN_FNAM) sssdatfile
729        character*(MAX_LEN_FNAM) topexmeanfile        character*(MAX_LEN_FNAM) topexmeanfile
# Line 499  c     driftfile     - reference data fil Line 732  c     driftfile     - reference data fil
732        character*(MAX_LEN_FNAM) ctdtfile        character*(MAX_LEN_FNAM) ctdtfile
733        character*(MAX_LEN_FNAM) ctdsfile        character*(MAX_LEN_FNAM) ctdsfile
734        character*(MAX_LEN_FNAM) xbtfile        character*(MAX_LEN_FNAM) xbtfile
735          character*(MAX_LEN_FNAM) argotfile
736          character*(MAX_LEN_FNAM) argosfile
737          character*(MAX_LEN_FNAM) argofile
738        character*(MAX_LEN_FNAM) udriftfile        character*(MAX_LEN_FNAM) udriftfile
739        character*(MAX_LEN_FNAM) vdriftfile              character*(MAX_LEN_FNAM) vdriftfile      
740    
741    
742    c     Flags used in the model-data comparison:
743    c     ========================================
744    c
745    c     using_ers - flag that indicates the use of ERS data
746    
747          common /cost_data_flags/
748         &                         using_topex,
749         &                         using_ers
750          logical using_topex
751          logical using_ers
752    
753  c     Calendar information for the observations:  c     Calendar information for the observations:
754  c     ==========================================  c     ==========================================
755  c  c
# Line 512  c     ersstartdate   - start date of the Line 759  c     ersstartdate   - start date of the
759  c     sshperiod      - sampling interval for the sea surface height data.  c     sshperiod      - sampling interval for the sea surface height data.
760    
761        common /cost_data_times_i/        common /cost_data_times_i/
762         &                           scatxstartdate,
763         &                           scatystartdate,
764       &                           sststartdate,       &                           sststartdate,
765       &                           sssstartdate,       &                           sssstartdate,
766       &                           topexstartdate,       &                           topexstartdate,
767       &                           ersstartdate       &                           ersstartdate
768          integer scatxstartdate(4)
769          integer scatystartdate(4)
770        integer sststartdate(4)        integer sststartdate(4)
771        integer sssstartdate(4)        integer sssstartdate(4)
772        integer topexstartdate(4)        integer topexstartdate(4)
# Line 524  c     sshperiod      - sampling interval Line 775  c     sshperiod      - sampling interval
775    
776        common /cost_data_times_r/        common /cost_data_times_r/
777       &                           topexperiod,       &                           topexperiod,
778       &                           ersperiod       &                           ersperiod,
779         &                           scatperiod
780        _RL topexperiod        _RL topexperiod
781        _RL ersperiod        _RL ersperiod
782          _RL scatperiod
783    
784  c     ==================================================================  c     ==================================================================
785  c     END OF HEADER COST  c     END OF HEADER COST
Legend:
Removed from v.1.2.6.2  
changed lines
  Added in v.1.2.6.3
  ViewVC Help
Powered by ViewVC 1.1.22