/[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.2 by heimbach, Fri Nov 15 04:03:25 2002 UTC revision 1.6 by heimbach, Wed Sep 7 14:59:00 2005 UTC
# Line 29  c     the original code. Christian Ecker Line 29  c     the original code. Christian Ecker
29  c     to the ECCO release of the MITgcmUV in order to get the off-  c     to the ECCO release of the MITgcmUV in order to get the off-
30  c     line version going. The on-line optimisation uses a simple  c     line version going. The on-line optimisation uses a simple
31  c     example, whereas the off-line version deals with the ECCO  c     example, whereas the off-line version deals with the ECCO
32  c     release's output. The off-line version can, of course, only  c     release output. The off-line version can, of course, only
33  c     do one optimization step at a time.  c     do one optimization step at a time.
34  c  c
35  c     started:  Christian Eckert  eckert@mit.edu  15-Feb-2000  c     started:  Christian Eckert  eckert@mit.edu  15-Feb-2000
# 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 121  c     == end of interface == Line 120  c     == end of interface ==
120       &  xx_tauu_file,       &  xx_tauu_file,
121       &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,       &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,
122       &  xx_tauv_file,       &  xx_tauv_file,
123       &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod       &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod,
124         &  xx_atemp_file,
125         &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
126         &  xx_aqh_file,
127         &  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,
135         &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
136         &  xx_vwind_file,
137         &  xx_vwindstartdate1, xx_vwindstartdate2, xx_vwindperiod,
138         &  xx_obcsn_file,
139         &  xx_obcsnstartdate1,  xx_obcsnstartdate2,  xx_obcsnperiod,
140         &  xx_obcss_file,
141         &  xx_obcssstartdate1,  xx_obcssstartdate2,  xx_obcssperiod,
142         &  xx_obcsw_file,
143         &  xx_obcswstartdate1,  xx_obcswstartdate2,  xx_obcswperiod,
144         &  xx_obcse_file,
145         &  xx_obcsestartdate1,  xx_obcsestartdate2,  xx_obcseperiod,
146         &  xx_diffkr_file,
147         &  xx_kapgm_file,
148         &  xx_tr1_file,
149         &  xx_sst_file,
150         &  xx_sststartdate1, xx_sststartdate2, xx_sstperiod,
151         &  xx_sss_file,
152         &  xx_sssstartdate1, xx_sssstartdate2, xx_sssperiod,
153         &  xx_hfacc_file,
154         &  xx_efluxy_file,
155         &  xx_efluxp_file,
156         &  xx_bottomdrag_file,
157         &  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 148  c--   Preset the optimization parameters Line 186  c--   Preset the optimization parameters
186        eps        = -1.e-6        eps        = -1.e-6
187        nupdate    =  1        nupdate    =  1
188        ff         = 0.        ff         = 0.
189          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 199  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 218  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.2  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22