/[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.4 by mlosch, Fri Dec 3 01:06:33 2004 UTC revision 1.5 by heimbach, Tue Jan 25 04:42:21 2005 UTC
# Line 147  c     == end of interface == Line 147  c     == end of interface ==
147       &  xx_efluxy_file,       &  xx_efluxy_file,
148       &  xx_efluxp_file,       &  xx_efluxp_file,
149       &  xx_bottomdrag_file,       &  xx_bottomdrag_file,
150       &  doInitXX       &  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 175  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 226  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 245  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.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22