/[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.16 by heimbach, Thu Jul 28 13:47:49 2005 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 87  c--   Read the namelist input. Line 84  c--   Read the namelist input.
84       &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,       &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
85       &  xx_aqh_file,       &  xx_aqh_file,
86       &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,       &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,
87         &  xx_precip_file,
88         &  xx_precipstartdate1, xx_precipstartdate2, xx_precipperiod,
89       &  xx_uwind_file,       &  xx_uwind_file,
90       &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,       &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
91       &  xx_vwind_file,       &  xx_vwind_file,
# Line 108  c--   Read the namelist input. Line 107  c--   Read the namelist input.
107       &  xx_efluxy_file,       &  xx_efluxy_file,
108       &  xx_efluxp_file,       &  xx_efluxp_file,
109       &  xx_bottomdrag_file,       &  xx_bottomdrag_file,
110       &  doInitXX       &  xx_edtaux_file,
111         &  xx_edtauy_file,
112         &  xx_uvel_file,
113         &  xx_vvel_file,
114         &  xx_etan_file,
115         &  doInitXX,
116         &  doPackDiag,
117         &  doZscaleUnpack,
118         &  doZscalePack
119    
120    
121        namelist /ctrl_packnames/        namelist /ctrl_packnames/
122       &  yadmark, yctrlid, yctrlposunpack, yctrlpospack,       &  yadmark, yctrlid, yctrlposunpack, yctrlpospack,
# Line 117  c--   Read the namelist input. Line 125  c--   Read the namelist input.
125        _BEGIN_MASTER( myThid )        _BEGIN_MASTER( myThid )
126    
127  c--     Set default values.  c--     Set default values.
128          xx_theta_file      = ' '          doInitXX            = .TRUE.
129          xx_salt_file       = ' '  #ifdef ALLOW_ADMTLM
130            doAdmTlm            = .TRUE.
131    #else
132            doAdmTlm            = .FALSE.
133    #endif
134            doPackDiag          = .FALSE.
135            doZscaleUnpack      = .FALSE.
136            doZscalePack        = .FALSE.
137    
138            xx_theta_file      = 'xx_theta'
139            xx_salt_file       = 'xx_salt'
140          xx_hfluxstartdate1 =   0          xx_hfluxstartdate1 =   0
141          xx_hfluxstartdate2 =   0          xx_hfluxstartdate2 =   0
142          xx_hfluxperiod     =   0. _d 0          xx_hfluxperiod     =   0. _d 0
143          xx_hflux_file      = ' '          xx_hflux_file      = 'xx_hfl'
144          xx_sfluxstartdate1 =   0          xx_sfluxstartdate1 =   0
145          xx_sfluxstartdate2 =   0          xx_sfluxstartdate2 =   0
146          xx_sfluxperiod     =   0. _d 0          xx_sfluxperiod     =   0. _d 0
147          xx_sflux_file      = ' '          xx_sflux_file      = 'xx_sfl'
148          xx_tauustartdate1  =   0          xx_tauustartdate1  =   0
149          xx_tauustartdate2  =   0          xx_tauustartdate2  =   0
150          xx_tauuperiod      =   0. _d 0          xx_tauuperiod      =   0. _d 0
151          xx_tauu_file       = ' '          xx_tauu_file       = 'xx_tauu'
152          xx_tauvstartdate1  =   0          xx_tauvstartdate1  =   0
153          xx_tauvstartdate2  =   0          xx_tauvstartdate2  =   0
154          xx_tauvperiod      =   0. _d 0          xx_tauvperiod      =   0. _d 0
155          xx_tauv_file       = ' '          xx_tauv_file       = 'xx_tauv'
156          xx_atempstartdate1  =   0          xx_atempstartdate1  =   0
157          xx_atempstartdate2  =   0          xx_atempstartdate2  =   0
158          xx_atempperiod      =   0. _d 0          xx_atempperiod      =   0. _d 0
159          xx_atemp_file       = ' '          xx_atemp_file       = 'xx_atemp'
160          xx_aqhstartdate1    =   0          xx_aqhstartdate1    =   0
161          xx_aqhstartdate2    =   0          xx_aqhstartdate2    =   0
162          xx_aqhperiod        =   0. _d 0          xx_aqhperiod        =   0. _d 0
163          xx_aqh_file         = ' '          xx_aqh_file         = 'xx_aqh'
164            xx_precipstartdate1 =   0
165            xx_precipstartdate2 =   0
166            xx_precipperiod     =   0. _d 0
167            xx_precip_file      = 'xx_precip'
168          xx_uwindstartdate1  =   0          xx_uwindstartdate1  =   0
169          xx_uwindstartdate2  =   0          xx_uwindstartdate2  =   0
170          xx_uwindperiod      =   0. _d 0          xx_uwindperiod      =   0. _d 0
171          xx_uwind_file       = ' '          xx_uwind_file       = 'xx_uwind'
172          xx_vwindstartdate1  =   0          xx_vwindstartdate1  =   0
173          xx_vwindstartdate2  =   0          xx_vwindstartdate2  =   0
174          xx_vwindperiod      =   0. _d 0          xx_vwindperiod      =   0. _d 0
175          xx_vwind_file       = ' '          xx_vwind_file       = 'xx_vwind'
176          xx_obcsnstartdate1  =   0          xx_obcsnstartdate1  =   0
177          xx_obcsnstartdate2  =   0          xx_obcsnstartdate2  =   0
178          xx_obcsnperiod      =   0. _d 0          xx_obcsnperiod      =   0. _d 0
179          xx_obcsn_file       = ' '          xx_obcsn_file       = 'xx_obcsn'
180          xx_obcssstartdate1  =   0          xx_obcssstartdate1  =   0
181          xx_obcssstartdate2  =   0          xx_obcssstartdate2  =   0
182          xx_obcssperiod      =   0. _d 0          xx_obcssperiod      =   0. _d 0
183          xx_obcss_file       = ' '          xx_obcss_file       = 'xx_obcss'
184          xx_obcswstartdate1  =   0          xx_obcswstartdate1  =   0
185          xx_obcswstartdate2  =   0          xx_obcswstartdate2  =   0
186          xx_obcswperiod      =   0. _d 0          xx_obcswperiod      =   0. _d 0
187          xx_obcsw_file       = ' '          xx_obcsw_file       = 'xx_obcsw'
188          xx_obcsestartdate1  =   0          xx_obcsestartdate1  =   0
189          xx_obcsestartdate2  =   0          xx_obcsestartdate2  =   0
190          xx_obcseperiod      =   0. _d 0          xx_obcseperiod      =   0. _d 0
191          xx_obcse_file       = ' '          xx_obcse_file       = 'xx_obcse'
192          xx_diffkr_file      = ' '          xx_diffkr_file      = 'xx_diffkr'
193          xx_kapgm_file       = ' '          xx_kapgm_file       = 'xx_kapgm'
194          xx_tr1_file         = ' '          xx_tr1_file         = 'xx_ptr'
195          xx_sst_file         = ' '          xx_sst_file         = 'xx_sst'
196          xx_sss_file         = ' '          xx_sss_file         = 'xx_sss'
197          xx_hfacc_file       = ' '          xx_hfacc_file       = 'xx_hfacc'
198          xx_efluxy_file      = ' '          xx_efluxy_file      = 'xx_efluxy'
199          xx_efluxp_file      = ' '          xx_efluxp_file      = 'xx_efluxp'
200          xx_bottomdrag_file  = ' '          xx_bottomdrag_file  = 'xx_bottomdrag'
201            xx_edtaux_file      = 'xx_edtaux'
202          doInitXX            = .TRUE.          xx_edtauy_file      = 'xx_edtauy'
203            xx_uvel_file        = 'xx_uvel'
204            xx_vvel_file        = 'xx_vvel'
205            xx_etan_file        = 'xx_etan'
206    
207  #ifdef ALLOW_TANGENTLINEAR_RUN  #ifdef ALLOW_TANGENTLINEAR_RUN
208          yadprefix           = 'g_'          yadprefix           = 'g_'
# Line 195  c--     Set default values. Line 220  c--     Set default values.
220          maskname            = ' '          maskname            = ' '
221          metaname            = ' '          metaname            = ' '
222    
223  c--     Check versions.  c--     Next, read the cost data file.
224            WRITE(msgBuf,'(A)') 'CTRL_READPARMS: opening data.ctrl'
225          open(unit=scrunit1,status='scratch')          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
226         &                    SQUEEZE_RIGHT , 1)
227    
228            CALL OPEN_COPY_DATA_FILE(
229         I                          'data.ctrl', 'CTRL_READPARMS',
230         O                          iUnit,
231         I                          myThid )
232    
233            READ(unit = iUnit, nml = ctrl_nml)
234            READ(unit = iUnit, nml = ctrl_packnames)
235    
236            WRITE(msgBuf,'(A)')
237         &       'CTRL_READPARMS: finished reading data.ctrl'
238            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
239         &                  SQUEEZE_RIGHT , 1)
240    
241  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 )  
242    
243        _END_MASTER( myThid )        _END_MASTER( myThid )
244    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22