clm_table_module.f90

来自「CLM集合卡曼滤波数据同化算法」· F90 代码 · 共 151 行

  MODULE clm_table_module

! ======================================================================
!      Source file: clm_table_module.f90
! Original version: Yongjiu Dai, September 15, 1999
!
! Vegetation : static, derived parameters
!              Soil albedo for different colored soils
! Note: user could replace them with the site data 
!
!=======================================================================

!    IGBP land cover types

!  0 ocean
!
!    natural vegetation
!  1 evergreen needleleaf forests
!  2 evergreen broadleaf forests
!  3 deciduous needleleaf forests
!  4 deciduous broadleaf forests
!  5 mixed forests
!  6 closed shrublands
!  7 open shrublands
!  8 woody savannas
!  9 svannas
! 10 grasslands
! 11 permanent wetlands

!    developed and mosaic lands
! 12 croplands
! 13 urban and built-up lands
! 14 cropland / natural vegetation mosaics

!    non-vegetation lands
! 15 snow and ice
! 16 barren or sparsely vegetated
! 17 water bodies
! 18 bare soil
!
!-----------------------------------------------------------------------
      IMPLICIT NONE

! the following four parameters(" vegc_table, seasf_table, lai_table, lai0_table")
! are from the BATS definition. We suggest the users please pay attention when you
! do off-line run, compare them with the observational data. if you have good 
! observational information on them, please replace them.
!------------------------------
! 1) Maximum fractional cover of vegetation [-]
!
      real, dimension(18) :: &
!     vegc_table =   (/0.80,     0.90,     0.80,     0.80,     0.80,&
!            0.80,     0.30,     0.80,     0.80,     0.80,     0.80,&
!            0.80,     0.60,     0.85,     0.00,     0.08,     0.00,	0.00/)

      vegc_table =   (/1.00,     1.00,     1.00,     1.00,     1.00,&
             1.00,     1.00,     1.00,     1.00,     1.00,     1.00,&
             1.00,     1.00,     1.00,     0.00,     1.00,     0.00,	0.00/)
!
! 2) Diff. between vegc and fractional cover at 269k [-]
!
      real, dimension(18) :: &
!     seasf_table =  (/0.10,     0.50,     0.30,     0.30,     0.20,&
!            0.30,     0.10,     0.20,     0.30,     0.10,     0.40,&
!            0.60,     0.30,     0.30,     0.00,     0.03,     0.00,	0.00/)

      seasf_table =  (/0.00,     0.00,     0.00,     0.00,     0.00,&
             0.00,     0.00,     0.00,     0.00,     0.00,     0.00,&
             0.00,     0.00,     0.00,     0.00,     0.00,     0.00,	0.00/)
!
! 3) Maximum leaf area index [-]
!
      real, dimension(18) :: &
      lai_table =    (/6.00,     6.00,     6.00,     6.00,     6.00,&
             6.00,     6.00,     6.00,     6.00,     2.00,     6.00,&
             6.00,     5.00,     6.00,     0.00,     6.00,     0.00,	0.00/)
!
! 4) Minimum leaf area index  [-]
!
      real, dimension(18) :: &
      lai0_table =   (/5.00,     5.00,     1.00,     1.00,     3.00,&
             2.00,     1.00,     2.00,     1.00,     0.50,     0.50,&
             0.50,     1.00,     2.00,     0.00,     0.50,     0.00,	0.00/)
!------------------------------
!
! 5) Stem area index [-]
!
      real, dimension(18) :: &
      sai_table =    (/2.00,     2.00,     2.00,     2.00,     2.00,&
             2.00,     2.00,     2.00,     2.00,     4.00,     2.00,&
             0.50,     2.00,     2.00,     2.00,     2.00,     2.00,	0.00/)
!
! 6) Aerodynamic roughness length [m]
!    overwritten by the data provided by xubin 2/8/00
      real, dimension(18) :: &
      z0m_table =    (/1.00,     2.20,     1.00,     0.80,     0.80,&
             0.10,     0.10,     0.70,     0.10,     0.03,     0.03,&
             0.06,     0.50,     0.06,     0.01,     0.05,     0.002,	0.05/)
!
! 7) Displacement height [m]
! overwritten by the data provided by xubin 2/8/00
      real, dimension(18) :: &
      displa_table = (/11.0,    23.00,     11.0,     13.0,     13.0,&
             0.30,     0.30,     6.50,     0.70,     0.30,     0.30,&
             0.30,     3.00,     0.30,     0.00,     0.10,     0.00,	0.00/)
!
! 8) Inverse sqrt of leaf dimension [m**-0.5]
!
      real, dimension(18) :: &
      sqrtdi_table = (/5.00,     5.00,     5.00,     5.00,     5.00,&
             5.00,     5.00,     5.00,     5.00,     5.00,     5.00,&
             5.00,     5.00,     5.00,     5.00,     5.00,     5.00,	0.00/)
!
! 9) Vegetation albvedo for wavelengths < 0.7 microns [-]
!    overwritten by the data provided by xubin 6/30/00
!
      real, dimension(18) :: &
      albvgs_table = (/0.04,     0.04,     0.05,     0.07,     0.06,&
             0.07,     0.14,     0.06,     0.07,     0.07,     0.06,&
             0.06,     0.06,     0.06,     0.80,     0.19,     0.07,	0.20/)
!
! 10) Vegetation albedo for wavelengths > 0.7 microns [-]
!     overwritten by the data provided by xubin 6/30/00
!
      real, dimension(18) :: &
      albvgl_table = (/0.20,     0.20,     0.23,     0.24,     0.24,&
             0.26,     0.32,     0.21,     0.26,     0.25,     0.18,&
             0.24,     0.22,     0.22,     0.60,     0.38,     0.20,	0.40/)
!
! 11a) Fitted numerical index of rooting distribution
!     overwritten by the data provided by xubin 2/8/00

      real, dimension(18) :: &
      roota =       (/6.706,     7.344,    7.066,    5.990,    4.453,&
             6.326,   7.718,     7.604,    8.235,    10.74,    0.000,&
             5.558,   5.558,     5.558,    0.000,    4.372,    0.000,	0.00/)     
!
! 11b) Fitted numerical index of rooting distribution
!
      real, dimension(18) :: &
      rootb =	    (/2.175,     1.303,    1.953,    1.955,    1.631,&
             1.567,   1.262,     2.300,    1.627,    2.608,    0.000,&
             2.614,   2.614,     2.614,    0.000,    0.978,    0.000,	0.00/)
!
! 12) Soil albedo for different colored soils (saturated soil, visible beam) [-]
!
      real, dimension(8)  :: &
      solour = (/0.12, 0.11, 0.10, 0.09, 0.08, 0.07, 0.06, 0.05/)

  END MODULE clm_table_module