Changeset 3902


Ignore:
Timestamp:
May 17, 2021, 6:45:51 PM (4 years ago)
Author:
jyg
Message:

Bug fix in wx_pbl_var_mod.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/wx_pbl_var_mod.F90

    r3897 r3902  
    401401      qsat = MIN(0.5,qsat)
    402402      zcor = 1./(1.-RETV*qsat)
    403       qsat0_x(j) = qsat*zcor
    404       dqsatdT0_x(j) = FOEDE(Ts0_x(j),zdelta,zcvm5,qsat0_x(j),zcor)
     403      qsat0_x(j) = fqsat*qsat*zcor
     404      dqsatdT0_x(j) = fqsat*FOEDE(Ts0_x(j),zdelta,zcvm5,qsat0_x(j),zcor)
    405405!!      dqsatdT0_x(j) = (RLvtt*(1.-zdelta)+RLSTT*zdelta)*qsat0_x(j)/(Rv*Ts0_x(j)*Ts0_x(j))
    406406   ENDDO
     
    413413      qsat = MIN(0.5,qsat)
    414414      zcor = 1./(1.-RETV*qsat)
    415       qsat0_w(j) = qsat*zcor
    416       dqsatdT0_w(j) = FOEDE(Ts0_w(j),zdelta,zcvm5,qsat0_w(j),zcor)
     415      qsat0_w(j) = fqsat*qsat*zcor
     416      dqsatdT0_w(j) = fqsat*FOEDE(Ts0_w(j),zdelta,zcvm5,qsat0_w(j),zcor)
    417417!!      dqsatdT0_w(j) = (RLvtt*(1.-zdelta)+RLSTT*zdelta)*qsat0_w(j)/(Rv*Ts0_w(j)*Ts0_w(j))
    418418   ENDDO
Note: See TracChangeset for help on using the changeset viewer.