Index: LMDZ6/trunk/libf/phylmd/wx_pbl_var_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/wx_pbl_var_mod.F90	(revision 3901)
+++ LMDZ6/trunk/libf/phylmd/wx_pbl_var_mod.F90	(revision 3902)
@@ -401,6 +401,6 @@
       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,6 +413,6 @@
       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