Opened 9 years ago

Closed 7 years ago

#31 closed defect (fixed)

Problem with small frac_nobio

Reported by: jgipsl Owned by: aducharne
Priority: major Milestone: ORCHIDEE 2.0
Component: Physical processes Version:
Keywords: Cc:

Description

The modele stops with CWRR while coupling LMDZ-ORCHIDEE(trunk sept 2012), with ipslerr(3) everywhere in hydrol.F90.

During hydrol_split_soil there is a problem with mass conservation in grid cells where frac_nobio is very small.

The bug disapears when preventiong frac_nobio from being too small in slowproc. See attached file but why ?

Correction done everywhere frac_nobio is modified :

 DO ji = 1, kjpindex
      IF ( ( SUM(frac_nobio(ji,:)) .LT. 0.001 ) ) THEN
        frac_nobio(ji,:) = zero
        SUMveg = SUM(frac_nobio(ji,:))+SUM(veget_max(ji,:))
        frac_nobio(ji,:) = frac_nobio(ji,:)/SUMveg
        veget_max(ji,:) = veget_max(ji,:)/SUMveg
      ENDIF
ENDDO

Change History (4)

comment:1 Changed 9 years ago by jpolcher

Where does the model crash ? In which ipslerr call ?

It does not make sense that limiting frac_nobio solves the problem as frac_nobio=0 is a perfectly legitimate value.

Could it be that some operations within hydrol_split_soil have forgotten that it is only vegtot+tot_frac_nobio=1 ?

comment:2 Changed 8 years ago by peylin

  • Owner changed from somebody to aducharne
  • Status changed from new to assigned

comment:3 Changed 8 years ago by aducharne

  • Status changed from assigned to accepted

comment:4 Changed 7 years ago by aducharne

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

Commit done in rev [1118]

Note: See TracTickets for help on using tickets.