pkg/exf/exf_readparms.F

c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_readparms.F,v 1.1.6.6 2003/03/04 00:17:40 heimbach Exp $

#include "EXF_CPPOPTIONS.h"
#ifdef ALLOW_OBCS
# include "OBCS_OPTIONS.h"
#endif

      subroutine exf_readparms( mythid )

c     ==================================================================
c     SUBROUTINE exf_readparms
c     ==================================================================
c
c     o This routine initialises the package that calculates external
c       forcing fields for a given timestep of the MITgcmUV. Parameters
c       for this package are set in "data.externalforcing". Some additional
c       precompiler switches have to be specified in "EXF_CPPOPTIONS.h".
c
c     started: Christian Eckert eckert@mit.edu  30-Jun-1999
c
c     changed: Christian Eckert eckert@mit.edu  11-Jan-2000
c              - Restructured the code in order to create a package
c                for the MITgcmUV.
c              Christian Eckert eckert@mit.edu  12-Feb-2000
c              - Changed Routine names (package prefix: exf_)
c     changed: Patrick Heimbach, heimbach@mit.edu  04-May-2000
c              - changed the handling of precip and sflux with respect
c                to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP
c     changed: Ralf.Giering@FastOpt.de 25-Mai-20000
c              - moved relaxation and climatology to extra routines
c              Patrick Heimbach, heimbach@mit.edu  04-May-2000
c              - added obcs parameters
c     changed: Virginie Thierry, vthierry@ucsd.edu 04-June-2001
c              - added new obcs parameters (for each boundaries) 
c     included runoff D. Stammer, Nov. 25, 2001
c     included pressure forcing. heimbach@mit.edu 05-Nov-2002
c     added "repeatPeriod" for cycling of forcing datasets 19-Dec-2002
c     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
c
c     ==================================================================
c     SUBROUTINE exf_readparms
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"
#include "cal.h"
#include "exf.h"
#include "exf_param.h"
#include "exf_constants.h"

c     == routine arguments ==

      integer mythid

c     == local variables ==

      integer i

c     == external ==

      integer  ilnblnk
      external ilnblnk

c     == end of interface ==

c     Surface flux data.
      namelist /exf_nml/
     &    repeatPeriod,
     &    hfluxstartdate1,    hfluxstartdate2,   hfluxperiod,
     &    atempstartdate1,    atempstartdate2,   atempperiod,
     &      aqhstartdate1,      aqhstartdate2,     aqhperiod,
     &    sfluxstartdate1,    sfluxstartdate2,   sfluxperiod,
     &     evapstartdate1,     evapstartdate2,    evapperiod,
     &   precipstartdate1,   precipstartdate2,  precipperiod,
     &   runoffstartdate1,   runoffstartdate2,  runoffperiod,
     &  ustressstartdate1,  ustressstartdate2, ustressperiod,
     &  vstressstartdate1,  vstressstartdate2, vstressperiod,
     &    uwindstartdate1,    uwindstartdate2,   uwindperiod,
     &    vwindstartdate1,    vwindstartdate2,   vwindperiod,
     &   swfluxstartdate1,   swfluxstartdate2,  swfluxperiod,
     &   lwfluxstartdate1,   lwfluxstartdate2,  lwfluxperiod,
     &   swdownstartdate1,   swdownstartdate2,  swdownperiod,
     &   lwdownstartdate1,   lwdownstartdate2,  lwdownperiod,
     &    obcsNstartdate1,    obcsNstartdate2,   obcsNperiod,
     &    obcsSstartdate1,    obcsSstartdate2,   obcsSperiod,
     &    obcsEstartdate1,    obcsEstartdate2,   obcsEperiod,
     &    obcsWstartdate1,    obcsWstartdate2,   obcsWperiod,
     &apressurestartdate1,apressurestartdate2,apressureperiod,
     &          hfluxfile,          atempfile,       aqhfile,
     &          sfluxfile,         precipfile,    runofffile,
     &        ustressfile,        vstressfile,      evapfile,
     &          uwindfile,          vwindfile,
     &         swfluxfile,         lwfluxfile, apressurefile,
     &         swdownfile,         lwdownfile,
     &          exf_iprec,         exf_yftype,
     &     exf_inscal_hfl,     exf_inscal_ust, exf_inscal_vst,
     &     exf_inscal_swf,     exf_inscal_sst, exf_inscal_sss,
     &     exf_inscal_sfl,  exf_inscal_runoff,exf_inscal_apressure,
     &    exf_outscal_hfl,    exf_outscal_ust,exf_outscal_vst,
     &    exf_outscal_swf,    exf_outscal_sst,exf_outscal_sss,
     &    exf_outscal_sfl,  exf_outscal_apressure,
     &         hfluxconst,
     &         atempconst,
     &         aqhconst,
     &         sfluxconst,
     &         evapconst,
     &         precipconst,
     &         runoffconst,
     &         ustressconst,
     &         vstressconst,
     &         uwindconst,
     &         vwindconst,
     &         swfluxconst,
     &         lwfluxconst,
     &         swdownconst,
     &         lwdownconst,
     &         apressureconst,
     &         EXFwindOnBgrid

      _BEGIN_MASTER(mythid)

c     Set default values.

c     Calendar data.
      hfluxstartdate1    = 0
      hfluxstartdate2    = 0
      hfluxperiod        = 0.0 _d 0
      hfluxconst         = 0.0 _d 0

      atempstartdate1    = 0
      atempstartdate2    = 0
      atempperiod        = 0.0 _d 0
      atempconst         = 0.0 _d 0

      aqhstartdate1      = 0
      aqhstartdate2      = 0
      aqhperiod          = 0.0 _d 0
      aqhconst           = 0.0 _d 0

      sfluxstartdate1    = 0
      sfluxstartdate2    = 0
      sfluxperiod        = 0.0 _d 0
      sfluxconst         = 0.0 _d 0

      evapstartdate1   = 0
      evapstartdate2   = 0
      evapperiod       = 0.0 _d 0
      evapconst        = 0.0 _d 0

      precipstartdate1   = 0
      precipstartdate2   = 0
      precipperiod       = 0.0 _d 0
      precipconst        = 0.0 _d 0

      runoffstartdate1   = 0
      runoffstartdate2   = 0
      runoffperiod       = 0.0 _d 0
      runoffconst        = 0.0 _d 0

      ustressstartdate1  = 0
      ustressstartdate2  = 0
      ustressperiod      = 0.0 _d 0
      ustressconst       = 0.0 _d 0

      vstressstartdate1  = 0
      vstressstartdate2  = 0
      vstressperiod      = 0.0 _d 0
      vstressconst       = 0.0 _d 0

      uwindstartdate1    = 0
      uwindstartdate2    = 0
      uwindperiod        = 0.0 _d 0
      uwindconst         = 0.0 _d 0

      vwindstartdate1    = 0
      vwindstartdate2    = 0
      vwindperiod        = 0.0 _d 0
      vwindconst         = 0.0 _d 0

      swfluxstartdate1   = 0
      swfluxstartdate2   = 0
      swfluxperiod       = 0.0 _d 0
      swfluxconst        = 0.0 _d 0

      lwfluxstartdate1   = 0
      lwfluxstartdate2   = 0
      lwfluxperiod       = 0.0 _d 0
      lwfluxconst        = 0.0 _d 0

      swdownstartdate1   = 0
      swdownstartdate2   = 0
      swdownperiod       = 0.0 _d 0
      swdownconst        = 0.0 _d 0

      lwdownstartdate1   = 0
      lwdownstartdate2   = 0
      lwdownperiod       = 0.0 _d 0
      lwdownconst        = 0.0 _d 0

      obcsNstartdate1    = 0
      obcsNstartdate2    = 0
      obcsNperiod        = 0.0 _d 0

      obcsSstartdate1    = 0
      obcsSstartdate2    = 0
      obcsSperiod        = 0.0 _d 0

      obcsEstartdate1    = 0
      obcsEstartdate2    = 0
      obcsEperiod        = 0.0 _d 0

      obcsWstartdate1    = 0
      obcsWstartdate2    = 0
      obcsWperiod        = 0.0 _d 0

      apressurestartdate1    = 0
      apressurestartdate2    = 0
      apressureperiod        = 0.0 _d 0
      apressureconst         = 0.0 _d 0

      repeatPeriod           = 0.0 _d 0

c     Data files.
      hfluxfile          = ' '
      atempfile          = ' '
      aqhfile            = ' '
      evapfile           = ' '
      precipfile         = ' '
      sfluxfile          = ' '
      runofffile         = ' '
      ustressfile        = ' '
      vstressfile        = ' '
      uwindfile          = ' '
      vwindfile          = ' '
      swfluxfile         = ' '
      lwfluxfile         = ' '
      swdownfile         = ' '
      lwdownfile         = ' '
      apressurefile      = ' '

c     Initialise the date arrays.
      do i = 1,4
         hfluxstartdate(i)    = 0
         atempstartdate(i)    = 0
         aqhstartdate(i)      = 0
         evapstartdate(i)     = 0
         precipstartdate(i)   = 0
         sfluxstartdate(i)    = 0
         runoffstartdate(i)   = 0
         ustressstartdate(i)  = 0
         vstressstartdate(i)  = 0
         uwindstartdate(i)    = 0
         vwindstartdate(i)    = 0
         swfluxstartdate(i)   = 0
         lwfluxstartdate(i)   = 0
         swdownstartdate(i)   = 0
         lwdownstartdate(i)   = 0
         obcsNstartdate(i)    = 0
         obcsSstartdate(i)    = 0
         obcsEstartdate(i)    = 0
         obcsWstartdate(i)    = 0
         apressurestartdate(i)= 0
      enddo

c     Initialise file type and field precision
      exf_iprec       = 32
      exf_yftype      = 'RL'

c     Input scaling factors.
      exf_inscal_hfl       =  1. _d 0
      exf_inscal_ust       =  1. _d 0
      exf_inscal_vst       =  1. _d 0
      exf_inscal_swf       =  1. _d 0
      exf_inscal_sst       =  1. _d 0
      exf_inscal_sss       =  1. _d 0
      exf_inscal_apressure =  1. _d 0
      exf_inscal_sfl       =  1. _d 0
cds   convert runoff from m/yr to m/s and change sign.
      exf_inscal_runoff    = -1.0/(86400.*365.0)

c     Output scaling factors.
      exf_outscal_hfl      =  1. _d 0
      exf_outscal_ust      =  1. _d 0
      exf_outscal_vst      =  1. _d 0
      exf_outscal_swf      =  1. _d 0
      exf_outscal_sst      =  1. _d 0
      exf_outscal_sss      =  1. _d 0
      exf_outscal_apressure=  1. _d 0
      exf_outscal_sfl      =  1. _d 0
c
      EXFwindOnBgrid = .FALSE.

c     Check for the availability of the right calendar version.
      if ( calendarversion .ne. usescalendarversion ) then
         print*,' exf_readparms: You are not using the appropriate'
         print*,'           version of the calendar package.'
         print*
         print*,' You are using Calendar version: ', calendarversion
         print*,' Please use    Calendar version: ', usescalendarversion
         stop ' stopped in exf_readparms.'
      endif

c     Next, read the forcing data file.
      call nml_filter( 'data.exf', scrunit1, myThid )
      if (scrunit1 .eq. 0) then
         stop 'exf_readparms: reading namelist failed'
      end if
      read(  scrunit1, nml = exf_nml )
      close( scrunit1 )

c     Complete the start date specifications for the forcing
c     fields to get a complete calendar date array.

c     check for consistency

      if (.NOT. 
     &     (exf_iprec .EQ. 32 .OR. exf_iprec .EQ. 64)
     &     ) then
         stop 'stop in exf_readparms: value of exf_iprec not allowed'
      else if (.NOT. 
     &        (exf_yftype .EQ. 'RS' .OR. 
     &        exf_yftype .EQ. 'RL')
     &        ) then
         stop 'stop in exf_readparms: value of exf_yftype not allowed'
      end if

#ifdef ALLOW_ATM_WIND
      call cal_FullDate(    uwindstartdate1,    uwindstartdate2,
     &                      uwindstartdate ,             mythid )
      call cal_FullDate(    vwindstartdate1,    vwindstartdate2,
     &                      vwindstartdate ,             mythid )
#else
      call cal_FullDate(  ustressstartdate1,  ustressstartdate2,
     &                    ustressstartdate ,             mythid )
      call cal_FullDate(  vstressstartdate1,  vstressstartdate2,
     &                    vstressstartdate ,             mythid )
#endif

#ifdef ALLOW_ATM_TEMP
      call cal_FullDate(    atempstartdate1,    atempstartdate2,
     &                      atempstartdate ,             mythid )
      call cal_FullDate(      aqhstartdate1,      aqhstartdate2,
     &                        aqhstartdate ,             mythid )
      call cal_FullDate(   lwfluxstartdate1,   lwfluxstartdate2,
     &                     lwfluxstartdate ,             mythid )
      call cal_FullDate(   precipstartdate1,   precipstartdate2,
     &                     precipstartdate ,             mythid )
#else
      call cal_FullDate(    hfluxstartdate1,    hfluxstartdate2,
     &                      hfluxstartdate ,             mythid )
      call cal_FullDate(    sfluxstartdate1,    sfluxstartdate2,
     &                      sfluxstartdate ,             mythid )
#endif

#if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
      call cal_FullDate(   swfluxstartdate1,   swfluxstartdate2,
     &                     swfluxstartdate ,             mythid )
#endif

#ifdef EXF_READ_EVAP
      call cal_FullDate(     evapstartdate1,     evapstartdate2,
     &                       evapstartdate ,             mythid )
#endif

#ifdef ALLOW_RUNOFF
      call cal_FullDate(   runoffstartdate1,   runoffstartdate2,
     &                     runoffstartdate ,             mythid )
#endif

#ifdef ALLOW_DOWNWARD_RADIATION
      call cal_FullDate(   swdownstartdate1,   swdownstartdate2,
     &                     swdownstartdate ,             mythid )
      call cal_FullDate(   lwdownstartdate1,   lwdownstartdate2,
     &                     lwdownstartdate ,             mythid )
#endif

#ifdef ALLOW_OBCS
#ifdef ALLOW_OBCS_NORTH
      call cal_FullDate(    obcsNstartdate1,    obcsNstartdate2,
     &                      obcsNstartdate ,             mythid )
#endif
#ifdef ALLOW_OBCS_SOUTH
      call cal_FullDate(    obcsSstartdate1,    obcsSstartdate2,
     &                      obcsSstartdate ,             mythid )
#endif
#ifdef ALLOW_OBCS_EAST
      call cal_FullDate(    obcsEstartdate1,    obcsEstartdate2,
     &                      obcsEstartdate ,             mythid )
#endif
#ifdef ALLOW_OBCS_WEST
      call cal_FullDate(    obcsWstartdate1,    obcsWstartdate2,
     &                      obcsWstartdate ,             mythid )
#endif
#endif /* ALLOW_OBCS */

#ifdef ATMOSPHERIC_LOADING
      call cal_FullDate(apressurestartdate1,apressurestartdate2,
     &                  apressurestartdate ,             mythid )
#endif

      _END_MASTER( mythid )

      _BARRIER

c--   Summarize the External forcing's setup.
      call exf_summary( mythid )

c--   set climatology parameters
      call exf_clim_readparms( mythid )

c--   summarize climatologic forcing configuration
      call exf_clim_summary( mythid )

      end
1	heimbach	1.7	c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_readparms.F,v 1.1.6.6 2003/03/04 00:17:40 heimbach Exp $
2	heimbach	1.1
3			#include "EXF_CPPOPTIONS.h"
4	heimbach	1.3	#ifdef ALLOW_OBCS
5			# include "OBCS_OPTIONS.h"
6			#endif
7	heimbach	1.1
8	heimbach	1.3	subroutine exf_readparms( mythid )
9	heimbach	1.1
10			c ==================================================================
11			c SUBROUTINE exf_readparms
12			c ==================================================================
13			c
14			c o This routine initialises the package that calculates external
15			c forcing fields for a given timestep of the MITgcmUV. Parameters
16			c for this package are set in "data.externalforcing". Some additional
17	heimbach	1.2	c precompiler switches have to be specified in "EXF_CPPOPTIONS.h".
18	heimbach	1.1	c
19			c started: Christian Eckert eckert@mit.edu 30-Jun-1999
20			c
21			c changed: Christian Eckert eckert@mit.edu 11-Jan-2000
22			c - Restructured the code in order to create a package
23			c for the MITgcmUV.
24			c Christian Eckert eckert@mit.edu 12-Feb-2000
25			c - Changed Routine names (package prefix: exf_)
26	heimbach	1.2	c changed: Patrick Heimbach, heimbach@mit.edu 04-May-2000
27	heimbach	1.1	c - changed the handling of precip and sflux with respect
28			c to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP
29			c changed: Ralf.Giering@FastOpt.de 25-Mai-20000
30	heimbach	1.2	c - moved relaxation and climatology to extra routines
31	heimbach	1.1	c Patrick Heimbach, heimbach@mit.edu 04-May-2000
32			c - added obcs parameters
33	heimbach	1.2	c changed: Virginie Thierry, vthierry@ucsd.edu 04-June-2001
34			c - added new obcs parameters (for each boundaries)
35	heimbach	1.3	c included runoff D. Stammer, Nov. 25, 2001
36			c included pressure forcing. heimbach@mit.edu 05-Nov-2002
37	dimitri	1.5	c added "repeatPeriod" for cycling of forcing datasets 19-Dec-2002
38	dimitri	1.6	c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
39	heimbach	1.1	c
40			c ==================================================================
41			c SUBROUTINE exf_readparms
42			c ==================================================================
43
44			implicit none
45
46			c == global variables ==
47
48			#include "EEPARAMS.h"
49			#include "SIZE.h"
50			#include "cal.h"
51			#include "exf.h"
52			#include "exf_param.h"
53			#include "exf_constants.h"
54
55			c == routine arguments ==
56
57			integer mythid
58
59			c == local variables ==
60
61			integer i
62
63			c == external ==
64
65			integer ilnblnk
66			external ilnblnk
67
68			c == end of interface ==
69
70			c Surface flux data.
71			namelist /exf_nml/
72	cheisey	1.4	& repeatPeriod,
73	heimbach	1.1	& hfluxstartdate1, hfluxstartdate2, hfluxperiod,
74			& atempstartdate1, atempstartdate2, atempperiod,
75			& aqhstartdate1, aqhstartdate2, aqhperiod,
76			& sfluxstartdate1, sfluxstartdate2, sfluxperiod,
77	dimitri	1.5	& evapstartdate1, evapstartdate2, evapperiod,
78	heimbach	1.1	& precipstartdate1, precipstartdate2, precipperiod,
79	heimbach	1.3	& runoffstartdate1, runoffstartdate2, runoffperiod,
80	heimbach	1.1	& ustressstartdate1, ustressstartdate2, ustressperiod,
81			& vstressstartdate1, vstressstartdate2, vstressperiod,
82			& uwindstartdate1, uwindstartdate2, uwindperiod,
83			& vwindstartdate1, vwindstartdate2, vwindperiod,
84			& swfluxstartdate1, swfluxstartdate2, swfluxperiod,
85			& lwfluxstartdate1, lwfluxstartdate2, lwfluxperiod,
86	dimitri	1.6	& swdownstartdate1, swdownstartdate2, swdownperiod,
87			& lwdownstartdate1, lwdownstartdate2, lwdownperiod,
88	heimbach	1.2	& obcsNstartdate1, obcsNstartdate2, obcsNperiod,
89			& obcsSstartdate1, obcsSstartdate2, obcsSperiod,
90			& obcsEstartdate1, obcsEstartdate2, obcsEperiod,
91			& obcsWstartdate1, obcsWstartdate2, obcsWperiod,
92	heimbach	1.3	&apressurestartdate1,apressurestartdate2,apressureperiod,
93	heimbach	1.1	& hfluxfile, atempfile, aqhfile,
94	heimbach	1.3	& sfluxfile, precipfile, runofffile,
95	dimitri	1.5	& ustressfile, vstressfile, evapfile,
96	heimbach	1.3	& uwindfile, vwindfile,
97			& swfluxfile, lwfluxfile, apressurefile,
98	dimitri	1.6	& swdownfile, lwdownfile,
99	heimbach	1.3	& exf_iprec, exf_yftype,
100	dimitri	1.6	& exf_inscal_hfl, exf_inscal_ust, exf_inscal_vst,
101			& exf_inscal_swf, exf_inscal_sst, exf_inscal_sss,
102			& exf_inscal_sfl, exf_inscal_runoff,exf_inscal_apressure,
103			& exf_outscal_hfl, exf_outscal_ust,exf_outscal_vst,
104			& exf_outscal_swf, exf_outscal_sst,exf_outscal_sss,
105			& exf_outscal_sfl, exf_outscal_apressure,
106	heimbach	1.7	& hfluxconst,
107			& atempconst,
108			& aqhconst,
109			& sfluxconst,
110			& evapconst,
111			& precipconst,
112			& runoffconst,
113			& ustressconst,
114			& vstressconst,
115			& uwindconst,
116			& vwindconst,
117			& swfluxconst,
118			& lwfluxconst,
119			& swdownconst,
120			& lwdownconst,
121			& apressureconst,
122			& EXFwindOnBgrid
123	heimbach	1.1
124			_BEGIN_MASTER(mythid)
125
126			c Set default values.
127
128			c Calendar data.
129			hfluxstartdate1 = 0
130			hfluxstartdate2 = 0
131			hfluxperiod = 0.0 _d 0
132	heimbach	1.7	hfluxconst = 0.0 _d 0
133	heimbach	1.1
134			atempstartdate1 = 0
135			atempstartdate2 = 0
136			atempperiod = 0.0 _d 0
137	heimbach	1.7	atempconst = 0.0 _d 0
138	heimbach	1.1
139			aqhstartdate1 = 0
140			aqhstartdate2 = 0
141			aqhperiod = 0.0 _d 0
142	heimbach	1.7	aqhconst = 0.0 _d 0
143	heimbach	1.1
144			sfluxstartdate1 = 0
145			sfluxstartdate2 = 0
146			sfluxperiod = 0.0 _d 0
147	heimbach	1.7	sfluxconst = 0.0 _d 0
148	heimbach	1.1
149	dimitri	1.5	evapstartdate1 = 0
150			evapstartdate2 = 0
151			evapperiod = 0.0 _d 0
152	heimbach	1.7	evapconst = 0.0 _d 0
153	dimitri	1.5
154	heimbach	1.1	precipstartdate1 = 0
155			precipstartdate2 = 0
156			precipperiod = 0.0 _d 0
157	heimbach	1.7	precipconst = 0.0 _d 0
158	heimbach	1.1
159	heimbach	1.3	runoffstartdate1 = 0
160			runoffstartdate2 = 0
161			runoffperiod = 0.0 _d 0
162	heimbach	1.7	runoffconst = 0.0 _d 0
163	heimbach	1.3
164	heimbach	1.1	ustressstartdate1 = 0
165			ustressstartdate2 = 0
166			ustressperiod = 0.0 _d 0
167	heimbach	1.7	ustressconst = 0.0 _d 0
168	heimbach	1.1
169			vstressstartdate1 = 0
170			vstressstartdate2 = 0
171			vstressperiod = 0.0 _d 0
172	heimbach	1.7	vstressconst = 0.0 _d 0
173	heimbach	1.1
174			uwindstartdate1 = 0
175			uwindstartdate2 = 0
176			uwindperiod = 0.0 _d 0
177	heimbach	1.7	uwindconst = 0.0 _d 0
178	heimbach	1.1
179			vwindstartdate1 = 0
180			vwindstartdate2 = 0
181			vwindperiod = 0.0 _d 0
182	heimbach	1.7	vwindconst = 0.0 _d 0
183	heimbach	1.1
184			swfluxstartdate1 = 0
185			swfluxstartdate2 = 0
186			swfluxperiod = 0.0 _d 0
187	heimbach	1.7	swfluxconst = 0.0 _d 0
188	heimbach	1.1
189			lwfluxstartdate1 = 0
190			lwfluxstartdate2 = 0
191			lwfluxperiod = 0.0 _d 0
192	heimbach	1.7	lwfluxconst = 0.0 _d 0
193	heimbach	1.1
194	dimitri	1.6	swdownstartdate1 = 0
195			swdownstartdate2 = 0
196			swdownperiod = 0.0 _d 0
197	heimbach	1.7	swdownconst = 0.0 _d 0
198	dimitri	1.6
199			lwdownstartdate1 = 0
200			lwdownstartdate2 = 0
201			lwdownperiod = 0.0 _d 0
202	heimbach	1.7	lwdownconst = 0.0 _d 0
203	dimitri	1.6
204	heimbach	1.2	obcsNstartdate1 = 0
205			obcsNstartdate2 = 0
206			obcsNperiod = 0.0 _d 0
207
208			obcsSstartdate1 = 0
209			obcsSstartdate2 = 0
210			obcsSperiod = 0.0 _d 0
211
212			obcsEstartdate1 = 0
213			obcsEstartdate2 = 0
214			obcsEperiod = 0.0 _d 0
215
216			obcsWstartdate1 = 0
217			obcsWstartdate2 = 0
218			obcsWperiod = 0.0 _d 0
219	heimbach	1.1
220	heimbach	1.3	apressurestartdate1 = 0
221			apressurestartdate2 = 0
222			apressureperiod = 0.0 _d 0
223	heimbach	1.7	apressureconst = 0.0 _d 0
224	cheisey	1.4
225			repeatPeriod = 0.0 _d 0
226
227	heimbach	1.1	c Data files.
228			hfluxfile = ' '
229			atempfile = ' '
230			aqhfile = ' '
231	dimitri	1.5	evapfile = ' '
232	heimbach	1.1	precipfile = ' '
233			sfluxfile = ' '
234	heimbach	1.3	runofffile = ' '
235	heimbach	1.1	ustressfile = ' '
236			vstressfile = ' '
237			uwindfile = ' '
238			vwindfile = ' '
239			swfluxfile = ' '
240			lwfluxfile = ' '
241	dimitri	1.6	swdownfile = ' '
242			lwdownfile = ' '
243	heimbach	1.3	apressurefile = ' '
244	heimbach	1.1
245			c Initialise the date arrays.
246			do i = 1,4
247			hfluxstartdate(i) = 0
248			atempstartdate(i) = 0
249			aqhstartdate(i) = 0
250	dimitri	1.5	evapstartdate(i) = 0
251	heimbach	1.1	precipstartdate(i) = 0
252			sfluxstartdate(i) = 0
253	heimbach	1.3	runoffstartdate(i) = 0
254	heimbach	1.1	ustressstartdate(i) = 0
255			vstressstartdate(i) = 0
256			uwindstartdate(i) = 0
257			vwindstartdate(i) = 0
258			swfluxstartdate(i) = 0
259			lwfluxstartdate(i) = 0
260	dimitri	1.6	swdownstartdate(i) = 0
261			lwdownstartdate(i) = 0
262	heimbach	1.2	obcsNstartdate(i) = 0
263			obcsSstartdate(i) = 0
264			obcsEstartdate(i) = 0
265			obcsWstartdate(i) = 0
266	heimbach	1.3	apressurestartdate(i)= 0
267	heimbach	1.1	enddo
268
269			c Initialise file type and field precision
270			exf_iprec = 32
271			exf_yftype = 'RL'
272
273	dimitri	1.6	c Input scaling factors.
274			exf_inscal_hfl = 1. _d 0
275			exf_inscal_ust = 1. _d 0
276			exf_inscal_vst = 1. _d 0
277			exf_inscal_swf = 1. _d 0
278			exf_inscal_sst = 1. _d 0
279			exf_inscal_sss = 1. _d 0
280			exf_inscal_apressure = 1. _d 0
281			exf_inscal_sfl = 1. _d 0
282			cds convert runoff from m/yr to m/s and change sign.
283			exf_inscal_runoff = -1.0/(86400.*365.0)
284
285			c Output scaling factors.
286			exf_outscal_hfl = 1. _d 0
287			exf_outscal_ust = 1. _d 0
288			exf_outscal_vst = 1. _d 0
289			exf_outscal_swf = 1. _d 0
290			exf_outscal_sst = 1. _d 0
291			exf_outscal_sss = 1. _d 0
292			exf_outscal_apressure= 1. _d 0
293			exf_outscal_sfl = 1. _d 0
294			c
295			EXFwindOnBgrid = .FALSE.
296	heimbach	1.3
297	heimbach	1.1	c Check for the availability of the right calendar version.
298			if ( calendarversion .ne. usescalendarversion ) then
299			print*,' exf_readparms: You are not using the appropriate'
300			print*,' version of the calendar package.'
301			print*
302			print*,' You are using Calendar version: ', calendarversion
303			print*,' Please use Calendar version: ', usescalendarversion
304			stop ' stopped in exf_readparms.'
305			endif
306
307			c Next, read the forcing data file.
308			call nml_filter( 'data.exf', scrunit1, myThid )
309			if (scrunit1 .eq. 0) then
310			stop 'exf_readparms: reading namelist failed'
311			end if
312			read( scrunit1, nml = exf_nml )
313			close( scrunit1 )
314
315			c Complete the start date specifications for the forcing
316			c fields to get a complete calendar date array.
317
318			c check for consistency
319
320			if (.NOT.
321			& (exf_iprec .EQ. 32 .OR. exf_iprec .EQ. 64)
322			& ) then
323			stop 'stop in exf_readparms: value of exf_iprec not allowed'
324			else if (.NOT.
325			& (exf_yftype .EQ. 'RS' .OR.
326			& exf_yftype .EQ. 'RL')
327			& ) then
328			stop 'stop in exf_readparms: value of exf_yftype not allowed'
329			end if
330	dimitri	1.5
331	dimitri	1.6	#ifdef ALLOW_ATM_WIND
332			call cal_FullDate( uwindstartdate1, uwindstartdate2,
333			& uwindstartdate , mythid )
334			call cal_FullDate( vwindstartdate1, vwindstartdate2,
335			& vwindstartdate , mythid )
336			#else
337			call cal_FullDate( ustressstartdate1, ustressstartdate2,
338			& ustressstartdate , mythid )
339			call cal_FullDate( vstressstartdate1, vstressstartdate2,
340			& vstressstartdate , mythid )
341	heimbach	1.3	#endif
342
343	heimbach	1.1	#ifdef ALLOW_ATM_TEMP
344	dimitri	1.6	call cal_FullDate( atempstartdate1, atempstartdate2,
345			& atempstartdate , mythid )
346			call cal_FullDate( aqhstartdate1, aqhstartdate2,
347			& aqhstartdate , mythid )
348			call cal_FullDate( lwfluxstartdate1, lwfluxstartdate2,
349			& lwfluxstartdate , mythid )
350			call cal_FullDate( precipstartdate1, precipstartdate2,
351			& precipstartdate , mythid )
352	heimbach	1.1	#else
353	dimitri	1.6	call cal_FullDate( hfluxstartdate1, hfluxstartdate2,
354			& hfluxstartdate , mythid )
355			call cal_FullDate( sfluxstartdate1, sfluxstartdate2,
356			& sfluxstartdate , mythid )
357	heimbach	1.1	#endif
358
359	dimitri	1.6	#if defined(ALLOW_ATM_TEMP) \|\| defined(SHORTWAVE_HEATING)
360			call cal_FullDate( swfluxstartdate1, swfluxstartdate2,
361			& swfluxstartdate , mythid )
362	heimbach	1.1	#endif
363
364	dimitri	1.6	#ifdef EXF_READ_EVAP
365			call cal_FullDate( evapstartdate1, evapstartdate2,
366			& evapstartdate , mythid )
367	heimbach	1.1	#endif
368
369	dimitri	1.6	#ifdef ALLOW_RUNOFF
370			call cal_FullDate( runoffstartdate1, runoffstartdate2,
371			& runoffstartdate , mythid )
372	heimbach	1.1	#endif
373
374	dimitri	1.6	#ifdef ALLOW_DOWNWARD_RADIATION
375			call cal_FullDate( swdownstartdate1, swdownstartdate2,
376			& swdownstartdate , mythid )
377			call cal_FullDate( lwdownstartdate1, lwdownstartdate2,
378			& lwdownstartdate , mythid )
379	heimbach	1.1	#endif
380
381			#ifdef ALLOW_OBCS
382	heimbach	1.2	#ifdef ALLOW_OBCS_NORTH
383	dimitri	1.6	call cal_FullDate( obcsNstartdate1, obcsNstartdate2,
384			& obcsNstartdate , mythid )
385	heimbach	1.2	#endif
386			#ifdef ALLOW_OBCS_SOUTH
387	dimitri	1.6	call cal_FullDate( obcsSstartdate1, obcsSstartdate2,
388			& obcsSstartdate , mythid )
389	heimbach	1.2	#endif
390			#ifdef ALLOW_OBCS_EAST
391	dimitri	1.6	call cal_FullDate( obcsEstartdate1, obcsEstartdate2,
392			& obcsEstartdate , mythid )
393	heimbach	1.2	#endif
394			#ifdef ALLOW_OBCS_WEST
395	dimitri	1.6	call cal_FullDate( obcsWstartdate1, obcsWstartdate2,
396			& obcsWstartdate , mythid )
397	heimbach	1.3	#endif
398	dimitri	1.6	#endif /* ALLOW_OBCS */
399	heimbach	1.3
400			#ifdef ATMOSPHERIC_LOADING
401	dimitri	1.6	call cal_FullDate(apressurestartdate1,apressurestartdate2,
402			& apressurestartdate , mythid )
403	heimbach	1.1	#endif
404
405			_END_MASTER( mythid )
406
407			_BARRIER
408
409			c-- Summarize the External forcing's setup.
410			call exf_summary( mythid )
411
412			c-- set climatology parameters
413			call exf_clim_readparms( mythid )
414
415			c-- summarize climatologic forcing configuration
416			call exf_clim_summary( mythid )
417
418			end