clmtype.f90

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

     real(r8) :: u10                       !10-m wind (m/s)
     real(r8) :: fv                        !friction velocity (m/s)
     real(r8) :: fm                        !used in u10 calculation

! Temperatures

     real(r8) :: t_veg                        !vegetation temperature (Kelvin)
     real(r8) :: t_grnd                       !ground temperature (Kelvin)
     real(r8) :: t_rad                        !radiative temperature (Kelvin)
     real(r8) :: t_ref2m                      !2 m height surface air temperature (Kelvin)
     real(r8) :: t_soisno(-nlevsno+1:nlevsoi) !soil temperature (Kelvin)
     real(r8) :: t_lake(1:nlevlak)            !lake temperature (Kelvin)
     real(r8) :: t_snow                       !vertically averaged snow temperature
     real(r8) :: dt_veg                       !change in t_veg, last iteration (Kelvin)
     real(r8) :: dt_grnd                      !change in t_grnd, last iteration (Kelvin)

! Soil properties

     real(r8) :: btran          !transpiration wetness factor (0 to 1) 

! Photosynthesis

     real(r8) :: fpsn           !photosynthesis (umol CO2 /m**2 /s)

!*************************************************************************
! hydrology
!*************************************************************************

     logical  :: do_capsnow                      !true => do snow capping  

     real(r8) :: qflx_infl                       !infiltration (mm H2O /s) 
     real(r8) :: qflx_surf                       !surface runoff (mm H2O /s) 
     real(r8) :: qflx_drain                      !sub-surface runoff (mm H2O /s) 
     real(r8) :: qflx_top_soil                   !net water input into soil from top (mm/s)
     real(r8) :: qflx_evap_soi                   !soil evaporation (mm H2O/s) (+ = to atm)
     real(r8) :: qflx_evap_veg                   !vegetation evaporation (mm H2O/s) (+ = to atm)
     real(r8) :: qflx_tran_veg                   !vegetation transpiration (mm H2O/s) (+ = to atm)
     real(r8) :: qflx_snomelt                    !snow melt (mm H2O /s)
     real(r8) :: qflx_evap_tot                   !qflx_evap_soi + qflx_evap_veg + qflx_tran_veg
     real(r8) :: qflx_prec_intr                  !interception of precipitation [mm/s]
     real(r8) :: qflx_prec_grnd                  !water onto ground including canopy runoff [kg/(m2 s)]
     real(r8) :: qflx_rain_grnd                  !rain on ground after interception (mm H2O/s) [+]
     real(r8) :: qflx_snow_grnd                  !snow on ground after interception (mm H2O/s) [+]
     real(r8) :: qflx_evap_grnd                  !ground surface evaporation rate (mm H2O/s) [+]
     real(r8) :: qflx_dew_grnd                   !ground surface dew formation (mm H2O /s) [+]
     real(r8) :: qflx_sub_snow                   !sublimation rate from snow pack (mm H2O /s) [+]
     real(r8) :: qflx_dew_snow                   !surface dew added to snow pack (mm H2O /s) [+]
     real(r8) :: qflx_snowcap                    !excess precipitation due to snow capping (mm H2O /s) [+]
     real(r8) :: qflx_qrgwl                      !qflx_surf at glaciers, wetlands, lakes
     real(r8) :: h2osno                          !snow water (mm H2O)
     real(r8) :: h2ocan                          !canopy water (mm H2O)
     real(r8) :: h2osoi_liq(-nlevsno+1:nlevsoi)  !liquid water (kg/m2) (new)
     real(r8) :: h2osoi_ice(-nlevsno+1:nlevsoi)  !ice lens (kg/m2) (new)
     real(r8) :: h2osoi_vol(nlevsoi)             !volumetric soil water (0<=h2osoi_vol<=watsat) [m3/m3]
     real(r8) :: snowdp                          !snow height (m) 
     real(r8) :: snowage                         !non dimensional snow age [-] (new)
     real(r8) :: snowice                         !average snow ice lens
     real(r8) :: snowliq                         !average snow liquid water
     real(r8) :: h2osno_old                      !snow mass for previous time step (kg/m2) (new)
     integer  :: frac_veg_nosno                  !fraction of vegetation not covered by snow (0 OR 1 now) [-] (new)
     integer  :: frac_veg_nosno_alb              !fraction of vegetation not covered by snow (0 OR 1 now) [-] (new)
     real(r8) :: frac_sno                        !fraction of ground covered by snow (0 to 1) 
     real(r8) :: frac_iceold(-nlevsno+1:nlevsoi) !fraction of ice relative to the total water (new)
     real(r8) :: rsw                             !soil water content for root zone
     real(r8) :: eff_porosity(nlevsoi)           !effective porosity = porosity - vol_ice
     real(r8) :: sfact                           !term for implicit correction to evaporation
     real(r8) :: sfactmax                        !maximim of "sfact"

     integer  :: imelt(-nlevsno+1:nlevsoi)       !flag for melting (=1), freezing (=2), Not=0 (new)        

!*************************************************************************
! surfacealbedo (for next time step)
!*************************************************************************

     real(r8) :: parsun         !average absorbed PAR for sunlit leaves (W/m**2)
     real(r8) :: parsha         !average absorbed PAR for shaded leaves (W/m**2)
     real(r8) :: albd(numrad)   !surface albedo (direct)                     
     real(r8) :: albi(numrad)   !surface albedo (diffuse)                    
     real(r8) :: albgrd(numrad) !ground albedo (direct)                      
     real(r8) :: albgri(numrad) !ground albedo (diffuse)                     
     real(r8) :: fabd(numrad)   !flux absorbed by veg per unit direct flux   
     real(r8) :: fabi(numrad)   !flux absorbed by veg per unit diffuse flux  
     real(r8) :: ftdd(numrad)   !down direct flux below veg per unit dir flx 
     real(r8) :: ftid(numrad)   !down diffuse flux below veg per unit dir flx
     real(r8) :: ftii(numrad)   !down diffuse flux below veg per unit dif flx
     real(r8) :: fsun           !sunlit fraction of canopy                   

!*************************************************************************
! ecosysdynamics
!*************************************************************************

     real(r8) :: hbot           !canopy bottom (m)
     real(r8) :: htop           !canopy top (m)
     real(r8) :: tlai           !one-sided leaf area index, no burying by snow
     real(r8) :: tsai           !one-sided stem area index, no burying by snow
     real(r8) :: elai           !one-sided leaf area index with burying by snow
     real(r8) :: esai           !one-sided stem area index with burying by snow
     real(r8) :: fwet           !fraction of canopy that is wet (0 to 1)
     real(r8) :: fdry           !fraction of foliage that is green and dry [-] (new)

!*************************************************************************
! terms from pft_varcon - to avoid indirect indexing
!*************************************************************************

     real(r8) :: z0mr           ! ratio of momentum roughness length to canopy top height [-]
     real(r8) :: z0m            ! momentum roughness length [m]
     real(r8) :: displar        ! ratio of displacement height to canopy top height [-]
     real(r8) :: displa         ! displacement height [m]
     real(r8) :: dleaf          ! leaf dimension [m]
     real(r8) :: xl             ! pft_varcon leaf/stem orientation index
     real(r8) :: rhol(numrad)   ! pft_varcon leaf reflectance  : 1=vis, 2=nir 
     real(r8) :: rhos(numrad)   ! pft_varcon stem reflectance  : 1=vis, 2=nir 
     real(r8) :: taul(numrad)   ! pft_varcon leaf transmittance: 1=vis, 2=nir 
     real(r8) :: taus(numrad)   ! pft_varcon stem transmittance: 1=vis, 2=nir 
     real(r8) :: qe25           ! quantum efficiency at 25c (umol co2 / umol photon)
     real(r8) :: vcmx25         ! maximum rate of carboxylation at 25c (umol co2/m**2/s)
     real(r8) :: mp             ! slope for conductance-to-photosynthesis relationship
     real(r8) :: c3psn          ! photosynthetic pathway: 0. = c4, 1. = c3

!*************************************************************************
! terms due to splitting the code into Biogeophys1 and Biogeophys2
!*************************************************************************

     real(r8) cgrnd  ! deriv. of soil energy flux wrt to soil temp [w/m2/k]
     real(r8) cgrndl ! deriv, of soil sensible heat flux wrt soil temp [w/m2/k]
     real(r8) cgrnds ! deriv of soil latent heat flux wrt soil temp [w/m**2/k]
     real(r8) tg     ! ground surface temperature [K]
     real(r8) tssbef(-nlevsno:nlevsoi)  ! soil/snow temperature before update
     real(r8) qg     ! ground specific humidity [kg/kg]
     real(r8) dqgdT  ! d(qg)/dT
     real(r8) emg    ! ground emissivity
     real(r8) emv    ! vegetation emissivity
     real(r8) htvp   ! latent heat of vapor of water (or sublimation) [j/kg]
     real(r8) z0mg   ! roughness length over ground, momentum [m]
     real(r8) z0hg   ! roughness length over ground, sensible heat [m]
     real(r8) z0qg   ! roughness length over ground, latent heat [m]
     real(r8) z0mv   ! roughness length over vegetation, momentum [m]
     real(r8) z0hv   ! roughness length over vegetation, sensible heat [m]
     real(r8) z0qv   ! roughness length over vegetation, latent heat [m]
     real(r8) beta   ! coefficient of convective velocity [-]
     real(r8) zii    ! convective boundary height [m]
     real(r8) thm    ! intermediate variable (forc_t+0.0098*forc_hgt_t)
     real(r8) thv    ! virtual potential temperature (kelvin)
     real(r8) ur     ! wind speed at reference height [m/s] ***DO WE NEED THIS???
     real(r8) dlrad  ! downward longwave radiation below the canopy [W/m2]
     real(r8) ulrad  ! upward longwave radiation above the canopy [W/m2]
     real(r8) qmelt  ! snow melt [mm/s]

! -----------------------------------------------------------------




  end type clm1d

  SAVE

end module clmtype