pkg/ctrl/ctrl_readparms.F

C
C $Header: /u/gcmpack/MITgcm/pkg/ctrl/ctrl_readparms.F,v 1.25 2006/10/27 05:16:54 heimbach Exp $
C $Name:  $

#include "AD_CONFIG.h"
#include "CTRL_CPPOPTIONS.h"


      subroutine ctrl_readparms( mythid )

c     ==================================================================
c     SUBROUTINE ctrl_readparms
c     ==================================================================
c
c     o read ctrl parameters
c       split from ctrl_init
c
c     started: heimbach@mit.edu 12-Jun-2003
c
c     ==================================================================
c     SUBROUTINE ctrl_readparms
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"
#include "GRID.h"
#include "ctrl.h"

#ifdef ALLOW_OBCS_CONTROL
# include "OBCS.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

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_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_tr1_file,
     &  xx_sst_file, xx_sss_file,
     &  xx_sststartdate1, xx_sststartdate2, xx_sstperiod,
     &  xx_sssstartdate1, xx_sssstartdate2, xx_sssperiod,
     &  xx_depth_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,
     &  doInitXX,
     &  doPackDiag,
     &  doZscaleUnpack, doZscalePack,
     &  doMainUnpack, doMainPack,
     &  doAdmtlmBypassAD,
     &  delZexp, forcingPrecond

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

      _BEGIN_MASTER( myThid )

c--     Set default values.
        doInitXX            = .TRUE.
#ifdef ALLOW_ADMTLM
        doAdmTlm            = .TRUE.
#else
        doAdmTlm            = .FALSE.
#endif
        doPackDiag          = .FALSE.
        doZscaleUnpack      = .FALSE.
        doZscalePack        = .FALSE.
        doMainUnpack        = .TRUE.
        doMainPack          = .TRUE.
        doAdmtlmBypassAD    = .FALSE.

        delZexp = 0.
        forcingPrecond = 1. _d 0

        xx_theta_file      = 'xx_theta'
        xx_salt_file       = 'xx_salt'
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_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_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_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_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_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_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'
c
        xx_sssstartdate1    =   0
        xx_sssstartdate2    =   0
        xx_sssperiod        =   0. _d 0
        xx_sss_file         = 'xx_sss'
c
        xx_diffkr_file      = 'xx_diffkr'
        xx_kapgm_file       = 'xx_kapgm'
        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'

#ifdef ALLOW_TANGENTLINEAR_RUN
        yadprefix           = 'g_'
        yadmark             = 'g_'
#else
        yadprefix           = 'ad'
        yadmark             = 'ad'
#endif
        yctrlid             = 'MIT_CE_000'
        yctrlposunpack      = '.opt'
        yctrlpospack        = '.opt'
        ctrlname            = ' '
        costname            = ' '
        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)

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

        CLOSE( iUnit )

      _END_MASTER( myThid )

      _BARRIER

      return
      end

1	edhill	1.5	C
2	gforget	1.26	C $Header: /u/gcmpack/MITgcm/pkg/ctrl/ctrl_readparms.F,v 1.25 2006/10/27 05:16:54 heimbach Exp $
3	heimbach	1.10	C $Name: $
4	heimbach	1.2
5	edhill	1.5	#include "AD_CONFIG.h"
6	heimbach	1.2	#include "CTRL_CPPOPTIONS.h"
7
8
9			subroutine ctrl_readparms( mythid )
10
11			c ==================================================================
12			c SUBROUTINE ctrl_readparms
13			c ==================================================================
14			c
15			c o read ctrl parameters
16			c split from ctrl_init
17			c
18			c started: heimbach@mit.edu 12-Jun-2003
19			c
20			c ==================================================================
21			c SUBROUTINE ctrl_readparms
22			c ==================================================================
23
24			implicit none
25
26			c == global variables ==
27
28			#include "EEPARAMS.h"
29			#include "SIZE.h"
30			#include "PARAMS.h"
31			#include "GRID.h"
32			#include "ctrl.h"
33
34			#ifdef ALLOW_OBCS_CONTROL
35			# include "OBCS.h"
36			#endif
37
38			c == routine arguments ==
39
40			integer mythid
41
42			c == local variables ==
43
44			integer bi,bj
45			integer i,j,k
46			integer ntmp
47			integer ivarindex
48	heimbach	1.15	integer iUnit
49	heimbach	1.2	integer iobcs
50			integer il
51			integer errio
52			integer startrec
53			integer endrec
54			integer difftime(4)
55			_RL diffsecs
56			_RL dummy
57
58			character*(80) ymaskobcs
59			character*(max_len_prec) record
60			character*(max_len_mbuf) msgbuf
61
62			integer nwetc3d
63
64			c == external ==
65
66			integer ilnblnk
67			external ilnblnk
68
69			c == end of interface ==
70
71			c-- Read the namelist input.
72			namelist /ctrl_nml/
73	heimbach	1.25	& xx_theta_file, xx_salt_file,
74			& xx_hflux_file, xx_hflux_remo_intercept, xx_hflux_remo_slope,
75	heimbach	1.2	& xx_hfluxstartdate1, xx_hfluxstartdate2, xx_hfluxperiod,
76	heimbach	1.25	& xx_sflux_file, xx_sflux_remo_intercept, xx_sflux_remo_slope,
77	heimbach	1.2	& xx_sfluxstartdate1, xx_sfluxstartdate2, xx_sfluxperiod,
78	heimbach	1.25	& xx_tauu_file, xx_tauu_remo_intercept, xx_tauu_remo_slope,
79	heimbach	1.2	& xx_tauustartdate1, xx_tauustartdate2, xx_tauuperiod,
80	heimbach	1.25	& xx_tauv_file, xx_tauv_remo_intercept, xx_tauv_remo_slope,
81	heimbach	1.2	& xx_tauvstartdate1, xx_tauvstartdate2, xx_tauvperiod,
82	heimbach	1.25	& xx_atemp_file, xx_atemp_remo_intercept, xx_atemp_remo_slope,
83	heimbach	1.2	& xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
84	heimbach	1.25	& xx_aqh_file, xx_aqh_remo_intercept, xx_aqh_remo_slope,
85	heimbach	1.2	& xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,
86	heimbach	1.25	& xx_precip_file, xx_precip_remo_intercept, xx_precip_remo_slope,
87	heimbach	1.16	& xx_precipstartdate1, xx_precipstartdate2, xx_precipperiod,
88	heimbach	1.25	& xx_swflux_file, xx_swflux_remo_intercept, xx_swflux_remo_slope,
89	heimbach	1.17	& xx_swfluxstartdate1, xx_swfluxstartdate2, xx_swfluxperiod,
90	heimbach	1.25	& xx_swdown_file, xx_swdown_remo_intercept, xx_swdown_remo_slope,
91	heimbach	1.18	& xx_swdownstartdate1, xx_swdownstartdate2, xx_swdownperiod,
92	heimbach	1.25	& xx_lwflux_file, xx_lwflux_remo_intercept, xx_lwflux_remo_slope,
93			& xx_lwfluxstartdate1, xx_lwfluxstartdate2, xx_lwfluxperiod,
94			& xx_lwdown_file, xx_lwdown_remo_intercept, xx_lwdown_remo_slope,
95			& xx_lwdownstartdate1, xx_lwdownstartdate2, xx_lwdownperiod,
96			& xx_evap_file, xx_evap_remo_intercept, xx_evap_remo_slope,
97			& xx_evapstartdate1, xx_evapstartdate2, xx_evapperiod,
98			& xx_snowprecip_file, xx_snowprecip_remo_intercept,
99			& xx_snowprecip_remo_slope, xx_snowprecipperiod,
100			& xx_snowprecipstartdate1, xx_snowprecipstartdate2,
101			& xx_apressure_file, xx_apressure_remo_intercept,
102			& xx_apressure_remo_slope, xx_apressureperiod,
103			& xx_apressurestartdate1, xx_apressurestartdate2,
104			& xx_runoff_file, xx_runoff_remo_intercept, xx_runoff_remo_slope,
105			& xx_runoffstartdate1, xx_runoffstartdate2, xx_runoffperiod,
106			& xx_uwind_file, xx_uwind_remo_intercept, xx_uwind_remo_slope,
107	heimbach	1.2	& xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
108	heimbach	1.25	& xx_vwind_file, xx_vwind_remo_intercept, xx_vwind_remo_slope,
109	heimbach	1.2	& xx_vwindstartdate1, xx_vwindstartdate2, xx_vwindperiod,
110			& xx_obcsn_file,
111			& xx_obcsnstartdate1, xx_obcsnstartdate2, xx_obcsnperiod,
112			& xx_obcss_file,
113			& xx_obcssstartdate1, xx_obcssstartdate2, xx_obcssperiod,
114			& xx_obcsw_file,
115			& xx_obcswstartdate1, xx_obcswstartdate2, xx_obcswperiod,
116			& xx_obcse_file,
117			& xx_obcsestartdate1, xx_obcsestartdate2, xx_obcseperiod,
118	heimbach	1.25	& xx_diffkr_file, xx_kapgm_file, xx_tr1_file,
119			& xx_sst_file, xx_sss_file,
120	heimbach	1.19	& xx_sststartdate1, xx_sststartdate2, xx_sstperiod,
121			& xx_sssstartdate1, xx_sssstartdate2, xx_sssperiod,
122	heimbach	1.24	& xx_depth_file,
123	heimbach	1.25	& xx_efluxy_file, xx_efluxp_file,
124	heimbach	1.7	& xx_bottomdrag_file,
125	heimbach	1.25	& xx_edtaux_file, xx_edtauy_file,
126			& xx_uvel_file, xx_vvel_file, xx_etan_file,
127	heimbach	1.10	& doInitXX,
128	heimbach	1.11	& doPackDiag,
129	heimbach	1.25	& doZscaleUnpack, doZscalePack,
130			& doMainUnpack, doMainPack,
131	heimbach	1.23	& doAdmtlmBypassAD,
132	gforget	1.26	& delZexp, forcingPrecond
133	heimbach	1.2
134			namelist /ctrl_packnames/
135	heimbach	1.6	& yadmark, yctrlid, yctrlposunpack, yctrlpospack,
136	heimbach	1.2	& ctrlname, costname, scalname, maskname, metaname
137
138			_BEGIN_MASTER( myThid )
139
140			c-- Set default values.
141	heimbach	1.8	doInitXX = .TRUE.
142			#ifdef ALLOW_ADMTLM
143			doAdmTlm = .TRUE.
144			#else
145			doAdmTlm = .FALSE.
146			#endif
147	heimbach	1.10	doPackDiag = .FALSE.
148	heimbach	1.11	doZscaleUnpack = .FALSE.
149			doZscalePack = .FALSE.
150	heimbach	1.20	doMainUnpack = .TRUE.
151			doMainPack = .TRUE.
152	heimbach	1.21	doAdmtlmBypassAD = .FALSE.
153	heimbach	1.8
154	heimbach	1.23	delZexp = 0.
155	gforget	1.26	forcingPrecond = 1. _d 0
156	heimbach	1.23
157	heimbach	1.13	xx_theta_file = 'xx_theta'
158			xx_salt_file = 'xx_salt'
159	heimbach	1.22	c
160	heimbach	1.2	xx_hfluxstartdate1 = 0
161			xx_hfluxstartdate2 = 0
162			xx_hfluxperiod = 0. _d 0
163	heimbach	1.13	xx_hflux_file = 'xx_hfl'
164	heimbach	1.22	xx_hflux_remo_intercept = 0. _d 0
165			xx_hflux_remo_slope = 0. _d 0
166			c
167	heimbach	1.2	xx_sfluxstartdate1 = 0
168			xx_sfluxstartdate2 = 0
169			xx_sfluxperiod = 0. _d 0
170	heimbach	1.13	xx_sflux_file = 'xx_sfl'
171	heimbach	1.22	xx_sflux_remo_intercept = 0. _d 0
172			xx_sflux_remo_slope = 0. _d 0
173			c
174	heimbach	1.2	xx_tauustartdate1 = 0
175			xx_tauustartdate2 = 0
176			xx_tauuperiod = 0. _d 0
177	heimbach	1.13	xx_tauu_file = 'xx_tauu'
178	heimbach	1.22	xx_tauu_remo_intercept = 0. _d 0
179			xx_tauu_remo_slope = 0. _d 0
180			c
181	heimbach	1.2	xx_tauvstartdate1 = 0
182			xx_tauvstartdate2 = 0
183			xx_tauvperiod = 0. _d 0
184	heimbach	1.13	xx_tauv_file = 'xx_tauv'
185	heimbach	1.22	xx_tauv_remo_intercept = 0. _d 0
186			xx_tauv_remo_slope = 0. _d 0
187			c
188	heimbach	1.2	xx_atempstartdate1 = 0
189			xx_atempstartdate2 = 0
190			xx_atempperiod = 0. _d 0
191	heimbach	1.13	xx_atemp_file = 'xx_atemp'
192	heimbach	1.22	xx_atemp_remo_intercept = 0. _d 0
193			xx_atemp_remo_slope = 0. _d 0
194			c
195	heimbach	1.2	xx_aqhstartdate1 = 0
196			xx_aqhstartdate2 = 0
197			xx_aqhperiod = 0. _d 0
198	heimbach	1.13	xx_aqh_file = 'xx_aqh'
199	heimbach	1.22	xx_aqh_remo_intercept = 0. _d 0
200			xx_aqh_remo_slope = 0. _d 0
201			c
202	heimbach	1.16	xx_precipstartdate1 = 0
203			xx_precipstartdate2 = 0
204			xx_precipperiod = 0. _d 0
205			xx_precip_file = 'xx_precip'
206	heimbach	1.22	xx_precip_remo_intercept = 0. _d 0
207			xx_precip_remo_slope = 0. _d 0
208			c
209	heimbach	1.17	xx_swfluxstartdate1 = 0
210			xx_swfluxstartdate2 = 0
211			xx_swfluxperiod = 0. _d 0
212			xx_swflux_file = 'xx_swflux'
213	heimbach	1.22	xx_swflux_remo_intercept = 0. _d 0
214			xx_swflux_remo_slope = 0. _d 0
215			c
216	heimbach	1.18	xx_swdownstartdate1 = 0
217			xx_swdownstartdate2 = 0
218			xx_swdownperiod = 0. _d 0
219			xx_swdown_file = 'xx_swdown'
220	heimbach	1.22	xx_swdown_remo_intercept = 0. _d 0
221			xx_swdown_remo_slope = 0. _d 0
222			c
223	heimbach	1.25	xx_lwfluxstartdate1 = 0
224			xx_lwfluxstartdate2 = 0
225			xx_lwfluxperiod = 0. _d 0
226			xx_lwflux_file = 'xx_lwflux'
227			xx_lwflux_remo_intercept = 0. _d 0
228			xx_lwflux_remo_slope = 0. _d 0
229			c
230			xx_lwdownstartdate1 = 0
231			xx_lwdownstartdate2 = 0
232			xx_lwdownperiod = 0. _d 0
233			xx_lwdown_file = 'xx_lwdown'
234			xx_lwdown_remo_intercept = 0. _d 0
235			xx_lwdown_remo_slope = 0. _d 0
236			c
237			xx_evapstartdate1 = 0
238			xx_evapstartdate2 = 0
239			xx_evapperiod = 0. _d 0
240			xx_evap_file = 'xx_evap'
241			xx_evap_remo_intercept = 0. _d 0
242			xx_evap_remo_slope = 0. _d 0
243			c
244			xx_snowprecipstartdate1 = 0
245			xx_snowprecipstartdate2 = 0
246			xx_snowprecipperiod = 0. _d 0
247			xx_snowprecip_file = 'xx_snowprecip'
248			xx_snowprecip_remo_intercept = 0. _d 0
249			xx_snowprecip_remo_slope = 0. _d 0
250			c
251			xx_apressurestartdate1 = 0
252			xx_apressurestartdate2 = 0
253			xx_apressureperiod = 0. _d 0
254			xx_apressure_file = 'xx_apressure'
255			xx_apressure_remo_intercept = 0. _d 0
256			xx_apressure_remo_slope = 0. _d 0
257			c
258			xx_runoffstartdate1 = 0
259			xx_runoffstartdate2 = 0
260			xx_runoffperiod = 0. _d 0
261			xx_runoff_file = 'xx_runoff'
262			xx_runoff_remo_intercept = 0. _d 0
263			xx_runoff_remo_slope = 0. _d 0
264			c
265	heimbach	1.2	xx_uwindstartdate1 = 0
266			xx_uwindstartdate2 = 0
267			xx_uwindperiod = 0. _d 0
268	heimbach	1.13	xx_uwind_file = 'xx_uwind'
269	heimbach	1.22	xx_uwind_remo_intercept = 0. _d 0
270			xx_uwind_remo_slope = 0. _d 0
271			c
272	heimbach	1.2	xx_vwindstartdate1 = 0
273			xx_vwindstartdate2 = 0
274			xx_vwindperiod = 0. _d 0
275	heimbach	1.13	xx_vwind_file = 'xx_vwind'
276	heimbach	1.22	xx_vwind_remo_intercept = 0. _d 0
277			xx_vwind_remo_slope = 0. _d 0
278			c
279	heimbach	1.2	xx_obcsnstartdate1 = 0
280			xx_obcsnstartdate2 = 0
281			xx_obcsnperiod = 0. _d 0
282	heimbach	1.13	xx_obcsn_file = 'xx_obcsn'
283	heimbach	1.22	c
284	heimbach	1.2	xx_obcssstartdate1 = 0
285			xx_obcssstartdate2 = 0
286			xx_obcssperiod = 0. _d 0
287	heimbach	1.13	xx_obcss_file = 'xx_obcss'
288	heimbach	1.22	c
289	heimbach	1.2	xx_obcswstartdate1 = 0
290			xx_obcswstartdate2 = 0
291			xx_obcswperiod = 0. _d 0
292	heimbach	1.13	xx_obcsw_file = 'xx_obcsw'
293	heimbach	1.22	c
294	heimbach	1.2	xx_obcsestartdate1 = 0
295			xx_obcsestartdate2 = 0
296			xx_obcseperiod = 0. _d 0
297	heimbach	1.13	xx_obcse_file = 'xx_obcse'
298	heimbach	1.22	c
299	heimbach	1.19	xx_sststartdate1 = 0
300			xx_sststartdate2 = 0
301			xx_sstperiod = 0. _d 0
302	heimbach	1.13	xx_sst_file = 'xx_sst'
303	heimbach	1.22	c
304	heimbach	1.19	xx_sssstartdate1 = 0
305			xx_sssstartdate2 = 0
306			xx_sssperiod = 0. _d 0
307	heimbach	1.13	xx_sss_file = 'xx_sss'
308	heimbach	1.22	c
309			xx_diffkr_file = 'xx_diffkr'
310			xx_kapgm_file = 'xx_kapgm'
311			xx_tr1_file = 'xx_ptr'
312	heimbach	1.24	xx_depth_file = 'xx_depth'
313	heimbach	1.13	xx_efluxy_file = 'xx_efluxy'
314			xx_efluxp_file = 'xx_efluxp'
315			xx_bottomdrag_file = 'xx_bottomdrag'
316			xx_edtaux_file = 'xx_edtaux'
317			xx_edtauy_file = 'xx_edtauy'
318	heimbach	1.14	xx_uvel_file = 'xx_uvel'
319			xx_vvel_file = 'xx_vvel'
320			xx_etan_file = 'xx_etan'
321	heimbach	1.2
322	heimbach	1.4	#ifdef ALLOW_TANGENTLINEAR_RUN
323			yadprefix = 'g_'
324			yadmark = 'g_'
325			#else
326			yadprefix = 'ad'
327	heimbach	1.2	yadmark = 'ad'
328			#endif
329			yctrlid = 'MIT_CE_000'
330	heimbach	1.6	yctrlposunpack = '.opt'
331			yctrlpospack = '.opt'
332	heimbach	1.2	ctrlname = ' '
333			costname = ' '
334			scalname = ' '
335			maskname = ' '
336			metaname = ' '
337
338	heimbach	1.15	c-- Next, read the cost data file.
339			WRITE(msgBuf,'(A)') 'CTRL_READPARMS: opening data.ctrl'
340			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
341			& SQUEEZE_RIGHT , 1)
342
343			CALL OPEN_COPY_DATA_FILE(
344			I 'data.ctrl', 'CTRL_READPARMS',
345			O iUnit,
346			I myThid )
347
348			READ(unit = iUnit, nml = ctrl_nml)
349			READ(unit = iUnit, nml = ctrl_packnames)
350
351			WRITE(msgBuf,'(A)')
352			& 'CTRL_READPARMS: finished reading data.ctrl'
353			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
354			& SQUEEZE_RIGHT , 1)
355	heimbach	1.2
356	heimbach	1.15	CLOSE( iUnit )
357	heimbach	1.2
358			_END_MASTER( myThid )
359
360			_BARRIER
361
362			return
363			end
364