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