pkg/exf/EXF_PARAM.h

C $Header: /u/gcmpack/MITgcm/pkg/exf/EXF_PARAM.h,v 1.5 2007/05/10 22:17:53 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 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     readStressOnAgrid  :: read wind-streess located on model-grid, A-grid position
C     readStressOnCgrid  :: read wind-streess located on model-grid, C-grid position
C     stressIsOnCgrid    :: ustress & vstress are positioned on Arakawa C-grid
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 readStressOnAgrid
      logical readStressOnCgrid
      logical stressIsOnCgrid
      logical useStabilityFct_overIce

      common /exf_param_l/
     &                     useExfYearlyFields, twoDigitYear,
     &                     useExfCheckRange,
     &                     readStressOnAgrid, readStressOnCgrid,
     &                     stressIsOnCgrid, useStabilityFct_overIce
      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_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 65000 (allows to read-in a 1x1 global data set)
      INTEGER    exf_interp_bufferSize
      PARAMETER( exf_interp_bufferSize = 65000 )
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.6	C $Header: /u/gcmpack/MITgcm/pkg/exf/EXF_PARAM.h,v 1.5 2007/05/10 22:17:53 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	jmc	1.1	integer obcsNstartdate1
257			integer obcsNstartdate2
258			_RL obcsNstartdate
259			_RL obcsNperiod
260
261			integer obcsSstartdate1
262			integer obcsSstartdate2
263			_RL obcsSstartdate
264			_RL obcsSperiod
265
266			integer obcsEstartdate1
267			integer obcsEstartdate2
268			_RL obcsEstartdate
269			_RL obcsEperiod
270
271			integer obcsWstartdate1
272			integer obcsWstartdate2
273			_RL obcsWstartdate
274			_RL obcsWperiod
275
276			c File names.
277			character*(128) hfluxfile
278			character*(128) atempfile
279			character*(128) aqhfile
280			character*(128) evapfile
281			character*(128) precipfile
282			character*(128) snowprecipfile
283			character*(128) sfluxfile
284			character*(128) runofffile
285			character*(128) ustressfile
286			character*(128) vstressfile
287			character*(128) uwindfile
288			character*(128) vwindfile
289			character*(128) wspeedfile
290			character*(128) swfluxfile
291			character*(128) lwfluxfile
292			character*(128) swdownfile
293			character*(128) lwdownfile
294			character*(128) apressurefile
295	heimbach	1.2	character*(128) climsstfile
296			character*(128) climsssfile
297	jmc	1.1
298			C useExfYearlyFields :: when set, automatically add extension
299			C _YEAR to input file names
300			C twoDigitYear :: when set, use 2-digit year extension YR
301			C instead of _YEAR for useExfYearlyFields
302	jmc	1.6	C readStressOnAgrid :: read wind-streess located on model-grid, A-grid position
303			C readStressOnCgrid :: read wind-streess located on model-grid, C-grid position
304			C stressIsOnCgrid :: ustress & vstress are positioned on Arakawa C-grid
305	jmc	1.4	C useStabilityFct_overIce :: over sea-ice, compute turbulent transfert
306			C coeff. function of stability (like over
307			C open ocean) rather than using fixed Coeff.
308	jmc	1.1	logical useExfYearlyFields, twoDigitYear
309			logical useExfCheckRange
310	jmc	1.6	logical readStressOnAgrid
311			logical readStressOnCgrid
312			logical stressIsOnCgrid
313	jmc	1.4	logical useStabilityFct_overIce
314	jmc	1.1
315			common /exf_param_l/
316			& useExfYearlyFields, twoDigitYear,
317	jmc	1.4	& useExfCheckRange,
318	jmc	1.6	& readStressOnAgrid, readStressOnCgrid,
319			& stressIsOnCgrid, useStabilityFct_overIce
320			common /exf_param_i/
321	jmc	1.1	& hfluxstartdate1, hfluxstartdate2,
322			& atempstartdate1, atempstartdate2,
323			& aqhstartdate1, aqhstartdate2,
324			& sfluxstartdate1, sfluxstartdate2,
325			& evapstartdate1, evapstartdate2,
326			& runoffstartdate1, runoffstartdate2,
327			& precipstartdate1, precipstartdate2,
328			& snowprecipstartdate1, snowprecipstartdate2,
329			& ustressstartdate1, ustressstartdate2,
330			& vstressstartdate1, vstressstartdate2,
331			& uwindstartdate1, uwindstartdate2,
332			& vwindstartdate1, vwindstartdate2,
333			& wspeedstartdate1, wspeedstartdate2,
334			& swfluxstartdate1, swfluxstartdate2,
335			& lwfluxstartdate1, lwfluxstartdate2,
336			& swdownstartdate1, swdownstartdate2,
337			& lwdownstartdate1, lwdownstartdate2,
338			& obcsNstartdate1, obcsNstartdate2,
339			& obcsSstartdate1, obcsSstartdate2,
340			& obcsEstartdate1, obcsEstartdate2,
341			& obcsWstartdate1, obcsWstartdate2,
342			& apressurestartdate1,apressurestartdate2
343
344			common /exf_param_r/
345			& year2sec, windstressmax,
346			& repeatPeriod, exf_monFreq,
347			& exf_scal_BulkCdn,
348			& hfluxperiod, hfluxstartdate,
349			& atempperiod, atempstartdate,
350			& aqhperiod, aqhstartdate,
351			& sfluxperiod, sfluxstartdate,
352			& evapperiod, evapstartdate,
353			& precipperiod, precipstartdate,
354			& snowprecipperiod, snowprecipstartdate,
355			& runoffperiod, runoffstartdate,
356			& ustressperiod, ustressstartdate,
357			& vstressperiod, vstressstartdate,
358			& uwindperiod, uwindstartdate,
359			& vwindperiod, vwindstartdate,
360			& wspeedperiod, wspeedstartdate,
361			& swfluxperiod, swfluxstartdate,
362			& lwfluxperiod, lwfluxstartdate,
363			& swdownperiod, swdownstartdate,
364			& lwdownperiod, lwdownstartdate,
365			& obcsNperiod, obcsNstartdate,
366			& obcsSperiod, obcsSstartdate,
367			& obcsEperiod, obcsEstartdate,
368			& obcsWperiod, obcsWstartdate,
369			& apressureperiod, apressurestartdate,
370			& hfluxconst,
371			& atempconst,
372			& aqhconst,
373			& sfluxconst,
374			& evapconst,
375			& precipconst,
376			& snowprecipconst,
377			& runoffconst,
378			& ustressconst,
379			& vstressconst,
380			& uwindconst,
381			& vwindconst,
382			& wspeedconst,
383			& swfluxconst,
384			& lwfluxconst,
385			& swdownconst,
386			& lwdownconst,
387			& apressureconst
388
389			common /exf_param_trend_removal/
390			& hflux_exfremo_intercept,
391			& atemp_exfremo_intercept,
392			& aqh_exfremo_intercept,
393			& sflux_exfremo_intercept,
394			& evap_exfremo_intercept,
395			& precip_exfremo_intercept,
396			& snowprecip_exfremo_intercept,
397			& runoff_exfremo_intercept,
398			& ustress_exfremo_intercept,
399			& vstress_exfremo_intercept,
400			& uwind_exfremo_intercept,
401			& vwind_exfremo_intercept,
402			& wspeed_exfremo_intercept,
403			& swflux_exfremo_intercept,
404			& lwflux_exfremo_intercept,
405			& swdown_exfremo_intercept,
406			& lwdown_exfremo_intercept,
407			& apressure_exfremo_intercept,
408			& hflux_exfremo_slope,
409			& atemp_exfremo_slope,
410			& aqh_exfremo_slope,
411			& sflux_exfremo_slope,
412			& evap_exfremo_slope,
413			& precip_exfremo_slope,
414			& snowprecip_exfremo_slope,
415			& runoff_exfremo_slope,
416			& ustress_exfremo_slope,
417			& vstress_exfremo_slope,
418			& uwind_exfremo_slope,
419			& vwind_exfremo_slope,
420			& wspeed_exfremo_slope,
421			& swflux_exfremo_slope,
422			& lwflux_exfremo_slope,
423			& swdown_exfremo_slope,
424			& lwdown_exfremo_slope,
425			& apressure_exfremo_slope
426
427			common /exf_param_c/
428			& hfluxfile,
429			& atempfile,
430			& aqhfile,
431			& sfluxfile,
432			& evapfile,
433			& precipfile,
434			& snowprecipfile,
435			& runofffile,
436			& ustressfile,
437			& vstressfile,
438			& uwindfile,
439			& vwindfile,
440			& wspeedfile,
441			& swfluxfile,
442			& lwfluxfile,
443			& swdownfile,
444			& lwdownfile,
445			& apressurefile
446
447	heimbach	1.2	common /exf_clim_i/
448			& climsststartdate1, climsststartdate2,
449			& climsssstartdate1, climsssstartdate2
450
451			common /exf_clim_c/
452			& climsstfile,
453			& climsssfile
454
455			common /exf_clim_r/
456			& climtempfreeze,
457			& climsstperiod, climsststartdate,
458			& climsssperiod, climsssstartdate,
459			& climsstconst, climsssconst,
460			& climsst_exfremo_intercept, climsst_exfremo_slope,
461			& climsss_exfremo_intercept, climsss_exfremo_slope,
462			& exf_inscal_climsst, exf_inscal_climsss
463
464	jmc	1.1	c file precision and field type
465
466	jmc	1.6	common /exf_param_type/
467	jmc	1.1	& exf_iprec,
468			& exf_yftype
469
470			integer exf_iprec
471			character*(2) exf_yftype
472
473			c exf_inscal_* input scaling factors
474			c exf_offset_atemp input air temperature offset
475			c (for conversion from C to K, if needed)
476			c exf_outscale_* output scaling factors
477
478			_RL exf_inscal_hflux
479			_RL exf_inscal_sflux
480			_RL exf_inscal_ustress
481			_RL exf_inscal_vstress
482			_RL exf_inscal_uwind
483			_RL exf_inscal_vwind
484			_RL exf_inscal_wspeed
485			_RL exf_inscal_swflux
486			_RL exf_inscal_lwflux
487			_RL exf_inscal_precip
488			_RL exf_inscal_snowprecip
489			_RL exf_inscal_sst
490			_RL exf_inscal_sss
491			_RL exf_inscal_atemp
492			_RL exf_offset_atemp
493			_RL exf_inscal_aqh
494			_RL exf_inscal_evap
495			_RL exf_inscal_apressure
496			_RL exf_inscal_runoff
497			_RL exf_inscal_swdown
498			_RL exf_inscal_lwdown
499	heimbach	1.2	_RL exf_inscal_climsst
500			_RL exf_inscal_climsss
501	jmc	1.1
502			_RL exf_outscal_hflux
503			_RL exf_outscal_sflux
504			_RL exf_outscal_ustress
505			_RL exf_outscal_vstress
506			_RL exf_outscal_swflux
507			_RL exf_outscal_sst
508			_RL exf_outscal_sss
509			_RL exf_outscal_apressure
510
511			common /exf_param_scal/
512			& exf_inscal_hflux
513			& , exf_inscal_sflux
514			& , exf_inscal_ustress
515			& , exf_inscal_vstress
516			& , exf_inscal_uwind
517			& , exf_inscal_vwind
518			& , exf_inscal_wspeed
519			& , exf_inscal_swflux
520			& , exf_inscal_lwflux
521			& , exf_inscal_precip
522			& , exf_inscal_snowprecip
523			& , exf_inscal_sst
524			& , exf_inscal_sss
525			& , exf_inscal_atemp
526			& , exf_offset_atemp
527			& , exf_inscal_aqh
528			& , exf_inscal_evap
529			& , exf_inscal_apressure
530			& , exf_inscal_runoff
531			& , exf_inscal_swdown
532			& , exf_inscal_lwdown
533			& , exf_outscal_hflux
534			& , exf_outscal_sflux
535			& , exf_outscal_ustress
536			& , exf_outscal_vstress
537			& , exf_outscal_swflux
538			& , exf_outscal_sst
539			& , exf_outscal_sss
540			& , exf_outscal_apressure
541
542			#ifndef USE_EXF_INTERPOLATION
543			c-- set dummy dimension 1
544	jmc	1.5	INTEGER exf_interp_bufferSize
545			PARAMETER( exf_interp_bufferSize = 1 )
546	jmc	1.1	integer MAX_LAT_INC
547			parameter(MAX_LAT_INC = 1)
548			#else
549	jmc	1.5	C To read input data without dynamical allocation (EXF_INTERP_USE_DYNALLOC undef),
550	jmc	1.6	C buffer size currently set to 65000 (allows to read-in a 1x1 global data set)
551	jmc	1.5	INTEGER exf_interp_bufferSize
552	jmc	1.6	PARAMETER( exf_interp_bufferSize = 65000 )
553	jmc	1.1	c for lat interpolation, arraysize currently set to 2176 max data values
554			integer MAX_LAT_INC
555			parameter(MAX_LAT_INC = 2176)
556			_RL ustress_lon0, ustress_lon_inc
557			_RL ustress_lat0, ustress_lat_inc(MAX_LAT_INC)
558			INTEGER ustress_nlon, ustress_nlat
559			_RL vstress_lon0, vstress_lon_inc
560			_RL vstress_lat0, vstress_lat_inc(MAX_LAT_INC)
561			INTEGER vstress_nlon, vstress_nlat
562			_RL hflux_lon0, hflux_lon_inc
563			_RL hflux_lat0, hflux_lat_inc(MAX_LAT_INC)
564			INTEGER hflux_nlon, hflux_nlat
565			_RL sflux_lon0, sflux_lon_inc
566			_RL sflux_lat0, sflux_lat_inc(MAX_LAT_INC)
567			INTEGER sflux_nlon, sflux_nlat
568			_RL swflux_lon0, swflux_lon_inc
569			_RL swflux_lat0, swflux_lat_inc(MAX_LAT_INC)
570			INTEGER swflux_nlon, swflux_nlat
571			_RL runoff_lon0, runoff_lon_inc
572			_RL runoff_lat0, runoff_lat_inc(MAX_LAT_INC)
573			INTEGER runoff_nlon, runoff_nlat
574			_RL atemp_lon0, atemp_lon_inc
575			_RL atemp_lat0, atemp_lat_inc(MAX_LAT_INC)
576			INTEGER atemp_nlon, atemp_nlat
577			_RL aqh_lon0, aqh_lon_inc
578			_RL aqh_lat0, aqh_lat_inc(MAX_LAT_INC)
579			INTEGER aqh_nlon, aqh_nlat
580			_RL evap_lon0, evap_lon_inc
581			_RL evap_lat0, evap_lat_inc(MAX_LAT_INC)
582			INTEGER evap_nlon, evap_nlat
583			_RL precip_lon0, precip_lon_inc
584			_RL precip_lat0, precip_lat_inc(MAX_LAT_INC)
585			INTEGER precip_nlon, precip_nlat
586			_RL snowprecip_lon0, snowprecip_lon_inc
587			_RL snowprecip_lat0, snowprecip_lat_inc(MAX_LAT_INC)
588			INTEGER snowprecip_nlon, snowprecip_nlat
589			_RL uwind_lon0, uwind_lon_inc
590			_RL uwind_lat0, uwind_lat_inc(MAX_LAT_INC)
591			INTEGER uwind_nlon, uwind_nlat
592			_RL vwind_lon0, vwind_lon_inc
593			_RL vwind_lat0, vwind_lat_inc(MAX_LAT_INC)
594			INTEGER vwind_nlon, vwind_nlat
595			_RL wspeed_lon0, wspeed_lon_inc
596			_RL wspeed_lat0, wspeed_lat_inc(MAX_LAT_INC)
597			INTEGER wspeed_nlon, wspeed_nlat
598			_RL lwflux_lon0, lwflux_lon_inc
599			_RL lwflux_lat0, lwflux_lat_inc(MAX_LAT_INC)
600			INTEGER lwflux_nlon, lwflux_nlat
601			_RL swdown_lon0, swdown_lon_inc
602			_RL swdown_lat0, swdown_lat_inc(MAX_LAT_INC)
603			INTEGER swdown_nlon, swdown_nlat
604			_RL lwdown_lon0, lwdown_lon_inc
605			_RL lwdown_lat0, lwdown_lat_inc(MAX_LAT_INC)
606			INTEGER lwdown_nlon, lwdown_nlat
607			_RL apressure_lon0,apressure_lon_inc
608			_RL apressure_lat0,apressure_lat_inc(MAX_LAT_INC)
609			INTEGER apressure_nlon,apressure_nlat
610
611			common /exf_interpolation/
612			& ustress_lon0, ustress_lon_inc,
613			& ustress_lat0, ustress_lat_inc,
614			& ustress_nlon, ustress_nlat,
615			& vstress_lon0, vstress_lon_inc,
616			& vstress_lat0, vstress_lat_inc,
617			& vstress_nlon, vstress_nlat,
618			& hflux_lon0, hflux_lon_inc,
619			& hflux_lat0, hflux_lat_inc,
620			& hflux_nlon, hflux_nlat,
621			& sflux_lon0, sflux_lon_inc,
622			& sflux_lat0, sflux_lat_inc,
623			& sflux_nlon, sflux_nlat,
624			& swflux_lon0, swflux_lon_inc,
625			& swflux_lat0, swflux_lat_inc,
626			& swflux_nlon, swflux_nlat,
627			& runoff_lon0, runoff_lon_inc,
628			& runoff_lat0, runoff_lat_inc,
629			& runoff_nlon, runoff_nlat,
630			& atemp_lon0, atemp_lon_inc,
631			& atemp_lat0, atemp_lat_inc,
632			& atemp_nlon, atemp_nlat,
633			& aqh_lon0, aqh_lon_inc,
634			& aqh_lat0, aqh_lat_inc,
635			& aqh_nlon, aqh_nlat,
636			& evap_lon0, evap_lon_inc,
637			& evap_lat0, evap_lat_inc,
638			& evap_nlon, evap_nlat,
639			& precip_lon0, precip_lon_inc,
640			& precip_lat0, precip_lat_inc,
641			& precip_nlon, precip_nlat,
642			& snowprecip_lon0, snowprecip_lon_inc,
643			& snowprecip_lat0, snowprecip_lat_inc,
644			& snowprecip_nlon, snowprecip_nlat,
645			& uwind_lon0, uwind_lon_inc,
646			& uwind_lat0, uwind_lat_inc,
647			& uwind_nlon, uwind_nlat,
648			& vwind_lon0, vwind_lon_inc,
649			& vwind_lat0, vwind_lat_inc,
650			& vwind_nlon, vwind_nlat,
651			& wspeed_lon0, wspeed_lon_inc,
652			& wspeed_lat0, wspeed_lat_inc,
653			& wspeed_nlon, wspeed_nlat,
654			& lwflux_lon0, lwflux_lon_inc,
655			& lwflux_lat0, lwflux_lat_inc,
656			& lwflux_nlon, lwflux_nlat,
657			& swdown_lon0, swdown_lon_inc,
658			& swdown_lat0, swdown_lat_inc,
659			& swdown_nlon, swdown_nlat,
660			& lwdown_lon0, lwdown_lon_inc,
661			& lwdown_lat0, lwdown_lat_inc,
662			& lwdown_nlon, lwdown_nlat,
663			& apressure_lon0,apressure_lon_inc,
664			& apressure_lat0,apressure_lat_inc,
665			& apressure_nlon,apressure_nlat
666	heimbach	1.2
667			_RL climsst_lon0, climsst_lon_inc
668			_RL climsst_lat0, climsst_lat_inc(MAX_LAT_INC)
669			INTEGER climsst_nlon, climsst_nlat
670			_RL climsss_lon0, climsss_lon_inc
671			_RL climsss_lat0, climsss_lat_inc(MAX_LAT_INC)
672			INTEGER climsss_nlon, climsss_nlat
673			common /exf_clim_interpolation/
674			& climsst_lon0, climsst_lon_inc,
675			& climsst_lat0, climsst_lat_inc,
676			& climsst_nlon, climsst_nlat,
677			& climsss_lon0, climsss_lon_inc,
678			& climsss_lat0, climsss_lat_inc,
679			& climsss_nlon, climsss_nlat
680
681	jmc	1.1	#endif