Opened 6 years ago

Closed 5 years ago

#375 closed defect (fixed)

Error in computation of kfact_root

Reported by: nvuilsce Owned by: somebody
Priority: major Milestone: IPSLCM6.v1
Component: Physical processes Version:
Keywords: Cc:

Description

kfact_root aims at accounting for root density profile in order to compute the conductivity (https://forge.ipsl.jussieu.fr/orchidee/browser/trunk/ORCHIDEE/src_sechiba/hydrol.f90#L683).
It accounts for veget_max in order to modulate the effect depending of the size of PFT fraction. If veget_max equals 0, there is no effect on the computation of kfact_root (kfact_root =1), and the higher the value of veget_max, the lower the value of kfact_root.

However, kfact_root is a parameter related to each soiltile. In this respect, the relevant information needed here is the fraction of each PFT within the soiltile (and not the fraction of each PFT within the grid cell).
For a same PFT distribution within the soiltile, we want to get the same kfact_root value for the soiltile, irrespective of the size of the soiltile (and consequently the size of the PFT that belong to this soiltile).

Case 1 : veget_max_case1= 0.5, 0, 0, 0, 0, 0.25, 0.25, 0, 0, 0, 0, 0, 0

           soiltile_case1= 0.5, 0.5, 0

Case 2 : veget_max_case2= 0, 0, 0, 0, 0, 0.5, 0.5, 0, 0, 0, 0, 0, 0

           soiltile_case2= 0, 1, 0
   
kfact_root_case1 should equal kfact_root_case2 for soiltile 2.

In order to correct this, I propose to replace veget_max by vegetmax_soil (=veget_max/soiltile).

Change History (2)

comment:1 Changed 6 years ago by nvuilsce

ORCHIDEE trunk has been modified in revision r4431.

comment:2 Changed 5 years ago by nvuilsce

  • Resolution set to fixed
  • Status changed from new to closed
Note: See TracTickets for help on using tickets.