Changeset 7519 for branches/ORCHIDEE_2_2

Timestamp:

2022-03-14T13:06:55+01:00 (3 years ago)

Author:

josefine.ghattas

Message:

Correct comments and moved capa_ice as proposed by Catherine Ottle, see ticket #809

Location:

branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters

Files:

• constantes_soil_var.f90 (modified) (2 diffs)
• constantes_var.f90 (modified) (1 diff)

branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes_soil_var.f90

@@ -81 +81 @@
                                                         !! (kg/m3)
 !$OMP THREADPRIVATE(sn_dens)
-  REAL(r_std), SAVE :: sn_capa                          !! Heat capacity for snow 
+  REAL(r_std), SAVE :: sn_capa                          !! Volumetric heat capacity for snow 
                                                         !! @tex $(J.m^{-3}.K^{-1})$ @endtex
 !$OMP THREADPRIVATE(sn_capa)
-  REAL(r_std), SAVE :: water_capa = 4.18e+6             !! Water heat capacity 
+  REAL(r_std), PARAMETER :: capa_ice = 2.228*1.E3       !! Specific heat capacity of ice (J/kg/K)
+
+  REAL(r_std), SAVE :: water_capa = 4.18e+6             !! Volumetric water heat capacity 
                                                         !! @tex $(J.m^{-3}.K^{-1})$ @endtex
 !$OMP THREADPRIVATE(water_capa)
-  REAL(r_std), SAVE :: brk_capa = 2.0e+6                !! Heat capacity of generic rock
+  REAL(r_std), SAVE :: brk_capa = 2.0e+6                !! Volumetric heat capacity of generic rock
                                                         !! @tex $(J.m^{-3}.K^{-1})$ @endtex
 !$OMP THREADPRIVATE(brk_capa)
@@ -199 +201 @@
 &     0.25_r_std /)  ! oxisols                 
 
-  REAL(r_std),PARAMETER,DIMENSION(nscm_usda) :: so_capa_dry_usda = &     !! Dry soil Heat capacity of soils,J.m^{-3}.K^{-1}
+  REAL(r_std),PARAMETER,DIMENSION(nscm_usda) :: so_capa_dry_usda = &     !! Dry soil Volumetric Heat capacity of soils,J.m^{-3}.K^{-1}
  & (/ 1.47e+6_r_std, 1.41e+6_r_std, 1.34e+6_r_std, 1.27e+6_r_std, &      !! Pielke [2002, 2013]
  &    1.21e+6_r_std, 1.21e+6_r_std, 1.18e+6_r_std, 1.32e+6_r_std, &

branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes_var.f90

@@ -379 +379 @@
 
   !! Heat capacity
-  REAL(r_std), PARAMETER :: capa_ice = 2.228*1.E3       !! Heat capacity of ice (J/kg/K)
   REAL(r_std), PARAMETER :: rho_water = 1000.           !! Density of water (kg/m3)
   REAL(r_std), PARAMETER :: rho_ice = 920.              !! Density of ice (kg/m3)