/[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.6 by heimbach, Wed Sep 7 14:59:00 2005 UTC
# Line 107  c     == local variables == Line 107  c     == local variables ==
107  c     == external ==  c     == external ==
108    
109        integer  ilnblnk        integer  ilnblnk
       external ilnblnk  
110    
111  c     == end of interface ==  c     == end of interface ==
112    
# Line 126  c     == end of interface == Line 125  c     == end of interface ==
125       &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,       &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
126       &  xx_aqh_file,       &  xx_aqh_file,
127       &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,       &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,
128         &  xx_precip_file,
129         &  xx_precipstartdate1, xx_precipstartdate2, xx_precipperiod,
130         &  xx_swflux_file,
131         &  xx_swfluxstartdate1, xx_swfluxstartdate2, xx_swfluxperiod,
132         &  xx_swdown_file,
133         &  xx_swdownstartdate1, xx_swdownstartdate2, xx_swdownperiod,
134       &  xx_uwind_file,       &  xx_uwind_file,
135       &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,       &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
136       &  xx_vwind_file,       &  xx_vwind_file,
# Line 142  c     == end of interface == Line 147  c     == end of interface ==
147       &  xx_kapgm_file,       &  xx_kapgm_file,
148       &  xx_tr1_file,       &  xx_tr1_file,
149       &  xx_sst_file,       &  xx_sst_file,
150         &  xx_sststartdate1, xx_sststartdate2, xx_sstperiod,
151       &  xx_sss_file,       &  xx_sss_file,
152         &  xx_sssstartdate1, xx_sssstartdate2, xx_sssperiod,
153       &  xx_hfacc_file,       &  xx_hfacc_file,
154       &  xx_efluxy_file,       &  xx_efluxy_file,
155       &  xx_efluxp_file,       &  xx_efluxp_file,
156       &  xx_bottomdrag_file,       &  xx_bottomdrag_file,
157       &  doInitXX       &  xx_edtaux_file,
158         &  xx_edtauy_file,
159         &  xx_uvel_file,
160         &  xx_vvel_file,
161         &  xx_etan_file,
162         &  doInitXX,
163         &  doPackDiag,
164         &  doZscaleUnpack,
165         &  doZscalePack
166    
167        namelist /CTRL_PACKNAMES/        namelist /CTRL_PACKNAMES/
168       &  yadmark,       &  yadmark, yctrlid, yctrlposunpack, yctrlpospack,
169       &  ctrlname, costname, scalname, maskname, metaname       &  ctrlname, costname, scalname, maskname, metaname
170    
171        namelist /ECCO_PARMS/        namelist /OPTIM/
      &                      expId  
       
       namelist /ECCO_OPTIM/  
172       &                      optimcycle,       &                      optimcycle,
173       &                      numiter, nfunc, fmin, iprint,       &                      numiter, nfunc, fmin, iprint,
174       &                      epsf, epsx, epsg,       &                      epsf, epsx, epsg,
# Line 175  c--   Preset the optimization parameters Line 187  c--   Preset the optimization parameters
187        nupdate    =  1        nupdate    =  1
188        ff         = 0.        ff         = 0.
189        expId      = 'MIT_CE_000'        expId      = 'MIT_CE_000'
190          yctrlid    = expId
191    
192        modeldataunit = 14        modeldataunit = 14
193        scrunit1      = 11        scrunit1      = 11
194    
 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.'  
   
195  c--   Read control parameters from file.  c--   Read control parameters from file.
196        open(unit=scrunit1,status='scratch')        open(unit=scrunit1,status='scratch')
197    
# Line 226  c--   Read control parameters from file. Line 216  c--   Read control parameters from file.
216        close( scrunit1 )        close( scrunit1 )
217        print*, ' OPTIM_NUMBMOD: Control options have been read.'        print*, ' OPTIM_NUMBMOD: Control options have been read.'
218    
219    cph(
220          expId = yctrlid
221    cph)
222    
223  c--   Read optimization parameters from file.  c--   Read optimization parameters from file.
224        open(unit=scrunit1,status='scratch')        open(unit=scrunit1,status='scratch')
225    
# Line 245  c--   Read optimization parameters from Line 239  c--   Read optimization parameters from
239        close( modeldataunit )        close( modeldataunit )
240    
241        rewind( scrunit1 )        rewind( scrunit1 )
242        read(unit = scrunit1, nml = ecco_optim)        read(unit = scrunit1, nml = optim)
243        close( scrunit1 )        close( scrunit1 )
244        print*, ' OPTIM_NUMBMOD: Minimization options have been read.'        print*, ' OPTIM_NUMBMOD: Minimization options have been read.'
245    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22