/[MITgcm]/MITgcm/optim/optim_numbmod.F
ViewVC logotype

Diff of /MITgcm/optim/optim_numbmod.F

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

revision 1.3 by heimbach, Thu Nov 6 22:02:33 2003 UTC revision 1.5 by heimbach, Tue Jan 25 04:42:21 2005 UTC
# Line 121  c     == end of interface == Line 121  c     == end of interface ==
121       &  xx_tauu_file,       &  xx_tauu_file,
122       &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,       &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,
123       &  xx_tauv_file,       &  xx_tauv_file,
124       &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod       &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod,
125         &  xx_atemp_file,
126         &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
127         &  xx_aqh_file,
128         &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,
129         &  xx_uwind_file,
130         &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
131         &  xx_vwind_file,
132         &  xx_vwindstartdate1, xx_vwindstartdate2, xx_vwindperiod,
133         &  xx_obcsn_file,
134         &  xx_obcsnstartdate1,  xx_obcsnstartdate2,  xx_obcsnperiod,
135         &  xx_obcss_file,
136         &  xx_obcssstartdate1,  xx_obcssstartdate2,  xx_obcssperiod,
137         &  xx_obcsw_file,
138         &  xx_obcswstartdate1,  xx_obcswstartdate2,  xx_obcswperiod,
139         &  xx_obcse_file,
140         &  xx_obcsestartdate1,  xx_obcsestartdate2,  xx_obcseperiod,
141         &  xx_diffkr_file,
142         &  xx_kapgm_file,
143         &  xx_tr1_file,
144         &  xx_sst_file,
145         &  xx_sss_file,
146         &  xx_hfacc_file,
147         &  xx_efluxy_file,
148         &  xx_efluxp_file,
149         &  xx_bottomdrag_file,
150         &  doInitXX,
151         &  doPackDiag,
152         &  doZscaleUnpack,
153         &  doZscalePack
154    
155        namelist /CTRL_PACKNAMES/        namelist /CTRL_PACKNAMES/
156       &  yadmark,       &  yadmark, yctrlid, yctrlposunpack, yctrlpospack,
157       &  ctrlname, costname, scalname, maskname, metaname       &  ctrlname, costname, scalname, maskname, metaname
158    
159        namelist /ECCO_PARMS/        namelist /OPTIM/
      &                      expId  
       
       namelist /ECCO_OPTIM/  
160       &                      optimcycle,       &                      optimcycle,
161       &                      numiter, nfunc, fmin, iprint,       &                      numiter, nfunc, fmin, iprint,
162       &                      epsf, epsx, epsg,       &                      epsf, epsx, epsg,
# Line 149  c--   Preset the optimization parameters Line 175  c--   Preset the optimization parameters
175        nupdate    =  1        nupdate    =  1
176        ff         = 0.        ff         = 0.
177        expId      = 'MIT_CE_000'        expId      = 'MIT_CE_000'
178          yctrlid    = expId
179    
180        modeldataunit = 14        modeldataunit = 14
181        scrunit1      = 11        scrunit1      = 11
182    
 c--   Read ecco parameters from file.  
       open(unit=scrunit1,status='scratch')  
   
       open(unit = modeldataunit,file = 'data.ecco',  
      &     status = 'old', iostat = errio)  
       if ( errio .lt. 0 ) then  
          stop ' stopped in optim_numbmod'  
       endif  
   
       do while ( .true. )  
          read(modeldataunit, fmt='(a)', end=20) record  
          il = max(ilnblnk(record),1)  
          if ( record(1:1) .ne. commentcharacter )  
      &        write(unit=scrunit1, fmt='(a)') record(:il)  
       enddo  
    20 continue  
       close( modeldataunit )  
   
       rewind( scrunit1 )  
       read(unit = scrunit1, nml = ecco_parms)  
       close( scrunit1 )  
       print*, ' OPTIM_NUMBMOD: ECCO options have been read.'  
   
183  c--   Read control parameters from file.  c--   Read control parameters from file.
184        open(unit=scrunit1,status='scratch')        open(unit=scrunit1,status='scratch')
185    
# Line 200  c--   Read control parameters from file. Line 204  c--   Read control parameters from file.
204        close( scrunit1 )        close( scrunit1 )
205        print*, ' OPTIM_NUMBMOD: Control options have been read.'        print*, ' OPTIM_NUMBMOD: Control options have been read.'
206    
207    cph(
208          expId = yctrlid
209    cph)
210    
211  c--   Read optimization parameters from file.  c--   Read optimization parameters from file.
212        open(unit=scrunit1,status='scratch')        open(unit=scrunit1,status='scratch')
213    
# Line 219  c--   Read optimization parameters from Line 227  c--   Read optimization parameters from
227        close( modeldataunit )        close( modeldataunit )
228    
229        rewind( scrunit1 )        rewind( scrunit1 )
230        read(unit = scrunit1, nml = ecco_optim)        read(unit = scrunit1, nml = optim)
231        close( scrunit1 )        close( scrunit1 )
232        print*, ' OPTIM_NUMBMOD: Minimization options have been read.'        print*, ' OPTIM_NUMBMOD: Minimization options have been read.'
233    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22