Ignore:
Timestamp:
Jan 18, 2019, 2:18:01 PM (6 years ago)
Author:
aboissinot
Message:

Watersat subroutine is replaced by Psat_water in the thermal plume model.
cf. revision 1993 or 29/08/2018 JL comment.

Location:
trunk/LMDZ.GENERIC/libf/phystd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/phystd/thermcell_env.F90

    r2060 r2071  
    1212      USE print_control_mod, ONLY: prt_level
    1313      USE thermcell_mod, ONLY: RKAPPA
    14       USE watercommon_h, ONLY: RLvCp, RETV, watersat
     14      USE watercommon_h, ONLY: RLvCp, RETV, Psat_water
    1515      USE planete_mod, ONLY: preff
    1616     
     
    5454      INTEGER ig, ll
    5555     
     56      REAL dummy
     57     
    5658!==============================================================================
    5759! Initialization
     
    7577      DO ll=1,nlay
    7678         DO ig=1,ngrid
    77             CALL watersat(pt(ig,ll), pplev(ig,ll), pqsat(ig,ll))
     79            CALL Psat_water(pt(ig,ll), pplev(ig,ll), dummy, pqsat(ig,ll))
    7880         ENDDO
    7981      ENDDO
  • trunk/LMDZ.GENERIC/libf/phystd/thermcell_plume.F90

    r2069 r2071  
    2323      USE ioipsl_getin_p_mod, ONLY : getin_p
    2424      USE print_control_mod, ONLY: prt_level
    25       USE watercommon_h, ONLY: RLvCp, RETV, watersat
     25      USE watercommon_h, ONLY: RLvCp, RETV, Psat_water
    2626      USE thermcell_mod
    2727     
     
    119119      REAL zltup                                ! useless here
    120120     
     121      REAL dummy
     122     
    121123      LOGICAL active(ngrid)                     ! if the plume is active at ig,l (speed and incoming mass flux > 0 or l=lmin)
    122124      LOGICAL activetmp(ngrid)                  ! if the plus is active at ig,l (active=true and outgoing mass flux > 0)
     
    239241         
    240242         DO ig=1,ngrid
    241             CALL watersat(ztemp(ig), pplev(ig,l), zqsat(ig))
     243            CALL Psat_water(ztemp(ig), pplev(ig,l), dummy, zqsat(ig))
    242244         ENDDO
    243245         
     
    431433         DO ig=1,ngrid
    432434            IF (activetmp(ig)) THEN
    433                CALL watersat(ztemp(ig), pplev(ig,l), zqsatth(ig,l))
     435               CALL Psat_water(ztemp(ig), pplev(ig,l), dummy, zqsatth(ig,l))
    434436            ENDIF
    435437         ENDDO
Note: See TracChangeset for help on using the changeset viewer.