pkg/exf/exf_param.h

c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_param.h,v 1.15 2004/03/24 21:46:20 heimbach Exp $
c
c
c     ==================================================================
c     HEADER exf_param
c     ==================================================================
c
c     o Header file for the surface flux data. Used by the external
c       forcing package.
c
c     started: Christian Eckert eckert@mit.edu  30-Jun-1999
c
c     changed: Christian Eckert eckert@mit.edu  14-Jan-2000
c              - Restructured the original version in order to have a
c                better interface to the MITgcmUV.
c
c              Christian Eckert eckert@mit.edu  12-Feb-2000
c              - Changed some variables names (package prefix: exf_)
c
c              Patrick Heimbach, heimbach@mit.edu  04-May-2000
c              - included exf_iprec, exf_yftype to enable easy
c                switch between 32bit/64 bit data format
c
c              Patrick Heimbach, heimbach@mit.edu  01-May-2001
c              - added obcs parameters
c
c     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
c
c     ==================================================================
c     HEADER exf_param
c     ==================================================================

c     Calendar data.
      _RL     repeatPeriod

c     Sea-water albedo
      _RL     exf_albedo

c     Maximum absolute windstress, used to reset unreastically high
c     data values
      _RL     windstressmax

      integer hfluxstartdate1
      integer hfluxstartdate2
      _RL     hfluxstartdate
      _RL     hfluxperiod
      _RL     hfluxconst
      character*1 hfluxmask
      parameter(  hfluxmask = 's' )

      integer atempstartdate1
      integer atempstartdate2
      _RL     atempstartdate
      _RL     atempperiod
      _RL     atempconst
      character*1 atempmask
      parameter(  atempmask = 's' )

      integer aqhstartdate1
      integer aqhstartdate2
      _RL     aqhstartdate
      _RL     aqhperiod
      _RL     aqhconst
      character*1 aqhmask
      parameter(  aqhmask = 's' )

      integer sfluxstartdate1
      integer sfluxstartdate2
      _RL     sfluxstartdate
      _RL     sfluxperiod
      _RL     sfluxconst
      character*1 sfluxmask
      parameter(  sfluxmask = 's' )

      integer evapstartdate1
      integer evapstartdate2
      _RL     evapstartdate
      _RL     evapperiod
      _RL     evapconst
      character*1 evapmask
      parameter(  evapmask = 's' )

      integer precipstartdate1
      integer precipstartdate2
      _RL     precipstartdate
      _RL     precipperiod
      _RL     precipconst
      character*1 precipmask
      parameter(  precipmask = 's' )

      integer runoffstartdate1
      integer runoffstartdate2
      _RL     runoffstartdate
      _RL     runoffperiod
      _RL     runoffconst
      character*1 runoffmask
      parameter(  runoffmask = 's' )

      integer ustressstartdate1
      integer ustressstartdate2
      _RL     ustressstartdate
      _RL     ustressperiod
      _RL     ustressconst
      character*1 ustressmask
      parameter(  ustressmask = 'u' )

      integer vstressstartdate1
      integer vstressstartdate2
      _RL     vstressstartdate
      _RL     vstressperiod
      _RL     vstressconst
      character*1 vstressmask
      parameter(  vstressmask = 'v' )

      integer uwindstartdate1
      integer uwindstartdate2
      _RL     uwindstartdate
      _RL     uwindperiod
      _RL     uwindconst
      character*1 uwindmask
      parameter(  uwindmask = 'u' )

      integer vwindstartdate1
      integer vwindstartdate2
      _RL     vwindstartdate
      _RL     vwindperiod
      _RL     vwindconst
      character*1 vwindmask
      parameter(  vwindmask = 'v' )

      integer swfluxstartdate1
      integer swfluxstartdate2
      _RL     swfluxstartdate
      _RL     swfluxperiod
      _RL     swfluxconst
      character*1 swfluxmask
      parameter(  swfluxmask = 's' )

      integer lwfluxstartdate1
      integer lwfluxstartdate2
      _RL     lwfluxstartdate
      _RL     lwfluxperiod
      _RL     lwfluxconst
      character*1 lwfluxmask
      parameter(  lwfluxmask = 's' )

      integer swdownstartdate1
      integer swdownstartdate2
      _RL     swdownstartdate
      _RL     swdownperiod
      _RL     swdownconst
      character*1 swdownmask
      parameter(  swdownmask = 's' )

      integer lwdownstartdate1
      integer lwdownstartdate2
      _RL     lwdownstartdate
      _RL     lwdownperiod
      _RL     lwdownconst
      character*1 lwdownmask
      parameter(  lwdownmask = 's' )

      integer obcsNstartdate1
      integer obcsNstartdate2
      _RL     obcsNstartdate
      _RL     obcsNperiod

      integer obcsSstartdate1
      integer obcsSstartdate2
      _RL     obcsSstartdate
      _RL     obcsSperiod

      integer obcsEstartdate1
      integer obcsEstartdate2
      _RL     obcsEstartdate
      _RL     obcsEperiod

      integer obcsWstartdate1
      integer obcsWstartdate2
      _RL     obcsWstartdate
      _RL     obcsWperiod

      integer apressurestartdate1
      integer apressurestartdate2
      _RL     apressurestartdate
      _RL     apressureperiod
      _RL     apressureconst
      character*1 apressuremask
      parameter(  apressuremask = 's' )

c     File names.
      character*(128) hfluxfile
      character*(128) atempfile
      character*(128) aqhfile
      character*(128) evapfile
      character*(128) precipfile
      character*(128) sfluxfile
      character*(128) runofffile
      character*(128) ustressfile
      character*(128) vstressfile
      character*(128) uwindfile
      character*(128) vwindfile
      character*(128) swfluxfile
      character*(128) lwfluxfile
      character*(128) swdownfile
      character*(128) lwdownfile
      character*(128) apressurefile

      common /exf_param_i/
     &                          hfluxstartdate1,   hfluxstartdate2,
     &                          atempstartdate1,   atempstartdate2,
     &                          aqhstartdate1,     aqhstartdate2,
     &                          sfluxstartdate1,   sfluxstartdate2,
     &                          evapstartdate1,    evapstartdate2,
     &                          runoffstartdate1,  runoffstartdate2,
     &                          precipstartdate1,  precipstartdate2,
     &                          ustressstartdate1, ustressstartdate2,
     &                          vstressstartdate1, vstressstartdate2,
     &                          uwindstartdate1,   uwindstartdate2,
     &                          vwindstartdate1,   vwindstartdate2,
     &                          swfluxstartdate1,  swfluxstartdate2,
     &                          lwfluxstartdate1,  lwfluxstartdate2,
     &                          swdownstartdate1,  swdownstartdate2,
     &                          lwdownstartdate1,  lwdownstartdate2,
     &                          obcsNstartdate1,   obcsNstartdate2,
     &                          obcsSstartdate1,   obcsSstartdate2,
     &                          obcsEstartdate1,   obcsEstartdate2,
     &                          obcsWstartdate1,   obcsWstartdate2,
     &                          apressurestartdate1,apressurestartdate2

      common /exf_param_r/
     &                          windstressmax,
     &                          repeatPeriod,      exf_albedo,
     &                          hfluxperiod,       hfluxstartdate,
     &                          atempperiod,       atempstartdate,
     &                          aqhperiod,         aqhstartdate,
     &                          sfluxperiod,       sfluxstartdate,
     &                          evapperiod,        evapstartdate,
     &                          precipperiod,      precipstartdate,
     &                          runoffperiod,      runoffstartdate,
     &                          ustressperiod,     ustressstartdate,
     &                          vstressperiod,     vstressstartdate,
     &                          uwindperiod,       uwindstartdate,
     &                          vwindperiod,       vwindstartdate,
     &                          swfluxperiod,      swfluxstartdate,
     &                          lwfluxperiod,      lwfluxstartdate,
     &                          swdownperiod,      swdownstartdate,
     &                          lwdownperiod,      lwdownstartdate,
     &                          obcsNperiod,       obcsNstartdate,
     &                          obcsSperiod,       obcsSstartdate,
     &                          obcsEperiod,       obcsEstartdate,
     &                          obcsWperiod,       obcsWstartdate,
     &                          apressureperiod,   apressurestartdate,
     &                          hfluxconst,
     &                          atempconst,
     &                          aqhconst,
     &                          sfluxconst,
     &                          evapconst,
     &                          precipconst,
     &                          runoffconst,
     &                          ustressconst,
     &                          vstressconst,
     &                          uwindconst,
     &                          vwindconst,
     &                          swfluxconst,
     &                          lwfluxconst,
     &                          swdownconst,
     &                          lwdownconst,
     &                          apressureconst

      common /exf_param_c/
     &                          hfluxfile,
     &                          atempfile,
     &                          aqhfile,
     &                          sfluxfile,
     &                          evapfile,
     &                          precipfile,
     &                          runofffile,
     &                          ustressfile,
     &                          vstressfile,
     &                          uwindfile,
     &                          vwindfile,
     &                          swfluxfile,
     &                          lwfluxfile,
     &                          swdownfile,
     &                          lwdownfile,
     &                          apressurefile

c     file precision and field type

      common /exf_param_type/ 
     &                        exf_iprec,
     &                        exf_yftype

      integer exf_iprec
      character*(2) exf_yftype

c     exf_inscal_*      input scaling factors
c     exf_offset_atemp  input air temperature offset
c                       (for conversion from C to K, if needed)
c     exf_outscale_*    output scaling factors

      _RL     exf_inscal_hflux
      _RL     exf_inscal_sflux
      _RL     exf_inscal_ustress
      _RL     exf_inscal_vstress
      _RL     exf_inscal_uwind
      _RL     exf_inscal_vwind
      _RL     exf_inscal_swflux
      _RL     exf_inscal_lwflux
      _RL     exf_inscal_precip
      _RL     exf_inscal_sst
      _RL     exf_inscal_sss
      _RL     exf_inscal_atemp
      _RL     exf_offset_atemp
      _RL     exf_inscal_aqh
      _RL     exf_inscal_evap
      _RL     exf_inscal_apressure
      _RL     exf_inscal_runoff
      _RL     exf_inscal_swdown
      _RL     exf_inscal_lwdown

      _RL     exf_outscal_hflux
      _RL     exf_outscal_sflux
      _RL     exf_outscal_ustress
      _RL     exf_outscal_vstress
      _RL     exf_outscal_swflux
      _RL     exf_outscal_sst
      _RL     exf_outscal_sss
      _RL     exf_outscal_apressure

      common /exf_param_scal/
     &                      exf_inscal_hflux
     &                    , exf_inscal_sflux
     &                    , exf_inscal_ustress
     &                    , exf_inscal_vstress
     &                    , exf_inscal_uwind
     &                    , exf_inscal_vwind
     &                    , exf_inscal_swflux
     &                    , exf_inscal_lwflux
     &                    , exf_inscal_precip
     &                    , exf_inscal_sst
     &                    , exf_inscal_sss
     &                    , exf_inscal_atemp
     &                    , exf_offset_atemp
     &                    , exf_inscal_aqh
     &                    , exf_inscal_evap
     &                    , exf_inscal_apressure
     &                    , exf_inscal_runoff
     &                    , exf_inscal_swdown
     &                    , exf_inscal_lwdown
     &                    , exf_outscal_hflux
     &                    , exf_outscal_sflux
     &                    , exf_outscal_ustress
     &                    , exf_outscal_vstress
     &                    , exf_outscal_swflux
     &                    , exf_outscal_sst
     &                    , exf_outscal_sss
     &                    , exf_outscal_apressure

#ifndef USE_EXF_INTERPOLATION
c-- set dummy dimension 1
       integer MAX_LAT_INC
       parameter(MAX_LAT_INC = 1)
#else
c for lat interpolation, arraysize currently set to 250 max data values
       integer MAX_LAT_INC
       parameter(MAX_LAT_INC = 250)
      _RL ustress_lon0, ustress_lon_inc
      _RL ustress_lat0, ustress_lat_inc(MAX_LAT_INC)
      INTEGER ustress_nlon, ustress_nlat
      _RL vstress_lon0, vstress_lon_inc
      _RL vstress_lat0, vstress_lat_inc(MAX_LAT_INC)
      INTEGER vstress_nlon, vstress_nlat
      _RL hflux_lon0, hflux_lon_inc
      _RL hflux_lat0, hflux_lat_inc(MAX_LAT_INC)
      INTEGER hflux_nlon, hflux_nlat
      _RL sflux_lon0, sflux_lon_inc
      _RL sflux_lat0, sflux_lat_inc(MAX_LAT_INC)
      INTEGER sflux_nlon, sflux_nlat
      _RL swflux_lon0, swflux_lon_inc
      _RL swflux_lat0, swflux_lat_inc(MAX_LAT_INC)
      INTEGER swflux_nlon, swflux_nlat
      _RL runoff_lon0, runoff_lon_inc
      _RL runoff_lat0, runoff_lat_inc(MAX_LAT_INC)
      INTEGER runoff_nlon, runoff_nlat
      _RL atemp_lon0, atemp_lon_inc
      _RL atemp_lat0, atemp_lat_inc(MAX_LAT_INC)
      INTEGER atemp_nlon, atemp_nlat
      _RL aqh_lon0, aqh_lon_inc
      _RL aqh_lat0, aqh_lat_inc(MAX_LAT_INC)
      INTEGER aqh_nlon, aqh_nlat
      _RL evap_lon0, evap_lon_inc
      _RL evap_lat0, evap_lat_inc(MAX_LAT_INC)
      INTEGER evap_nlon, evap_nlat
      _RL precip_lon0, precip_lon_inc
      _RL precip_lat0, precip_lat_inc(MAX_LAT_INC)
      INTEGER precip_nlon, precip_nlat
      _RL uwind_lon0, uwind_lon_inc
      _RL uwind_lat0, uwind_lat_inc(MAX_LAT_INC)
      INTEGER uwind_nlon, uwind_nlat
      _RL vwind_lon0, vwind_lon_inc
      _RL vwind_lat0, vwind_lat_inc(MAX_LAT_INC)
      INTEGER vwind_nlon, vwind_nlat
      _RL lwflux_lon0, lwflux_lon_inc
      _RL lwflux_lat0, lwflux_lat_inc(MAX_LAT_INC)
      INTEGER lwflux_nlon, lwflux_nlat
      _RL swdown_lon0, swdown_lon_inc
      _RL swdown_lat0, swdown_lat_inc(MAX_LAT_INC)
      INTEGER swdown_nlon, swdown_nlat
      _RL lwdown_lon0, lwdown_lon_inc
      _RL lwdown_lat0, lwdown_lat_inc(MAX_LAT_INC)
      INTEGER lwdown_nlon, lwdown_nlat
      _RL apressure_lon0,apressure_lon_inc
      _RL apressure_lat0,apressure_lat_inc(MAX_LAT_INC)
      INTEGER apressure_nlon,apressure_nlat

      common /exf_interpolation/
     & ustress_lon0, ustress_lon_inc,
     & ustress_lat0, ustress_lat_inc,
     & ustress_nlon, ustress_nlat,
     & vstress_lon0, vstress_lon_inc,
     & vstress_lat0, vstress_lat_inc,
     & vstress_nlon, vstress_nlat,
     & hflux_lon0, hflux_lon_inc,
     & hflux_lat0, hflux_lat_inc,
     & hflux_nlon, hflux_nlat,
     & sflux_lon0, sflux_lon_inc,
     & sflux_lat0, sflux_lat_inc,
     & sflux_nlon, sflux_nlat,
     & swflux_lon0, swflux_lon_inc,
     & swflux_lat0, swflux_lat_inc,
     & swflux_nlon, swflux_nlat,
     & runoff_lon0, runoff_lon_inc,
     & runoff_lat0, runoff_lat_inc,
     & runoff_nlon, runoff_nlat,
     & atemp_lon0, atemp_lon_inc,
     & atemp_lat0, atemp_lat_inc,
     & atemp_nlon, atemp_nlat,
     & aqh_lon0, aqh_lon_inc,
     & aqh_lat0, aqh_lat_inc,
     & aqh_nlon, aqh_nlat,
     & evap_lon0, evap_lon_inc,
     & evap_lat0, evap_lat_inc,
     & evap_nlon, evap_nlat,
     & precip_lon0, precip_lon_inc,
     & precip_lat0, precip_lat_inc,
     & precip_nlon, precip_nlat,
     & uwind_lon0, uwind_lon_inc,
     & uwind_lat0, uwind_lat_inc,
     & uwind_nlon, uwind_nlat,
     & vwind_lon0, vwind_lon_inc,
     & vwind_lat0, vwind_lat_inc,
     & vwind_nlon, vwind_nlat,
     & lwflux_lon0, lwflux_lon_inc,
     & lwflux_lat0, lwflux_lat_inc,
     & lwflux_nlon, lwflux_nlat,
     & swdown_lon0, swdown_lon_inc,
     & swdown_lat0, swdown_lat_inc,
     & swdown_nlon, swdown_nlat,
     & lwdown_lon0, lwdown_lon_inc,
     & lwdown_lat0, lwdown_lat_inc,
     & lwdown_nlon, lwdown_nlat,
     & apressure_lon0,apressure_lon_inc,
     & apressure_lat0,apressure_lat_inc,
     & apressure_nlon,apressure_nlat
#endif
1	mlosch	1.16	c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_param.h,v 1.15 2004/03/24 21:46:20 heimbach Exp $
2	heimbach	1.1	c
3			c
4			c ==================================================================
5			c HEADER exf_param
6			c ==================================================================
7			c
8			c o Header file for the surface flux data. Used by the external
9			c forcing package.
10			c
11			c started: Christian Eckert eckert@mit.edu 30-Jun-1999
12			c
13			c changed: Christian Eckert eckert@mit.edu 14-Jan-2000
14			c - Restructured the original version in order to have a
15			c better interface to the MITgcmUV.
16			c
17			c Christian Eckert eckert@mit.edu 12-Feb-2000
18			c - Changed some variables names (package prefix: exf_)
19			c
20			c Patrick Heimbach, heimbach@mit.edu 04-May-2000
21			c - included exf_iprec, exf_yftype to enable easy
22			c switch between 32bit/64 bit data format
23			c
24			c Patrick Heimbach, heimbach@mit.edu 01-May-2001
25			c - added obcs parameters
26			c
27	dimitri	1.6	c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
28	dimitri	1.5	c
29	heimbach	1.1	c ==================================================================
30			c HEADER exf_param
31			c ==================================================================
32
33			c Calendar data.
34	cheisey	1.4	_RL repeatPeriod
35
36	dimitri	1.14	c Sea-water albedo
37			_RL exf_albedo
38
39	mlosch	1.16	c Maximum absolute windstress, used to reset unreastically high
40			c data values
41			_RL windstressmax
42
43	heimbach	1.1	integer hfluxstartdate1
44			integer hfluxstartdate2
45	dimitri	1.11	_RL hfluxstartdate
46	heimbach	1.1	_RL hfluxperiod
47	heimbach	1.7	_RL hfluxconst
48	heimbach	1.1	character*1 hfluxmask
49			parameter( hfluxmask = 's' )
50
51			integer atempstartdate1
52			integer atempstartdate2
53	dimitri	1.11	_RL atempstartdate
54	heimbach	1.1	_RL atempperiod
55	heimbach	1.7	_RL atempconst
56	heimbach	1.1	character*1 atempmask
57			parameter( atempmask = 's' )
58
59			integer aqhstartdate1
60			integer aqhstartdate2
61	dimitri	1.11	_RL aqhstartdate
62	heimbach	1.1	_RL aqhperiod
63	heimbach	1.7	_RL aqhconst
64	heimbach	1.1	character*1 aqhmask
65			parameter( aqhmask = 's' )
66
67			integer sfluxstartdate1
68			integer sfluxstartdate2
69	dimitri	1.11	_RL sfluxstartdate
70	heimbach	1.1	_RL sfluxperiod
71	heimbach	1.7	_RL sfluxconst
72	heimbach	1.1	character*1 sfluxmask
73			parameter( sfluxmask = 's' )
74
75	dimitri	1.5	integer evapstartdate1
76			integer evapstartdate2
77	dimitri	1.11	_RL evapstartdate
78	dimitri	1.5	_RL evapperiod
79	heimbach	1.7	_RL evapconst
80	dimitri	1.5	character*1 evapmask
81			parameter( evapmask = 's' )
82
83	heimbach	1.1	integer precipstartdate1
84			integer precipstartdate2
85	dimitri	1.11	_RL precipstartdate
86	heimbach	1.1	_RL precipperiod
87	heimbach	1.7	_RL precipconst
88	heimbach	1.1	character*1 precipmask
89			parameter( precipmask = 's' )
90
91	heimbach	1.3	integer runoffstartdate1
92			integer runoffstartdate2
93	dimitri	1.11	_RL runoffstartdate
94	heimbach	1.3	_RL runoffperiod
95	heimbach	1.7	_RL runoffconst
96	heimbach	1.3	character*1 runoffmask
97			parameter( runoffmask = 's' )
98
99	heimbach	1.1	integer ustressstartdate1
100			integer ustressstartdate2
101	dimitri	1.11	_RL ustressstartdate
102	heimbach	1.1	_RL ustressperiod
103	heimbach	1.7	_RL ustressconst
104	heimbach	1.1	character*1 ustressmask
105			parameter( ustressmask = 'u' )
106
107			integer vstressstartdate1
108			integer vstressstartdate2
109	dimitri	1.11	_RL vstressstartdate
110	heimbach	1.1	_RL vstressperiod
111	heimbach	1.7	_RL vstressconst
112	heimbach	1.1	character*1 vstressmask
113			parameter( vstressmask = 'v' )
114
115			integer uwindstartdate1
116			integer uwindstartdate2
117	dimitri	1.11	_RL uwindstartdate
118	heimbach	1.1	_RL uwindperiod
119	heimbach	1.7	_RL uwindconst
120	heimbach	1.1	character*1 uwindmask
121			parameter( uwindmask = 'u' )
122
123			integer vwindstartdate1
124			integer vwindstartdate2
125	dimitri	1.11	_RL vwindstartdate
126	heimbach	1.1	_RL vwindperiod
127	heimbach	1.7	_RL vwindconst
128	heimbach	1.1	character*1 vwindmask
129			parameter( vwindmask = 'v' )
130
131			integer swfluxstartdate1
132			integer swfluxstartdate2
133	dimitri	1.11	_RL swfluxstartdate
134	heimbach	1.1	_RL swfluxperiod
135	heimbach	1.7	_RL swfluxconst
136	heimbach	1.1	character*1 swfluxmask
137			parameter( swfluxmask = 's' )
138
139			integer lwfluxstartdate1
140			integer lwfluxstartdate2
141	dimitri	1.11	_RL lwfluxstartdate
142	heimbach	1.1	_RL lwfluxperiod
143	heimbach	1.7	_RL lwfluxconst
144	heimbach	1.1	character*1 lwfluxmask
145			parameter( lwfluxmask = 's' )
146
147	dimitri	1.6	integer swdownstartdate1
148			integer swdownstartdate2
149	dimitri	1.11	_RL swdownstartdate
150	dimitri	1.6	_RL swdownperiod
151	heimbach	1.7	_RL swdownconst
152	dimitri	1.6	character*1 swdownmask
153			parameter( swdownmask = 's' )
154
155			integer lwdownstartdate1
156			integer lwdownstartdate2
157	dimitri	1.11	_RL lwdownstartdate
158	dimitri	1.6	_RL lwdownperiod
159	heimbach	1.7	_RL lwdownconst
160	dimitri	1.6	character*1 lwdownmask
161			parameter( lwdownmask = 's' )
162
163	heimbach	1.2	integer obcsNstartdate1
164			integer obcsNstartdate2
165	dimitri	1.11	_RL obcsNstartdate
166	heimbach	1.2	_RL obcsNperiod
167
168			integer obcsSstartdate1
169			integer obcsSstartdate2
170	dimitri	1.11	_RL obcsSstartdate
171	heimbach	1.2	_RL obcsSperiod
172
173			integer obcsEstartdate1
174			integer obcsEstartdate2
175	dimitri	1.11	_RL obcsEstartdate
176	heimbach	1.2	_RL obcsEperiod
177
178			integer obcsWstartdate1
179			integer obcsWstartdate2
180	dimitri	1.11	_RL obcsWstartdate
181	heimbach	1.2	_RL obcsWperiod
182	heimbach	1.1
183	heimbach	1.3	integer apressurestartdate1
184			integer apressurestartdate2
185	dimitri	1.11	_RL apressurestartdate
186	heimbach	1.3	_RL apressureperiod
187	heimbach	1.7	_RL apressureconst
188	heimbach	1.3	character*1 apressuremask
189			parameter( apressuremask = 's' )
190
191	heimbach	1.1	c File names.
192			character*(128) hfluxfile
193			character*(128) atempfile
194			character*(128) aqhfile
195	dimitri	1.5	character*(128) evapfile
196	heimbach	1.1	character*(128) precipfile
197			character*(128) sfluxfile
198	heimbach	1.3	character*(128) runofffile
199	heimbach	1.1	character*(128) ustressfile
200			character*(128) vstressfile
201			character*(128) uwindfile
202			character*(128) vwindfile
203			character*(128) swfluxfile
204			character*(128) lwfluxfile
205	dimitri	1.6	character*(128) swdownfile
206			character*(128) lwdownfile
207	heimbach	1.3	character*(128) apressurefile
208	heimbach	1.1
209			common /exf_param_i/
210	heimbach	1.2	& hfluxstartdate1, hfluxstartdate2,
211			& atempstartdate1, atempstartdate2,
212			& aqhstartdate1, aqhstartdate2,
213			& sfluxstartdate1, sfluxstartdate2,
214	dimitri	1.5	& evapstartdate1, evapstartdate2,
215	heimbach	1.3	& runoffstartdate1, runoffstartdate2,
216	heimbach	1.2	& precipstartdate1, precipstartdate2,
217			& ustressstartdate1, ustressstartdate2,
218			& vstressstartdate1, vstressstartdate2,
219			& uwindstartdate1, uwindstartdate2,
220			& vwindstartdate1, vwindstartdate2,
221			& swfluxstartdate1, swfluxstartdate2,
222			& lwfluxstartdate1, lwfluxstartdate2,
223	dimitri	1.6	& swdownstartdate1, swdownstartdate2,
224			& lwdownstartdate1, lwdownstartdate2,
225	heimbach	1.2	& obcsNstartdate1, obcsNstartdate2,
226			& obcsSstartdate1, obcsSstartdate2,
227			& obcsEstartdate1, obcsEstartdate2,
228			& obcsWstartdate1, obcsWstartdate2,
229	dimitri	1.11	& apressurestartdate1,apressurestartdate2
230	heimbach	1.1
231			common /exf_param_r/
232	mlosch	1.16	& windstressmax,
233	dimitri	1.14	& repeatPeriod, exf_albedo,
234	dimitri	1.11	& hfluxperiod, hfluxstartdate,
235			& atempperiod, atempstartdate,
236			& aqhperiod, aqhstartdate,
237			& sfluxperiod, sfluxstartdate,
238			& evapperiod, evapstartdate,
239			& precipperiod, precipstartdate,
240			& runoffperiod, runoffstartdate,
241			& ustressperiod, ustressstartdate,
242			& vstressperiod, vstressstartdate,
243			& uwindperiod, uwindstartdate,
244			& vwindperiod, vwindstartdate,
245			& swfluxperiod, swfluxstartdate,
246			& lwfluxperiod, lwfluxstartdate,
247			& swdownperiod, swdownstartdate,
248			& lwdownperiod, lwdownstartdate,
249			& obcsNperiod, obcsNstartdate,
250			& obcsSperiod, obcsSstartdate,
251			& obcsEperiod, obcsEstartdate,
252			& obcsWperiod, obcsWstartdate,
253			& apressureperiod, apressurestartdate,
254	heimbach	1.7	& hfluxconst,
255			& atempconst,
256			& aqhconst,
257			& sfluxconst,
258			& evapconst,
259			& precipconst,
260			& runoffconst,
261			& ustressconst,
262			& vstressconst,
263			& uwindconst,
264			& vwindconst,
265			& swfluxconst,
266			& lwfluxconst,
267			& swdownconst,
268			& lwdownconst,
269			& apressureconst
270	heimbach	1.1
271			common /exf_param_c/
272			& hfluxfile,
273			& atempfile,
274			& aqhfile,
275			& sfluxfile,
276	dimitri	1.5	& evapfile,
277	heimbach	1.1	& precipfile,
278	heimbach	1.3	& runofffile,
279	heimbach	1.1	& ustressfile,
280			& vstressfile,
281			& uwindfile,
282			& vwindfile,
283			& swfluxfile,
284	heimbach	1.3	& lwfluxfile,
285	dimitri	1.6	& swdownfile,
286			& lwdownfile,
287	heimbach	1.3	& apressurefile
288	heimbach	1.1
289			c file precision and field type
290
291			common /exf_param_type/
292			& exf_iprec,
293			& exf_yftype
294
295			integer exf_iprec
296			character*(2) exf_yftype
297	heimbach	1.3
298	dimitri	1.13	c exf_inscal_* input scaling factors
299			c exf_offset_atemp input air temperature offset
300			c (for conversion from C to K, if needed)
301			c exf_outscale_* output scaling factors
302	heimbach	1.3
303	heimbach	1.8	_RL exf_inscal_hflux
304			_RL exf_inscal_sflux
305			_RL exf_inscal_ustress
306			_RL exf_inscal_vstress
307			_RL exf_inscal_uwind
308			_RL exf_inscal_vwind
309			_RL exf_inscal_swflux
310			_RL exf_inscal_lwflux
311			_RL exf_inscal_precip
312	dimitri	1.6	_RL exf_inscal_sst
313			_RL exf_inscal_sss
314	heimbach	1.8	_RL exf_inscal_atemp
315	dimitri	1.13	_RL exf_offset_atemp
316	heimbach	1.8	_RL exf_inscal_aqh
317			_RL exf_inscal_evap
318	dimitri	1.6	_RL exf_inscal_apressure
319			_RL exf_inscal_runoff
320	heimbach	1.8	_RL exf_inscal_swdown
321			_RL exf_inscal_lwdown
322
323			_RL exf_outscal_hflux
324			_RL exf_outscal_sflux
325			_RL exf_outscal_ustress
326			_RL exf_outscal_vstress
327			_RL exf_outscal_swflux
328	dimitri	1.6	_RL exf_outscal_sst
329			_RL exf_outscal_sss
330			_RL exf_outscal_apressure
331	heimbach	1.3
332			common /exf_param_scal/
333	heimbach	1.8	& exf_inscal_hflux
334			& , exf_inscal_sflux
335			& , exf_inscal_ustress
336			& , exf_inscal_vstress
337			& , exf_inscal_uwind
338			& , exf_inscal_vwind
339			& , exf_inscal_swflux
340			& , exf_inscal_lwflux
341			& , exf_inscal_precip
342	dimitri	1.6	& , exf_inscal_sst
343			& , exf_inscal_sss
344	heimbach	1.8	& , exf_inscal_atemp
345	dimitri	1.13	& , exf_offset_atemp
346	heimbach	1.8	& , exf_inscal_aqh
347			& , exf_inscal_evap
348	dimitri	1.6	& , exf_inscal_apressure
349			& , exf_inscal_runoff
350	heimbach	1.8	& , exf_inscal_swdown
351			& , exf_inscal_lwdown
352			& , exf_outscal_hflux
353			& , exf_outscal_sflux
354			& , exf_outscal_ustress
355			& , exf_outscal_vstress
356			& , exf_outscal_swflux
357	dimitri	1.6	& , exf_outscal_sst
358			& , exf_outscal_sss
359			& , exf_outscal_apressure
360
361	heimbach	1.15	#ifndef USE_EXF_INTERPOLATION
362			c-- set dummy dimension 1
363			integer MAX_LAT_INC
364			parameter(MAX_LAT_INC = 1)
365			#else
366	dimitri	1.9	c for lat interpolation, arraysize currently set to 250 max data values
367			integer MAX_LAT_INC
368			parameter(MAX_LAT_INC = 250)
369			_RL ustress_lon0, ustress_lon_inc
370			_RL ustress_lat0, ustress_lat_inc(MAX_LAT_INC)
371			INTEGER ustress_nlon, ustress_nlat
372			_RL vstress_lon0, vstress_lon_inc
373			_RL vstress_lat0, vstress_lat_inc(MAX_LAT_INC)
374			INTEGER vstress_nlon, vstress_nlat
375			_RL hflux_lon0, hflux_lon_inc
376			_RL hflux_lat0, hflux_lat_inc(MAX_LAT_INC)
377			INTEGER hflux_nlon, hflux_nlat
378			_RL sflux_lon0, sflux_lon_inc
379			_RL sflux_lat0, sflux_lat_inc(MAX_LAT_INC)
380			INTEGER sflux_nlon, sflux_nlat
381			_RL swflux_lon0, swflux_lon_inc
382			_RL swflux_lat0, swflux_lat_inc(MAX_LAT_INC)
383			INTEGER swflux_nlon, swflux_nlat
384			_RL runoff_lon0, runoff_lon_inc
385			_RL runoff_lat0, runoff_lat_inc(MAX_LAT_INC)
386			INTEGER runoff_nlon, runoff_nlat
387			_RL atemp_lon0, atemp_lon_inc
388	dimitri	1.10	_RL atemp_lat0, atemp_lat_inc(MAX_LAT_INC)
389	dimitri	1.9	INTEGER atemp_nlon, atemp_nlat
390			_RL aqh_lon0, aqh_lon_inc
391	dimitri	1.10	_RL aqh_lat0, aqh_lat_inc(MAX_LAT_INC)
392	dimitri	1.9	INTEGER aqh_nlon, aqh_nlat
393			_RL evap_lon0, evap_lon_inc
394	dimitri	1.10	_RL evap_lat0, evap_lat_inc(MAX_LAT_INC)
395	dimitri	1.9	INTEGER evap_nlon, evap_nlat
396			_RL precip_lon0, precip_lon_inc
397	dimitri	1.10	_RL precip_lat0, precip_lat_inc(MAX_LAT_INC)
398	dimitri	1.9	INTEGER precip_nlon, precip_nlat
399			_RL uwind_lon0, uwind_lon_inc
400	dimitri	1.10	_RL uwind_lat0, uwind_lat_inc(MAX_LAT_INC)
401	dimitri	1.9	INTEGER uwind_nlon, uwind_nlat
402			_RL vwind_lon0, vwind_lon_inc
403	dimitri	1.10	_RL vwind_lat0, vwind_lat_inc(MAX_LAT_INC)
404	dimitri	1.9	INTEGER vwind_nlon, vwind_nlat
405			_RL lwflux_lon0, lwflux_lon_inc
406	dimitri	1.10	_RL lwflux_lat0, lwflux_lat_inc(MAX_LAT_INC)
407	dimitri	1.9	INTEGER lwflux_nlon, lwflux_nlat
408			_RL swdown_lon0, swdown_lon_inc
409	dimitri	1.10	_RL swdown_lat0, swdown_lat_inc(MAX_LAT_INC)
410	dimitri	1.9	INTEGER swdown_nlon, swdown_nlat
411			_RL lwdown_lon0, lwdown_lon_inc
412	dimitri	1.10	_RL lwdown_lat0, lwdown_lat_inc(MAX_LAT_INC)
413	dimitri	1.9	INTEGER lwdown_nlon, lwdown_nlat
414			_RL apressure_lon0,apressure_lon_inc
415	dimitri	1.10	_RL apressure_lat0,apressure_lat_inc(MAX_LAT_INC)
416	dimitri	1.9	INTEGER apressure_nlon,apressure_nlat
417
418	dimitri	1.10	common /exf_interpolation/
419	dimitri	1.9	& ustress_lon0, ustress_lon_inc,
420			& ustress_lat0, ustress_lat_inc,
421			& ustress_nlon, ustress_nlat,
422			& vstress_lon0, vstress_lon_inc,
423			& vstress_lat0, vstress_lat_inc,
424			& vstress_nlon, vstress_nlat,
425			& hflux_lon0, hflux_lon_inc,
426			& hflux_lat0, hflux_lat_inc,
427			& hflux_nlon, hflux_nlat,
428			& sflux_lon0, sflux_lon_inc,
429			& sflux_lat0, sflux_lat_inc,
430			& sflux_nlon, sflux_nlat,
431			& swflux_lon0, swflux_lon_inc,
432			& swflux_lat0, swflux_lat_inc,
433			& swflux_nlon, swflux_nlat,
434			& runoff_lon0, runoff_lon_inc,
435			& runoff_lat0, runoff_lat_inc,
436			& runoff_nlon, runoff_nlat,
437	dimitri	1.12	& atemp_lon0, atemp_lon_inc,
438	dimitri	1.9	& atemp_lat0, atemp_lat_inc,
439			& atemp_nlon, atemp_nlat,
440			& aqh_lon0, aqh_lon_inc,
441			& aqh_lat0, aqh_lat_inc,
442			& aqh_nlon, aqh_nlat,
443			& evap_lon0, evap_lon_inc,
444			& evap_lat0, evap_lat_inc,
445			& evap_nlon, evap_nlat,
446			& precip_lon0, precip_lon_inc,
447			& precip_lat0, precip_lat_inc,
448			& precip_nlon, precip_nlat,
449			& uwind_lon0, uwind_lon_inc,
450			& uwind_lat0, uwind_lat_inc,
451			& uwind_nlon, uwind_nlat,
452			& vwind_lon0, vwind_lon_inc,
453			& vwind_lat0, vwind_lat_inc,
454			& vwind_nlon, vwind_nlat,
455			& lwflux_lon0, lwflux_lon_inc,
456			& lwflux_lat0, lwflux_lat_inc,
457			& lwflux_nlon, lwflux_nlat,
458			& swdown_lon0, swdown_lon_inc,
459			& swdown_lat0, swdown_lat_inc,
460			& swdown_nlon, swdown_nlat,
461			& lwdown_lon0, lwdown_lon_inc,
462			& lwdown_lat0, lwdown_lat_inc,
463			& lwdown_nlon, lwdown_nlat,
464			& apressure_lon0,apressure_lon_inc,
465			& apressure_lat0,apressure_lat_inc,
466			& apressure_nlon,apressure_nlat
467			#endif