clmtc_module.f90

CLM集合卡曼滤波数据同化算法

  MODULE clmtc_module

!=======================================================================
!      Source file: clmtc_module.f90
! Original version: Yongjiu Dai, September 15, 1999
!
!=======================================================================

        IMPLICIT NONE

! -----------------------------------------------------------------
! [1] the land surface model works by gathering all the land points on a
!     [lon] x [lat] grid into a vector of [lpt] land points. this
!     is then expanded into a vector of [kpt] subgrid points, allowing
!     for up to [msub] subgrid points per land point. [ixy], [jxy], [kvec],
!     and [klnd] are indices for the mapping: [lon] x [lat] grid <->
!     [lpt] vector of land points <-> [kpt] vector of subgrid points. [wsg2g]
!     are the weights to obtain the land average from the subgrid points.
! -----------------------------------------------------------------

  integer, dimension(:), allocatable :: &
        ixy,        &! longitude index
        jxy,        &! latitude index
        ivt,        &! vegetation type
        isc,        &! color classes for soil albedos
        ktile        ! tile number for land grid

  integer, dimension(:), allocatable :: &
        klnd         ! land index for subgrid point

  integer, dimension(:,:), allocatable :: &
        kvec         ! subgrid index for land point

  real, dimension(:), allocatable :: &
        dlat,       &! latitude in radians
        dlon,       &! longitude in radians
        sand,       &! percent sand
        clay         ! percent clay

  real, dimension(:,:), allocatable :: &
        wsg2g        ! subgrid weights

! -----------------------------------------------------------------
! [2] the "big" vectors of [kpt] subgrid points are processed as [numlv]-1
!    "little" vectors of [lvec] points and one "little" vector of [mpt]
!     points. for each "little" vector: [beg] (<= [kpt]) is the first
!     point in the "big" [kpt] vector to process. [numpnt] is the number of
!     points in the "big" [kpt] vector to process
! -----------------------------------------------------------------
!*integer, dimension(:), allocatable :: &   ! (numlv)
!*      beg         &! first point in "big" vector to process
!*      numpnt       ! number of points in "big" kpt vector to process

! -----------------------------------------------------------------
! [3] time-invariant boundary data for the [kpt] subgrid points
! -----------------------------------------------------------------

! Soil parameters
  real, dimension(:,:), allocatable :: &
        csol,       &! heat capacity of soil solids [J/(m3 K)]
        porsl,      &! fraction of soil that is voids [-]
        phi0,       &! minimum soil suction [mm]
        bsw,        &! clapp and hornbereger "b" parameter [-]
        dkmg,       &! thermal conductivity of soil minerals [W/m-K]
        dksatu,     &! thermal conductivity of saturated soil [W/m-K]
        dkdry,      &! thermal conductivity for dry soil  [W/(m-K)]
        hksati       ! hydraulic conductivity at saturation [mm h2o/s]

  real, dimension(:), allocatable :: &
        albsol       ! soil albedo for different coloured soils [-]

! Vegetation static parameters
  real, dimension(:), allocatable :: &
        z0m,        &! aerodynamic roughness length [m]
        displa,     &! displacement height [m]
        sqrtdi,     &! inverse sqrt of leaf dimension [m**-0.5]
        effcon,     &! quantum efficiency of RuBP regeneration (mol CO2 / mol quanta)
        vmax25,     &! maximum carboxylation rate at 25 C at canopy top
        slti,       &! slope of low temperature inhibition function      (0.2)
        hlti,       &! 1/2 point of low temperature inhibition function  (288.16)
        shti,       &! slope of high temperature inhibition function     (0.3)
        hhti,       &! 1/2 point of high temperature inhibition function (313.16)
        trda,       &! temperature coefficient in gs-a model             (1.3)
        trdm,       &! temperature coefficient in gs-a model             (328.16)
        trop,       &! temperature coefficient in gs-a model             (298.16)
        gradm,      &! conductance-photosynthesis slope parameter
        binter,     &! conductance-photosynthesis intercep
        extkn,      &!

        albvgs,     &! veg. albedo for wavelengths < 0.7 microns [-]
        albvgl,     &! veg. albedo for wavelengths > 0.7 microns [-]
        chil         ! leaf angle distribution factor

  real, dimension(:,:,:), allocatable :: &
        ref,        &! leaf reflectance (iw=iband, il=life and dead)
        tran         ! leaf transmittance (iw=iband, il=life and dead)

  real, dimension(:,:), allocatable :: &
        rootfr       ! fraction of roots in each soil layer

! -----------------------------------------------------------------
! [4] CLM time step
! -----------------------------------------------------------------

  real  dtime        ! CLM time step [seconds]

  real  glai(18,18,12) 
  real  gsai(18,18,12)

  END MODULE clmtc_module