/[MITgcm]/MITgcm/pkg/ctrl/ctrl_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_readparms.F

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

revision 1.7 by heimbach, Thu Nov 11 06:28:33 2004 UTC revision 1.26 by gforget, Wed Jun 20 19:31:03 2007 UTC
# Line 34  c     == global variables == Line 34  c     == global variables ==
34  #ifdef ALLOW_OBCS_CONTROL  #ifdef ALLOW_OBCS_CONTROL
35  # include "OBCS.h"  # include "OBCS.h"
36  #endif  #endif
 #ifdef ALLOW_ECCO_OPTIMIZATION  
 # include "optim.h"  
 #endif  
37    
38  c     == routine arguments ==  c     == routine arguments ==
39    
# Line 48  c     == local variables == Line 45  c     == local variables ==
45        integer i,j,k        integer i,j,k
46        integer ntmp        integer ntmp
47        integer ivarindex        integer ivarindex
48          integer iUnit
49        integer iobcs        integer iobcs
50        integer il        integer il
51        integer errio        integer errio
# Line 73  c     == end of interface == Line 70  c     == end of interface ==
70    
71  c--   Read the namelist input.  c--   Read the namelist input.
72        namelist /ctrl_nml/        namelist /ctrl_nml/
73       &  xx_theta_file,       &  xx_theta_file, xx_salt_file,
74       &  xx_salt_file,       &  xx_hflux_file, xx_hflux_remo_intercept, xx_hflux_remo_slope,
      &  xx_hflux_file,  
75       &  xx_hfluxstartdate1, xx_hfluxstartdate2, xx_hfluxperiod,       &  xx_hfluxstartdate1, xx_hfluxstartdate2, xx_hfluxperiod,
76       &  xx_sflux_file,       &  xx_sflux_file, xx_sflux_remo_intercept, xx_sflux_remo_slope,
77       &  xx_sfluxstartdate1, xx_sfluxstartdate2, xx_sfluxperiod,       &  xx_sfluxstartdate1, xx_sfluxstartdate2, xx_sfluxperiod,
78       &  xx_tauu_file,       &  xx_tauu_file, xx_tauu_remo_intercept, xx_tauu_remo_slope,
79       &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,       &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,
80       &  xx_tauv_file,       &  xx_tauv_file, xx_tauv_remo_intercept, xx_tauv_remo_slope,
81       &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod,       &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod,
82       &  xx_atemp_file,       &  xx_atemp_file, xx_atemp_remo_intercept, xx_atemp_remo_slope,
83       &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,       &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
84       &  xx_aqh_file,       &  xx_aqh_file, xx_aqh_remo_intercept, xx_aqh_remo_slope,
85       &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,       &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,
86       &  xx_uwind_file,       &  xx_precip_file, xx_precip_remo_intercept, xx_precip_remo_slope,
87         &  xx_precipstartdate1, xx_precipstartdate2, xx_precipperiod,
88         &  xx_swflux_file, xx_swflux_remo_intercept, xx_swflux_remo_slope,
89         &  xx_swfluxstartdate1, xx_swfluxstartdate2, xx_swfluxperiod,
90         &  xx_swdown_file, xx_swdown_remo_intercept, xx_swdown_remo_slope,
91         &  xx_swdownstartdate1, xx_swdownstartdate2, xx_swdownperiod,
92         &  xx_lwflux_file, xx_lwflux_remo_intercept, xx_lwflux_remo_slope,
93         &  xx_lwfluxstartdate1, xx_lwfluxstartdate2, xx_lwfluxperiod,
94         &  xx_lwdown_file, xx_lwdown_remo_intercept, xx_lwdown_remo_slope,
95         &  xx_lwdownstartdate1, xx_lwdownstartdate2, xx_lwdownperiod,
96         &  xx_evap_file, xx_evap_remo_intercept, xx_evap_remo_slope,
97         &  xx_evapstartdate1, xx_evapstartdate2, xx_evapperiod,
98         &  xx_snowprecip_file, xx_snowprecip_remo_intercept,
99         &  xx_snowprecip_remo_slope,  xx_snowprecipperiod,
100         &  xx_snowprecipstartdate1, xx_snowprecipstartdate2,
101         &  xx_apressure_file, xx_apressure_remo_intercept,
102         &  xx_apressure_remo_slope, xx_apressureperiod,
103         &  xx_apressurestartdate1, xx_apressurestartdate2,
104         &  xx_runoff_file, xx_runoff_remo_intercept, xx_runoff_remo_slope,
105         &  xx_runoffstartdate1, xx_runoffstartdate2, xx_runoffperiod,
106         &  xx_uwind_file, xx_uwind_remo_intercept, xx_uwind_remo_slope,
107       &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,       &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
108       &  xx_vwind_file,       &  xx_vwind_file, xx_vwind_remo_intercept, xx_vwind_remo_slope,
109       &  xx_vwindstartdate1, xx_vwindstartdate2, xx_vwindperiod,       &  xx_vwindstartdate1, xx_vwindstartdate2, xx_vwindperiod,
110       &  xx_obcsn_file,       &  xx_obcsn_file,
111       &  xx_obcsnstartdate1,  xx_obcsnstartdate2,  xx_obcsnperiod,       &  xx_obcsnstartdate1,  xx_obcsnstartdate2,  xx_obcsnperiod,
# Line 99  c--   Read the namelist input. Line 115  c--   Read the namelist input.
115       &  xx_obcswstartdate1,  xx_obcswstartdate2,  xx_obcswperiod,       &  xx_obcswstartdate1,  xx_obcswstartdate2,  xx_obcswperiod,
116       &  xx_obcse_file,       &  xx_obcse_file,
117       &  xx_obcsestartdate1,  xx_obcsestartdate2,  xx_obcseperiod,       &  xx_obcsestartdate1,  xx_obcsestartdate2,  xx_obcseperiod,
118       &  xx_diffkr_file,       &  xx_diffkr_file, xx_kapgm_file, xx_tr1_file,
119       &  xx_kapgm_file,       &  xx_sst_file, xx_sss_file,
120       &  xx_tr1_file,       &  xx_sststartdate1, xx_sststartdate2, xx_sstperiod,
121       &  xx_sst_file,       &  xx_sssstartdate1, xx_sssstartdate2, xx_sssperiod,
122       &  xx_sss_file,       &  xx_depth_file,
123       &  xx_hfacc_file,       &  xx_efluxy_file, xx_efluxp_file,
      &  xx_efluxy_file,  
      &  xx_efluxp_file,  
124       &  xx_bottomdrag_file,       &  xx_bottomdrag_file,
125       &  doInitXX       &  xx_edtaux_file, xx_edtauy_file,
126         &  xx_uvel_file, xx_vvel_file, xx_etan_file,
127         &  doInitXX,
128         &  doPackDiag,
129         &  doZscaleUnpack, doZscalePack,
130         &  doMainUnpack, doMainPack,
131         &  doAdmtlmBypassAD,
132         &  delZexp, forcingPrecond
133    
134        namelist /ctrl_packnames/        namelist /ctrl_packnames/
135       &  yadmark, yctrlid, yctrlposunpack, yctrlpospack,       &  yadmark, yctrlid, yctrlposunpack, yctrlpospack,
# Line 117  c--   Read the namelist input. Line 138  c--   Read the namelist input.
138        _BEGIN_MASTER( myThid )        _BEGIN_MASTER( myThid )
139    
140  c--     Set default values.  c--     Set default values.
141          xx_theta_file      = ' '          doInitXX            = .TRUE.
142          xx_salt_file       = ' '  #ifdef ALLOW_ADMTLM
143            doAdmTlm            = .TRUE.
144    #else
145            doAdmTlm            = .FALSE.
146    #endif
147            doPackDiag          = .FALSE.
148            doZscaleUnpack      = .FALSE.
149            doZscalePack        = .FALSE.
150            doMainUnpack        = .TRUE.
151            doMainPack          = .TRUE.
152            doAdmtlmBypassAD    = .FALSE.
153    
154            delZexp = 0.
155            forcingPrecond = 1. _d 0
156    
157            xx_theta_file      = 'xx_theta'
158            xx_salt_file       = 'xx_salt'
159    c
160          xx_hfluxstartdate1 =   0          xx_hfluxstartdate1 =   0
161          xx_hfluxstartdate2 =   0          xx_hfluxstartdate2 =   0
162          xx_hfluxperiod     =   0. _d 0          xx_hfluxperiod     =   0. _d 0
163          xx_hflux_file      = ' '          xx_hflux_file      = 'xx_hfl'
164            xx_hflux_remo_intercept =   0. _d 0
165            xx_hflux_remo_slope =   0. _d 0
166    c
167          xx_sfluxstartdate1 =   0          xx_sfluxstartdate1 =   0
168          xx_sfluxstartdate2 =   0          xx_sfluxstartdate2 =   0
169          xx_sfluxperiod     =   0. _d 0          xx_sfluxperiod     =   0. _d 0
170          xx_sflux_file      = ' '          xx_sflux_file      = 'xx_sfl'
171            xx_sflux_remo_intercept =   0. _d 0
172            xx_sflux_remo_slope =   0. _d 0
173    c
174          xx_tauustartdate1  =   0          xx_tauustartdate1  =   0
175          xx_tauustartdate2  =   0          xx_tauustartdate2  =   0
176          xx_tauuperiod      =   0. _d 0          xx_tauuperiod      =   0. _d 0
177          xx_tauu_file       = ' '          xx_tauu_file       = 'xx_tauu'
178            xx_tauu_remo_intercept =   0. _d 0
179            xx_tauu_remo_slope =   0. _d 0
180    c
181          xx_tauvstartdate1  =   0          xx_tauvstartdate1  =   0
182          xx_tauvstartdate2  =   0          xx_tauvstartdate2  =   0
183          xx_tauvperiod      =   0. _d 0          xx_tauvperiod      =   0. _d 0
184          xx_tauv_file       = ' '          xx_tauv_file       = 'xx_tauv'
185            xx_tauv_remo_intercept =   0. _d 0
186            xx_tauv_remo_slope =   0. _d 0
187    c
188          xx_atempstartdate1  =   0          xx_atempstartdate1  =   0
189          xx_atempstartdate2  =   0          xx_atempstartdate2  =   0
190          xx_atempperiod      =   0. _d 0          xx_atempperiod      =   0. _d 0
191          xx_atemp_file       = ' '          xx_atemp_file       = 'xx_atemp'
192            xx_atemp_remo_intercept =   0. _d 0
193            xx_atemp_remo_slope =   0. _d 0
194    c
195          xx_aqhstartdate1    =   0          xx_aqhstartdate1    =   0
196          xx_aqhstartdate2    =   0          xx_aqhstartdate2    =   0
197          xx_aqhperiod        =   0. _d 0          xx_aqhperiod        =   0. _d 0
198          xx_aqh_file         = ' '          xx_aqh_file         = 'xx_aqh'
199            xx_aqh_remo_intercept =   0. _d 0
200            xx_aqh_remo_slope =   0. _d 0
201    c
202            xx_precipstartdate1 =   0
203            xx_precipstartdate2 =   0
204            xx_precipperiod     =   0. _d 0
205            xx_precip_file      = 'xx_precip'
206            xx_precip_remo_intercept =   0. _d 0
207            xx_precip_remo_slope =   0. _d 0
208    c
209            xx_swfluxstartdate1 =   0
210            xx_swfluxstartdate2 =   0
211            xx_swfluxperiod     =   0. _d 0
212            xx_swflux_file      = 'xx_swflux'
213            xx_swflux_remo_intercept =   0. _d 0
214            xx_swflux_remo_slope =   0. _d 0
215    c
216            xx_swdownstartdate1 =   0
217            xx_swdownstartdate2 =   0
218            xx_swdownperiod     =   0. _d 0
219            xx_swdown_file      = 'xx_swdown'
220            xx_swdown_remo_intercept =   0. _d 0
221            xx_swdown_remo_slope =   0. _d 0
222    c
223            xx_lwfluxstartdate1 =   0
224            xx_lwfluxstartdate2 =   0
225            xx_lwfluxperiod     =   0. _d 0
226            xx_lwflux_file      = 'xx_lwflux'
227            xx_lwflux_remo_intercept =   0. _d 0
228            xx_lwflux_remo_slope =   0. _d 0
229    c
230            xx_lwdownstartdate1 =   0
231            xx_lwdownstartdate2 =   0
232            xx_lwdownperiod     =   0. _d 0
233            xx_lwdown_file      = 'xx_lwdown'
234            xx_lwdown_remo_intercept =   0. _d 0
235            xx_lwdown_remo_slope =   0. _d 0
236    c
237            xx_evapstartdate1 =   0
238            xx_evapstartdate2 =   0
239            xx_evapperiod     =   0. _d 0
240            xx_evap_file      = 'xx_evap'
241            xx_evap_remo_intercept =   0. _d 0
242            xx_evap_remo_slope =   0. _d 0
243    c
244            xx_snowprecipstartdate1 =   0
245            xx_snowprecipstartdate2 =   0
246            xx_snowprecipperiod     =   0. _d 0
247            xx_snowprecip_file      = 'xx_snowprecip'
248            xx_snowprecip_remo_intercept =   0. _d 0
249            xx_snowprecip_remo_slope =   0. _d 0
250    c
251            xx_apressurestartdate1 =   0
252            xx_apressurestartdate2 =   0
253            xx_apressureperiod     =   0. _d 0
254            xx_apressure_file      = 'xx_apressure'
255            xx_apressure_remo_intercept =   0. _d 0
256            xx_apressure_remo_slope =   0. _d 0
257    c
258            xx_runoffstartdate1 =   0
259            xx_runoffstartdate2 =   0
260            xx_runoffperiod     =   0. _d 0
261            xx_runoff_file      = 'xx_runoff'
262            xx_runoff_remo_intercept =   0. _d 0
263            xx_runoff_remo_slope =   0. _d 0
264    c
265          xx_uwindstartdate1  =   0          xx_uwindstartdate1  =   0
266          xx_uwindstartdate2  =   0          xx_uwindstartdate2  =   0
267          xx_uwindperiod      =   0. _d 0          xx_uwindperiod      =   0. _d 0
268          xx_uwind_file       = ' '          xx_uwind_file       = 'xx_uwind'
269            xx_uwind_remo_intercept =   0. _d 0
270            xx_uwind_remo_slope =   0. _d 0
271    c
272          xx_vwindstartdate1  =   0          xx_vwindstartdate1  =   0
273          xx_vwindstartdate2  =   0          xx_vwindstartdate2  =   0
274          xx_vwindperiod      =   0. _d 0          xx_vwindperiod      =   0. _d 0
275          xx_vwind_file       = ' '          xx_vwind_file       = 'xx_vwind'
276            xx_vwind_remo_intercept =   0. _d 0
277            xx_vwind_remo_slope =   0. _d 0
278    c
279          xx_obcsnstartdate1  =   0          xx_obcsnstartdate1  =   0
280          xx_obcsnstartdate2  =   0          xx_obcsnstartdate2  =   0
281          xx_obcsnperiod      =   0. _d 0          xx_obcsnperiod      =   0. _d 0
282          xx_obcsn_file       = ' '          xx_obcsn_file       = 'xx_obcsn'
283    c
284          xx_obcssstartdate1  =   0          xx_obcssstartdate1  =   0
285          xx_obcssstartdate2  =   0          xx_obcssstartdate2  =   0
286          xx_obcssperiod      =   0. _d 0          xx_obcssperiod      =   0. _d 0
287          xx_obcss_file       = ' '          xx_obcss_file       = 'xx_obcss'
288    c
289          xx_obcswstartdate1  =   0          xx_obcswstartdate1  =   0
290          xx_obcswstartdate2  =   0          xx_obcswstartdate2  =   0
291          xx_obcswperiod      =   0. _d 0          xx_obcswperiod      =   0. _d 0
292          xx_obcsw_file       = ' '          xx_obcsw_file       = 'xx_obcsw'
293    c
294          xx_obcsestartdate1  =   0          xx_obcsestartdate1  =   0
295          xx_obcsestartdate2  =   0          xx_obcsestartdate2  =   0
296          xx_obcseperiod      =   0. _d 0          xx_obcseperiod      =   0. _d 0
297          xx_obcse_file       = ' '          xx_obcse_file       = 'xx_obcse'
298          xx_diffkr_file      = ' '  c
299          xx_kapgm_file       = ' '          xx_sststartdate1    =   0
300          xx_tr1_file         = ' '          xx_sststartdate2    =   0
301          xx_sst_file         = ' '          xx_sstperiod        =   0. _d 0
302          xx_sss_file         = ' '          xx_sst_file         = 'xx_sst'
303          xx_hfacc_file       = ' '  c
304          xx_efluxy_file      = ' '          xx_sssstartdate1    =   0
305          xx_efluxp_file      = ' '          xx_sssstartdate2    =   0
306          xx_bottomdrag_file  = ' '          xx_sssperiod        =   0. _d 0
307            xx_sss_file         = 'xx_sss'
308          doInitXX            = .TRUE.  c
309            xx_diffkr_file      = 'xx_diffkr'
310            xx_kapgm_file       = 'xx_kapgm'
311            xx_tr1_file         = 'xx_ptr'
312            xx_depth_file       = 'xx_depth'
313            xx_efluxy_file      = 'xx_efluxy'
314            xx_efluxp_file      = 'xx_efluxp'
315            xx_bottomdrag_file  = 'xx_bottomdrag'
316            xx_edtaux_file      = 'xx_edtaux'
317            xx_edtauy_file      = 'xx_edtauy'
318            xx_uvel_file        = 'xx_uvel'
319            xx_vvel_file        = 'xx_vvel'
320            xx_etan_file        = 'xx_etan'
321    
322  #ifdef ALLOW_TANGENTLINEAR_RUN  #ifdef ALLOW_TANGENTLINEAR_RUN
323          yadprefix           = 'g_'          yadprefix           = 'g_'
# Line 195  c--     Set default values. Line 335  c--     Set default values.
335          maskname            = ' '          maskname            = ' '
336          metaname            = ' '          metaname            = ' '
337    
338  c--     Check versions.  c--     Next, read the cost data file.
339            WRITE(msgBuf,'(A)') 'CTRL_READPARMS: opening data.ctrl'
340          open(unit=scrunit1,status='scratch')          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
341         &                    SQUEEZE_RIGHT , 1)
342    
343            CALL OPEN_COPY_DATA_FILE(
344         I                          'data.ctrl', 'CTRL_READPARMS',
345         O                          iUnit,
346         I                          myThid )
347    
348            READ(unit = iUnit, nml = ctrl_nml)
349            READ(unit = iUnit, nml = ctrl_packnames)
350    
351            WRITE(msgBuf,'(A)')
352         &       'CTRL_READPARMS: finished reading data.ctrl'
353            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
354         &                  SQUEEZE_RIGHT , 1)
355    
356  c--     Next, read the ecco data file.          CLOSE( iUnit )
         open(unit = modeldataunit,file = 'data.ctrl',  
      &       status = 'old', iostat = errio)  
         if ( errio .lt. 0 ) then  
           stop ' stopped in ctrl_readparms'  
         endif  
   
         do while ( .true. )  
           read(modeldataunit, fmt='(a)', end=1001) record  
           il = max(ilnblnk(record),1)  
           if ( record(1:1) .ne. commentcharacter )  
      &        write(unit=scrunit1, fmt='(a)') record(:il)  
         enddo  
  1001   continue  
         close( modeldataunit )  
   
         rewind( scrunit1 )  
         read(unit = scrunit1, nml = ctrl_nml)  
         read(unit = scrunit1, nml = ctrl_packnames)  
         close( scrunit1 )  
357    
358        _END_MASTER( myThid )        _END_MASTER( myThid )
359    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22