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