pkg/ctrl/ctrl_readparms.F

C
C $Header: /u/gcmpack/MITgcm/pkg/ctrl/ctrl_readparms.F,v 1.9 2004/12/04 17:12:00 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 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_hfluxstartdate1, xx_hfluxstartdate2, xx_hfluxperiod,
     &  xx_sflux_file,
     &  xx_sfluxstartdate1, xx_sfluxstartdate2, xx_sfluxperiod,
     &  xx_tauu_file,
     &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,
     &  xx_tauv_file,
     &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod,
     &  xx_atemp_file,
     &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
     &  xx_aqh_file,
     &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,
     &  xx_uwind_file,
     &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
     &  xx_vwind_file,
     &  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_hfacc_file,
     &  xx_efluxy_file,
     &  xx_efluxp_file,
     &  xx_bottomdrag_file,
     &  doInitXX,
     &  doPackDiag

      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.

        xx_theta_file      = ' '
        xx_salt_file       = ' '
        xx_hfluxstartdate1 =   0
        xx_hfluxstartdate2 =   0
        xx_hfluxperiod     =   0. _d 0
        xx_hflux_file      = ' '
        xx_sfluxstartdate1 =   0
        xx_sfluxstartdate2 =   0
        xx_sfluxperiod     =   0. _d 0
        xx_sflux_file      = ' '
        xx_tauustartdate1  =   0
        xx_tauustartdate2  =   0
        xx_tauuperiod      =   0. _d 0
        xx_tauu_file       = ' '
        xx_tauvstartdate1  =   0
        xx_tauvstartdate2  =   0
        xx_tauvperiod      =   0. _d 0
        xx_tauv_file       = ' '
        xx_atempstartdate1  =   0
        xx_atempstartdate2  =   0
        xx_atempperiod      =   0. _d 0
        xx_atemp_file       = ' '
        xx_aqhstartdate1    =   0
        xx_aqhstartdate2    =   0
        xx_aqhperiod        =   0. _d 0
        xx_aqh_file         = ' '
        xx_uwindstartdate1  =   0
        xx_uwindstartdate2  =   0
        xx_uwindperiod      =   0. _d 0
        xx_uwind_file       = ' '
        xx_vwindstartdate1  =   0
        xx_vwindstartdate2  =   0
        xx_vwindperiod      =   0. _d 0
        xx_vwind_file       = ' '
        xx_obcsnstartdate1  =   0
        xx_obcsnstartdate2  =   0
        xx_obcsnperiod      =   0. _d 0
        xx_obcsn_file       = ' '
        xx_obcssstartdate1  =   0
        xx_obcssstartdate2  =   0
        xx_obcssperiod      =   0. _d 0
        xx_obcss_file       = ' '
        xx_obcswstartdate1  =   0
        xx_obcswstartdate2  =   0
        xx_obcswperiod      =   0. _d 0
        xx_obcsw_file       = ' '
        xx_obcsestartdate1  =   0
        xx_obcsestartdate2  =   0
        xx_obcseperiod      =   0. _d 0
        xx_obcse_file       = ' '
        xx_diffkr_file      = ' '
        xx_kapgm_file       = ' '
        xx_tr1_file         = ' '
        xx_sst_file         = ' '
        xx_sss_file         = ' '
        xx_hfacc_file       = ' '
        xx_efluxy_file      = ' '
        xx_efluxp_file      = ' '
        xx_bottomdrag_file  = ' '

#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--     Check versions.

        open(unit=scrunit1,status='scratch')

c--     Next, read the ecco data file.
        open(unit = modeldataunit,file = 'data.ctrl',
     &       status = 'old', iostat = errio)
        if ( errio .lt. 0 ) then
          stop ' stopped in ctrl_readparms'
        endif

        do while ( .true. )
          read(modeldataunit, fmt='(a)', end=1001) record
          il = max(ilnblnk(record),1)
          if ( record(1:1) .ne. commentcharacter )
     &        write(unit=scrunit1, fmt='(a)') record(:il)
        enddo
 1001   continue
        close( modeldataunit )

        rewind( scrunit1 )
        read(unit = scrunit1, nml = ctrl_nml)
        read(unit = scrunit1, nml = ctrl_packnames)
        close( scrunit1 )

      _END_MASTER( myThid )

      _BARRIER

      return
      end

1	edhill	1.5	C
2	heimbach	1.10	C $Header: /u/gcmpack/MITgcm/pkg/ctrl/ctrl_readparms.F,v 1.9 2004/12/04 17:12:00 heimbach Exp $
3			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
49			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			& xx_theta_file,
74			& xx_salt_file,
75			& xx_hflux_file,
76			& xx_hfluxstartdate1, xx_hfluxstartdate2, xx_hfluxperiod,
77			& xx_sflux_file,
78			& xx_sfluxstartdate1, xx_sfluxstartdate2, xx_sfluxperiod,
79			& xx_tauu_file,
80			& xx_tauustartdate1, xx_tauustartdate2, xx_tauuperiod,
81			& xx_tauv_file,
82			& xx_tauvstartdate1, xx_tauvstartdate2, xx_tauvperiod,
83			& xx_atemp_file,
84			& xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,
85			& xx_aqh_file,
86			& xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,
87			& xx_uwind_file,
88			& xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,
89			& xx_vwind_file,
90			& xx_vwindstartdate1, xx_vwindstartdate2, xx_vwindperiod,
91			& xx_obcsn_file,
92			& xx_obcsnstartdate1, xx_obcsnstartdate2, xx_obcsnperiod,
93			& xx_obcss_file,
94			& xx_obcssstartdate1, xx_obcssstartdate2, xx_obcssperiod,
95			& xx_obcsw_file,
96			& xx_obcswstartdate1, xx_obcswstartdate2, xx_obcswperiod,
97			& xx_obcse_file,
98			& xx_obcsestartdate1, xx_obcsestartdate2, xx_obcseperiod,
99			& xx_diffkr_file,
100			& xx_kapgm_file,
101			& xx_tr1_file,
102			& xx_sst_file,
103			& xx_sss_file,
104			& xx_hfacc_file,
105			& xx_efluxy_file,
106			& xx_efluxp_file,
107	heimbach	1.7	& xx_bottomdrag_file,
108	heimbach	1.10	& doInitXX,
109			& doPackDiag
110	heimbach	1.2
111			namelist /ctrl_packnames/
112	heimbach	1.6	& yadmark, yctrlid, yctrlposunpack, yctrlpospack,
113	heimbach	1.2	& ctrlname, costname, scalname, maskname, metaname
114
115			_BEGIN_MASTER( myThid )
116
117			c-- Set default values.
118	heimbach	1.8	doInitXX = .TRUE.
119			#ifdef ALLOW_ADMTLM
120			doAdmTlm = .TRUE.
121			#else
122			doAdmTlm = .FALSE.
123			#endif
124	heimbach	1.10	doPackDiag = .FALSE.
125	heimbach	1.8
126	heimbach	1.2	xx_theta_file = ' '
127			xx_salt_file = ' '
128			xx_hfluxstartdate1 = 0
129			xx_hfluxstartdate2 = 0
130			xx_hfluxperiod = 0. _d 0
131			xx_hflux_file = ' '
132			xx_sfluxstartdate1 = 0
133			xx_sfluxstartdate2 = 0
134			xx_sfluxperiod = 0. _d 0
135			xx_sflux_file = ' '
136			xx_tauustartdate1 = 0
137			xx_tauustartdate2 = 0
138			xx_tauuperiod = 0. _d 0
139			xx_tauu_file = ' '
140			xx_tauvstartdate1 = 0
141			xx_tauvstartdate2 = 0
142			xx_tauvperiod = 0. _d 0
143			xx_tauv_file = ' '
144			xx_atempstartdate1 = 0
145			xx_atempstartdate2 = 0
146			xx_atempperiod = 0. _d 0
147			xx_atemp_file = ' '
148			xx_aqhstartdate1 = 0
149			xx_aqhstartdate2 = 0
150			xx_aqhperiod = 0. _d 0
151			xx_aqh_file = ' '
152			xx_uwindstartdate1 = 0
153			xx_uwindstartdate2 = 0
154			xx_uwindperiod = 0. _d 0
155			xx_uwind_file = ' '
156			xx_vwindstartdate1 = 0
157			xx_vwindstartdate2 = 0
158			xx_vwindperiod = 0. _d 0
159			xx_vwind_file = ' '
160			xx_obcsnstartdate1 = 0
161			xx_obcsnstartdate2 = 0
162			xx_obcsnperiod = 0. _d 0
163			xx_obcsn_file = ' '
164			xx_obcssstartdate1 = 0
165			xx_obcssstartdate2 = 0
166			xx_obcssperiod = 0. _d 0
167			xx_obcss_file = ' '
168			xx_obcswstartdate1 = 0
169			xx_obcswstartdate2 = 0
170			xx_obcswperiod = 0. _d 0
171			xx_obcsw_file = ' '
172			xx_obcsestartdate1 = 0
173			xx_obcsestartdate2 = 0
174			xx_obcseperiod = 0. _d 0
175			xx_obcse_file = ' '
176			xx_diffkr_file = ' '
177			xx_kapgm_file = ' '
178			xx_tr1_file = ' '
179			xx_sst_file = ' '
180			xx_sss_file = ' '
181			xx_hfacc_file = ' '
182			xx_efluxy_file = ' '
183			xx_efluxp_file = ' '
184			xx_bottomdrag_file = ' '
185
186	heimbach	1.4	#ifdef ALLOW_TANGENTLINEAR_RUN
187			yadprefix = 'g_'
188			yadmark = 'g_'
189			#else
190			yadprefix = 'ad'
191	heimbach	1.2	yadmark = 'ad'
192			#endif
193			yctrlid = 'MIT_CE_000'
194	heimbach	1.6	yctrlposunpack = '.opt'
195			yctrlpospack = '.opt'
196	heimbach	1.2	ctrlname = ' '
197			costname = ' '
198			scalname = ' '
199			maskname = ' '
200			metaname = ' '
201
202			c-- Check versions.
203
204			open(unit=scrunit1,status='scratch')
205
206			c-- Next, read the ecco data file.
207			open(unit = modeldataunit,file = 'data.ctrl',
208			& status = 'old', iostat = errio)
209			if ( errio .lt. 0 ) then
210			stop ' stopped in ctrl_readparms'
211			endif
212
213			do while ( .true. )
214			read(modeldataunit, fmt='(a)', end=1001) record
215			il = max(ilnblnk(record),1)
216			if ( record(1:1) .ne. commentcharacter )
217			& write(unit=scrunit1, fmt='(a)') record(:il)
218			enddo
219			1001 continue
220			close( modeldataunit )
221
222			rewind( scrunit1 )
223			read(unit = scrunit1, nml = ctrl_nml)
224			read(unit = scrunit1, nml = ctrl_packnames)
225			close( scrunit1 )
226
227			_END_MASTER( myThid )
228
229			_BARRIER
230
231			return
232			end
233