Opened 7 years ago

Closed 7 years ago

#215 closed defect (fixed)

bug report on DGVM

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


At Line 1429 in slowproc.f90, it says "IF (ok_stomate) CALL stomate_veget_update(kjpindex,veget_max)", and stomate_veget_update (defined in stomate.f90) says, IF veget_max < min_vegfrac (0.001), everything will be reset to zero.

It's very simple to debug, just add "IF (ok_stomate .AND. .NOT. ok_dgvm)" before the CALL. Below is a point test after this change, and veget_max evolves normally.

Attachments (1)

zd.gif (6.3 KB) - added by dzhu 7 years ago.
a point test

Download all attachments as: .zip

Change History (5)

Changed 7 years ago by dzhu

a point test

comment:1 Changed 7 years ago by jgipsl

The correction has been done in the trunk revision [3038].

The bug was introduced in revision [2718], see the related ticket #107

comment:2 Changed 7 years ago by jgipsl

Nicolas Vuichard found that it is not only a bug for the DGVM. The subroutine stomate_veget_update must not be called before lcchange.

A correction calling lcchange in the beginning of the sechiba time step will be tested. The lines in stomate_veget_update will be put back in the end of lcchange subroutine.

comment:3 Changed 7 years ago by jgipsl

  • Owner changed from somebody to jgipsl
  • Status changed from new to accepted

comment:4 Changed 7 years ago by jgipsl

  • Resolution set to fixed
  • Status changed from accepted to closed

This has finaly been corrected also for the case DGVM=n in rev [3094], see ticket #107

Note: See TracTickets for help on using tickets.