/[MITgcm]/MITgcm/pkg/ecco/ecco_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/ecco/ecco_readparms.F

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

revision 1.15 by heimbach, Fri Oct 14 19:54:49 2005 UTC revision 1.45 by heimbach, Sat Feb 6 11:30:16 2010 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "COST_CPPOPTIONS.h"  #include "COST_CPPOPTIONS.h"
5    
# Line 34  c     == local variables == Line 35  c     == local variables ==
35  C     msgBuf      - Informational/error meesage buffer  C     msgBuf      - Informational/error meesage buffer
36  C     iUnit       - Work variable for IO unit number  C     iUnit       - Work variable for IO unit number
37        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
38        INTEGER iUnit        INTEGER k, iUnit, num_file, num_var
39    
40  c     == end of interface ==  c     == end of interface ==
41    
42  c--   Read the namelist input.  c--   Read the namelist input.
43        namelist /ecco_cost_nml/        namelist /ecco_cost_nml/
44       &                   data_errfile,       &                   data_errfile,
45       &                   tbarfile, sbarfile, psbarfile,       &                   tbarfile, sbarfile,sstbarfile,
46         &                   psbarfile, bpbarfile,
47       &                   ubarfile, vbarfile, wbarfile,       &                   ubarfile, vbarfile, wbarfile,
48       &                   tauxbarfile, tauybarfile,       &                   tauxbarfile, tauybarfile,
49       &                   hfluxbarfile, sfluxbarfile,       &                   hfluxmeanbarfile, sfluxmeanbarfile,
50       &                   using_topex,       &                   costTranspDataFile,  topexmeanfile,
51         &                   using_topex, using_ers, using_gfo,
52       &                   topexstartdate1, topexstartdate2,       &                   topexstartdate1, topexstartdate2,
53       &                   topexperiod, topexfile, topexmeanfile,       &                   topexperiod, topexfile, tpTimeMaskFile,
54       &                   using_ers,       &                   topexintercept, topexslope, ersTimeMaskFile,
55       &                   ersstartdate1, ersstartdate2,       &                   ersstartdate1, ersstartdate2, gfoTimeMaskFile,
56       &                   ersperiod, ersfile,       &                   ersperiod, ersfile,
57       &                   scatstartdate1, scatstartdate2,       &                   ersintercept, ersslope,
58         &                   gfostartdate1, gfostartdate2,
59         &                   gfoperiod, gfofile,
60         &                   gfointercept, gfoslope,
61         &                   scatstartdate1, scatstartdate2,
62       &                   scatperiod, scatxdatfile, scatydatfile,       &                   scatperiod, scatxdatfile, scatydatfile,
63       &                   scatx_errfile, scaty_errfile,       &                   scatx_errfile, scaty_errfile,
64       &                   ssh_errfile,       &                   ssh_errfile,
65         &                   sshv4cost_errfile, sshv4cost_errfactor,
66         &                   tp_errfile, ers_errfile, gfo_errfile,
67       &                   tmistartdate1, tmistartdate2, tmidatfile,       &                   tmistartdate1, tmistartdate2, tmidatfile,
68       &                   sststartdate1, sststartdate2, sstdatfile,       &                   sststartdate1, sststartdate2, sstdatfile,
69       &                   sssstartdate1, sssstartdate2, sssdatfile,       &                   sssstartdate1, sssstartdate2, sssdatfile,
70       &                   tdatfile, sdatfile,       &                   bpstartdate1, bpstartdate2, bpdatfile,
71       &                   ctdtfile, ctdsfile,       &                   tdatfile, sdatfile, ctdtfile, ctdsfile,
72       &                   ctdtclimfile, ctdsclimfile,       &                   ctdtclimfile, ctdsclimfile, xbtfile,
      &                   xbtfile,  
73       &                   argotstartdate1, argotstartdate2, argotfile,       &                   argotstartdate1, argotstartdate2, argotfile,
74       &                   argosstartdate1, argosstartdate2, argosfile,       &                   argosstartdate1, argosstartdate2, argosfile,
75       &                   udriftfile, vdriftfile,       &                   udriftfile, vdriftfile,
# Line 69  c--   Read the namelist input. Line 77  c--   Read the namelist input.
77       &                   curmtrufile,curmtrvfile,       &                   curmtrufile,curmtrvfile,
78       &                   salterrfile, temperrfile, velerrfile,       &                   salterrfile, temperrfile, velerrfile,
79       &                   salt0errfile, temp0errfile, vel0errfile,       &                   salt0errfile, temp0errfile, vel0errfile,
80         &                   ssterrfile, ssserrfile, bperrfile,
81       &                   hflux_errfile, sflux_errfile,       &                   hflux_errfile, sflux_errfile,
82       &                   tauu_errfile, tauv_errfile,       &                   tauu_errfile, tauv_errfile,
83       &                   atemp_errfile, aqh_errfile,       &                   atemp_errfile, aqh_errfile,
84       &                   precip_errfile, swflux_errfile, swdown_errfile,       &                   precip_errfile, swflux_errfile, swdown_errfile,
85         &                   lwflux_errfile, lwdown_errfile, evap_errfile,
86         &                   snowprecip_errfile, apressure_errfile,
87       &                   uwind_errfile, vwind_errfile,       &                   uwind_errfile, vwind_errfile,
88       &                   geoid_errfile, geoid_covariancefile,       &                   geoid_errfile, geoid_covariancefile,
89         &                   bottomdrag_errfile,edtau_errfile,
90         &                   kapgm_errfile,kapredi_errfile,diffkr_errfile,
91       &                   mult_hflux, mult_sflux, mult_tauu, mult_tauv,       &                   mult_hflux, mult_sflux, mult_tauu, mult_tauv,
92       &                   mult_atemp, mult_aqh,       &                   mult_hfluxmm, mult_sfluxmm,
93         &                   mult_atemp, mult_aqh,
94       &                   mult_precip, mult_swflux, mult_swdown,       &                   mult_precip, mult_swflux, mult_swdown,
95         &                   mult_lwflux, mult_lwdown, mult_evap,
96         &                   mult_snowprecip, mult_apressure,
97       &                   mult_uwind, mult_vwind,       &                   mult_uwind, mult_vwind,
98       &                   mult_hmean, mult_h, mult_temp, mult_salt,       &                   mult_hmean, mult_h, mult_sshv4cost,
99       &                   mult_tmi, mult_sst, mult_sss,       &                   mult_tp, mult_ers, mult_gfo,
100       &                   mult_ctdt, mult_ctds,       &                   mult_temp, mult_salt,
101       &                   mult_ctdtclim, mult_ctdsclim,       &                   mult_tmi, mult_sst, mult_sss, mult_bp,
102         &                   mult_ctdt, mult_ctds,
103         &                   mult_ctdtclim, mult_ctdsclim,
104       &                   mult_xbt, mult_drift,       &                   mult_xbt, mult_drift,
105       &                   mult_argot, mult_argos,       &                   mult_argot, mult_argos,
106       &                   mult_scatx, mult_scaty,       &                   mult_scatx, mult_scaty,
# Line 91  c--   Read the namelist input. Line 109  c--   Read the namelist input.
109       &                   mult_obcsn, mult_obcss,       &                   mult_obcsn, mult_obcss,
110       &                   mult_obcsw, mult_obcse, mult_obcsvol,       &                   mult_obcsw, mult_obcse, mult_obcsvol,
111       &                   mult_curmtr, mult_ageos,       &                   mult_curmtr, mult_ageos,
112       &                   mult_kapgm, mult_diffkr,       &                   mult_kapgm, mult_kapredi, mult_diffkr,
113         &                   mult_edtau, mult_bottomdrag,
114         &                   mult_transp,
115         &                   mult_smooth_ic, mult_smooth_bc,
116       &                   whflux0, wsflux0, wtau0,       &                   whflux0, wsflux0, wtau0,
117       &                   watemp0, waqh0,       &                   wbottomdrag0,
118         &                   watemp0, waqh0,
119       &                   wprecip0, wswflux0, wswdown0, wwind0,       &                   wprecip0, wswflux0, wswdown0, wwind0,
120       &                   wmean_hflux, wmean_sflux, wmean_tau,       &                   wsnowprecip0, wlwflux0, wlwdown0, wevap0,
121       &                   wmean_atemp, wmean_aqh, wmean_precip,       &                   wapressure0, wdiffkr0, wkapgm0, wedtau0,
122         &                   wkapredi0, wmean_hflux, wmean_sflux, wmean_tau,
123         &                   wmean_atemp, wmean_aqh, wmean_precip,
124       &                   wmean_swflux, wmean_swdown, wmean_wind,       &                   wmean_swflux, wmean_swdown, wmean_wind,
125         &                   wmean_lwflux, wmean_lwdown, wmean_evap,
126         &                   wmean_snowprecip, wmean_apressure,
127       &                   cost_iprec, cost_yftype       &                   cost_iprec, cost_yftype
128    
129    #ifdef ALLOW_GENCOST_CONTRIBUTION
130          namelist /ecco_gencost_nml/
131         &         gencost_datafile,
132         &         gencost_errfile,
133         &         gencost_mask,
134         &         mult_gencost
135    #endif
136    
137        _BEGIN_MASTER( myThid )        _BEGIN_MASTER( myThid )
138    
139  c--     Set default values.  c--     Set default values.
140          data_errfile          =    ' '          data_errfile          =    ' '
141          tbarfile              =    'tbar'          tbarfile              =    'tbar'
142          sbarfile              =    'sbar'          sbarfile              =    'sbar'
143            sstbarfile            =    'sstbar'
144          psbarfile             =    'psbar'          psbarfile             =    'psbar'
145            bpbarfile             =    'bpbar'
146          ubarfile              =    'ubar'          ubarfile              =    'ubar'
147          vbarfile              =    'vbar'          vbarfile              =    'vbar'
148          wbarfile              =    'wbar'          wbarfile              =    'wbar'
149          tauxbarfile           =    'tauxbar'          tauxbarfile           =    'tauxbar'
150          tauybarfile           =    'tauybar'          tauybarfile           =    'tauybar'
151          hfluxbarfile          =    'hfluxbar'          hfluxmeanbarfile      =    'hfluxmeanbar'
152          sfluxbarfile          =    'sfluxbar'          sfluxmeanbarfile      =    'sfluxmeanbar'
153            costTranspDataFile    =    ' '
154          using_topex           = .false.          using_topex           = .false.
155          topexstartdate1       =      0          topexstartdate1       =      0
156          topexstartdate2       =      0          topexstartdate2       =      0
157          topexperiod           =      0. _d 0          topexperiod           =      0. _d 0
158            topexintercept        =      0. _d 0
159            topexslope            =      0. _d 0
160          topexfile             =    ' '          topexfile             =    ' '
161          topexmeanfile         =    ' '          topexmeanfile         =    ' '
162          using_ers             = .false.          using_ers             = .false.
163          ersstartdate1         =      0          ersstartdate1         =      0
164          ersstartdate2         =      0          ersstartdate2         =      0
165          ersperiod             =      0. _d 0          ersperiod             =      0. _d 0
166            ersintercept          =      0. _d 0
167            ersslope              =      0. _d 0
168          ersfile               =    ' '          ersfile               =    ' '
169            using_gfo             = .false.
170            gfostartdate1         =      0
171            gfostartdate2         =      0
172            gfoperiod             =      0. _d 0
173            gfointercept          =      0. _d 0
174            gfoslope              =      0. _d 0
175            gfofile               =    ' '
176          scatstartdate1        =      0          scatstartdate1        =      0
177          scatstartdate2        =      0          scatstartdate2        =      0
178          scatperiod            =      0. _d 0          scatperiod            =      0. _d 0
179          scatxdatfile          =    ' '          scatxdatfile          =    ' '
180          scatydatfile          =    ' '          scatydatfile          =    ' '
181          ssh_errfile           =    ' '          ssh_errfile           =    ' '
182            tp_errfile            =    ' '
183            ers_errfile           =    ' '
184            gfo_errfile           =    ' '
185          scatx_errfile         =    ' '          scatx_errfile         =    ' '
186          scaty_errfile         =    ' '          scaty_errfile         =    ' '
187          tmistartdate1         =      0          tmistartdate1         =      0
# Line 142  c--     Set default values. Line 193  c--     Set default values.
193          sssstartdate1         =      0          sssstartdate1         =      0
194          sssstartdate2         =      0          sssstartdate2         =      0
195          sssdatfile            =    ' '          sssdatfile            =    ' '
196            bpstartdate1          =      0
197            bpstartdate2          =      0
198            bpdatfile             =    ' '
199          tdatfile              =    ' '          tdatfile              =    ' '
200          sdatfile              =    ' '          sdatfile              =    ' '
201          ctdtfile              =    ' '          ctdtfile              =    ' '
# Line 164  c--     Set default values. Line 218  c--     Set default values.
218          salterrfile           =    ' '          salterrfile           =    ' '
219          temperrfile           =    ' '          temperrfile           =    ' '
220          velerrfile            =    ' '          velerrfile            =    ' '
221            salt0errfile          =    ' '
222            temp0errfile          =    ' '
223            vel0errfile           =    ' '
224            ssterrfile            =    ' '
225            ssserrfile            =    ' '
226            bperrfile             =    ' '
227          hflux_errfile         =    ' '          hflux_errfile         =    ' '
228          sflux_errfile         =    ' '          sflux_errfile         =    ' '
229          tauu_errfile          =    ' '          tauu_errfile          =    ' '
# Line 173  c--     Set default values. Line 233  c--     Set default values.
233          precip_errfile        =    ' '          precip_errfile        =    ' '
234          swflux_errfile        =    ' '          swflux_errfile        =    ' '
235          swdown_errfile        =    ' '          swdown_errfile        =    ' '
236            lwflux_errfile        =    ' '
237            lwdown_errfile        =    ' '
238            evap_errfile          =    ' '
239            snowprecip_errfile    =    ' '
240            apressure_errfile     =    ' '
241          uwind_errfile         =    ' '          uwind_errfile         =    ' '
242          vwind_errfile         =    ' '          vwind_errfile         =    ' '
243          geoid_errfile         =    ' '          geoid_errfile         =    ' '
244          geoid_covariancefile  =    ' '          geoid_covariancefile  =    ' '
245            bottomdrag_errfile    =    ' '
246            edtau_errfile         =    ' '
247            kapgm_errfile         =    ' '
248            kapredi_errfile       =    ' '
249            diffkr_errfile        =    ' '
250            do k=1,NSSHV4COST
251             sshv4cost_errfile(k) =    ' '
252             sshv4cost_errfactor(k)=     1. _d 0
253            enddo
254            tpTimeMaskFile        =    ' '
255            ersTimeMaskFile       =    ' '
256            gfoTimeMaskFile       =    ' '
257    
258          mult_hflux            =      0. _d 0          mult_hflux            =      0. _d 0
259          mult_sflux            =      0. _d 0          mult_sflux            =      0. _d 0
260            mult_hfluxmm          =      0. _d 0
261            mult_sfluxmm          =      0. _d 0
262          mult_tauu             =      0. _d 0          mult_tauu             =      0. _d 0
263          mult_tauv             =      0. _d 0          mult_tauv             =      0. _d 0
264          mult_atemp            =      0. _d 0          mult_atemp            =      0. _d 0
# Line 186  c--     Set default values. Line 266  c--     Set default values.
266          mult_precip           =      0. _d 0          mult_precip           =      0. _d 0
267          mult_swflux           =      0. _d 0          mult_swflux           =      0. _d 0
268          mult_swdown           =      0. _d 0          mult_swdown           =      0. _d 0
269            mult_lwflux           =      0. _d 0
270            mult_lwdown           =      0. _d 0
271            mult_evap             =      0. _d 0
272            mult_snowprecip       =      0. _d 0
273            mult_apressure        =      0. _d 0
274          mult_uwind            =      0. _d 0          mult_uwind            =      0. _d 0
275          mult_vwind            =      0. _d 0          mult_vwind            =      0. _d 0
276          mult_hmean            =      0. _d 0          mult_hmean            =      0. _d 0
277          mult_h                =      0. _d 0          mult_h                =      1. _d 0
278            mult_tp               =      0. _d 0
279            mult_ers              =      0. _d 0
280            mult_gfo              =      0. _d 0
281          mult_temp             =      0. _d 0          mult_temp             =      0. _d 0
282          mult_salt             =      0. _d 0          mult_salt             =      0. _d 0
283          mult_temp0            =      0. _d 0          mult_temp0            =      0. _d 0
# Line 197  c--     Set default values. Line 285  c--     Set default values.
285          mult_tmi              =      0. _d 0          mult_tmi              =      0. _d 0
286          mult_sst              =      0. _d 0          mult_sst              =      0. _d 0
287          mult_sss              =      0. _d 0          mult_sss              =      0. _d 0
288            mult_bp               =      0. _d 0
289          mult_ctdt             =      0. _d 0          mult_ctdt             =      0. _d 0
290          mult_ctds             =      0. _d 0          mult_ctds             =      0. _d 0
291          mult_ctdtclim         =      0. _d 0          mult_ctdtclim         =      0. _d 0
# Line 218  c--     Set default values. Line 307  c--     Set default values.
307          mult_ageos            =      0. _d 0          mult_ageos            =      0. _d 0
308          mult_curmtr           =      0. _d 0          mult_curmtr           =      0. _d 0
309          mult_kapgm            =      0. _d 0          mult_kapgm            =      0. _d 0
310            mult_kapredi            =      0. _d 0
311          mult_diffkr           =      0. _d 0          mult_diffkr           =      0. _d 0
312            mult_edtau            =      0. _d 0
313            mult_bottomdrag       =      0. _d 0
314            do k=1,NSSHV4COST
315            mult_sshv4cost(k)       =      0. _d 0
316            enddo
317            mult_smooth_ic        =      0. _d 0
318            mult_smooth_bc        =      0. _d 0
319            mult_transp           =      0. _d 0
320          watemp0               =      0. _d 0          watemp0               =      0. _d 0
321          waqh0                 =      0. _d 0          waqh0                 =      0. _d 0
322          wprecip0              =      0. _d 0          wprecip0              =      0. _d 0
323          wswflux0              =      0. _d 0          wswflux0              =      0. _d 0
324          wswdown0              =      0. _d 0          wswdown0              =      0. _d 0
325            wlwflux0              =      0. _d 0
326            wlwdown0              =      0. _d 0
327            wevap0                =      0. _d 0
328            wsnowprecip0          =      0. _d 0
329            wapressure0           =      0. _d 0
330          wwind0                =      0. _d 0          wwind0                =      0. _d 0
331            wdiffkr0              =      1. _d 0
332            wkapgm0               =      1. _d 0
333            wkapredi0               =      1. _d 0
334            wedtau0               =      1. _d 0
335          whflux0               =      0. _d 0          whflux0               =      0. _d 0
336          wsflux0               =      0. _d 0          wsflux0               =      0. _d 0
337          wtau0                 =      0. _d 0          wtau0                 =      0. _d 0
338            wbottomdrag0          =      1. _d 0
339          wmean_atemp           =      0. _d 0          wmean_atemp           =      0. _d 0
340          wmean_aqh             =      0. _d 0          wmean_aqh             =      0. _d 0
341          wmean_precip          =      0. _d 0          wmean_precip          =      0. _d 0
342          wmean_swflux          =      0. _d 0          wmean_swflux          =      0. _d 0
343          wmean_swdown          =      0. _d 0          wmean_swdown          =      0. _d 0
344            wmean_lwflux          =      0. _d 0
345            wmean_lwdown          =      0. _d 0
346            wmean_evap            =      0. _d 0
347            wmean_snowprecip      =      0. _d 0
348            wmean_apressure       =      0. _d 0
349          wmean_wind            =      0. _d 0          wmean_wind            =      0. _d 0
350          wmean_hflux           =      0. _d 0          wmean_hflux           =      0. _d 0
351          wmean_sflux           =      0. _d 0          wmean_sflux           =      0. _d 0
352          wmean_tau             =      0. _d 0          wmean_tau             =      0. _d 0
353            cost_iprec  = precFloat32
354            cost_yftype = 'RL'
355    
356  c--     Next, read the cost data file.  c--     Next, read the cost data file.
357          WRITE(msgBuf,'(A)') 'ECCO_READPARMS: opening data.ecco'          WRITE(msgBuf,'(A)') 'ECCO_READPARMS: opening data.ecco'
# Line 250  c--     Next, read the cost data file. Line 365  c--     Next, read the cost data file.
365    
366          READ(  iUnit, nml = ecco_cost_nml )          READ(  iUnit, nml = ecco_cost_nml )
367    
368          WRITE(msgBuf,'(A)') 'ECCO_READPARMS: finished reading data.ecco'          WRITE(msgBuf,'(A)')
369         &       'ECCO_READPARMS: finished reading #1: ecco_cost_nml'
370          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
371       &                  SQUEEZE_RIGHT , 1)       &       SQUEEZE_RIGHT , 1)
372    
373          CLOSE( iUnit )  C--   Check parameters:
374            IF ( cost_yftype.NE.'RL' ) THEN
375              STOP 'ABNORMAL END: S/R ECCO_READPARMS: invalid cost_yftype'
376            ENDIF
377    
378  #ifdef ALLOW_CAL  #ifdef ALLOW_CAL
379  c--     Get the complete dates of the ...  c--     Get the complete dates of the ...
380  c--     ... SST data.  c--     ... SST data.
381          call cal_FullDate( sststartdate1,   sststartdate2,          WRITE(msgBuf,'(A)') 'ECCO_READPARMS: start assigning cost dates'
382       &                     sststartdate,    mythid )          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
383         &                  SQUEEZE_RIGHT , 1)
384    c--     ... SST data.
385            if ( sstdatfile .ne. ' ' )
386         &   call cal_FullDate( sststartdate1,   sststartdate2,
387         &                      sststartdate,    mythid )
388  c--     ... TMI data.  c--     ... TMI data.
389          call cal_FullDate( tmistartdate1,   tmistartdate2,          if ( tmidatfile .ne. ' ' )
390       &                     tmistartdate,    mythid )       &   call cal_FullDate( tmistartdate1,   tmistartdate2,
391         &                      tmistartdate,    mythid )
392  c--     ... SSS data.  c--     ... SSS data.
393          call cal_FullDate( sssstartdate1,   sssstartdate2,          if ( sssdatfile .ne. ' ' )
394       &                     sssstartdate,    mythid )           &   call cal_FullDate( sssstartdate1,   sssstartdate2,
395         &                      sssstartdate,    mythid )
396    c--     ... BP data.
397            if ( bpdatfile .ne. ' ' )
398         &   call cal_FullDate( bpstartdate1,   bpstartdate2,
399         &                      bpstartdate,    mythid )
400  c--     ... T/P data.  c--     ... T/P data.
401          call cal_FullDate( topexstartdate1, topexstartdate2,          if ( topexfile .ne. ' ' )
402         &  call cal_FullDate( topexstartdate1, topexstartdate2,
403       &                     topexstartdate,  mythid )       &                     topexstartdate,  mythid )
404  c--     ... ERS data.  c--     ... ERS data.
405          call cal_FullDate( ersstartdate1,   ersstartdate2,          if ( ersfile .ne. ' ' )
406       &                     ersstartdate,    mythid )       &   call cal_FullDate( ersstartdate1,   ersstartdate2,
407         &                      ersstartdate,    mythid )
408    c--     ... GFO data.
409            if ( gfofile .ne. ' ' )
410         &   call cal_FullDate( gfostartdate1,   gfostartdate2,
411         &                      gfostartdate,    mythid )
412  c--     ... SCAT data.  c--     ... SCAT data.
413          call cal_FullDate( scatstartdate1,   scatstartdate2,          if ( scatxdatfile .ne. ' ' )
414       &                     scatxstartdate,    mythid )       &   call cal_FullDate( scatstartdate1,   scatstartdate2,
415          call cal_FullDate( scatstartdate1,   scatstartdate2,         &                      scatxstartdate,    mythid )
416       &                     scatystartdate,    mythid )          if ( scatydatfile .ne. ' ' )
417         &   call cal_FullDate( scatstartdate1,   scatstartdate2,
418         &                      scatystartdate,    mythid )
419    
420            WRITE(msgBuf,'(A)') 'ECCO_READPARMS: end assigning cost dates'
421            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
422         &                  SQUEEZE_RIGHT , 1)
423  #endif  #endif
424    
425        _END_MASTER( mythid )        _END_MASTER( myThid )
426    
427    #ifdef ALLOW_GENCOST_CONTRIBUTION
428    
429          _BEGIN_MASTER( myThid )
430    
431    c--   Set default values.
432          do k=1,NGENCOST
433             gencost_avgperiod(k)  = '     '
434             gencost_datafile(k)   = ' '
435             gencost_errfile(k)    = ' '
436             gencost_mask(k)       = ' '
437             write(gencost_barfile(k)(1:8),'(a,i2.2)'),
438         &        'genbar', k
439             mult_gencost(k)      = 0. _d 0
440          enddo
441    
442          READ(  iUnit, nml = ecco_gencost_nml )
443    
444          WRITE(msgBuf,'(A)')
445         &     'ECCO_READPARMS: finished reading #2: ecco_gencost_nml'
446          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
447         &     SQUEEZE_RIGHT , 1)
448    
449    C--   Check parameters:
450    
451          do k=1,NGENCOST
452             if ( gencost_avgperiod(k) .NE. 'day' .OR.
453         &        gencost_avgperiod(k) .NE. 'DAY' .OR.
454         &        gencost_avgperiod(k) .NE. 'month' .OR.
455         &        gencost_avgperiod(k) .NE. 'MONTH' .OR.
456         &        gencost_avgperiod(k) .NE. 'year' .OR.
457         &        gencost_avgperiod(k) .NE. 'YEAR' ) then
458                STOP
459         &           'ECCO_READPARMS: gencost_avgperiod not properly set'
460             end if
461          end do
462    
463          _END_MASTER( myThid )
464    
465    #endif /* ALLOW_GENCOST_CONTRIBUTION */
466    
467          _BEGIN_MASTER( myThid )
468    
469          CLOSE( iUnit )
470    
471          _END_MASTER( myThid )
472    
473          WRITE(msgBuf,'(A)')
474         &     'ECCO_READPARMS: done'
475          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
476         &     SQUEEZE_RIGHT , 1)
477    
478        _BARRIER        _BARRIER
479    
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.45
  ViewVC Help
Powered by ViewVC 1.1.22