Set ice density to 917 kg/m3

CABLE-LSM · Jul 3, 2024 · a2e177b · a2e177b
1 parent efe332b
commit a2e177b
Show file tree

Hide file tree

Showing 2 changed files with 21 additions and 21 deletions.
diff --git a/src/offline/cable_soilparm.nml b/src/offline/cable_soilparm.nml
@@ -13,7 +13,7 @@ soilin%bch=4.2,7.1,11.4,5.15,10.4,10.4,7.12,5.83,7.1,
 soilin%clay=0.09,0.3,0.67,0.2,0.42,0.48,0.27,0.17,0.3,
 soilin%css=7*850,1920,2100,
 soilin%hyds=0.000166,0.000004,0.000001,0.000021,0.000002,0.000001,0.000006,0.0008,0.000001,
-soilin%rhosoil=1600,1600,1381,1373,1476,1521,1373,1537,910,
+soilin%rhosoil=1600,1600,1381,1373,1476,1521,1373,1537,917,
 soilin%sand=0.83,0.37,0.16,0.6,0.52,0.27,0.58,0.13,0.37,
 soilin%sfc=0.143,0.301,0.367,0.218,0.31,0.37,0.255,0.45,0.301,
 soilin%silt=0.08,0.33,0.17,0.2,0.06,0.25,0.15,0.7,0.33,

diff --git a/src/params/cable_phys_constants_mod.F90 b/src/params/cable_phys_constants_mod.F90
@@ -21,27 +21,27 @@ MODULE cable_phys_constants_mod
 
 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/s/m2)
+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
-REAL, PARAMETER :: hls = 2.8350e6         ! latent heatOFsublimation (J/kg)
-REAL, PARAMETER :: dheat  = 21.5e-6       ! molecular diffusivity for heat
-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    ! denisty of ice
+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.9257e6  ! 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 = 917.0    ! density of ice                            
 
 ! Teten coefficients
 REAL, PARAMETER :: tetena = 6.106         ! Magnus Tetans (Murray 1967)