cable_phys_constants_mod.F90 Source File


This source code is part of the Community Atmosphere Biosphere Land Exchange (CABLE) model. This work is licensed under the CSIRO Open Source Software License Agreement (variation of the BSD / MIT License).You may not use this this file except in compliance with this License. A copy of the License is available at https://trac.nci.org.au/trac/cable/wiki/license.



Files dependent on this one

sourcefile~~cable_phys_constants_mod.f90~~AfferentGraph sourcefile~cable_phys_constants_mod.f90 cable_phys_constants_mod.F90 sourcefile~cable_air.f90 cable_air.F90 sourcefile~cable_air.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cable_canopy.f90 cable_canopy.F90 sourcefile~cable_canopy.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cable_canopy.f90->sourcefile~cable_air.f90 sourcefile~cable_roughness.f90 cable_roughness.F90 sourcefile~cable_canopy.f90->sourcefile~cable_roughness.f90 sourcefile~cbl_dryleaf.f90 cbl_dryLeaf.F90 sourcefile~cable_canopy.f90->sourcefile~cbl_dryleaf.f90 sourcefile~cbl_latent_heat.f90 cbl_latent_heat.F90 sourcefile~cable_canopy.f90->sourcefile~cbl_latent_heat.f90 sourcefile~cbl_surfacewetness.f90 cbl_SurfaceWetness.F90 sourcefile~cable_canopy.f90->sourcefile~cbl_surfacewetness.f90 sourcefile~cable_carbon.f90 cable_carbon.F90 sourcefile~cable_carbon.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cable_driver_common.f90 cable_driver_common.F90 sourcefile~cable_driver_common.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cable_driver_common.f90->sourcefile~cable_carbon.f90 sourcefile~cable_input.f90 cable_input.F90 sourcefile~cable_driver_common.f90->sourcefile~cable_input.f90 sourcefile~cable_mpimaster.f90 cable_mpimaster.F90 sourcefile~cable_mpimaster.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cable_mpimaster.f90->sourcefile~cable_carbon.f90 sourcefile~cable_mpimaster.f90->sourcefile~cable_driver_common.f90 sourcefile~cbl_model_driver_offline.f90 cbl_model_driver_offline.F90 sourcefile~cable_mpimaster.f90->sourcefile~cbl_model_driver_offline.f90 sourcefile~cable_mpimaster.f90->sourcefile~cable_input.f90 sourcefile~casa_cable.f90 casa_cable.F90 sourcefile~cable_mpimaster.f90->sourcefile~casa_cable.f90 sourcefile~cable_parameters.f90 cable_parameters.F90 sourcefile~cable_parameters.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cable_roughness.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cable_serial.f90 cable_serial.F90 sourcefile~cable_serial.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cable_serial.f90->sourcefile~cable_driver_common.f90 sourcefile~cable_serial.f90->sourcefile~cbl_model_driver_offline.f90 sourcefile~cable_serial.f90->sourcefile~cable_input.f90 sourcefile~cable_serial.f90->sourcefile~casa_cable.f90 sourcefile~cbl_soilsnow_init_special.f90 cbl_soilsnow_init_special.F90 sourcefile~cable_serial.f90->sourcefile~cbl_soilsnow_init_special.f90 sourcefile~cbl_dryleaf.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cbl_latent_heat.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cbl_model_driver_offline.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cbl_model_driver_offline.f90->sourcefile~cable_air.f90 sourcefile~cbl_model_driver_offline.f90->sourcefile~cable_canopy.f90 sourcefile~cbl_model_driver_offline.f90->sourcefile~cable_carbon.f90 sourcefile~cbl_model_driver_offline.f90->sourcefile~cable_roughness.f90 sourcefile~cbl_snow_aging.f90 cbl_snow_aging.F90 sourcefile~cbl_model_driver_offline.f90->sourcefile~cbl_snow_aging.f90 sourcefile~cbl_albedo.f90 cbl_albedo.F90 sourcefile~cbl_model_driver_offline.f90->sourcefile~cbl_albedo.f90 sourcefile~cbl_soilsnow_main.f90 cbl_soilsnow_main.F90 sourcefile~cbl_model_driver_offline.f90->sourcefile~cbl_soilsnow_main.f90 sourcefile~cbl_snow_aging.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cbl_snow_albedo.f90 cbl_snow_albedo.F90 sourcefile~cbl_snow_albedo.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cbl_soilsnow_data.f90 cbl_soilsnow_data.F90 sourcefile~cbl_soilsnow_data.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cbl_surfacewetness.f90->sourcefile~cable_phys_constants_mod.f90 sourcefile~cable_input.f90->sourcefile~cable_parameters.f90 sourcefile~cable_mpiworker.f90 cable_mpiworker.F90 sourcefile~cable_mpiworker.f90->sourcefile~cable_carbon.f90 sourcefile~cable_mpiworker.f90->sourcefile~cable_driver_common.f90 sourcefile~cable_mpiworker.f90->sourcefile~cbl_model_driver_offline.f90 sourcefile~cable_mpiworker.f90->sourcefile~cable_input.f90 sourcefile~cable_mpiworker.f90->sourcefile~casa_cable.f90 sourcefile~cable_mpiworker.f90->sourcefile~cbl_soilsnow_init_special.f90 sourcefile~cable_offline_driver.f90 cable_offline_driver.F90 sourcefile~cable_offline_driver.f90->sourcefile~cable_driver_common.f90 sourcefile~cable_offline_driver.f90->sourcefile~cable_serial.f90 sourcefile~casa_cable.f90->sourcefile~cable_carbon.f90 sourcefile~casa_sumcflux.f90 casa_sumcflux.F90 sourcefile~casa_sumcflux.f90->sourcefile~cable_carbon.f90 sourcefile~casaonly_luc.f90 CASAONLY_LUC.F90 sourcefile~casaonly_luc.f90->sourcefile~cable_carbon.f90 sourcefile~casaonly_luc.f90->sourcefile~casa_cable.f90 sourcefile~cbl_albedo.f90->sourcefile~cbl_snow_albedo.f90 sourcefile~cbl_conductivity.f90 cbl_conductivity.F90 sourcefile~cbl_conductivity.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_gw.f90 cbl_GW.F90 sourcefile~cbl_gw.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_gw.f90->sourcefile~cbl_conductivity.f90 sourcefile~cbl_oldconductivity.f90 cbl_Oldconductivity.F90 sourcefile~cbl_gw.f90->sourcefile~cbl_oldconductivity.f90 sourcefile~cbl_trimb.f90 cbl_trimb.F90 sourcefile~cbl_gw.f90->sourcefile~cbl_trimb.f90 sourcefile~cbl_hyd_redistrib.f90 cbl_hyd_redistrib.F90 sourcefile~cbl_hyd_redistrib.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_oldconductivity.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_remove_trans.f90 cbl_remove_trans.F90 sourcefile~cbl_remove_trans.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_smoisturev.f90 cbl_smoisturev.F90 sourcefile~cbl_smoisturev.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_smoisturev.f90->sourcefile~cbl_trimb.f90 sourcefile~cbl_snowaccum.f90 cbl_snowAccum.F90 sourcefile~cbl_snowaccum.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_snowcheck.f90 cbl_snowCheck.F90 sourcefile~cbl_snowcheck.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_snowdensity.f90 cbl_snowDensity.F90 sourcefile~cbl_snowdensity.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_snowl_adjust.f90 cbl_snowl_adjust.F90 sourcefile~cbl_snowl_adjust.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_snowmelt.f90 cbl_snowMelt.F90 sourcefile~cbl_snowmelt.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_soilfreeze.f90 cbl_soilfreeze.F90 sourcefile~cbl_soilfreeze.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_soilsnow_init_special.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_hyd_redistrib.f90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_remove_trans.f90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_snowaccum.f90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_snowcheck.f90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_snowdensity.f90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_snowl_adjust.f90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_snowmelt.f90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_soilfreeze.f90 sourcefile~cbl_stempv.f90 cbl_stempv.F90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_stempv.f90 sourcefile~cbl_surfbv.f90 cbl_surfbv.F90 sourcefile~cbl_soilsnow_main.f90->sourcefile~cbl_surfbv.f90 sourcefile~cbl_stempv.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_stempv.f90->sourcefile~cbl_conductivity.f90 sourcefile~cbl_stempv.f90->sourcefile~cbl_oldconductivity.f90 sourcefile~cbl_stempv.f90->sourcefile~cbl_trimb.f90 sourcefile~cbl_surfbv.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_surfbv.f90->sourcefile~cbl_smoisturev.f90 sourcefile~cbl_thermal.f90 cbl_thermal.F90 sourcefile~cbl_thermal.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~cbl_thermal.f90->sourcefile~cbl_gw.f90 sourcefile~cbl_thermal.f90->sourcefile~cbl_snowaccum.f90 sourcefile~cbl_thermal.f90->sourcefile~cbl_snowcheck.f90 sourcefile~cbl_thermal.f90->sourcefile~cbl_snowdensity.f90 sourcefile~cbl_thermal.f90->sourcefile~cbl_snowl_adjust.f90 sourcefile~cbl_thermal.f90->sourcefile~cbl_snowmelt.f90 sourcefile~cbl_trimb.f90->sourcefile~cbl_soilsnow_data.f90 sourcefile~spincasacnp.f90 spincasacnp.F90 sourcefile~spincasacnp.f90->sourcefile~cable_carbon.f90 sourcefile~spincasacnp.f90->sourcefile~casa_cable.f90

Source Code

!******************************************************************************
! This source code is part of the Community Atmosphere Biosphere Land Exchange
! (CABLE) model. This work is licensed under the CSIRO Open Source Software
! License Agreement (variation of the BSD / MIT License).You may not use this
! this file except in compliance with this License. A copy of the License is
! available at https://trac.nci.org.au/trac/cable/wiki/license.
!******************************************************************************
MODULE cable_phys_constants_mod

!-----------------------------------------------------------------------------
! Description:
!   CABLE physical constants
!
! This MODULE is USEd throughout CABLE.
!
! Module specific documentation:https://trac.nci.org.au/trac/cable/wiki/TBC
! Where it fits in the model flow:https://trac.nci.org.au/trac/cable/wiki/TBC
!******************************************************************************

IMPLICIT NONE

PUBLIC

REAL, PARAMETER :: tfrz   = 273.16        ! Temp (K) corresp. to 0 C
REAL, PARAMETER :: sboltz = 5.67e-8       ! Stefan-Boltz. const (W/m2/K4)
REAL, PARAMETER :: emsoil = 1.0           ! soil emissivity
REAL, PARAMETER :: emleaf = 1.0           ! leaf emissivity
REAL, PARAMETER :: capp   = 1004.64       ! air spec. heat (J/kg/K)
REAL, PARAMETER :: hl     = 2.5014e6      ! latent heat of vaporization (J/kg)
!Below are constants used in CABLE model which are not as yet used in JAC-6.2
REAL, PARAMETER :: hlf    = 0.334e6       ! latent heat of fusion (J/kg)
REAL, PARAMETER :: hls    = 2.8350e6      ! latent heat of sublimation (J/kg)
REAL, PARAMETER :: dheat  = 21.5e-6       ! molecular diffusivity for heat (cm2/s)
REAL, PARAMETER :: grav   = 9.8086        ! gravity acceleration (m/s2)
REAL, PARAMETER :: rgas   = 8.3143        ! universal gas const (J/mol/K)
REAL, PARAMETER :: rmair  = 0.02897       ! molecular wt: dry air (kg/mol)
REAL, PARAMETER :: rmh2o  = 0.018016      ! molecular wt: water (kg/mol)
REAL, PARAMETER :: cgsnow = 2090.0        ! specific heat for snow (J/kg/K)
REAL, PARAMETER :: cs_rho_ice = 1.9341e6  ! heat capacity * density ice
REAL, PARAMETER :: cs_rho_wat = 4.218e6   ! heat capacity * density  water
REAL, PARAMETER :: csice = 2.100e3        ! specific heat for ice (J/kg/K)
REAL, PARAMETER :: cswat = 4.218e3        ! specific heat for water at 0°C (J/kg/K)
REAL, PARAMETER :: density_liq = 1000.0   ! density of liquid water
REAL, PARAMETER :: density_ice = 921.0    ! density of ice                              

! Teten coefficients
REAL, PARAMETER :: tetena = 6.106         ! Magnus Tetans (Murray 1967)
REAL, PARAMETER :: tetenb = 17.27
REAL, PARAMETER :: tetenc = 237.3
! mrd561 the parameters for sat above ice
REAL, PARAMETER :: tetena_ice = 6.1078    ! ??? refs?
REAL, PARAMETER :: tetenb_ice = 21.875
REAL, PARAMETER :: tetenc_ice = 265.5

! Aerodynamic parameters, diffusivities, water density:
REAL, PARAMETER :: vonk   = 0.40          ! von Karman constant
REAL, PARAMETER :: a33    = 1.25          ! inertial sublayer sw/us
REAL, PARAMETER :: csw    = 0.50          ! canopy sw decay (Weil theory)
REAL, PARAMETER :: ctl    = 0.40          ! Wagga wheat (RDD 1992-Challenges)
REAL, PARAMETER :: apol   = 0.70          ! Polhausen coeff: one-sided plate
REAL, PARAMETER :: prandt = 0.71          ! Prandtl number: visc/diffh
REAL, PARAMETER :: schmid = 0.60          ! Schmidt number: visc/diffw
REAL, PARAMETER :: diffwc = 1.60          ! diffw/diffc = H2O/CO2 diffusivity
REAL, PARAMETER :: rhow   = 1000.0        ! liquid water density   [kg/m3]
REAL, PARAMETER :: crd = 0.3              ! element drag coefficient
REAL, PARAMETER :: csd = 0.003            ! substrate drag coefficient

!jhan:hardwire for now. note beta2 = crd/csd
REAL, PARAMETER ::  beta2 = 0.3/0.003     ! ratio cr/cs
REAL, PARAMETER ::  ccd   = 15.0          ! constant in d/h equation
REAL, PARAMETER ::  ccw_c = 2.0           ! ccw=(zw-d)/(h-d)
REAL, PARAMETER ::  usuhm = 0.3           ! (max of us/uh)

! Turbulence parameters:
REAL, PARAMETER :: zetmul = 0.4     ! if niter=2, final zeta=zetmul*zetar(2)
                                    ! niter=4 ATM see cable_define_types.F90
REAL, PARAMETER :: zeta0  = 0.0     ! initial value of za/L
REAL, PARAMETER :: zetneg = -15.0   ! negative limit on za/L when niter>=3
REAL, PARAMETER :: zetpos = 1.0     ! positive limit on za/L when niter>=3
REAL, PARAMETER :: zdlin  = 1.0     ! height frac of d below which TL linear
REAL, PARAMETER :: umin   = 0.1     ! guarantees convergence, was 0.01

!model parameter shared across subroutines -> cable_phys_constants
REAL, PARAMETER :: snow_depth_thresh = 1.0

END MODULE cable_phys_constants_mod