clmtype.f90

来自「CCSM Research Tools: Community Atmospher」· F90 代码 · 共 325 行 · 第 1/2 页

#include <misc.h>
#include <preproc.h>

module clmtype

!=========================================================================
! Module for defining 1-D (vertical) CLM variable specification.
!=========================================================================

  use precision
  use clm_varpar, only : nlevsoi, nlevsno, nlevlak, numrad, nvoc, ndst
  implicit none
  public clm1d

  type clm1d

!*************************************************************************
! time invariant patch types and grid position
!*************************************************************************

     integer  :: nstep          !time step 
     integer  :: kpatch         !patch index  
     integer  :: itypveg        !vegetation type
     integer  :: itypwat        !water type
     integer  :: itypprc        !precipitation type
     integer  :: isoicol        !color classes for soil albedos
     logical  :: lakpoi         !true = lake point
     real(r8) :: lat            !latitude  for subgrid patches (radians)
     real(r8) :: lon            !longitude for subgrid patches (radians)
     real(r8) :: latdeg         !latitude  for subgrid patches (degrees)
     real(r8) :: londeg         !longitude for subgrid patches (degrees)
     real(r8) :: dtime          !time step interval

!*************************************************************************
! time variant level information
!*************************************************************************

     integer :: snl              ! number of snow layers

!*************************************************************************
! time-invariant boundary data for each of the subgrid patches
!*************************************************************************

! level values

     real(r8) :: zi(-nlevsno+0:nlevsoi) !interface level below a "z" level (m)
     real(r8) :: dz(-nlevsno+1:nlevsoi) !layer depth (m)
     real(r8) :: z (-nlevsno+1:nlevsoi) !layer thickness (m)

! soil physical properties

     real(r8) :: bsw   (nlevsoi)       !Clapp and Hornberger "b"
     real(r8) :: watsat(nlevsoi)       !volumetric soil water at saturation (porosity)
     real(r8) :: hksat (nlevsoi)       !hydraulic conductivity at saturation (mm H2O /s)
     real(r8) :: sucsat(nlevsoi)       !minimum soil suction (mm)
     real(r8) :: csol  (nlevsoi)       !heat capacity, soil solids (J/m**3/Kelvin)
     real(r8) :: tkmg  (nlevsoi)       !thermal conductivity, soil minerals  [W/m-K]  (new)
     real(r8) :: tkdry (nlevsoi)       !thermal conductivity, dry soil       (W/m/Kelvin)
     real(r8) :: tksatu(nlevsoi)       !thermal conductivity, saturated soil [W/m-K]  (new)
     real(r8) :: rootfr(nlevsoi)       !fraction of roots in each soil layer
     real(r8) :: rootr(nlevsoi)        !effective fraction of roots in each soil layer

! leaf constants

     real(r8) :: dewmx                 !Maximum allowed dew [mm]

! hydraulic constants of soil 

     real(r8) :: wtfact                !Fraction of model area with high water table
     real(r8) :: trsmx0                !Max transpiration for moist soil+100% veg. [mm/s]

! roughness lengths

     real(r8) :: zlnd                  !Roughness length for soil [m] (new)             
     real(r8) :: zsno                  !Roughness length for snow [m] (new)             
     real(r8) :: csoilc                !Drag coefficient for soil under canopy [-] (new)

! numerical finite-difference

     real(r8) :: cnfac                 !Crank Nicholson factor (between 0 and 1) (new)
     real(r8) :: capr                  !Tuning factor to turn first layer T into surface T (new)  
     real(r8) :: ssi                   !Irreducible water saturation of snow (new)
     real(r8) :: wimp                  !Water impremeable if porosity less than wimp (new)
     real(r8) :: pondmx                !Ponding depth (mm) (new)
     real(r8) :: smpmax                !wilting point potential in mm (new)
     real(r8) :: smpmin                !restriction for min of soil potential (mm) (new)

! water and energy balance check

     real(r8) :: begwb                 !water mass begining of the time step
     real(r8) :: endwb                 !water mass end of the time step
     real(r8) :: errh2o                !water conservation error (mm H2O)
     real(r8) :: errsoi                !soil/lake energy conservation error (W/m**2)
     real(r8) :: errseb                !surface energy conservation error (W/m**2)
     real(r8) :: errsol                !solar radiation conservation error (W/m**2)
     real(r8) :: errlon                !longwave radiation conservation error (W/m**2)
     real(r8) :: acc_errh2o            !accumulation of water balance error
     real(r8) :: acc_errseb            !accumulation of surface energy balance error
     real(r8) :: acc_errsoi            !accumulation of energy balance error
     real(r8) :: acc_errsol            !accumulation of energy balance error

!*************************************************************************
! subgrid patch version of atm model input
!*************************************************************************

     real(r8) :: forc_t              !atmospheric temperature (Kelvin)
     real(r8) :: forc_u              !atmospheric wind speed in east direction (m/s)
     real(r8) :: forc_v              !atmospheric wind speed in north direction (m/s)
     real(r8) :: forc_q              !atmospheric specific humidity (kg/kg)
     real(r8) :: forc_hgt            !atmospheric reference height (m) 
     real(r8) :: forc_hgt_u          !observational height of wind [m] (new)
     real(r8) :: forc_hgt_t          !observational height of temperature [m] (new)
     real(r8) :: forc_hgt_q          !observational height of humidity [m] (new)
     real(r8) :: forc_pbot           !atmospheric pressure (Pa)
     real(r8) :: forc_th             !atmospheric potential temperature (Kelvin)
     real(r8) :: forc_vp             !atmospheric vapor pressure (Pa)
     real(r8) :: forc_rho            !density (kg/m**3)
     real(r8) :: forc_co2            !atmospheric CO2 concentration (Pa)
     real(r8) :: forc_o2             !atmospheric O2 concentration (Pa)
     real(r8) :: forc_lwrad          !downward infrared (longwave) radiation (W/m**2)
     real(r8) :: forc_psrf           !surface pressure (Pa)
     real(r8) :: forc_solad(numrad)  !direct beam radiation (vis=forc_sols , nir=forc_soll )
     real(r8) :: forc_solai(numrad)  !diffuse radiation     (vis=forc_solsd, nir=forc_solld)
     real(r8) :: forc_rain           !rain rate [mm/s]
     real(r8) :: forc_snow           !snow rate [mm/s]

!*************************************************************************
! biogeophys
!*************************************************************************

! Surface solar radiation 

     real(r8) :: rssun          !sunlit stomatal resistance (s/m)
     real(r8) :: rssha          !shaded stomatal resistance (s/m)
     real(r8) :: psnsun         !sunlit leaf photosynthesis (umol CO2 /m**2/ s) 
     real(r8) :: psnsha         !shaded leaf photosynthesis (umol CO2 /m**2/ s)
     real(r8) :: laisun         !sunlit leaf area
     real(r8) :: laisha         !shaded leaf area
     real(r8) :: ndvi           !Normalized Difference Vegetation Index
     real(r8) :: sabg           !solar radiation absorbed by ground (W/m**2)
     real(r8) :: sabv           !solar radiation absorbed by vegetation (W/m**2)
     real(r8) :: fsa            !solar radiation absorbed (total) (W/m**2)
     real(r8) :: fsr            !solar radiation reflected (W/m**2)

! Surface energy fluxes

     real(r8) :: taux           !wind stress: e-w (kg/m/s**2)
     real(r8) :: tauy           !wind stress: n-s (kg/m/s**2)
     real(r8) :: eflx_lwrad_out !emitted infrared (longwave) radiation (W/m**2) 
     real(r8) :: eflx_lwrad_net !net infrared (longwave) rad (W/m**2) [+ = to atm]
     real(r8) :: eflx_sh_tot    !total sensible heat flux (W/m**2) [+ to atm]
     real(r8) :: eflx_sh_veg    !sensible heat flux from leaves (W/m**2) [+ to atm]
     real(r8) :: eflx_sh_grnd   !sensible heat flux from ground (W/m**2) [+ to atm]
     real(r8) :: eflx_lh_tot    !total latent heat flux (W/m8*2)  [+ to atm] 
     real(r8) :: eflx_lh_vege   !veg evaporation heat flux (W/m**2) [+ to atm]
     real(r8) :: eflx_lh_vegt   !veg transpiration heat flux (W/m**2) [+ to atm]
     real(r8) :: eflx_lh_grnd   !ground evaporation heat flux (W/m**2) [+ to atm]   
     real(r8) :: eflx_soil_grnd !soil heat flux (W/m**2) [+ = into soil]
     real(r8) :: eflx_snomelt   !snow melt heat flux (W/m**2)
     real(r8) :: eflx_impsoil   !implicit evaporation for soil temperature equation (W/m**2)

! velocities