Opened 7 years ago

Closed 6 years ago

#226 closed defect (fixed)

Inconsistencies between biomass/lai/leaf_frac/leaf_age/vcmax/assim_param

Reported by: maignan Owned by: somebody
Priority: major Milestone: IPSLCM6.v1
Component: Biogeochemical processes Version: trunc
Keywords: Cc:

Description

We have encountered problems in diffuco_trans_co2 with this part of the code:

                   IF ( (A_1 .EQ. 9999.) .OR. ( A_1 .LT. (-Rd(icinic)) ) ) THEN
                      WRITE(*,*) 'We have a problem in diffuco_trans_co2'
                      WRITE(*,*) 'We have a problem in diffuco_trans_co2'
                      WRITE(*,*) 'no real positive solution found for pft:',jv
                      WRITE(*,*) 't2m:',t2m(icinic)
                      WRITE(*,*) 'vpd:',vpd(icinic)
                      A_1 = -Rd(icinic)
                   ENDIF
                   assimi(icinic) = A_1
                   ! Eq. 18 of Yin et al. (2009)
                   Cc = ( gamma_star(icinic) * x1 + ( assimi(icinic) + Rd(icinic) ) * x2 )  &
                        / ( x1 - ( assimi(icinic) + Rd(icinic) ) )

when some times x1 - ( assimi(icinic) + Rd(icinic) ) is zero.

We have traced one such pixel causing the problem:

"scheme general"

1 - sechiba_main
1.1 - diffuco_trans_co2
2 - stomate_main
2.1 - stomatelpj


-------------------------------------------------------------------
toutes les variables à  0
-------------------------------------------------------------------
dans le dt 48

sechiba
assim_param = 0


stomatelpj
prescribe: variables output (PFT modifiés)
- leaf_frac (2,5,10,12)
- biomass (12)

vmax: variables output (PFT modifiés)
-vcmax (2, 5 ,10,12)
-leaf_frac (2, 5 ,10,12)
-leaf_age (2, 5 ,10,12)

turn: variables output (PFT modifiés)
-leaf_frac ( tout à 0 )
 
establish: variables output (PFT modifiés)
-biomass (3, 5 ,10)
-leaf_frac (3, 5, 10)
-leaf_age (2, 5 10, 12) -> age 2

-------------------------------------------------------------------
dans le dt = 49

sechiba
- assim_param (2, 5 ,10, 12)   <== c'est defini par vcmax dans le dt 48

stomate_main
maint_respiration: variables output (PFT modifiés)
- lai (3, 5, 10)
-------------------------------------------------------------------
dans le dt = 50

sechiba
- assim_param(2, 5, 10 ,12)
- biomass (3, 5, 10)

erreur dans diffuco_trans_co2

lai (3,5,10) != vcmax(2,5,10,12)

So we have to think of a way to maintain the coherency between all these linked variables: biomass/lai/leaf_frac/leaf_age/vcmax/assim_param
All suggestions are welcome.

Albert & Fabienne.

Change History (3)

comment:1 Changed 7 years ago by ajornet

comment:2 Changed 6 years ago by nvuilsce

This bug has been fixed in revision r4189. There is an error in the number of the tcket mentionned in the commit?

comment:3 Changed 6 years ago by nvuilsce

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