Opened 6 years ago

Last modified 4 years ago

#388 assigned enhancement

Incoherence in qsatcalc when temperature is too high

Reported by: maignan Owned by: peylin
Priority: minor Milestone: Not scheduled yet
Component: Physical processes Version: trunc
Keywords: Cc:


Sujet : 	Re: [orchidee-help] soil temperature exceed threshold (of 1M K) in grid w/o vegetation: model stops
Date : 	Thu, 24 Aug 2017 11:38:12 +0200
De : 	Fabienne Maignan <>
Pour : 	Daniel Goll <>,


Le 11/08/2017 à 12:03, Daniel Goll a écrit :
> Hi
> my simulation stop because the soil temperature in a gridbox without vegetation growing since 1yr of the simulation exceeds the temperature of the core of the sun (1M K; threshold 'missing insrc_sechiba/enerbil.f90) after several hundred of years. (The gridbox is a coastal point (where there is actually ocean in real life); so no vegetation is okay.)
> 1. Are you serious about the threshold? What is the rational behind this threshold? 
The logic is the following one:
In qsat_moisture.f90, the temperature is tested against the threshold max_temp defined in the module as 370K.

In the old version (subroutine qsat_old) if the temperature exceeded the threshold and the flag diag_qsat was activated (TRUE is the default value), the return value for qsat was set to 999999.
This mechanism has disappeared in the new version (qsatcalc).

In enerbil_begin, qsatcalc is called, if the flag diag_qsat is activated qsol_sat is tested against missing = 999998 and then stops if qsol_sat is higher.

> 2. Any suggestion what to do?
To be coherent, the return value for qsat should be set again to 999999 if the temperature is too high, or another more elegant solution should be implemented.
However I'm not aware of the reasons why the qsat subroutines were modified.
I'll open a ticket to track this.

Attachments (1)

FCTTRE.h (2.3 KB) - added by fcheruy 4 years ago.
calcul de qsat dans LMDZ

Download all attachments as: .zip

Change History (7)

comment:1 Changed 5 years ago by luyssaert

  • Milestone changed from IPSLCM6.v2 to IPSLCM6.v1

comment:2 Changed 5 years ago by jgipsl

  • Milestone IPSLCM6.v1 deleted

comment:3 Changed 5 years ago by aducharne

  • Milestone set to ORCHIDEE 3.0
  • Owner changed from somebody to peylin
  • Status changed from new to assigned

Réunion tickets du 18/12/18:
soit on met un stop qd temp_sol devient trop grande (on veut que le code plante ici et pas plus loin), soit on prend une valeur plus raisonnable à la valeur max de qsol_sat

Changed 4 years ago by fcheruy

calcul de qsat dans LMDZ

comment:4 Changed 4 years ago by aducharne

Proposition le 20/01/2020:

  • prendre la fonction qsatcalc de LMDZ.
  • mettre un warning si temp_sol > 100°C
  • mettre un stop qd on dépasse une valeur plus grande, cf h_gardefou de LMDZ

=> PP will do it and commit will be done with JG

Last edited 4 years ago by jgipsl (previous) (diff)

comment:5 Changed 4 years ago by peylin

  • Milestone changed from ORCHIDEE 3.0 to ORCHIDEE 4.0

comment:6 Changed 4 years ago by luyssaert

  • Milestone changed from ORCHIDEE 4.0 to Not scheduled yet
Note: See TracTickets for help on using tickets.