Changeset 3902 for LMDZ6/trunk/libf

Timestamp:

May 17, 2021, 6:45:51 PM (3 years ago)

Author:

jyg

Message:

Bug fix in wx_pbl_var_mod.F90

File:

• LMDZ6/trunk/libf/phylmd/wx_pbl_var_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/wx_pbl_var_mod.F90

@@ -401 +401 @@
       qsat = MIN(0.5,qsat)
       zcor = 1./(1.-RETV*qsat)
-      qsat0_x(j) = qsat*zcor
-      dqsatdT0_x(j) = FOEDE(Ts0_x(j),zdelta,zcvm5,qsat0_x(j),zcor)
+      qsat0_x(j) = fqsat*qsat*zcor
+      dqsatdT0_x(j) = fqsat*FOEDE(Ts0_x(j),zdelta,zcvm5,qsat0_x(j),zcor)
 !!      dqsatdT0_x(j) = (RLvtt*(1.-zdelta)+RLSTT*zdelta)*qsat0_x(j)/(Rv*Ts0_x(j)*Ts0_x(j))
    ENDDO
@@ -413 +413 @@
       qsat = MIN(0.5,qsat)
       zcor = 1./(1.-RETV*qsat)
-      qsat0_w(j) = qsat*zcor
-      dqsatdT0_w(j) = FOEDE(Ts0_w(j),zdelta,zcvm5,qsat0_w(j),zcor)
+      qsat0_w(j) = fqsat*qsat*zcor
+      dqsatdT0_w(j) = fqsat*FOEDE(Ts0_w(j),zdelta,zcvm5,qsat0_w(j),zcor)
 !!      dqsatdT0_w(j) = (RLvtt*(1.-zdelta)+RLSTT*zdelta)*qsat0_w(j)/(Rv*Ts0_w(j)*Ts0_w(j))
    ENDDO