pkg/exf/exf_param.h

c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_param.h,v 1.28 2007/01/18 23:44:43 dimitri 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     year in seconds
      _RL     year2sec

c     Calendar data.
      _RL     repeatPeriod

c     Monitor Frequency (s)
      _RL     exf_monFreq

c     Sea-water albedo
      _RL     exf_albedo
c     longwave surface emissivities (ice and snow emissivities are used
c     in conjunction with thsice/seaice)
      _RL     ocean_emissivity
      _RL     ice_emissivity
      _RL     snow_emissivity

c     Drag coefficient scaling factor
      _RL     exf_scal_BulkCdn

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

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

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

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

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

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

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

      integer snowprecipstartdate1
      integer snowprecipstartdate2
      _RL     snowprecipstartdate
      _RL     snowprecipperiod
      _RL     snowprecipconst
      _RL     snowprecip_exfremo_intercept 
      _RL     snowprecip_exfremo_slope
      character*1 snowprecipmask
      parameter(  snowprecipmask = 's' )

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

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

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

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

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

      integer wspeedstartdate1
      integer wspeedstartdate2
      _RL     wspeedstartdate
      _RL     wspeedperiod
      _RL     wspeedconst
      _RL     wspeed_exfremo_intercept 
      _RL     wspeed_exfremo_slope
      character*1 wspeedmask
      parameter(  wspeedmask = 's' )

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

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

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

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

      integer apressurestartdate1
      integer apressurestartdate2
      _RL     apressurestartdate
      _RL     apressureperiod
      _RL     apressureconst
      _RL     apressure_exfremo_intercept 
      _RL     apressure_exfremo_slope
      character*1 apressuremask
      parameter(  apressuremask = '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

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

C     useExfYearlyFields :: when set, automatically add extension
C                           _YEAR to input file names
C     twoDigitYear       :: when set, use 2-digit year extension YR
C                           instead of _YEAR for useExfYearlyFields
      logical useExfYearlyFields, twoDigitYear
      logical useExfCheckRange

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

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

      common /exf_param_trend_removal/
     &                     hflux_exfremo_intercept,
     &                     atemp_exfremo_intercept,
     &                     aqh_exfremo_intercept,
     &                     sflux_exfremo_intercept,
     &                     evap_exfremo_intercept,
     &                     precip_exfremo_intercept,
     &                     snowprecip_exfremo_intercept,
     &                     runoff_exfremo_intercept,
     &                     ustress_exfremo_intercept,
     &                     vstress_exfremo_intercept,
     &                     uwind_exfremo_intercept,
     &                     vwind_exfremo_intercept,
     &                     wspeed_exfremo_intercept,
     &                     swflux_exfremo_intercept,
     &                     lwflux_exfremo_intercept,
     &                     swdown_exfremo_intercept,
     &                     lwdown_exfremo_intercept,
     &                     apressure_exfremo_intercept,
     &                     hflux_exfremo_slope,
     &                     atemp_exfremo_slope,
     &                     aqh_exfremo_slope,
     &                     sflux_exfremo_slope,
     &                     evap_exfremo_slope,
     &                     precip_exfremo_slope,
     &                     snowprecip_exfremo_slope,
     &                     runoff_exfremo_slope,
     &                     ustress_exfremo_slope,
     &                     vstress_exfremo_slope,
     &                     uwind_exfremo_slope,
     &                     vwind_exfremo_slope,
     &                     wspeed_exfremo_slope,
     &                     swflux_exfremo_slope,
     &                     lwflux_exfremo_slope,
     &                     swdown_exfremo_slope,
     &                     lwdown_exfremo_slope,
     &                     apressure_exfremo_slope

      common /exf_param_c/
     &                     hfluxfile,
     &                     atempfile,
     &                     aqhfile,
     &                     sfluxfile,
     &                     evapfile,
     &                     precipfile,
     &                     snowprecipfile,
     &                     runofffile,
     &                     ustressfile,
     &                     vstressfile,
     &                     uwindfile,
     &                     vwindfile,
     &                     wspeedfile,
     &                     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_wspeed
      _RL     exf_inscal_swflux
      _RL     exf_inscal_lwflux
      _RL     exf_inscal_precip
      _RL     exf_inscal_snowprecip
      _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_wspeed
     &                    , exf_inscal_swflux
     &                    , exf_inscal_lwflux
     &                    , exf_inscal_precip
     &                    , exf_inscal_snowprecip
     &                    , 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 2176 max data values
       integer MAX_LAT_INC
       parameter(MAX_LAT_INC = 2176)
      _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 snowprecip_lon0, snowprecip_lon_inc
      _RL snowprecip_lat0, snowprecip_lat_inc(MAX_LAT_INC)
      INTEGER snowprecip_nlon, snowprecip_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 wspeed_lon0, wspeed_lon_inc
      _RL wspeed_lat0, wspeed_lat_inc(MAX_LAT_INC)
      INTEGER wspeed_nlon, wspeed_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,
     & snowprecip_lon0, snowprecip_lon_inc,
     & snowprecip_lat0, snowprecip_lat_inc,
     & snowprecip_nlon, snowprecip_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,
     & wspeed_lon0, wspeed_lon_inc,
     & wspeed_lat0, wspeed_lat_inc,
     & wspeed_nlon, wspeed_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	dimitri	1.29	c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_param.h,v 1.28 2007/01/18 23:44:43 dimitri 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	heimbach	1.19	c year in seconds
34			_RL year2sec
35
36	heimbach	1.1	c Calendar data.
37	cheisey	1.4	_RL repeatPeriod
38
39	jmc	1.26	c Monitor Frequency (s)
40			_RL exf_monFreq
41
42	dimitri	1.14	c Sea-water albedo
43			_RL exf_albedo
44	mlosch	1.25	c longwave surface emissivities (ice and snow emissivities are used
45			c in conjunction with thsice/seaice)
46			_RL ocean_emissivity
47			_RL ice_emissivity
48			_RL snow_emissivity
49	dimitri	1.14
50	dimitri	1.28	c Drag coefficient scaling factor
51			_RL exf_scal_BulkCdn
52
53	mlosch	1.16	c Maximum absolute windstress, used to reset unreastically high
54			c data values
55			_RL windstressmax
56
57	heimbach	1.1	integer hfluxstartdate1
58			integer hfluxstartdate2
59	dimitri	1.11	_RL hfluxstartdate
60	heimbach	1.1	_RL hfluxperiod
61	heimbach	1.7	_RL hfluxconst
62	heimbach	1.23	_RL hflux_exfremo_intercept
63			_RL hflux_exfremo_slope
64	heimbach	1.1	character*1 hfluxmask
65			parameter( hfluxmask = 's' )
66
67			integer atempstartdate1
68			integer atempstartdate2
69	dimitri	1.11	_RL atempstartdate
70	heimbach	1.1	_RL atempperiod
71	heimbach	1.7	_RL atempconst
72	heimbach	1.23	_RL atemp_exfremo_intercept
73			_RL atemp_exfremo_slope
74	heimbach	1.1	character*1 atempmask
75			parameter( atempmask = 's' )
76
77			integer aqhstartdate1
78			integer aqhstartdate2
79	dimitri	1.11	_RL aqhstartdate
80	heimbach	1.1	_RL aqhperiod
81	heimbach	1.7	_RL aqhconst
82	heimbach	1.23	_RL aqh_exfremo_intercept
83			_RL aqh_exfremo_slope
84	heimbach	1.1	character*1 aqhmask
85			parameter( aqhmask = 's' )
86
87			integer sfluxstartdate1
88			integer sfluxstartdate2
89	dimitri	1.11	_RL sfluxstartdate
90	heimbach	1.1	_RL sfluxperiod
91	heimbach	1.7	_RL sfluxconst
92	heimbach	1.23	_RL sflux_exfremo_intercept
93			_RL sflux_exfremo_slope
94	heimbach	1.1	character*1 sfluxmask
95			parameter( sfluxmask = 's' )
96
97	dimitri	1.5	integer evapstartdate1
98			integer evapstartdate2
99	dimitri	1.11	_RL evapstartdate
100	dimitri	1.5	_RL evapperiod
101	heimbach	1.7	_RL evapconst
102	heimbach	1.23	_RL evap_exfremo_intercept
103			_RL evap_exfremo_slope
104	dimitri	1.5	character*1 evapmask
105			parameter( evapmask = 's' )
106
107	heimbach	1.1	integer precipstartdate1
108			integer precipstartdate2
109	dimitri	1.11	_RL precipstartdate
110	heimbach	1.1	_RL precipperiod
111	heimbach	1.7	_RL precipconst
112	heimbach	1.23	_RL precip_exfremo_intercept
113			_RL precip_exfremo_slope
114	heimbach	1.1	character*1 precipmask
115			parameter( precipmask = 's' )
116
117	heimbach	1.24	integer snowprecipstartdate1
118			integer snowprecipstartdate2
119			_RL snowprecipstartdate
120			_RL snowprecipperiod
121			_RL snowprecipconst
122			_RL snowprecip_exfremo_intercept
123			_RL snowprecip_exfremo_slope
124			character*1 snowprecipmask
125			parameter( snowprecipmask = 's' )
126
127	heimbach	1.3	integer runoffstartdate1
128			integer runoffstartdate2
129	dimitri	1.11	_RL runoffstartdate
130	heimbach	1.3	_RL runoffperiod
131	heimbach	1.7	_RL runoffconst
132	heimbach	1.23	_RL runoff_exfremo_intercept
133			_RL runoff_exfremo_slope
134	heimbach	1.3	character*1 runoffmask
135			parameter( runoffmask = 's' )
136
137	heimbach	1.1	integer ustressstartdate1
138			integer ustressstartdate2
139	dimitri	1.11	_RL ustressstartdate
140	heimbach	1.1	_RL ustressperiod
141	heimbach	1.7	_RL ustressconst
142	heimbach	1.23	_RL ustress_exfremo_intercept
143			_RL ustress_exfremo_slope
144	heimbach	1.1	character*1 ustressmask
145			parameter( ustressmask = 'u' )
146
147			integer vstressstartdate1
148			integer vstressstartdate2
149	dimitri	1.11	_RL vstressstartdate
150	heimbach	1.1	_RL vstressperiod
151	heimbach	1.7	_RL vstressconst
152	heimbach	1.23	_RL vstress_exfremo_intercept
153			_RL vstress_exfremo_slope
154	heimbach	1.1	character*1 vstressmask
155			parameter( vstressmask = 'v' )
156
157			integer uwindstartdate1
158			integer uwindstartdate2
159	dimitri	1.11	_RL uwindstartdate
160	heimbach	1.1	_RL uwindperiod
161	heimbach	1.7	_RL uwindconst
162	heimbach	1.23	_RL uwind_exfremo_intercept
163			_RL uwind_exfremo_slope
164	heimbach	1.1	character*1 uwindmask
165	mlosch	1.17	parameter( uwindmask = 's' )
166	heimbach	1.1
167			integer vwindstartdate1
168			integer vwindstartdate2
169	dimitri	1.11	_RL vwindstartdate
170	heimbach	1.1	_RL vwindperiod
171	heimbach	1.7	_RL vwindconst
172	heimbach	1.23	_RL vwind_exfremo_intercept
173			_RL vwind_exfremo_slope
174	heimbach	1.1	character*1 vwindmask
175	mlosch	1.17	parameter( vwindmask = 's' )
176	heimbach	1.1
177	heimbach	1.24	integer wspeedstartdate1
178			integer wspeedstartdate2
179			_RL wspeedstartdate
180			_RL wspeedperiod
181			_RL wspeedconst
182			_RL wspeed_exfremo_intercept
183			_RL wspeed_exfremo_slope
184			character*1 wspeedmask
185			parameter( wspeedmask = 's' )
186
187	heimbach	1.1	integer swfluxstartdate1
188			integer swfluxstartdate2
189	dimitri	1.11	_RL swfluxstartdate
190	heimbach	1.1	_RL swfluxperiod
191	heimbach	1.7	_RL swfluxconst
192	heimbach	1.23	_RL swflux_exfremo_intercept
193			_RL swflux_exfremo_slope
194	heimbach	1.1	character*1 swfluxmask
195			parameter( swfluxmask = 's' )
196
197			integer lwfluxstartdate1
198			integer lwfluxstartdate2
199	dimitri	1.11	_RL lwfluxstartdate
200	heimbach	1.1	_RL lwfluxperiod
201	heimbach	1.7	_RL lwfluxconst
202	heimbach	1.23	_RL lwflux_exfremo_intercept
203			_RL lwflux_exfremo_slope
204	heimbach	1.1	character*1 lwfluxmask
205			parameter( lwfluxmask = 's' )
206
207	dimitri	1.6	integer swdownstartdate1
208			integer swdownstartdate2
209	dimitri	1.11	_RL swdownstartdate
210	dimitri	1.6	_RL swdownperiod
211	heimbach	1.7	_RL swdownconst
212	heimbach	1.23	_RL swdown_exfremo_intercept
213			_RL swdown_exfremo_slope
214	dimitri	1.6	character*1 swdownmask
215			parameter( swdownmask = 's' )
216
217			integer lwdownstartdate1
218			integer lwdownstartdate2
219	dimitri	1.11	_RL lwdownstartdate
220	dimitri	1.6	_RL lwdownperiod
221	heimbach	1.7	_RL lwdownconst
222	heimbach	1.23	_RL lwdown_exfremo_intercept
223			_RL lwdown_exfremo_slope
224	dimitri	1.6	character*1 lwdownmask
225			parameter( lwdownmask = 's' )
226
227	heimbach	1.23	integer apressurestartdate1
228			integer apressurestartdate2
229			_RL apressurestartdate
230			_RL apressureperiod
231			_RL apressureconst
232			_RL apressure_exfremo_intercept
233			_RL apressure_exfremo_slope
234			character*1 apressuremask
235			parameter( apressuremask = 's' )
236
237	heimbach	1.2	integer obcsNstartdate1
238			integer obcsNstartdate2
239	dimitri	1.11	_RL obcsNstartdate
240	heimbach	1.2	_RL obcsNperiod
241
242			integer obcsSstartdate1
243			integer obcsSstartdate2
244	dimitri	1.11	_RL obcsSstartdate
245	heimbach	1.2	_RL obcsSperiod
246
247			integer obcsEstartdate1
248			integer obcsEstartdate2
249	dimitri	1.11	_RL obcsEstartdate
250	heimbach	1.2	_RL obcsEperiod
251
252			integer obcsWstartdate1
253			integer obcsWstartdate2
254	dimitri	1.11	_RL obcsWstartdate
255	heimbach	1.2	_RL obcsWperiod
256	heimbach	1.1
257			c File names.
258			character*(128) hfluxfile
259			character*(128) atempfile
260			character*(128) aqhfile
261	dimitri	1.5	character*(128) evapfile
262	heimbach	1.1	character*(128) precipfile
263	heimbach	1.24	character*(128) snowprecipfile
264	heimbach	1.1	character*(128) sfluxfile
265	heimbach	1.3	character*(128) runofffile
266	heimbach	1.1	character*(128) ustressfile
267			character*(128) vstressfile
268			character*(128) uwindfile
269			character*(128) vwindfile
270	heimbach	1.24	character*(128) wspeedfile
271	heimbach	1.1	character*(128) swfluxfile
272			character*(128) lwfluxfile
273	dimitri	1.6	character*(128) swdownfile
274			character*(128) lwdownfile
275	heimbach	1.3	character*(128) apressurefile
276	heimbach	1.1
277	dimitri	1.22	C useExfYearlyFields :: when set, automatically add extension
278			C _YEAR to input file names
279			C twoDigitYear :: when set, use 2-digit year extension YR
280			C instead of _YEAR for useExfYearlyFields
281			logical useExfYearlyFields, twoDigitYear
282	heimbach	1.18	logical useExfCheckRange
283
284			common /exf_param_l/
285	heimbach	1.24	& useExfYearlyFields, twoDigitYear,
286			& useExfCheckRange
287	heimbach	1.1	common /exf_param_i/
288	heimbach	1.24	& hfluxstartdate1, hfluxstartdate2,
289			& atempstartdate1, atempstartdate2,
290			& aqhstartdate1, aqhstartdate2,
291			& sfluxstartdate1, sfluxstartdate2,
292			& evapstartdate1, evapstartdate2,
293			& runoffstartdate1, runoffstartdate2,
294			& precipstartdate1, precipstartdate2,
295			& snowprecipstartdate1, snowprecipstartdate2,
296			& ustressstartdate1, ustressstartdate2,
297			& vstressstartdate1, vstressstartdate2,
298			& uwindstartdate1, uwindstartdate2,
299			& vwindstartdate1, vwindstartdate2,
300			& wspeedstartdate1, wspeedstartdate2,
301			& swfluxstartdate1, swfluxstartdate2,
302			& lwfluxstartdate1, lwfluxstartdate2,
303			& swdownstartdate1, swdownstartdate2,
304			& lwdownstartdate1, lwdownstartdate2,
305			& obcsNstartdate1, obcsNstartdate2,
306			& obcsSstartdate1, obcsSstartdate2,
307			& obcsEstartdate1, obcsEstartdate2,
308			& obcsWstartdate1, obcsWstartdate2,
309			& apressurestartdate1,apressurestartdate2
310	heimbach	1.1
311			common /exf_param_r/
312	heimbach	1.24	& year2sec, windstressmax,
313	jmc	1.26	& repeatPeriod, exf_monFreq,
314			& exf_albedo, ocean_emissivity,
315	mlosch	1.25	& ice_emissivity, snow_emissivity,
316	dimitri	1.28	& exf_scal_BulkCdn,
317	heimbach	1.24	& hfluxperiod, hfluxstartdate,
318			& atempperiod, atempstartdate,
319			& aqhperiod, aqhstartdate,
320			& sfluxperiod, sfluxstartdate,
321			& evapperiod, evapstartdate,
322			& precipperiod, precipstartdate,
323			& snowprecipperiod, snowprecipstartdate,
324			& runoffperiod, runoffstartdate,
325			& ustressperiod, ustressstartdate,
326			& vstressperiod, vstressstartdate,
327			& uwindperiod, uwindstartdate,
328			& vwindperiod, vwindstartdate,
329			& wspeedperiod, wspeedstartdate,
330			& swfluxperiod, swfluxstartdate,
331			& lwfluxperiod, lwfluxstartdate,
332			& swdownperiod, swdownstartdate,
333			& lwdownperiod, lwdownstartdate,
334			& obcsNperiod, obcsNstartdate,
335			& obcsSperiod, obcsSstartdate,
336			& obcsEperiod, obcsEstartdate,
337			& obcsWperiod, obcsWstartdate,
338			& apressureperiod, apressurestartdate,
339			& hfluxconst,
340			& atempconst,
341			& aqhconst,
342			& sfluxconst,
343			& evapconst,
344			& precipconst,
345			& snowprecipconst,
346			& runoffconst,
347			& ustressconst,
348			& vstressconst,
349			& uwindconst,
350			& vwindconst,
351			& wspeedconst,
352			& swfluxconst,
353			& lwfluxconst,
354			& swdownconst,
355			& lwdownconst,
356			& apressureconst
357	heimbach	1.1
358	heimbach	1.23	common /exf_param_trend_removal/
359	heimbach	1.24	& hflux_exfremo_intercept,
360			& atemp_exfremo_intercept,
361			& aqh_exfremo_intercept,
362			& sflux_exfremo_intercept,
363			& evap_exfremo_intercept,
364			& precip_exfremo_intercept,
365			& snowprecip_exfremo_intercept,
366			& runoff_exfremo_intercept,
367			& ustress_exfremo_intercept,
368			& vstress_exfremo_intercept,
369			& uwind_exfremo_intercept,
370			& vwind_exfremo_intercept,
371			& wspeed_exfremo_intercept,
372			& swflux_exfremo_intercept,
373			& lwflux_exfremo_intercept,
374			& swdown_exfremo_intercept,
375			& lwdown_exfremo_intercept,
376			& apressure_exfremo_intercept,
377			& hflux_exfremo_slope,
378			& atemp_exfremo_slope,
379			& aqh_exfremo_slope,
380			& sflux_exfremo_slope,
381			& evap_exfremo_slope,
382			& precip_exfremo_slope,
383			& snowprecip_exfremo_slope,
384			& runoff_exfremo_slope,
385			& ustress_exfremo_slope,
386			& vstress_exfremo_slope,
387			& uwind_exfremo_slope,
388			& vwind_exfremo_slope,
389			& wspeed_exfremo_slope,
390			& swflux_exfremo_slope,
391			& lwflux_exfremo_slope,
392			& swdown_exfremo_slope,
393			& lwdown_exfremo_slope,
394			& apressure_exfremo_slope
395	heimbach	1.23
396	heimbach	1.1	common /exf_param_c/
397	heimbach	1.24	& hfluxfile,
398			& atempfile,
399			& aqhfile,
400			& sfluxfile,
401			& evapfile,
402			& precipfile,
403			& snowprecipfile,
404			& runofffile,
405			& ustressfile,
406			& vstressfile,
407			& uwindfile,
408			& vwindfile,
409			& wspeedfile,
410			& swfluxfile,
411			& lwfluxfile,
412			& swdownfile,
413			& lwdownfile,
414	dimitri	1.29	& apressurefile
415	heimbach	1.1
416			c file precision and field type
417
418			common /exf_param_type/
419	heimbach	1.24	& exf_iprec,
420			& exf_yftype
421	heimbach	1.1
422			integer exf_iprec
423			character*(2) exf_yftype
424	heimbach	1.3
425	dimitri	1.13	c exf_inscal_* input scaling factors
426			c exf_offset_atemp input air temperature offset
427			c (for conversion from C to K, if needed)
428			c exf_outscale_* output scaling factors
429	heimbach	1.3
430	heimbach	1.8	_RL exf_inscal_hflux
431			_RL exf_inscal_sflux
432			_RL exf_inscal_ustress
433			_RL exf_inscal_vstress
434			_RL exf_inscal_uwind
435			_RL exf_inscal_vwind
436	heimbach	1.24	_RL exf_inscal_wspeed
437	heimbach	1.8	_RL exf_inscal_swflux
438			_RL exf_inscal_lwflux
439			_RL exf_inscal_precip
440	heimbach	1.24	_RL exf_inscal_snowprecip
441	dimitri	1.6	_RL exf_inscal_sst
442			_RL exf_inscal_sss
443	heimbach	1.8	_RL exf_inscal_atemp
444	dimitri	1.13	_RL exf_offset_atemp
445	heimbach	1.8	_RL exf_inscal_aqh
446			_RL exf_inscal_evap
447	dimitri	1.6	_RL exf_inscal_apressure
448			_RL exf_inscal_runoff
449	heimbach	1.8	_RL exf_inscal_swdown
450			_RL exf_inscal_lwdown
451
452			_RL exf_outscal_hflux
453			_RL exf_outscal_sflux
454			_RL exf_outscal_ustress
455			_RL exf_outscal_vstress
456			_RL exf_outscal_swflux
457	dimitri	1.6	_RL exf_outscal_sst
458			_RL exf_outscal_sss
459			_RL exf_outscal_apressure
460	heimbach	1.3
461			common /exf_param_scal/
462	heimbach	1.8	& exf_inscal_hflux
463			& , exf_inscal_sflux
464			& , exf_inscal_ustress
465			& , exf_inscal_vstress
466			& , exf_inscal_uwind
467			& , exf_inscal_vwind
468	heimbach	1.24	& , exf_inscal_wspeed
469	heimbach	1.8	& , exf_inscal_swflux
470			& , exf_inscal_lwflux
471			& , exf_inscal_precip
472	heimbach	1.24	& , exf_inscal_snowprecip
473	dimitri	1.6	& , exf_inscal_sst
474			& , exf_inscal_sss
475	heimbach	1.8	& , exf_inscal_atemp
476	dimitri	1.13	& , exf_offset_atemp
477	heimbach	1.8	& , exf_inscal_aqh
478			& , exf_inscal_evap
479	dimitri	1.6	& , exf_inscal_apressure
480			& , exf_inscal_runoff
481	heimbach	1.8	& , exf_inscal_swdown
482			& , exf_inscal_lwdown
483			& , exf_outscal_hflux
484			& , exf_outscal_sflux
485			& , exf_outscal_ustress
486			& , exf_outscal_vstress
487			& , exf_outscal_swflux
488	dimitri	1.6	& , exf_outscal_sst
489			& , exf_outscal_sss
490			& , exf_outscal_apressure
491
492	heimbach	1.15	#ifndef USE_EXF_INTERPOLATION
493			c-- set dummy dimension 1
494			integer MAX_LAT_INC
495			parameter(MAX_LAT_INC = 1)
496			#else
497	dimitri	1.21	c for lat interpolation, arraysize currently set to 2176 max data values
498	dimitri	1.9	integer MAX_LAT_INC
499	dimitri	1.21	parameter(MAX_LAT_INC = 2176)
500	dimitri	1.9	_RL ustress_lon0, ustress_lon_inc
501			_RL ustress_lat0, ustress_lat_inc(MAX_LAT_INC)
502			INTEGER ustress_nlon, ustress_nlat
503			_RL vstress_lon0, vstress_lon_inc
504			_RL vstress_lat0, vstress_lat_inc(MAX_LAT_INC)
505			INTEGER vstress_nlon, vstress_nlat
506			_RL hflux_lon0, hflux_lon_inc
507			_RL hflux_lat0, hflux_lat_inc(MAX_LAT_INC)
508			INTEGER hflux_nlon, hflux_nlat
509			_RL sflux_lon0, sflux_lon_inc
510			_RL sflux_lat0, sflux_lat_inc(MAX_LAT_INC)
511			INTEGER sflux_nlon, sflux_nlat
512			_RL swflux_lon0, swflux_lon_inc
513			_RL swflux_lat0, swflux_lat_inc(MAX_LAT_INC)
514			INTEGER swflux_nlon, swflux_nlat
515			_RL runoff_lon0, runoff_lon_inc
516			_RL runoff_lat0, runoff_lat_inc(MAX_LAT_INC)
517			INTEGER runoff_nlon, runoff_nlat
518			_RL atemp_lon0, atemp_lon_inc
519	dimitri	1.10	_RL atemp_lat0, atemp_lat_inc(MAX_LAT_INC)
520	dimitri	1.9	INTEGER atemp_nlon, atemp_nlat
521			_RL aqh_lon0, aqh_lon_inc
522	dimitri	1.10	_RL aqh_lat0, aqh_lat_inc(MAX_LAT_INC)
523	dimitri	1.9	INTEGER aqh_nlon, aqh_nlat
524			_RL evap_lon0, evap_lon_inc
525	dimitri	1.10	_RL evap_lat0, evap_lat_inc(MAX_LAT_INC)
526	dimitri	1.9	INTEGER evap_nlon, evap_nlat
527			_RL precip_lon0, precip_lon_inc
528	dimitri	1.10	_RL precip_lat0, precip_lat_inc(MAX_LAT_INC)
529	dimitri	1.9	INTEGER precip_nlon, precip_nlat
530	heimbach	1.24	_RL snowprecip_lon0, snowprecip_lon_inc
531			_RL snowprecip_lat0, snowprecip_lat_inc(MAX_LAT_INC)
532			INTEGER snowprecip_nlon, snowprecip_nlat
533	dimitri	1.9	_RL uwind_lon0, uwind_lon_inc
534	dimitri	1.10	_RL uwind_lat0, uwind_lat_inc(MAX_LAT_INC)
535	dimitri	1.9	INTEGER uwind_nlon, uwind_nlat
536			_RL vwind_lon0, vwind_lon_inc
537	dimitri	1.10	_RL vwind_lat0, vwind_lat_inc(MAX_LAT_INC)
538	dimitri	1.9	INTEGER vwind_nlon, vwind_nlat
539	heimbach	1.24	_RL wspeed_lon0, wspeed_lon_inc
540			_RL wspeed_lat0, wspeed_lat_inc(MAX_LAT_INC)
541			INTEGER wspeed_nlon, wspeed_nlat
542	dimitri	1.9	_RL lwflux_lon0, lwflux_lon_inc
543	dimitri	1.10	_RL lwflux_lat0, lwflux_lat_inc(MAX_LAT_INC)
544	dimitri	1.9	INTEGER lwflux_nlon, lwflux_nlat
545			_RL swdown_lon0, swdown_lon_inc
546	dimitri	1.10	_RL swdown_lat0, swdown_lat_inc(MAX_LAT_INC)
547	dimitri	1.9	INTEGER swdown_nlon, swdown_nlat
548			_RL lwdown_lon0, lwdown_lon_inc
549	dimitri	1.10	_RL lwdown_lat0, lwdown_lat_inc(MAX_LAT_INC)
550	dimitri	1.9	INTEGER lwdown_nlon, lwdown_nlat
551			_RL apressure_lon0,apressure_lon_inc
552	dimitri	1.10	_RL apressure_lat0,apressure_lat_inc(MAX_LAT_INC)
553	dimitri	1.9	INTEGER apressure_nlon,apressure_nlat
554
555	dimitri	1.10	common /exf_interpolation/
556	dimitri	1.9	& ustress_lon0, ustress_lon_inc,
557			& ustress_lat0, ustress_lat_inc,
558			& ustress_nlon, ustress_nlat,
559			& vstress_lon0, vstress_lon_inc,
560			& vstress_lat0, vstress_lat_inc,
561			& vstress_nlon, vstress_nlat,
562			& hflux_lon0, hflux_lon_inc,
563			& hflux_lat0, hflux_lat_inc,
564			& hflux_nlon, hflux_nlat,
565			& sflux_lon0, sflux_lon_inc,
566			& sflux_lat0, sflux_lat_inc,
567			& sflux_nlon, sflux_nlat,
568			& swflux_lon0, swflux_lon_inc,
569			& swflux_lat0, swflux_lat_inc,
570			& swflux_nlon, swflux_nlat,
571			& runoff_lon0, runoff_lon_inc,
572			& runoff_lat0, runoff_lat_inc,
573			& runoff_nlon, runoff_nlat,
574	dimitri	1.12	& atemp_lon0, atemp_lon_inc,
575	dimitri	1.9	& atemp_lat0, atemp_lat_inc,
576			& atemp_nlon, atemp_nlat,
577			& aqh_lon0, aqh_lon_inc,
578			& aqh_lat0, aqh_lat_inc,
579			& aqh_nlon, aqh_nlat,
580			& evap_lon0, evap_lon_inc,
581			& evap_lat0, evap_lat_inc,
582			& evap_nlon, evap_nlat,
583			& precip_lon0, precip_lon_inc,
584			& precip_lat0, precip_lat_inc,
585			& precip_nlon, precip_nlat,
586	heimbach	1.24	& snowprecip_lon0, snowprecip_lon_inc,
587			& snowprecip_lat0, snowprecip_lat_inc,
588			& snowprecip_nlon, snowprecip_nlat,
589	dimitri	1.9	& uwind_lon0, uwind_lon_inc,
590			& uwind_lat0, uwind_lat_inc,
591			& uwind_nlon, uwind_nlat,
592			& vwind_lon0, vwind_lon_inc,
593			& vwind_lat0, vwind_lat_inc,
594			& vwind_nlon, vwind_nlat,
595	heimbach	1.24	& wspeed_lon0, wspeed_lon_inc,
596			& wspeed_lat0, wspeed_lat_inc,
597			& wspeed_nlon, wspeed_nlat,
598	dimitri	1.9	& lwflux_lon0, lwflux_lon_inc,
599			& lwflux_lat0, lwflux_lat_inc,
600			& lwflux_nlon, lwflux_nlat,
601			& swdown_lon0, swdown_lon_inc,
602			& swdown_lat0, swdown_lat_inc,
603			& swdown_nlon, swdown_nlat,
604			& lwdown_lon0, lwdown_lon_inc,
605			& lwdown_lat0, lwdown_lat_inc,
606			& lwdown_nlon, lwdown_nlat,
607			& apressure_lon0,apressure_lon_inc,
608			& apressure_lat0,apressure_lat_inc,
609			& apressure_nlon,apressure_nlat
610			#endif