Changeset 2188 for LMDZ5/trunk/libf/phylmd/surf_land_mod.F90

Timestamp:

Jan 30, 2015, 5:25:12 PM (10 years ago)

Author:

fhourdin

Message:

Modification of the numerical scheme for coupling of surface temperature
with Orchidee.

Modification du schéma numérique pour le couplage de la température
de surface avec Orchidee.
Cette modification concerne la gestion des pas de temps différents pour
l'appel au rayonnement et au modèle de surface.

Orchidee prend en entrée le flux LW descendant.
Dans la version actuelle, on recalculait ce flux descendant à partir
du flux net issu du rayonnement, comme si il ne variait pas dans
le temps.
pbl_surface recalculait ensuite un flux descendant en ajoutant
au flux net descendant l'emission du sol, simga Ts4.
Mais du coup, si la température du sol partait en dérive froide entre
deux appels au rayonnement, on ne tenait pas compte du rappel
exercé par le rayonnement : le flux descendant ne bouge
pas et réhauffe la surface tandis que celle ci se refroidit de
moins en moins par émission LW.

Dans la nouvelle version, on passe le LWdown directement depuis
le rayonnement.
Il s'agit de deux hypothèses extrêmes :
Avant on supposait que le LWnet ne bougeait pas d'un appel au rayonnement
sur l'autre. Maintenant on suppose que c'est le LWdown.
La réalité est entre les deux car le WLdown a tendance à
suivre la température de surface, quand la température de l'air
et de la surface sont fortement couplées.

La nouvelle façon de faire devrait limiter la dérive vers le froid
dans les conditions de découplage nocturne. Un découplage extrême
avait amené à limiter la fonction de Richardson dans le calcul du coefficient
de trraînée Cd f(Ri) à F(ri) > 0.6 plutôt que 0.1 sur le cas 1D Dice.
Les corrections pourraient également résoudre des plantages
aléatoires sur l'Himalaya.

Jean-Louis Dufresne, Fuxing Wang, Frédéric Hourdin

PS : les changements concernent l'interface entre physiq et
pbl_surface, ainsi que phyredem et phyetat0 pour garantir 1+1=2.
On en a profité pour faire un brin de toilette.

File:

• LMDZ5/trunk/libf/phylmd/surf_land_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/surf_land_mod.F90

@@ -85 +85 @@
     REAL, DIMENSION(klon) :: pref_tmp
     REAL, DIMENSION(klon) :: swdown     ! downwelling shortwave radiation at land surface
-    REAL, DIMENSION(klon) :: lwdown     ! downwelling longwave radiation at land surface
     REAL, DIMENSION(klon) :: epot_air           ! potential air temperature
     REAL, DIMENSION(klon) :: tsol_rad, emis_new ! output from interfsol not used
@@ -106 +105 @@
        pref_tmp(1:knon)  = pref(1:knon)/100.
 ! 
-!* Calculate incoming flux for SW and LW interval: swdown, lwdown
+!* Calculate incoming flux for SW and LW interval: swdown
 !
        swdown(:) = 0.0
-       lwdown(:) = 0.0
        DO i = 1, knon
           swdown(i) = swnet(i)/(1-albedo(i))
-          lwdown(i) = lwnet(i) + RSIGMA*tsurf(i)**4
        END DO
 !