/[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.1.2.2 by heimbach, Wed May 28 21:53:58 2003 UTC revision 1.13 by heimbach, Fri Mar 11 20:01:55 2016 UTC
# Line 1  Line 1 
1    
2    #include "CTRL_OPTIONS.h"
3    
4  cph #ifdef TARGET_CRAY_VECTOR  cph #ifdef TARGET_CRAY_VECTOR
5  cph #  define NML_OPTIM ECCO_OPTIM  cph #  define NML_OPTIM ECCO_OPTIM
6  cph #  define NML_ECCO  ECCO_PARMS  cph #  define NML_ECCO  ECCO_PARMS
# Line 73  c     == global variables == Line 75  c     == global variables ==
75  #include "EEPARAMS.h"  #include "EEPARAMS.h"
76  #include "SIZE.h"  #include "SIZE.h"
77    
 #include "ecco.h"  
78  #include "ctrl.h"  #include "ctrl.h"
79  #include "optim.h"  #include "optim.h"
80  #include "minimization.h"  #include "minimization.h"
81    #include "CTRL_OBCS.h"
82    #ifndef ECCO_CTRL_DEPRECATED
83    # include "ctrl_local_params.h"
84    #endif
85    
86  c     == routine arguments ==  c     == routine arguments ==
87    
# Line 107  c     == local variables == Line 112  c     == local variables ==
112  c     == external ==  c     == external ==
113    
114        integer  ilnblnk        integer  ilnblnk
       external ilnblnk  
115    
116  c     == end of interface ==  c     == end of interface ==
117    
118        namelist /CTRL_NML/        namelist /CTRL_NML/
119       &  xx_theta_file,       &  xx_theta_file, xx_salt_file,
120       &  xx_salt_file,       &  xx_hflux_file, xx_hflux_remo_intercept, xx_hflux_remo_slope,
      &  xx_hflux_file,  
121       &  xx_hfluxstartdate1, xx_hfluxstartdate2, xx_hfluxperiod,       &  xx_hfluxstartdate1, xx_hfluxstartdate2, xx_hfluxperiod,
122       &  xx_sflux_file,       &  xx_sflux_file, xx_sflux_remo_intercept, xx_sflux_remo_slope,
123       &  xx_sfluxstartdate1, xx_sfluxstartdate2, xx_sfluxperiod,       &  xx_sfluxstartdate1, xx_sfluxstartdate2, xx_sfluxperiod,
124       &  xx_tauu_file,       &  xx_tauu_file, xx_tauu_remo_intercept, xx_tauu_remo_slope,
125       &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,       &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,
126       &  xx_tauv_file,       &  xx_tauv_file, xx_tauv_remo_intercept, xx_tauv_remo_slope,
127       &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod       &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod,
128         &  xx_atemp_file, xx_atemp_remo_intercept, xx_atemp_remo_slope,
129         &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
130         &  xx_aqh_file, xx_aqh_remo_intercept, xx_aqh_remo_slope,
131         &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,
132         &  xx_precip_file, xx_precip_remo_intercept, xx_precip_remo_slope,
133         &  xx_precipstartdate1, xx_precipstartdate2, xx_precipperiod,
134         &  xx_swflux_file, xx_swflux_remo_intercept, xx_swflux_remo_slope,
135         &  xx_swfluxstartdate1, xx_swfluxstartdate2, xx_swfluxperiod,
136         &  xx_swdown_file, xx_swdown_remo_intercept, xx_swdown_remo_slope,
137         &  xx_swdownstartdate1, xx_swdownstartdate2, xx_swdownperiod,
138         &  xx_lwflux_file, xx_lwflux_remo_intercept, xx_lwflux_remo_slope,
139         &  xx_lwfluxstartdate1, xx_lwfluxstartdate2, xx_lwfluxperiod,
140         &  xx_lwdown_file, xx_lwdown_remo_intercept, xx_lwdown_remo_slope,
141         &  xx_lwdownstartdate1, xx_lwdownstartdate2, xx_lwdownperiod,
142         &  xx_evap_file, xx_evap_remo_intercept, xx_evap_remo_slope,
143         &  xx_evapstartdate1, xx_evapstartdate2, xx_evapperiod,
144         &  xx_snowprecip_file, xx_snowprecip_remo_intercept,
145         &  xx_snowprecip_remo_slope,  xx_snowprecipperiod,
146         &  xx_snowprecipstartdate1, xx_snowprecipstartdate2,
147         &  xx_apressure_file, xx_apressure_remo_intercept,
148         &  xx_apressure_remo_slope, xx_apressureperiod,
149         &  xx_apressurestartdate1, xx_apressurestartdate2,
150         &  xx_runoff_file, xx_runoff_remo_intercept, xx_runoff_remo_slope,
151         &  xx_runoffstartdate1, xx_runoffstartdate2, xx_runoffperiod,
152         &  xx_uwind_file, xx_uwind_remo_intercept, xx_uwind_remo_slope,
153         &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
154         &  xx_vwind_file, xx_vwind_remo_intercept, xx_vwind_remo_slope,
155         &  xx_vwindstartdate1, xx_vwindstartdate2, xx_vwindperiod,
156         &  xx_obcsn_file, xx_obcss_file, xx_obcsw_file, xx_obcse_file,
157         &  xx_obcsnstartdate1,  xx_obcsnstartdate2,  xx_obcsnperiod,
158         &  xx_obcssstartdate1,  xx_obcssstartdate2,  xx_obcssperiod,
159         &  xx_obcswstartdate1,  xx_obcswstartdate2,  xx_obcswperiod,
160         &  xx_obcsestartdate1,  xx_obcsestartdate2,  xx_obcseperiod,
161         &  xx_diffkr_file, xx_kapgm_file, xx_tr1_file,
162         &  xx_sst_file, xx_sss_file,
163         &  xx_sststartdate1, xx_sststartdate2, xx_sstperiod,
164         &  xx_sssstartdate1, xx_sssstartdate2, xx_sssperiod,
165         &  xx_depth_file, xx_efluxy_file, xx_efluxp_file,
166         &  xx_bottomdrag_file, xx_edtaux_file, xx_edtauy_file,
167         &  xx_uvel_file, xx_vvel_file, xx_etan_file,
168         &  xx_shifwflx_file,
169         &  xx_shifwflx_remo_intercept, xx_shifwflx_remo_slope,
170         &  xx_shifwflxstartdate1, xx_shifwflxstartdate2, xx_shifwflxperiod,
171         &  doInitXX, doPackDiag, doZscaleUnpack, doZscalePack,
172         &  doMainUnpack, doMainPack, doAdmtlmBypassAD, delZexp,
173         &  xx_hfluxm_file, doSinglePrecTapelev
174    
175        namelist /CTRL_PACKNAMES/        namelist /CTRL_PACKNAMES/
176       &  yadmark,       &  yadmark, yctrlid, yctrlposunpack, yctrlpospack,
177       &  ctrlname, costname, scalname, maskname, metaname       &  ctrlname, costname, scalname, maskname, metaname
178    
179        namelist /ECCO_PARMS/        namelist /OPTIM/
      &                      expId  
       
       namelist /ECCO_OPTIM/  
180       &                      optimcycle,       &                      optimcycle,
181       &                      numiter, nfunc, fmin, iprint,       &                      numiter, nfunc, fmin, iprint,
182       &                      epsf, epsx, epsg,       &                      epsf, epsx, epsg,
# Line 148  c--   Preset the optimization parameters Line 194  c--   Preset the optimization parameters
194        eps        = -1.e-6        eps        = -1.e-6
195        nupdate    =  1        nupdate    =  1
196        ff         = 0.        ff         = 0.
197        expId      = 'MIT_CE_000'  cdfer      expId      = 'MIT_CE_000'
198          yctrlid    = 'MIT_CE_000'
199    
200        modeldataunit = 14        modeldataunit = 14
201        scrunit1      = 11        scrunit1      = 11
202    
 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.'  
   
203  c--   Read control parameters from file.  c--   Read control parameters from file.
204        open(unit=scrunit1,status='scratch')        open(unit=scrunit1,status='scratch')
205    
# Line 200  c--   Read control parameters from file. Line 224  c--   Read control parameters from file.
224        close( scrunit1 )        close( scrunit1 )
225        print*, ' OPTIM_NUMBMOD: Control options have been read.'        print*, ' OPTIM_NUMBMOD: Control options have been read.'
226    
227    cph(
228    cdfer      expId = yctrlid
229    cph)
230    
231  c--   Read optimization parameters from file.  c--   Read optimization parameters from file.
232        open(unit=scrunit1,status='scratch')        open(unit=scrunit1,status='scratch')
233    
# Line 219  c--   Read optimization parameters from Line 247  c--   Read optimization parameters from
247        close( modeldataunit )        close( modeldataunit )
248    
249        rewind( scrunit1 )        rewind( scrunit1 )
250        read(unit = scrunit1, nml = ecco_optim)        read(unit = scrunit1, nml = optim)
251        close( scrunit1 )        close( scrunit1 )
252        print*, ' OPTIM_NUMBMOD: Minimization options have been read.'        print*, ' OPTIM_NUMBMOD: Minimization options have been read.'
253    
Legend:
Removed from v.1.1.2.2  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22