pkg/ctrl/ctrl_readparms.F

C $Header: /u/gcmpack/MITgcm/pkg/ctrl/ctrl_readparms.F,v 1.48 2014/04/04 21:19:04 jmc Exp $
C $Name:  $

#include "CTRL_OPTIONS.h"
#include "AD_CONFIG.h"

      subroutine ctrl_readparms( mythid )

c     ==================================================================
c     SUBROUTINE ctrl_readparms
c     ==================================================================
c
c     o read ctrl parameters
c
c     ==================================================================
c     SUBROUTINE ctrl_readparms
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"
#include "CTRL_SIZE.h"
#include "ctrl.h"
#include "CTRL_GENARR.h"
#ifdef ALLOW_AUTODIFF
#include "AUTODIFF_PARAMS.h"
#endif

c     == routine arguments ==

      integer mythid

c     == local variables ==

      integer bi,bj
      integer i,j,k
      integer ntmp
      integer ivarindex
      integer iUnit
      integer iobcs
      integer il
      integer errio
      integer startrec
      integer endrec
      integer difftime(4)
      _RL     diffsecs
      _RL     dummy

      character*(80)   ymaskobcs
      character*(max_len_prec) record
      character*(max_len_mbuf) msgbuf

      integer nwetc3d
      integer iarr

c     == external ==

      integer  ilnblnk
      external ilnblnk

c     == end of interface ==

c--   Read the namelist input.
      namelist /ctrl_nml/
     &  xx_theta_file, xx_salt_file,
     &  xx_hflux_file, xx_hflux_remo_intercept, xx_hflux_remo_slope,
     &  xx_hfluxstartdate1, xx_hfluxstartdate2, xx_hfluxperiod,
     &  xx_sflux_file, xx_sflux_remo_intercept, xx_sflux_remo_slope,
     &  xx_sfluxstartdate1, xx_sfluxstartdate2, xx_sfluxperiod,
     &  xx_tauu_file, xx_tauu_remo_intercept, xx_tauu_remo_slope,
     &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,
     &  xx_tauv_file, xx_tauv_remo_intercept, xx_tauv_remo_slope,
     &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod,
     &  xx_atemp_file, xx_atemp_remo_intercept, xx_atemp_remo_slope,
     &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
     &  xx_aqh_file, xx_aqh_remo_intercept, xx_aqh_remo_slope,
     &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,
     &  xx_precip_file, xx_precip_remo_intercept, xx_precip_remo_slope,
     &  xx_precipstartdate1, xx_precipstartdate2, xx_precipperiod,
     &  xx_swflux_file, xx_swflux_remo_intercept, xx_swflux_remo_slope,
     &  xx_swfluxstartdate1, xx_swfluxstartdate2, xx_swfluxperiod,
     &  xx_swdown_file, xx_swdown_remo_intercept, xx_swdown_remo_slope,
     &  xx_swdownstartdate1, xx_swdownstartdate2, xx_swdownperiod,
     &  xx_lwflux_file, xx_lwflux_remo_intercept, xx_lwflux_remo_slope,
     &  xx_lwfluxstartdate1, xx_lwfluxstartdate2, xx_lwfluxperiod,
     &  xx_lwdown_file, xx_lwdown_remo_intercept, xx_lwdown_remo_slope,
     &  xx_lwdownstartdate1, xx_lwdownstartdate2, xx_lwdownperiod,
     &  xx_evap_file, xx_evap_remo_intercept, xx_evap_remo_slope,
     &  xx_evapstartdate1, xx_evapstartdate2, xx_evapperiod,
     &  xx_snowprecip_file, xx_snowprecip_remo_intercept,
     &  xx_snowprecip_remo_slope,  xx_snowprecipperiod,
     &  xx_snowprecipstartdate1, xx_snowprecipstartdate2,
     &  xx_apressure_file, xx_apressure_remo_intercept,
     &  xx_apressure_remo_slope, xx_apressureperiod,
     &  xx_apressurestartdate1, xx_apressurestartdate2,
     &  xx_runoff_file, xx_runoff_remo_intercept, xx_runoff_remo_slope,
     &  xx_runoffstartdate1, xx_runoffstartdate2, xx_runoffperiod,
     &  xx_uwind_file, xx_uwind_remo_intercept, xx_uwind_remo_slope,
     &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
     &  xx_vwind_file, xx_vwind_remo_intercept, xx_vwind_remo_slope,
     &  xx_vwindstartdate1, xx_vwindstartdate2, xx_vwindperiod,
     &  xx_aqh_mean_file, xx_atemp_mean_file, xx_precip_mean_file,
     &  xx_swdown_mean_file, xx_uwind_mean_file, xx_vwind_mean_file,
     &  xx_obcsn_file,
     &  xx_obcsnstartdate1,  xx_obcsnstartdate2,  xx_obcsnperiod,
     &  xx_obcss_file,
     &  xx_obcssstartdate1,  xx_obcssstartdate2,  xx_obcssperiod,
     &  xx_obcsw_file,
     &  xx_obcswstartdate1,  xx_obcswstartdate2,  xx_obcswperiod,
     &  xx_obcse_file,
     &  xx_obcsestartdate1,  xx_obcsestartdate2,  xx_obcseperiod,
     &  xx_diffkr_file, xx_kapgm_file, xx_kapredi_file, xx_tr1_file,
     &  xx_sst_file, xx_sst_remo_intercept, xx_sst_remo_slope,
     &  xx_sss_file, xx_sss_remo_intercept, xx_sss_remo_slope,
     &  xx_sststartdate1, xx_sststartdate2, xx_sstperiod,
     &  xx_sssstartdate1, xx_sssstartdate2, xx_sssperiod,
     &  xx_depth_file, xx_gen2d_file, xx_gen3d_file,
     &  xx_efluxy_file, xx_efluxp_file, xx_bottomdrag_file,
     &  xx_edtaux_file, xx_edtauy_file,
     &  xx_uvel_file, xx_vvel_file, xx_etan_file,
     &  xx_siarea_file, xx_siheff_file, xx_sihsnow_file,
     &  xx_shifwflx_file,
     &  xx_shifwflx_remo_intercept, xx_shifwflx_remo_slope,
     &  xx_shifwflxstartdate1, xx_shifwflxstartdate2, xx_shifwflxperiod,
     &  doInitXX, ctrlSmoothCorrel3D, ctrlSmoothCorrel2D,
     &  doPackDiag,
     &  doZscaleUnpack, doZscalePack,
     &  doMainUnpack, doMainPack,
     &  doSinglePrecTapelev,
     &  doAdmtlmBypassAD,
     &  delZexp, forcingPrecond
cHFLUXM_CONTROL
     &  ,xx_hfluxm_file
cHFLUXM_CONTROL

      namelist /ctrl_packnames/
     &  yadmark, yctrlid, yctrlposunpack, yctrlpospack,
     &  ctrlname, costname, scalname, maskname, metaname

#if (defined (ALLOW_GENARR2D_CONTROL) || defined (ALLOW_GENARR3D_CONTROL) || defined (ALLOW_GENTIM2D_CONTROL))
      namelist /ctrl_nml_genarr/
     &  xx_genarr2d_file, xx_genarr3d_file, xx_gentim2d_file,
     &  xx_genarr2d_weight, xx_genarr3d_weight, xx_gentim2d_weight,
     &  xx_gentim2d_startdate1, xx_gentim2d_startdate2,
     &  xx_gentim2d_startdate, xx_gentim2d_period,
     &  xx_gentim2d_cumsum,xx_gentim2d_glosum,
     &  genarr2dPrecond, gentim2dPrecond,mult_gentim2d
#endif

      _BEGIN_MASTER( myThid )

c--     Set default values.
        doInitXX            = .TRUE.
#ifdef ALLOW_SMOOTH_CORREL3D
        ctrlSmoothCorrel3D  = .TRUE.
#else
        ctrlSmoothCorrel3D  = .FALSE.
#endif
#ifdef ALLOW_SMOOTH_CORREL2D
        ctrlSmoothCorrel2D  = .TRUE.
#else
        ctrlSmoothCorrel2D  = .FALSE.
#endif
#ifdef ALLOW_ADMTLM
        doAdmTlm            = .TRUE.
#else
        doAdmTlm            = .FALSE.
#endif
        doPackDiag          = .FALSE.
        doZscaleUnpack      = .FALSE.
        doZscalePack        = .FALSE.
        doMainUnpack        = .TRUE.
#if (defined (ALLOW_ADJOINT_RUN)||defined (ALLOW_TANGENTLINEAR_RUN))
        doMainPack          = .TRUE.
#else
        doMainPack          = .FALSE.
#endif
        doSinglePrecTapelev = .FALSE.
        doAdmtlmBypassAD    = .FALSE.

        delZexp = 0.
        forcingPrecond = 1. _d 0

        xx_theta_file      = 'xx_theta'
        xx_salt_file       = 'xx_salt'
c
        xx_gen2d_file      = 'xx_gen2d'
        xx_gen3d_file      = 'xx_gen3d'

#if (defined (ALLOW_GENARR2D_CONTROL) || defined (ALLOW_GENARR3D_CONTROL) || defined (ALLOW_GENTIM2D_CONTROL))
        do iarr = 1, maxCtrlArr2D
         write(xx_genarr2d_file(iarr)(1:14),'(a,i2.2)')
     &        'xx_genarr2d_', iarr
         write(xx_genarr2d_file(iarr)(15:MAX_LEN_FNAM),'(a)') ' '
         write(xx_genarr2d_weight(iarr)(1:MAX_LEN_FNAM),'(a)') ' '
         genarr2dPrecond(iarr) = 1. _d 0
        enddo
c
        do iarr = 1, maxCtrlArr3D
         write(xx_genarr3d_file(iarr)(1:14),'(a,i2.2)')
     &        'xx_genarr3d_', iarr
         write(xx_genarr3d_file(iarr)(15:MAX_LEN_FNAM),'(a)') ' '
         write(xx_genarr3d_weight(iarr)(1:MAX_LEN_FNAM),'(a)') ' '
        enddo
c
        do iarr = 1, maxCtrlTim2D
         write(xx_gentim2d_file(iarr)(1:14),'(a,i2.2)')
     &        'xx_gentim2d_', iarr
         write(xx_gentim2d_file(iarr)(15:MAX_LEN_FNAM),'(a)') ' '
         write(xx_gentim2d_weight(iarr)(1:MAX_LEN_FNAM),'(a)') ' '
         xx_gentim2d_startdate1(iarr) = 0
         xx_gentim2d_startdate2(iarr) = 0
         xx_gentim2d_period(iarr)     = 0. _d 0
         xx_gentim2d_cumsum(iarr)     = .FALSE.
         xx_gentim2d_glosum(iarr)     = .FALSE.
         gentim2dPrecond(iarr)        = 1. _d 0
         mult_gentim2d(iarr)          = 0. _d 0
        enddo
#endif
c
        xx_hfluxstartdate1 =   0
        xx_hfluxstartdate2 =   0
        xx_hfluxperiod     =   0. _d 0
        xx_hflux_file      = 'xx_hfl'
        xx_hflux_remo_intercept =   0. _d 0
        xx_hflux_remo_slope =   0. _d 0
c
        xx_sfluxstartdate1 =   0
        xx_sfluxstartdate2 =   0
        xx_sfluxperiod     =   0. _d 0
        xx_sflux_file      = 'xx_sfl'
        xx_sflux_remo_intercept =   0. _d 0
        xx_sflux_remo_slope =   0. _d 0
c
        xx_tauustartdate1  =   0
        xx_tauustartdate2  =   0
        xx_tauuperiod      =   0. _d 0
        xx_tauu_file       = 'xx_tauu'
        xx_tauu_remo_intercept =   0. _d 0
        xx_tauu_remo_slope =   0. _d 0
c
        xx_tauvstartdate1  =   0
        xx_tauvstartdate2  =   0
        xx_tauvperiod      =   0. _d 0
        xx_tauv_file       = 'xx_tauv'
        xx_tauv_remo_intercept =   0. _d 0
        xx_tauv_remo_slope =   0. _d 0
c
        xx_atempstartdate1  =   0
        xx_atempstartdate2  =   0
        xx_atempperiod      =   0. _d 0
        xx_atemp_file       = 'xx_atemp'
        xx_atemp_mean_file  = 'xx_atemp_mean'
        xx_atemp_remo_intercept =   0. _d 0
        xx_atemp_remo_slope =   0. _d 0
c
        xx_aqhstartdate1    =   0
        xx_aqhstartdate2    =   0
        xx_aqhperiod        =   0. _d 0
        xx_aqh_file         = 'xx_aqh'
        xx_aqh_mean_file    = 'xx_aqh_mean'
        xx_aqh_remo_intercept =   0. _d 0
        xx_aqh_remo_slope =   0. _d 0
c
        xx_precipstartdate1 =   0
        xx_precipstartdate2 =   0
        xx_precipperiod     =   0. _d 0
        xx_precip_file      = 'xx_precip'
        xx_precip_mean_file = 'xx_precip_mean'
        xx_precip_remo_intercept =   0. _d 0
        xx_precip_remo_slope =   0. _d 0
c
        xx_swfluxstartdate1 =   0
        xx_swfluxstartdate2 =   0
        xx_swfluxperiod     =   0. _d 0
        xx_swflux_file      = 'xx_swflux'
        xx_swflux_remo_intercept =   0. _d 0
        xx_swflux_remo_slope =   0. _d 0
c
        xx_swdownstartdate1 =   0
        xx_swdownstartdate2 =   0
        xx_swdownperiod     =   0. _d 0
        xx_swdown_file      = 'xx_swdown'
        xx_swdown_mean_file = 'xx_swdown_mean'
        xx_swdown_remo_intercept =   0. _d 0
        xx_swdown_remo_slope =   0. _d 0
c
        xx_lwfluxstartdate1 =   0
        xx_lwfluxstartdate2 =   0
        xx_lwfluxperiod     =   0. _d 0
        xx_lwflux_file      = 'xx_lwflux'
        xx_lwflux_remo_intercept =   0. _d 0
        xx_lwflux_remo_slope =   0. _d 0
c
        xx_lwdownstartdate1 =   0
        xx_lwdownstartdate2 =   0
        xx_lwdownperiod     =   0. _d 0
        xx_lwdown_file      = 'xx_lwdown'
        xx_lwdown_remo_intercept =   0. _d 0
        xx_lwdown_remo_slope =   0. _d 0
c
        xx_evapstartdate1 =   0
        xx_evapstartdate2 =   0
        xx_evapperiod     =   0. _d 0
        xx_evap_file      = 'xx_evap'
        xx_evap_remo_intercept =   0. _d 0
        xx_evap_remo_slope =   0. _d 0
c
        xx_snowprecipstartdate1 =   0
        xx_snowprecipstartdate2 =   0
        xx_snowprecipperiod     =   0. _d 0
        xx_snowprecip_file      = 'xx_snowprecip'
        xx_snowprecip_remo_intercept =   0. _d 0
        xx_snowprecip_remo_slope =   0. _d 0
c
        xx_apressurestartdate1 =   0
        xx_apressurestartdate2 =   0
        xx_apressureperiod     =   0. _d 0
        xx_apressure_file      = 'xx_apressure'
        xx_apressure_remo_intercept =   0. _d 0
        xx_apressure_remo_slope =   0. _d 0
c
        xx_runoffstartdate1 =   0
        xx_runoffstartdate2 =   0
        xx_runoffperiod     =   0. _d 0
        xx_runoff_file      = 'xx_runoff'
        xx_runoff_remo_intercept =   0. _d 0
        xx_runoff_remo_slope =   0. _d 0
c
        xx_uwindstartdate1  =   0
        xx_uwindstartdate2  =   0
        xx_uwindperiod      =   0. _d 0
        xx_uwind_file       = 'xx_uwind'
        xx_uwind_mean_file  = 'xx_uwind_mean'
        xx_uwind_remo_intercept =   0. _d 0
        xx_uwind_remo_slope =   0. _d 0
c
        xx_vwindstartdate1  =   0
        xx_vwindstartdate2  =   0
        xx_vwindperiod      =   0. _d 0
        xx_vwind_file       = 'xx_vwind'
        xx_vwind_mean_file  = 'xx_vwind_mean'
        xx_vwind_remo_intercept =   0. _d 0
        xx_vwind_remo_slope =   0. _d 0
c
        xx_obcsnstartdate1  =   0
        xx_obcsnstartdate2  =   0
        xx_obcsnperiod      =   0. _d 0
        xx_obcsn_file       = 'xx_obcsn'
c
        xx_obcssstartdate1  =   0
        xx_obcssstartdate2  =   0
        xx_obcssperiod      =   0. _d 0
        xx_obcss_file       = 'xx_obcss'
c
        xx_obcswstartdate1  =   0
        xx_obcswstartdate2  =   0
        xx_obcswperiod      =   0. _d 0
        xx_obcsw_file       = 'xx_obcsw'
c
        xx_obcsestartdate1  =   0
        xx_obcsestartdate2  =   0
        xx_obcseperiod      =   0. _d 0
        xx_obcse_file       = 'xx_obcse'
c
        xx_sststartdate1    =   0
        xx_sststartdate2    =   0
        xx_sstperiod        =   0. _d 0
        xx_sst_file         = 'xx_sst'
        xx_sst_remo_intercept =   0. _d 0
        xx_sst_remo_slope =   0. _d 0
c
        xx_sssstartdate1    =   0
        xx_sssstartdate2    =   0
        xx_sssperiod        =   0. _d 0
        xx_sss_file         = 'xx_sss'
        xx_sss_remo_intercept =   0. _d 0
        xx_sss_remo_slope =   0. _d 0
c
        xx_shifwflxstartdate1  =   0
        xx_shifwflxstartdate2  =   0
        xx_shifwflxperiod      =   0. _d 0
        xx_shifwflx_file       = 'xx_shifwflx'
        xx_shifwflx_remo_intercept=0. _d 0
        xx_shifwflx_remo_slope =   0. _d 0
c
        xx_diffkr_file      = 'xx_diffkr'
        xx_kapgm_file       = 'xx_kapgm'
        xx_kapredi_file     = 'xx_kapredi'
        xx_tr1_file         = 'xx_ptr'
        xx_depth_file       = 'xx_depth'
        xx_efluxy_file      = 'xx_efluxy'
        xx_efluxp_file      = 'xx_efluxp'
        xx_bottomdrag_file  = 'xx_bottomdrag'
        xx_edtaux_file      = 'xx_edtaux'
        xx_edtauy_file      = 'xx_edtauy'
        xx_uvel_file        = 'xx_uvel'
        xx_vvel_file        = 'xx_vvel'
        xx_etan_file        = 'xx_etan'
        xx_siarea_file      = 'xx_siarea'
        xx_siheff_file      = 'xx_siheff'
        xx_sihsnow_file     = 'xx_sihsnow'
cHFLUXM_CONTROL
        xx_hfluxm_file      = 'xx_hfluxm'
cHFLUXM_CONTROL

#ifdef ALLOW_TANGENTLINEAR_RUN
        yadprefix           = 'g_'
        yadmark             = 'g_'
#else
        yadprefix           = 'ad'
        yadmark             = 'ad'
#endif
        yctrlid             = 'MIT_CE_000'
        yctrlposunpack      = '.opt'
        yctrlpospack        = '.opt'
        ctrlname            = 'ecco_ctrl'
        costname            = 'ecco_cost'
        scalname            = ' '
        maskname            = ' '
        metaname            = ' '

c--     Next, read the cost data file.
        WRITE(msgBuf,'(A)') 'CTRL_READPARMS: opening data.ctrl'
        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT , 1)

        CALL OPEN_COPY_DATA_FILE(
     I                          'data.ctrl', 'CTRL_READPARMS',
     O                          iUnit,
     I                          myThid )

        READ(unit = iUnit, nml = ctrl_nml)
        READ(unit = iUnit, nml = ctrl_packnames)
#if (defined (ALLOW_GENARR2D_CONTROL) || defined (ALLOW_GENARR3D_CONTROL) || defined (ALLOW_GENTIM2D_CONTROL))
        READ(unit = iUnit, nml = ctrl_nml_genarr)
#endif

        WRITE(msgBuf,'(A)')
     &       'CTRL_READPARMS: finished reading data.ctrl'
        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                  SQUEEZE_RIGHT , 1)

        CLOSE( iUnit )

#ifdef ALLOW_AUTODIFF
c allow for switching off correl2d in adjoint
      useSmoothCorrel2DinFwdMode = ctrlSmoothCorrel2D
      useSmoothCorrel2DinAdMode =
     &         useSmoothCorrel2DinAdMode .AND. ctrlSmoothCorrel2D
       CALL WRITE_0D_L( useSmoothCorrel2DinAdMode, INDEX_NONE,
     &   'useSmoothCorrel2DinAdMode =',
     &   ' /* use ctrlSmoothCorrel2D in adjoint mode */')
#endif

      _END_MASTER( myThid )

      _BARRIER

      return
      end
1	gforget	1.49	C $Header: /u/gcmpack/MITgcm/pkg/ctrl/ctrl_readparms.F,v 1.48 2014/04/04 21:19:04 jmc Exp $
2	heimbach	1.10	C $Name: $
3	heimbach	1.2
4	jmc	1.40	#include "CTRL_OPTIONS.h"
5	jmc	1.37	#include "AD_CONFIG.h"
6	heimbach	1.2
7			subroutine ctrl_readparms( mythid )
8
9			c ==================================================================
10			c SUBROUTINE ctrl_readparms
11			c ==================================================================
12			c
13			c o read ctrl parameters
14			c
15			c ==================================================================
16			c SUBROUTINE ctrl_readparms
17			c ==================================================================
18
19			implicit none
20
21			c == global variables ==
22
23			#include "EEPARAMS.h"
24			#include "SIZE.h"
25			#include "PARAMS.h"
26	heimbach	1.39	#include "CTRL_SIZE.h"
27	heimbach	1.2	#include "ctrl.h"
28	heimbach	1.39	#include "CTRL_GENARR.h"
29	gforget	1.49	#ifdef ALLOW_AUTODIFF
30			#include "AUTODIFF_PARAMS.h"
31			#endif
32	heimbach	1.2
33			c == routine arguments ==
34
35			integer mythid
36
37			c == local variables ==
38
39			integer bi,bj
40			integer i,j,k
41			integer ntmp
42			integer ivarindex
43	heimbach	1.15	integer iUnit
44	heimbach	1.2	integer iobcs
45			integer il
46			integer errio
47			integer startrec
48			integer endrec
49			integer difftime(4)
50			_RL diffsecs
51			_RL dummy
52
53			character*(80) ymaskobcs
54			character*(max_len_prec) record
55			character*(max_len_mbuf) msgbuf
56
57			integer nwetc3d
58	heimbach	1.38	integer iarr
59	heimbach	1.2
60			c == external ==
61
62			integer ilnblnk
63			external ilnblnk
64
65			c == end of interface ==
66
67			c-- Read the namelist input.
68			namelist /ctrl_nml/
69	heimbach	1.25	& xx_theta_file, xx_salt_file,
70			& xx_hflux_file, xx_hflux_remo_intercept, xx_hflux_remo_slope,
71	heimbach	1.2	& xx_hfluxstartdate1, xx_hfluxstartdate2, xx_hfluxperiod,
72	heimbach	1.25	& xx_sflux_file, xx_sflux_remo_intercept, xx_sflux_remo_slope,
73	heimbach	1.2	& xx_sfluxstartdate1, xx_sfluxstartdate2, xx_sfluxperiod,
74	heimbach	1.25	& xx_tauu_file, xx_tauu_remo_intercept, xx_tauu_remo_slope,
75	heimbach	1.2	& xx_tauustartdate1, xx_tauustartdate2, xx_tauuperiod,
76	heimbach	1.25	& xx_tauv_file, xx_tauv_remo_intercept, xx_tauv_remo_slope,
77	heimbach	1.2	& xx_tauvstartdate1, xx_tauvstartdate2, xx_tauvperiod,
78	heimbach	1.25	& xx_atemp_file, xx_atemp_remo_intercept, xx_atemp_remo_slope,
79	heimbach	1.2	& xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
80	heimbach	1.25	& xx_aqh_file, xx_aqh_remo_intercept, xx_aqh_remo_slope,
81	heimbach	1.2	& xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,
82	heimbach	1.25	& xx_precip_file, xx_precip_remo_intercept, xx_precip_remo_slope,
83	heimbach	1.16	& xx_precipstartdate1, xx_precipstartdate2, xx_precipperiod,
84	heimbach	1.25	& xx_swflux_file, xx_swflux_remo_intercept, xx_swflux_remo_slope,
85	heimbach	1.17	& xx_swfluxstartdate1, xx_swfluxstartdate2, xx_swfluxperiod,
86	heimbach	1.25	& xx_swdown_file, xx_swdown_remo_intercept, xx_swdown_remo_slope,
87	heimbach	1.18	& xx_swdownstartdate1, xx_swdownstartdate2, xx_swdownperiod,
88	heimbach	1.25	& xx_lwflux_file, xx_lwflux_remo_intercept, xx_lwflux_remo_slope,
89			& xx_lwfluxstartdate1, xx_lwfluxstartdate2, xx_lwfluxperiod,
90			& xx_lwdown_file, xx_lwdown_remo_intercept, xx_lwdown_remo_slope,
91			& xx_lwdownstartdate1, xx_lwdownstartdate2, xx_lwdownperiod,
92			& xx_evap_file, xx_evap_remo_intercept, xx_evap_remo_slope,
93			& xx_evapstartdate1, xx_evapstartdate2, xx_evapperiod,
94	jmc	1.34	& xx_snowprecip_file, xx_snowprecip_remo_intercept,
95	heimbach	1.25	& xx_snowprecip_remo_slope, xx_snowprecipperiod,
96			& xx_snowprecipstartdate1, xx_snowprecipstartdate2,
97	jmc	1.34	& xx_apressure_file, xx_apressure_remo_intercept,
98	heimbach	1.25	& xx_apressure_remo_slope, xx_apressureperiod,
99			& xx_apressurestartdate1, xx_apressurestartdate2,
100			& xx_runoff_file, xx_runoff_remo_intercept, xx_runoff_remo_slope,
101			& xx_runoffstartdate1, xx_runoffstartdate2, xx_runoffperiod,
102			& xx_uwind_file, xx_uwind_remo_intercept, xx_uwind_remo_slope,
103	heimbach	1.2	& xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
104	heimbach	1.25	& xx_vwind_file, xx_vwind_remo_intercept, xx_vwind_remo_slope,
105	heimbach	1.2	& xx_vwindstartdate1, xx_vwindstartdate2, xx_vwindperiod,
106	heimbach	1.36	& xx_aqh_mean_file, xx_atemp_mean_file, xx_precip_mean_file,
107			& xx_swdown_mean_file, xx_uwind_mean_file, xx_vwind_mean_file,
108	heimbach	1.2	& xx_obcsn_file,
109			& xx_obcsnstartdate1, xx_obcsnstartdate2, xx_obcsnperiod,
110			& xx_obcss_file,
111			& xx_obcssstartdate1, xx_obcssstartdate2, xx_obcssperiod,
112			& xx_obcsw_file,
113			& xx_obcswstartdate1, xx_obcswstartdate2, xx_obcswperiod,
114			& xx_obcse_file,
115			& xx_obcsestartdate1, xx_obcsestartdate2, xx_obcseperiod,
116	gforget	1.29	& xx_diffkr_file, xx_kapgm_file, xx_kapredi_file, xx_tr1_file,
117	heimbach	1.35	& xx_sst_file, xx_sst_remo_intercept, xx_sst_remo_slope,
118			& xx_sss_file, xx_sss_remo_intercept, xx_sss_remo_slope,
119	heimbach	1.19	& xx_sststartdate1, xx_sststartdate2, xx_sstperiod,
120			& xx_sssstartdate1, xx_sssstartdate2, xx_sssperiod,
121	jmc	1.40	& xx_depth_file, xx_gen2d_file, xx_gen3d_file,
122	heimbach	1.38	& xx_efluxy_file, xx_efluxp_file, xx_bottomdrag_file,
123	heimbach	1.25	& xx_edtaux_file, xx_edtauy_file,
124			& xx_uvel_file, xx_vvel_file, xx_etan_file,
125	heimbach	1.27	& xx_siarea_file, xx_siheff_file, xx_sihsnow_file,
126	jmc	1.34	& xx_shifwflx_file,
127			& xx_shifwflx_remo_intercept, xx_shifwflx_remo_slope,
128	mlosch	1.33	& xx_shifwflxstartdate1, xx_shifwflxstartdate2, xx_shifwflxperiod,
129	gforget	1.41	& doInitXX, ctrlSmoothCorrel3D, ctrlSmoothCorrel2D,
130	heimbach	1.11	& doPackDiag,
131	heimbach	1.25	& doZscaleUnpack, doZscalePack,
132			& doMainUnpack, doMainPack,
133	heimbach	1.31	& doSinglePrecTapelev,
134	heimbach	1.23	& doAdmtlmBypassAD,
135	gforget	1.26	& delZexp, forcingPrecond
136	dfer	1.28	cHFLUXM_CONTROL
137			& ,xx_hfluxm_file
138			cHFLUXM_CONTROL
139	heimbach	1.2
140			namelist /ctrl_packnames/
141	heimbach	1.6	& yadmark, yctrlid, yctrlposunpack, yctrlpospack,
142	heimbach	1.2	& ctrlname, costname, scalname, maskname, metaname
143
144	heimbach	1.43	#if (defined (ALLOW_GENARR2D_CONTROL) \|\| defined (ALLOW_GENARR3D_CONTROL) \|\| defined (ALLOW_GENTIM2D_CONTROL))
145			namelist /ctrl_nml_genarr/
146			& xx_genarr2d_file, xx_genarr3d_file, xx_gentim2d_file,
147	heimbach	1.44	& xx_genarr2d_weight, xx_genarr3d_weight, xx_gentim2d_weight,
148	heimbach	1.43	& xx_gentim2d_startdate1, xx_gentim2d_startdate2,
149	heimbach	1.44	& xx_gentim2d_startdate, xx_gentim2d_period,
150	gforget	1.46	& xx_gentim2d_cumsum,xx_gentim2d_glosum,
151	gforget	1.45	& genarr2dPrecond, gentim2dPrecond,mult_gentim2d
152	heimbach	1.43	#endif
153
154	heimbach	1.2	_BEGIN_MASTER( myThid )
155
156			c-- Set default values.
157	heimbach	1.8	doInitXX = .TRUE.
158	jmc	1.42	#ifdef ALLOW_SMOOTH_CORREL3D
159	gforget	1.41	ctrlSmoothCorrel3D = .TRUE.
160			#else
161			ctrlSmoothCorrel3D = .FALSE.
162			#endif
163	jmc	1.42	#ifdef ALLOW_SMOOTH_CORREL2D
164	gforget	1.41	ctrlSmoothCorrel2D = .TRUE.
165			#else
166			ctrlSmoothCorrel2D = .FALSE.
167			#endif
168	heimbach	1.8	#ifdef ALLOW_ADMTLM
169			doAdmTlm = .TRUE.
170			#else
171			doAdmTlm = .FALSE.
172			#endif
173	heimbach	1.10	doPackDiag = .FALSE.
174	heimbach	1.11	doZscaleUnpack = .FALSE.
175			doZscalePack = .FALSE.
176	heimbach	1.20	doMainUnpack = .TRUE.
177	jmc	1.48	#if (defined (ALLOW_ADJOINT_RUN)\|\|defined (ALLOW_TANGENTLINEAR_RUN))
178	heimbach	1.20	doMainPack = .TRUE.
179	jmc	1.42	#else
180			doMainPack = .FALSE.
181			#endif
182	heimbach	1.31	doSinglePrecTapelev = .FALSE.
183	heimbach	1.21	doAdmtlmBypassAD = .FALSE.
184	heimbach	1.8
185	heimbach	1.23	delZexp = 0.
186	gforget	1.26	forcingPrecond = 1. _d 0
187	heimbach	1.23
188	heimbach	1.13	xx_theta_file = 'xx_theta'
189			xx_salt_file = 'xx_salt'
190	heimbach	1.22	c
191	heimbach	1.32	xx_gen2d_file = 'xx_gen2d'
192			xx_gen3d_file = 'xx_gen3d'
193	heimbach	1.43
194			#if (defined (ALLOW_GENARR2D_CONTROL) \|\| defined (ALLOW_GENARR3D_CONTROL) \|\| defined (ALLOW_GENTIM2D_CONTROL))
195	heimbach	1.38	do iarr = 1, maxCtrlArr2D
196			write(xx_genarr2d_file(iarr)(1:14),'(a,i2.2)')
197			& 'xx_genarr2d_', iarr
198			write(xx_genarr2d_file(iarr)(15:MAX_LEN_FNAM),'(a)') ' '
199	heimbach	1.44	write(xx_genarr2d_weight(iarr)(1:MAX_LEN_FNAM),'(a)') ' '
200			genarr2dPrecond(iarr) = 1. _d 0
201	heimbach	1.38	enddo
202			c
203			do iarr = 1, maxCtrlArr3D
204			write(xx_genarr3d_file(iarr)(1:14),'(a,i2.2)')
205			& 'xx_genarr3d_', iarr
206			write(xx_genarr3d_file(iarr)(15:MAX_LEN_FNAM),'(a)') ' '
207	heimbach	1.44	write(xx_genarr3d_weight(iarr)(1:MAX_LEN_FNAM),'(a)') ' '
208	heimbach	1.38	enddo
209	heimbach	1.32	c
210	heimbach	1.43	do iarr = 1, maxCtrlTim2D
211			write(xx_gentim2d_file(iarr)(1:14),'(a,i2.2)')
212			& 'xx_gentim2d_', iarr
213			write(xx_gentim2d_file(iarr)(15:MAX_LEN_FNAM),'(a)') ' '
214	heimbach	1.44	write(xx_gentim2d_weight(iarr)(1:MAX_LEN_FNAM),'(a)') ' '
215			xx_gentim2d_startdate1(iarr) = 0
216			xx_gentim2d_startdate2(iarr) = 0
217			xx_gentim2d_period(iarr) = 0. _d 0
218	gforget	1.46	xx_gentim2d_cumsum(iarr) = .FALSE.
219			xx_gentim2d_glosum(iarr) = .FALSE.
220	heimbach	1.44	gentim2dPrecond(iarr) = 1. _d 0
221	gforget	1.45	mult_gentim2d(iarr) = 0. _d 0
222	heimbach	1.43	enddo
223			#endif
224			c
225	heimbach	1.2	xx_hfluxstartdate1 = 0
226			xx_hfluxstartdate2 = 0
227			xx_hfluxperiod = 0. _d 0
228	heimbach	1.13	xx_hflux_file = 'xx_hfl'
229	heimbach	1.22	xx_hflux_remo_intercept = 0. _d 0
230			xx_hflux_remo_slope = 0. _d 0
231			c
232	heimbach	1.2	xx_sfluxstartdate1 = 0
233			xx_sfluxstartdate2 = 0
234			xx_sfluxperiod = 0. _d 0
235	heimbach	1.13	xx_sflux_file = 'xx_sfl'
236	heimbach	1.22	xx_sflux_remo_intercept = 0. _d 0
237			xx_sflux_remo_slope = 0. _d 0
238			c
239	heimbach	1.2	xx_tauustartdate1 = 0
240			xx_tauustartdate2 = 0
241			xx_tauuperiod = 0. _d 0
242	heimbach	1.13	xx_tauu_file = 'xx_tauu'
243	heimbach	1.22	xx_tauu_remo_intercept = 0. _d 0
244			xx_tauu_remo_slope = 0. _d 0
245			c
246	heimbach	1.2	xx_tauvstartdate1 = 0
247			xx_tauvstartdate2 = 0
248			xx_tauvperiod = 0. _d 0
249	heimbach	1.13	xx_tauv_file = 'xx_tauv'
250	heimbach	1.22	xx_tauv_remo_intercept = 0. _d 0
251			xx_tauv_remo_slope = 0. _d 0
252			c
253	heimbach	1.2	xx_atempstartdate1 = 0
254			xx_atempstartdate2 = 0
255			xx_atempperiod = 0. _d 0
256	heimbach	1.13	xx_atemp_file = 'xx_atemp'
257	heimbach	1.36	xx_atemp_mean_file = 'xx_atemp_mean'
258	heimbach	1.22	xx_atemp_remo_intercept = 0. _d 0
259			xx_atemp_remo_slope = 0. _d 0
260			c
261	heimbach	1.2	xx_aqhstartdate1 = 0
262			xx_aqhstartdate2 = 0
263			xx_aqhperiod = 0. _d 0
264	heimbach	1.13	xx_aqh_file = 'xx_aqh'
265	heimbach	1.36	xx_aqh_mean_file = 'xx_aqh_mean'
266	heimbach	1.22	xx_aqh_remo_intercept = 0. _d 0
267			xx_aqh_remo_slope = 0. _d 0
268			c
269	heimbach	1.16	xx_precipstartdate1 = 0
270			xx_precipstartdate2 = 0
271			xx_precipperiod = 0. _d 0
272			xx_precip_file = 'xx_precip'
273	heimbach	1.36	xx_precip_mean_file = 'xx_precip_mean'
274	heimbach	1.22	xx_precip_remo_intercept = 0. _d 0
275			xx_precip_remo_slope = 0. _d 0
276			c
277	heimbach	1.17	xx_swfluxstartdate1 = 0
278			xx_swfluxstartdate2 = 0
279			xx_swfluxperiod = 0. _d 0
280			xx_swflux_file = 'xx_swflux'
281	heimbach	1.22	xx_swflux_remo_intercept = 0. _d 0
282			xx_swflux_remo_slope = 0. _d 0
283			c
284	heimbach	1.18	xx_swdownstartdate1 = 0
285			xx_swdownstartdate2 = 0
286			xx_swdownperiod = 0. _d 0
287			xx_swdown_file = 'xx_swdown'
288	heimbach	1.36	xx_swdown_mean_file = 'xx_swdown_mean'
289	heimbach	1.22	xx_swdown_remo_intercept = 0. _d 0
290			xx_swdown_remo_slope = 0. _d 0
291			c
292	heimbach	1.25	xx_lwfluxstartdate1 = 0
293			xx_lwfluxstartdate2 = 0
294			xx_lwfluxperiod = 0. _d 0
295			xx_lwflux_file = 'xx_lwflux'
296			xx_lwflux_remo_intercept = 0. _d 0
297			xx_lwflux_remo_slope = 0. _d 0
298			c
299			xx_lwdownstartdate1 = 0
300			xx_lwdownstartdate2 = 0
301			xx_lwdownperiod = 0. _d 0
302			xx_lwdown_file = 'xx_lwdown'
303			xx_lwdown_remo_intercept = 0. _d 0
304			xx_lwdown_remo_slope = 0. _d 0
305			c
306			xx_evapstartdate1 = 0
307			xx_evapstartdate2 = 0
308			xx_evapperiod = 0. _d 0
309			xx_evap_file = 'xx_evap'
310			xx_evap_remo_intercept = 0. _d 0
311			xx_evap_remo_slope = 0. _d 0
312			c
313			xx_snowprecipstartdate1 = 0
314			xx_snowprecipstartdate2 = 0
315			xx_snowprecipperiod = 0. _d 0
316			xx_snowprecip_file = 'xx_snowprecip'
317			xx_snowprecip_remo_intercept = 0. _d 0
318			xx_snowprecip_remo_slope = 0. _d 0
319			c
320			xx_apressurestartdate1 = 0
321			xx_apressurestartdate2 = 0
322			xx_apressureperiod = 0. _d 0
323			xx_apressure_file = 'xx_apressure'
324			xx_apressure_remo_intercept = 0. _d 0
325			xx_apressure_remo_slope = 0. _d 0
326			c
327			xx_runoffstartdate1 = 0
328			xx_runoffstartdate2 = 0
329			xx_runoffperiod = 0. _d 0
330			xx_runoff_file = 'xx_runoff'
331			xx_runoff_remo_intercept = 0. _d 0
332			xx_runoff_remo_slope = 0. _d 0
333			c
334	heimbach	1.2	xx_uwindstartdate1 = 0
335			xx_uwindstartdate2 = 0
336			xx_uwindperiod = 0. _d 0
337	heimbach	1.13	xx_uwind_file = 'xx_uwind'
338	heimbach	1.36	xx_uwind_mean_file = 'xx_uwind_mean'
339	heimbach	1.22	xx_uwind_remo_intercept = 0. _d 0
340			xx_uwind_remo_slope = 0. _d 0
341			c
342	heimbach	1.2	xx_vwindstartdate1 = 0
343			xx_vwindstartdate2 = 0
344			xx_vwindperiod = 0. _d 0
345	heimbach	1.13	xx_vwind_file = 'xx_vwind'
346	heimbach	1.36	xx_vwind_mean_file = 'xx_vwind_mean'
347	heimbach	1.22	xx_vwind_remo_intercept = 0. _d 0
348			xx_vwind_remo_slope = 0. _d 0
349			c
350	heimbach	1.2	xx_obcsnstartdate1 = 0
351			xx_obcsnstartdate2 = 0
352			xx_obcsnperiod = 0. _d 0
353	heimbach	1.13	xx_obcsn_file = 'xx_obcsn'
354	heimbach	1.22	c
355	heimbach	1.2	xx_obcssstartdate1 = 0
356			xx_obcssstartdate2 = 0
357			xx_obcssperiod = 0. _d 0
358	heimbach	1.13	xx_obcss_file = 'xx_obcss'
359	heimbach	1.22	c
360	heimbach	1.2	xx_obcswstartdate1 = 0
361			xx_obcswstartdate2 = 0
362			xx_obcswperiod = 0. _d 0
363	heimbach	1.13	xx_obcsw_file = 'xx_obcsw'
364	heimbach	1.22	c
365	heimbach	1.2	xx_obcsestartdate1 = 0
366			xx_obcsestartdate2 = 0
367			xx_obcseperiod = 0. _d 0
368	heimbach	1.13	xx_obcse_file = 'xx_obcse'
369	heimbach	1.22	c
370	heimbach	1.19	xx_sststartdate1 = 0
371			xx_sststartdate2 = 0
372			xx_sstperiod = 0. _d 0
373	heimbach	1.13	xx_sst_file = 'xx_sst'
374	heimbach	1.35	xx_sst_remo_intercept = 0. _d 0
375			xx_sst_remo_slope = 0. _d 0
376	heimbach	1.22	c
377	heimbach	1.19	xx_sssstartdate1 = 0
378			xx_sssstartdate2 = 0
379			xx_sssperiod = 0. _d 0
380	heimbach	1.13	xx_sss_file = 'xx_sss'
381	heimbach	1.35	xx_sss_remo_intercept = 0. _d 0
382			xx_sss_remo_slope = 0. _d 0
383	heimbach	1.22	c
384	mlosch	1.33	xx_shifwflxstartdate1 = 0
385			xx_shifwflxstartdate2 = 0
386			xx_shifwflxperiod = 0. _d 0
387			xx_shifwflx_file = 'xx_shifwflx'
388			xx_shifwflx_remo_intercept=0. _d 0
389	jmc	1.34	xx_shifwflx_remo_slope = 0. _d 0
390	mlosch	1.33	c
391	heimbach	1.22	xx_diffkr_file = 'xx_diffkr'
392			xx_kapgm_file = 'xx_kapgm'
393	jmc	1.34	xx_kapredi_file = 'xx_kapredi'
394	heimbach	1.22	xx_tr1_file = 'xx_ptr'
395	heimbach	1.24	xx_depth_file = 'xx_depth'
396	heimbach	1.13	xx_efluxy_file = 'xx_efluxy'
397			xx_efluxp_file = 'xx_efluxp'
398			xx_bottomdrag_file = 'xx_bottomdrag'
399			xx_edtaux_file = 'xx_edtaux'
400			xx_edtauy_file = 'xx_edtauy'
401	heimbach	1.14	xx_uvel_file = 'xx_uvel'
402			xx_vvel_file = 'xx_vvel'
403			xx_etan_file = 'xx_etan'
404	jmc	1.34	xx_siarea_file = 'xx_siarea'
405			xx_siheff_file = 'xx_siheff'
406			xx_sihsnow_file = 'xx_sihsnow'
407	dfer	1.28	cHFLUXM_CONTROL
408			xx_hfluxm_file = 'xx_hfluxm'
409			cHFLUXM_CONTROL
410	heimbach	1.2
411	heimbach	1.4	#ifdef ALLOW_TANGENTLINEAR_RUN
412			yadprefix = 'g_'
413			yadmark = 'g_'
414			#else
415			yadprefix = 'ad'
416	heimbach	1.2	yadmark = 'ad'
417			#endif
418			yctrlid = 'MIT_CE_000'
419	heimbach	1.6	yctrlposunpack = '.opt'
420			yctrlpospack = '.opt'
421	heimbach	1.30	ctrlname = 'ecco_ctrl'
422			costname = 'ecco_cost'
423	heimbach	1.2	scalname = ' '
424			maskname = ' '
425			metaname = ' '
426
427	heimbach	1.15	c-- Next, read the cost data file.
428			WRITE(msgBuf,'(A)') 'CTRL_READPARMS: opening data.ctrl'
429			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
430			& SQUEEZE_RIGHT , 1)
431
432			CALL OPEN_COPY_DATA_FILE(
433			I 'data.ctrl', 'CTRL_READPARMS',
434			O iUnit,
435			I myThid )
436
437			READ(unit = iUnit, nml = ctrl_nml)
438			READ(unit = iUnit, nml = ctrl_packnames)
439	heimbach	1.43	#if (defined (ALLOW_GENARR2D_CONTROL) \|\| defined (ALLOW_GENARR3D_CONTROL) \|\| defined (ALLOW_GENTIM2D_CONTROL))
440			READ(unit = iUnit, nml = ctrl_nml_genarr)
441			#endif
442	heimbach	1.15
443	jmc	1.34	WRITE(msgBuf,'(A)')
444	heimbach	1.15	& 'CTRL_READPARMS: finished reading data.ctrl'
445			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
446			& SQUEEZE_RIGHT , 1)
447	heimbach	1.2
448	heimbach	1.15	CLOSE( iUnit )
449	heimbach	1.2
450	gforget	1.49	#ifdef ALLOW_AUTODIFF
451			c allow for switching off correl2d in adjoint
452			useSmoothCorrel2DinFwdMode = ctrlSmoothCorrel2D
453			useSmoothCorrel2DinAdMode =
454			& useSmoothCorrel2DinAdMode .AND. ctrlSmoothCorrel2D
455			CALL WRITE_0D_L( useSmoothCorrel2DinAdMode, INDEX_NONE,
456			& 'useSmoothCorrel2DinAdMode =',
457			& ' /* use ctrlSmoothCorrel2D in adjoint mode */')
458			#endif
459
460	heimbach	1.2	_END_MASTER( myThid )
461
462			_BARRIER
463
464			return
465			end