pkg/exf/EXF_PARAM.h

C $Header: /u/gcmpack/MITgcm/pkg/exf/EXF_PARAM.h,v 1.4 2007/05/08 03:48:08 jmc Exp $
C $Name:  $
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     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' )

c     Calendar data.
      integer climsststartdate1
      integer climsststartdate2
      _RL     climsststartdate
      _RL     climsstperiod
      _RL     climsstconst
      _RL     climsst_exfremo_intercept 
      _RL     climsst_exfremo_slope
      character*1 climsstmask
      parameter(  climsstmask = 's' )

      integer climsssstartdate1
      integer climsssstartdate2
      _RL     climsssstartdate
      _RL     climsssperiod
      _RL     climsssconst
      _RL     climsss_exfremo_intercept 
      _RL     climsss_exfremo_slope
      character*1 climsssmask
      parameter(  climsssmask = 's' )

c     freezing temperature is the minimum temperature allowed, used
c     to reset climatological temperatures fields where they have
c     values below climtempfreeze
      _RL climtempfreeze

      integer selectStressGridPosition

      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
      character*(128) climsstfile
      character*(128) climsssfile

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
C     useStabilityFct_overIce :: over sea-ice, compute turbulent transfert
C                                coeff. function of stability (like over
C                                open ocean) rather than using fixed Coeff.
      logical useExfYearlyFields, twoDigitYear
      logical useExfCheckRange
      logical useStabilityFct_overIce

      common /exf_param_l/
     &                     useExfYearlyFields, twoDigitYear,
     &                     useExfCheckRange,
     &                     useStabilityFct_overIce
      common /exf_param_i/ selectStressGridPosition,
     &                     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_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

      common /exf_clim_i/
     &                        climsststartdate1,  climsststartdate2,
     &                        climsssstartdate1,  climsssstartdate2

      common /exf_clim_c/
     &                        climsstfile,
     &                        climsssfile

      common /exf_clim_r/
     &                        climtempfreeze,
     &                        climsstperiod,      climsststartdate,
     &                        climsssperiod,      climsssstartdate,
     &                        climsstconst,       climsssconst,
     &     climsst_exfremo_intercept, climsst_exfremo_slope,
     &     climsss_exfremo_intercept, climsss_exfremo_slope,
     &     exf_inscal_climsst, exf_inscal_climsss

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_inscal_climsst
      _RL     exf_inscal_climsss

      _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    exf_interp_bufferSize
      PARAMETER( exf_interp_bufferSize = 1 )
       integer MAX_LAT_INC
       parameter(MAX_LAT_INC = 1)
#else
C  To read input data without dynamical allocation (EXF_INTERP_USE_DYNALLOC undef),
C  buffer size currently set to 20000 (allows to read-in a 2x2 global data set)
      INTEGER    exf_interp_bufferSize
      PARAMETER( exf_interp_bufferSize = 10000 )
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

      _RL climsst_lon0, climsst_lon_inc
      _RL climsst_lat0, climsst_lat_inc(MAX_LAT_INC)
      INTEGER climsst_nlon, climsst_nlat
      _RL climsss_lon0, climsss_lon_inc
      _RL climsss_lat0, climsss_lat_inc(MAX_LAT_INC)
      INTEGER climsss_nlon, climsss_nlat
      common /exf_clim_interpolation/
     & climsst_lon0, climsst_lon_inc,
     & climsst_lat0, climsst_lat_inc,
     & climsst_nlon, climsst_nlat,
     & climsss_lon0, climsss_lon_inc,
     & climsss_lat0, climsss_lat_inc,
     & climsss_nlon, climsss_nlat

#endif
1	jmc	1.5	C $Header: /u/gcmpack/MITgcm/pkg/exf/EXF_PARAM.h,v 1.4 2007/05/08 03:48:08 jmc Exp $
2	jmc	1.1	C $Name: $
3			c
4			c
5			c ==================================================================
6			c HEADER exf_param
7			c ==================================================================
8			c
9			c o Header file for the surface flux data. Used by the external
10			c forcing package.
11			c
12			c started: Christian Eckert eckert@mit.edu 30-Jun-1999
13			c
14			c changed: Christian Eckert eckert@mit.edu 14-Jan-2000
15			c - Restructured the original version in order to have a
16			c better interface to the MITgcmUV.
17			c
18			c Christian Eckert eckert@mit.edu 12-Feb-2000
19			c - Changed some variables names (package prefix: exf_)
20			c
21			c Patrick Heimbach, heimbach@mit.edu 04-May-2000
22			c - included exf_iprec, exf_yftype to enable easy
23			c switch between 32bit/64 bit data format
24			c
25			c Patrick Heimbach, heimbach@mit.edu 01-May-2001
26			c - added obcs parameters
27			c
28			c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
29			c
30			c ==================================================================
31			c HEADER exf_param
32			c ==================================================================
33
34			c year in seconds
35			_RL year2sec
36
37			c Calendar data.
38			_RL repeatPeriod
39
40			c Monitor Frequency (s)
41			_RL exf_monFreq
42
43			c Drag coefficient scaling factor
44			_RL exf_scal_BulkCdn
45
46			c Maximum absolute windstress, used to reset unreastically high
47			c data values
48			_RL windstressmax
49
50			integer hfluxstartdate1
51			integer hfluxstartdate2
52			_RL hfluxstartdate
53			_RL hfluxperiod
54			_RL hfluxconst
55			_RL hflux_exfremo_intercept
56			_RL hflux_exfremo_slope
57			character*1 hfluxmask
58			parameter( hfluxmask = 's' )
59
60			integer atempstartdate1
61			integer atempstartdate2
62			_RL atempstartdate
63			_RL atempperiod
64			_RL atempconst
65			_RL atemp_exfremo_intercept
66			_RL atemp_exfremo_slope
67			character*1 atempmask
68			parameter( atempmask = 's' )
69
70			integer aqhstartdate1
71			integer aqhstartdate2
72			_RL aqhstartdate
73			_RL aqhperiod
74			_RL aqhconst
75			_RL aqh_exfremo_intercept
76			_RL aqh_exfremo_slope
77			character*1 aqhmask
78			parameter( aqhmask = 's' )
79
80			integer sfluxstartdate1
81			integer sfluxstartdate2
82			_RL sfluxstartdate
83			_RL sfluxperiod
84			_RL sfluxconst
85			_RL sflux_exfremo_intercept
86			_RL sflux_exfremo_slope
87			character*1 sfluxmask
88			parameter( sfluxmask = 's' )
89
90			integer evapstartdate1
91			integer evapstartdate2
92			_RL evapstartdate
93			_RL evapperiod
94			_RL evapconst
95			_RL evap_exfremo_intercept
96			_RL evap_exfremo_slope
97			character*1 evapmask
98			parameter( evapmask = 's' )
99
100			integer precipstartdate1
101			integer precipstartdate2
102			_RL precipstartdate
103			_RL precipperiod
104			_RL precipconst
105			_RL precip_exfremo_intercept
106			_RL precip_exfremo_slope
107			character*1 precipmask
108			parameter( precipmask = 's' )
109
110			integer snowprecipstartdate1
111			integer snowprecipstartdate2
112			_RL snowprecipstartdate
113			_RL snowprecipperiod
114			_RL snowprecipconst
115			_RL snowprecip_exfremo_intercept
116			_RL snowprecip_exfremo_slope
117			character*1 snowprecipmask
118			parameter( snowprecipmask = 's' )
119
120			integer runoffstartdate1
121			integer runoffstartdate2
122			_RL runoffstartdate
123			_RL runoffperiod
124			_RL runoffconst
125			_RL runoff_exfremo_intercept
126			_RL runoff_exfremo_slope
127			character*1 runoffmask
128			parameter( runoffmask = 's' )
129
130			integer ustressstartdate1
131			integer ustressstartdate2
132			_RL ustressstartdate
133			_RL ustressperiod
134			_RL ustressconst
135			_RL ustress_exfremo_intercept
136			_RL ustress_exfremo_slope
137			character*1 ustressmask
138			parameter( ustressmask = 'u' )
139
140			integer vstressstartdate1
141			integer vstressstartdate2
142			_RL vstressstartdate
143			_RL vstressperiod
144			_RL vstressconst
145			_RL vstress_exfremo_intercept
146			_RL vstress_exfremo_slope
147			character*1 vstressmask
148			parameter( vstressmask = 'v' )
149
150			integer uwindstartdate1
151			integer uwindstartdate2
152			_RL uwindstartdate
153			_RL uwindperiod
154			_RL uwindconst
155			_RL uwind_exfremo_intercept
156			_RL uwind_exfremo_slope
157			character*1 uwindmask
158			parameter( uwindmask = 's' )
159
160			integer vwindstartdate1
161			integer vwindstartdate2
162			_RL vwindstartdate
163			_RL vwindperiod
164			_RL vwindconst
165			_RL vwind_exfremo_intercept
166			_RL vwind_exfremo_slope
167			character*1 vwindmask
168			parameter( vwindmask = 's' )
169
170			integer wspeedstartdate1
171			integer wspeedstartdate2
172			_RL wspeedstartdate
173			_RL wspeedperiod
174			_RL wspeedconst
175			_RL wspeed_exfremo_intercept
176			_RL wspeed_exfremo_slope
177			character*1 wspeedmask
178			parameter( wspeedmask = 's' )
179
180			integer swfluxstartdate1
181			integer swfluxstartdate2
182			_RL swfluxstartdate
183			_RL swfluxperiod
184			_RL swfluxconst
185			_RL swflux_exfremo_intercept
186			_RL swflux_exfremo_slope
187			character*1 swfluxmask
188			parameter( swfluxmask = 's' )
189
190			integer lwfluxstartdate1
191			integer lwfluxstartdate2
192			_RL lwfluxstartdate
193			_RL lwfluxperiod
194			_RL lwfluxconst
195			_RL lwflux_exfremo_intercept
196			_RL lwflux_exfremo_slope
197			character*1 lwfluxmask
198			parameter( lwfluxmask = 's' )
199
200			integer swdownstartdate1
201			integer swdownstartdate2
202			_RL swdownstartdate
203			_RL swdownperiod
204			_RL swdownconst
205			_RL swdown_exfremo_intercept
206			_RL swdown_exfremo_slope
207			character*1 swdownmask
208			parameter( swdownmask = 's' )
209
210			integer lwdownstartdate1
211			integer lwdownstartdate2
212			_RL lwdownstartdate
213			_RL lwdownperiod
214			_RL lwdownconst
215			_RL lwdown_exfremo_intercept
216			_RL lwdown_exfremo_slope
217			character*1 lwdownmask
218			parameter( lwdownmask = 's' )
219
220			integer apressurestartdate1
221			integer apressurestartdate2
222			_RL apressurestartdate
223			_RL apressureperiod
224			_RL apressureconst
225			_RL apressure_exfremo_intercept
226			_RL apressure_exfremo_slope
227			character*1 apressuremask
228			parameter( apressuremask = 's' )
229
230	heimbach	1.2	c Calendar data.
231			integer climsststartdate1
232			integer climsststartdate2
233			_RL climsststartdate
234			_RL climsstperiod
235			_RL climsstconst
236			_RL climsst_exfremo_intercept
237			_RL climsst_exfremo_slope
238			character*1 climsstmask
239			parameter( climsstmask = 's' )
240
241			integer climsssstartdate1
242			integer climsssstartdate2
243			_RL climsssstartdate
244			_RL climsssperiod
245			_RL climsssconst
246			_RL climsss_exfremo_intercept
247			_RL climsss_exfremo_slope
248			character*1 climsssmask
249			parameter( climsssmask = 's' )
250
251			c freezing temperature is the minimum temperature allowed, used
252			c to reset climatological temperatures fields where they have
253			c values below climtempfreeze
254			_RL climtempfreeze
255
256	heimbach	1.3	integer selectStressGridPosition
257
258	jmc	1.1	integer obcsNstartdate1
259			integer obcsNstartdate2
260			_RL obcsNstartdate
261			_RL obcsNperiod
262
263			integer obcsSstartdate1
264			integer obcsSstartdate2
265			_RL obcsSstartdate
266			_RL obcsSperiod
267
268			integer obcsEstartdate1
269			integer obcsEstartdate2
270			_RL obcsEstartdate
271			_RL obcsEperiod
272
273			integer obcsWstartdate1
274			integer obcsWstartdate2
275			_RL obcsWstartdate
276			_RL obcsWperiod
277
278			c File names.
279			character*(128) hfluxfile
280			character*(128) atempfile
281			character*(128) aqhfile
282			character*(128) evapfile
283			character*(128) precipfile
284			character*(128) snowprecipfile
285			character*(128) sfluxfile
286			character*(128) runofffile
287			character*(128) ustressfile
288			character*(128) vstressfile
289			character*(128) uwindfile
290			character*(128) vwindfile
291			character*(128) wspeedfile
292			character*(128) swfluxfile
293			character*(128) lwfluxfile
294			character*(128) swdownfile
295			character*(128) lwdownfile
296			character*(128) apressurefile
297	heimbach	1.2	character*(128) climsstfile
298			character*(128) climsssfile
299	jmc	1.1
300			C useExfYearlyFields :: when set, automatically add extension
301			C _YEAR to input file names
302			C twoDigitYear :: when set, use 2-digit year extension YR
303			C instead of _YEAR for useExfYearlyFields
304	jmc	1.4	C useStabilityFct_overIce :: over sea-ice, compute turbulent transfert
305			C coeff. function of stability (like over
306			C open ocean) rather than using fixed Coeff.
307	jmc	1.1	logical useExfYearlyFields, twoDigitYear
308			logical useExfCheckRange
309	jmc	1.4	logical useStabilityFct_overIce
310	jmc	1.1
311			common /exf_param_l/
312			& useExfYearlyFields, twoDigitYear,
313	jmc	1.4	& useExfCheckRange,
314			& useStabilityFct_overIce
315	heimbach	1.3	common /exf_param_i/ selectStressGridPosition,
316	jmc	1.1	& hfluxstartdate1, hfluxstartdate2,
317			& atempstartdate1, atempstartdate2,
318			& aqhstartdate1, aqhstartdate2,
319			& sfluxstartdate1, sfluxstartdate2,
320			& evapstartdate1, evapstartdate2,
321			& runoffstartdate1, runoffstartdate2,
322			& precipstartdate1, precipstartdate2,
323			& snowprecipstartdate1, snowprecipstartdate2,
324			& ustressstartdate1, ustressstartdate2,
325			& vstressstartdate1, vstressstartdate2,
326			& uwindstartdate1, uwindstartdate2,
327			& vwindstartdate1, vwindstartdate2,
328			& wspeedstartdate1, wspeedstartdate2,
329			& swfluxstartdate1, swfluxstartdate2,
330			& lwfluxstartdate1, lwfluxstartdate2,
331			& swdownstartdate1, swdownstartdate2,
332			& lwdownstartdate1, lwdownstartdate2,
333			& obcsNstartdate1, obcsNstartdate2,
334			& obcsSstartdate1, obcsSstartdate2,
335			& obcsEstartdate1, obcsEstartdate2,
336			& obcsWstartdate1, obcsWstartdate2,
337			& apressurestartdate1,apressurestartdate2
338
339			common /exf_param_r/
340			& year2sec, windstressmax,
341			& repeatPeriod, exf_monFreq,
342			& exf_scal_BulkCdn,
343			& hfluxperiod, hfluxstartdate,
344			& atempperiod, atempstartdate,
345			& aqhperiod, aqhstartdate,
346			& sfluxperiod, sfluxstartdate,
347			& evapperiod, evapstartdate,
348			& precipperiod, precipstartdate,
349			& snowprecipperiod, snowprecipstartdate,
350			& runoffperiod, runoffstartdate,
351			& ustressperiod, ustressstartdate,
352			& vstressperiod, vstressstartdate,
353			& uwindperiod, uwindstartdate,
354			& vwindperiod, vwindstartdate,
355			& wspeedperiod, wspeedstartdate,
356			& swfluxperiod, swfluxstartdate,
357			& lwfluxperiod, lwfluxstartdate,
358			& swdownperiod, swdownstartdate,
359			& lwdownperiod, lwdownstartdate,
360			& obcsNperiod, obcsNstartdate,
361			& obcsSperiod, obcsSstartdate,
362			& obcsEperiod, obcsEstartdate,
363			& obcsWperiod, obcsWstartdate,
364			& apressureperiod, apressurestartdate,
365			& hfluxconst,
366			& atempconst,
367			& aqhconst,
368			& sfluxconst,
369			& evapconst,
370			& precipconst,
371			& snowprecipconst,
372			& runoffconst,
373			& ustressconst,
374			& vstressconst,
375			& uwindconst,
376			& vwindconst,
377			& wspeedconst,
378			& swfluxconst,
379			& lwfluxconst,
380			& swdownconst,
381			& lwdownconst,
382			& apressureconst
383
384			common /exf_param_trend_removal/
385			& hflux_exfremo_intercept,
386			& atemp_exfremo_intercept,
387			& aqh_exfremo_intercept,
388			& sflux_exfremo_intercept,
389			& evap_exfremo_intercept,
390			& precip_exfremo_intercept,
391			& snowprecip_exfremo_intercept,
392			& runoff_exfremo_intercept,
393			& ustress_exfremo_intercept,
394			& vstress_exfremo_intercept,
395			& uwind_exfremo_intercept,
396			& vwind_exfremo_intercept,
397			& wspeed_exfremo_intercept,
398			& swflux_exfremo_intercept,
399			& lwflux_exfremo_intercept,
400			& swdown_exfremo_intercept,
401			& lwdown_exfremo_intercept,
402			& apressure_exfremo_intercept,
403			& hflux_exfremo_slope,
404			& atemp_exfremo_slope,
405			& aqh_exfremo_slope,
406			& sflux_exfremo_slope,
407			& evap_exfremo_slope,
408			& precip_exfremo_slope,
409			& snowprecip_exfremo_slope,
410			& runoff_exfremo_slope,
411			& ustress_exfremo_slope,
412			& vstress_exfremo_slope,
413			& uwind_exfremo_slope,
414			& vwind_exfremo_slope,
415			& wspeed_exfremo_slope,
416			& swflux_exfremo_slope,
417			& lwflux_exfremo_slope,
418			& swdown_exfremo_slope,
419			& lwdown_exfremo_slope,
420			& apressure_exfremo_slope
421
422			common /exf_param_c/
423			& hfluxfile,
424			& atempfile,
425			& aqhfile,
426			& sfluxfile,
427			& evapfile,
428			& precipfile,
429			& snowprecipfile,
430			& runofffile,
431			& ustressfile,
432			& vstressfile,
433			& uwindfile,
434			& vwindfile,
435			& wspeedfile,
436			& swfluxfile,
437			& lwfluxfile,
438			& swdownfile,
439			& lwdownfile,
440			& apressurefile
441
442	heimbach	1.2	common /exf_clim_i/
443			& climsststartdate1, climsststartdate2,
444			& climsssstartdate1, climsssstartdate2
445
446			common /exf_clim_c/
447			& climsstfile,
448			& climsssfile
449
450			common /exf_clim_r/
451			& climtempfreeze,
452			& climsstperiod, climsststartdate,
453			& climsssperiod, climsssstartdate,
454			& climsstconst, climsssconst,
455			& climsst_exfremo_intercept, climsst_exfremo_slope,
456			& climsss_exfremo_intercept, climsss_exfremo_slope,
457			& exf_inscal_climsst, exf_inscal_climsss
458
459	jmc	1.1	c file precision and field type
460
461			common /exf_param_type/
462			& exf_iprec,
463			& exf_yftype
464
465			integer exf_iprec
466			character*(2) exf_yftype
467
468			c exf_inscal_* input scaling factors
469			c exf_offset_atemp input air temperature offset
470			c (for conversion from C to K, if needed)
471			c exf_outscale_* output scaling factors
472
473			_RL exf_inscal_hflux
474			_RL exf_inscal_sflux
475			_RL exf_inscal_ustress
476			_RL exf_inscal_vstress
477			_RL exf_inscal_uwind
478			_RL exf_inscal_vwind
479			_RL exf_inscal_wspeed
480			_RL exf_inscal_swflux
481			_RL exf_inscal_lwflux
482			_RL exf_inscal_precip
483			_RL exf_inscal_snowprecip
484			_RL exf_inscal_sst
485			_RL exf_inscal_sss
486			_RL exf_inscal_atemp
487			_RL exf_offset_atemp
488			_RL exf_inscal_aqh
489			_RL exf_inscal_evap
490			_RL exf_inscal_apressure
491			_RL exf_inscal_runoff
492			_RL exf_inscal_swdown
493			_RL exf_inscal_lwdown
494	heimbach	1.2	_RL exf_inscal_climsst
495			_RL exf_inscal_climsss
496	jmc	1.1
497			_RL exf_outscal_hflux
498			_RL exf_outscal_sflux
499			_RL exf_outscal_ustress
500			_RL exf_outscal_vstress
501			_RL exf_outscal_swflux
502			_RL exf_outscal_sst
503			_RL exf_outscal_sss
504			_RL exf_outscal_apressure
505
506			common /exf_param_scal/
507			& exf_inscal_hflux
508			& , exf_inscal_sflux
509			& , exf_inscal_ustress
510			& , exf_inscal_vstress
511			& , exf_inscal_uwind
512			& , exf_inscal_vwind
513			& , exf_inscal_wspeed
514			& , exf_inscal_swflux
515			& , exf_inscal_lwflux
516			& , exf_inscal_precip
517			& , exf_inscal_snowprecip
518			& , exf_inscal_sst
519			& , exf_inscal_sss
520			& , exf_inscal_atemp
521			& , exf_offset_atemp
522			& , exf_inscal_aqh
523			& , exf_inscal_evap
524			& , exf_inscal_apressure
525			& , exf_inscal_runoff
526			& , exf_inscal_swdown
527			& , exf_inscal_lwdown
528			& , exf_outscal_hflux
529			& , exf_outscal_sflux
530			& , exf_outscal_ustress
531			& , exf_outscal_vstress
532			& , exf_outscal_swflux
533			& , exf_outscal_sst
534			& , exf_outscal_sss
535			& , exf_outscal_apressure
536
537			#ifndef USE_EXF_INTERPOLATION
538			c-- set dummy dimension 1
539	jmc	1.5	INTEGER exf_interp_bufferSize
540			PARAMETER( exf_interp_bufferSize = 1 )
541	jmc	1.1	integer MAX_LAT_INC
542			parameter(MAX_LAT_INC = 1)
543			#else
544	jmc	1.5	C To read input data without dynamical allocation (EXF_INTERP_USE_DYNALLOC undef),
545			C buffer size currently set to 20000 (allows to read-in a 2x2 global data set)
546			INTEGER exf_interp_bufferSize
547			PARAMETER( exf_interp_bufferSize = 10000 )
548	jmc	1.1	c for lat interpolation, arraysize currently set to 2176 max data values
549			integer MAX_LAT_INC
550			parameter(MAX_LAT_INC = 2176)
551			_RL ustress_lon0, ustress_lon_inc
552			_RL ustress_lat0, ustress_lat_inc(MAX_LAT_INC)
553			INTEGER ustress_nlon, ustress_nlat
554			_RL vstress_lon0, vstress_lon_inc
555			_RL vstress_lat0, vstress_lat_inc(MAX_LAT_INC)
556			INTEGER vstress_nlon, vstress_nlat
557			_RL hflux_lon0, hflux_lon_inc
558			_RL hflux_lat0, hflux_lat_inc(MAX_LAT_INC)
559			INTEGER hflux_nlon, hflux_nlat
560			_RL sflux_lon0, sflux_lon_inc
561			_RL sflux_lat0, sflux_lat_inc(MAX_LAT_INC)
562			INTEGER sflux_nlon, sflux_nlat
563			_RL swflux_lon0, swflux_lon_inc
564			_RL swflux_lat0, swflux_lat_inc(MAX_LAT_INC)
565			INTEGER swflux_nlon, swflux_nlat
566			_RL runoff_lon0, runoff_lon_inc
567			_RL runoff_lat0, runoff_lat_inc(MAX_LAT_INC)
568			INTEGER runoff_nlon, runoff_nlat
569			_RL atemp_lon0, atemp_lon_inc
570			_RL atemp_lat0, atemp_lat_inc(MAX_LAT_INC)
571			INTEGER atemp_nlon, atemp_nlat
572			_RL aqh_lon0, aqh_lon_inc
573			_RL aqh_lat0, aqh_lat_inc(MAX_LAT_INC)
574			INTEGER aqh_nlon, aqh_nlat
575			_RL evap_lon0, evap_lon_inc
576			_RL evap_lat0, evap_lat_inc(MAX_LAT_INC)
577			INTEGER evap_nlon, evap_nlat
578			_RL precip_lon0, precip_lon_inc
579			_RL precip_lat0, precip_lat_inc(MAX_LAT_INC)
580			INTEGER precip_nlon, precip_nlat
581			_RL snowprecip_lon0, snowprecip_lon_inc
582			_RL snowprecip_lat0, snowprecip_lat_inc(MAX_LAT_INC)
583			INTEGER snowprecip_nlon, snowprecip_nlat
584			_RL uwind_lon0, uwind_lon_inc
585			_RL uwind_lat0, uwind_lat_inc(MAX_LAT_INC)
586			INTEGER uwind_nlon, uwind_nlat
587			_RL vwind_lon0, vwind_lon_inc
588			_RL vwind_lat0, vwind_lat_inc(MAX_LAT_INC)
589			INTEGER vwind_nlon, vwind_nlat
590			_RL wspeed_lon0, wspeed_lon_inc
591			_RL wspeed_lat0, wspeed_lat_inc(MAX_LAT_INC)
592			INTEGER wspeed_nlon, wspeed_nlat
593			_RL lwflux_lon0, lwflux_lon_inc
594			_RL lwflux_lat0, lwflux_lat_inc(MAX_LAT_INC)
595			INTEGER lwflux_nlon, lwflux_nlat
596			_RL swdown_lon0, swdown_lon_inc
597			_RL swdown_lat0, swdown_lat_inc(MAX_LAT_INC)
598			INTEGER swdown_nlon, swdown_nlat
599			_RL lwdown_lon0, lwdown_lon_inc
600			_RL lwdown_lat0, lwdown_lat_inc(MAX_LAT_INC)
601			INTEGER lwdown_nlon, lwdown_nlat
602			_RL apressure_lon0,apressure_lon_inc
603			_RL apressure_lat0,apressure_lat_inc(MAX_LAT_INC)
604			INTEGER apressure_nlon,apressure_nlat
605
606			common /exf_interpolation/
607			& ustress_lon0, ustress_lon_inc,
608			& ustress_lat0, ustress_lat_inc,
609			& ustress_nlon, ustress_nlat,
610			& vstress_lon0, vstress_lon_inc,
611			& vstress_lat0, vstress_lat_inc,
612			& vstress_nlon, vstress_nlat,
613			& hflux_lon0, hflux_lon_inc,
614			& hflux_lat0, hflux_lat_inc,
615			& hflux_nlon, hflux_nlat,
616			& sflux_lon0, sflux_lon_inc,
617			& sflux_lat0, sflux_lat_inc,
618			& sflux_nlon, sflux_nlat,
619			& swflux_lon0, swflux_lon_inc,
620			& swflux_lat0, swflux_lat_inc,
621			& swflux_nlon, swflux_nlat,
622			& runoff_lon0, runoff_lon_inc,
623			& runoff_lat0, runoff_lat_inc,
624			& runoff_nlon, runoff_nlat,
625			& atemp_lon0, atemp_lon_inc,
626			& atemp_lat0, atemp_lat_inc,
627			& atemp_nlon, atemp_nlat,
628			& aqh_lon0, aqh_lon_inc,
629			& aqh_lat0, aqh_lat_inc,
630			& aqh_nlon, aqh_nlat,
631			& evap_lon0, evap_lon_inc,
632			& evap_lat0, evap_lat_inc,
633			& evap_nlon, evap_nlat,
634			& precip_lon0, precip_lon_inc,
635			& precip_lat0, precip_lat_inc,
636			& precip_nlon, precip_nlat,
637			& snowprecip_lon0, snowprecip_lon_inc,
638			& snowprecip_lat0, snowprecip_lat_inc,
639			& snowprecip_nlon, snowprecip_nlat,
640			& uwind_lon0, uwind_lon_inc,
641			& uwind_lat0, uwind_lat_inc,
642			& uwind_nlon, uwind_nlat,
643			& vwind_lon0, vwind_lon_inc,
644			& vwind_lat0, vwind_lat_inc,
645			& vwind_nlon, vwind_nlat,
646			& wspeed_lon0, wspeed_lon_inc,
647			& wspeed_lat0, wspeed_lat_inc,
648			& wspeed_nlon, wspeed_nlat,
649			& lwflux_lon0, lwflux_lon_inc,
650			& lwflux_lat0, lwflux_lat_inc,
651			& lwflux_nlon, lwflux_nlat,
652			& swdown_lon0, swdown_lon_inc,
653			& swdown_lat0, swdown_lat_inc,
654			& swdown_nlon, swdown_nlat,
655			& lwdown_lon0, lwdown_lon_inc,
656			& lwdown_lat0, lwdown_lat_inc,
657			& lwdown_nlon, lwdown_nlat,
658			& apressure_lon0,apressure_lon_inc,
659			& apressure_lat0,apressure_lat_inc,
660			& apressure_nlon,apressure_nlat
661	heimbach	1.2
662			_RL climsst_lon0, climsst_lon_inc
663			_RL climsst_lat0, climsst_lat_inc(MAX_LAT_INC)
664			INTEGER climsst_nlon, climsst_nlat
665			_RL climsss_lon0, climsss_lon_inc
666			_RL climsss_lat0, climsss_lat_inc(MAX_LAT_INC)
667			INTEGER climsss_nlon, climsss_nlat
668			common /exf_clim_interpolation/
669			& climsst_lon0, climsst_lon_inc,
670			& climsst_lat0, climsst_lat_inc,
671			& climsst_nlon, climsst_nlat,
672			& climsss_lon0, climsss_lon_inc,
673			& climsss_lat0, climsss_lat_inc,
674			& climsss_nlon, climsss_nlat
675
676	jmc	1.1	#endif