Opened 10 years ago

Closed 10 years ago

#39 closed enhancement (fixed)

so_cond and so_capa are deleted from the code and their only usage in the usage replaced by a more reasonable assumption

Reported by: jpolcher Owned by: somebody
Priority: major Milestone: ORCHIDEE 1.9.7
Component: Anthropogenic processes Version: orchidee_1_9_6
Keywords: Cc:

Description

To simplify the externalisation and reduce the number of useless parameters an effort was made to delete from the code so_capa and so_cond.

svn diff src_parameters/ src_sechiba/

Index: src_parameters/constantes_soil.f90
===================================================================
--- src_parameters/constantes_soil.f90 (revision 1051)
+++ src_parameters/constantes_soil.f90 (working copy)
@@ -90,12 +90,6 @@

INTEGER(i_std), SAVE
nscm=nscm_fao !! Default value for nscm !-
-!- Parameters for soil thermodynamics
-!-

  • REAL(r_std),PARAMETER :: so_cond = 1.5396 !! Average Thermal Conductivity of soils @tex $(W.m{-2}.K{-1})$ @endtex
  • REAL(r_std),PARAMETER :: so_capa = 2.0514e+6 !! Average Heat capacity of soils @tex $(J.m{-3}.K{-1})$ @endtex

-
-!-

! Values taken from : PIELKE,'MESOSCALE METEOROLOGICAL MODELING',P.384
! Dry soil heat capacity was decreased and conductivity increased.
!-

Index: src_sechiba/thermosoil.f90
===================================================================
--- src_sechiba/thermosoil.f90 (revision 1051)
+++ src_sechiba/thermosoil.f90 (working copy)
@@ -914,10 +914,14 @@

INTEGER(i_std)
ier, ji, jg
REAL(r_std)
sum
+ !! so_capa and so_cond are temporary variables which contain the average values of soil conductivity
+ !! and soil conductivity and are only needed in thermosoil_var_init to set the vertical layering.
+ REAL(r_std) :: so_capa, so_cond

!_ ================================================================================================================================


!! 1. Initialization of the parameters of the vertical discretization and of the attenuation depths

+ so_capa = (so_capa_wet + so_capa_dry)/deux
+ so_cond = (so_cond_wet + so_cond_dry)/deux

cstgrnd=SQRT(one_day / pi)
lskin = SQRT(so_cond / so_capa * one_day / pi)
fz1 = 0.3_r_std * cstgrnd

Change History (1)

comment:1 Changed 10 years ago by dsolyga

  • Milestone set to ORCHIDEE 1.9.7
  • Resolution set to fixed
  • Status changed from new to closed
  • Version set to orchidee_1_9_6

Change committed in the trunk. See [1082].

Note: See TracTickets for help on using tickets.