Changeset 2158 for LMDZ5/trunk/libf/phylmd

Timestamp:

Nov 27, 2014, 4:06:03 PM (10 years ago)

Author:

Laurent Fairhead

Message:

Modification pour éviter de faire des calculs avec un q négatif qui peut entraîner
des calculs de t2m irréalistes (raison du plantage avec ORCHIDEE en mode debug a priori)

---

Modifications to prevent calculations with a negative humidity that would lead to
unrealistic calculations for temperature

File:

• LMDZ5/trunk/libf/phylmd/coefcdrag.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/coefcdrag.F90

@@ -78 +78 @@
        ztsolv(i) = ts(i)
 !       ztvd(i) = t(i) * (psol(i)/pref(i))**RKAPPA
-       ztvd(i) = (t(i)+zdphi(i)/RCPD/(1.+RVTMP2*q(i))) &
-          *(1.+RETV*q(i))
+!       ztvd(i) = (t(i)+zdphi(i)/RCPD/(1.+RVTMP2*q(i))) &
+!          *(1.+RETV*q(i))
+       ztvd(i) = (t(i)+zdphi(i)/RCPD/(1.+RVTMP2*q(i)))
        trm0(i) = 1. + RETV * max(qsurf(i),0.0)
        trm1(i) = 1. + RETV * max(q(i),0.0)
        ztsolv(i) = ztsolv(i) * trm0(i)
-!       ztvd(i) = ztvd(i) * trm1(i)
+       ztvd(i) = ztvd(i) * trm1(i)
        zri1(i) = zdphi(i)*(ztvd(i)-ztsolv(i))/(zdu2(i)*ztvd(i))
 !