Changeset 1016 for trunk/LMDZ.GENERIC/libf/phystd/watercommon_h.F90

Timestamp:

Aug 7, 2013, 4:21:01 PM (12 years ago)

Author:

jleconte

Message:

07/08/2013 == JL

• Water cycle in double precision (largescale+moistadj)
• Improved wate rayleigh.
• First step for rayleigh with variable species. Now, just need to change optcv.
• changed some interpolation indices in callcorrk to limit dependency of OLR on the number of layers

File:

• trunk/LMDZ.GENERIC/libf/phystd/watercommon_h.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/watercommon_h.F90

@@ -175 +175 @@
       end subroutine Tsat_water
 
+!==================================================================
+      subroutine Psat_waterDP(T,p,psat,qsat)
+
+         implicit none
+
+!==================================================================
+!     Purpose
+!     -------
+!     Compute the saturation vapor pressure and mass mixing ratio at saturation (kg/kg)
+!     for a given pressure (Pa) and temperature (K)
+!     DOUBLE PRECISION
+!     Based on the Tetens formula from L.Li physical parametrization manual
+!
+!     Authors
+!     -------
+!     Jeremy Leconte (2012)
+!
+!==================================================================
+
+!        input
+         double precision, intent(in) :: T, p
+  
+!        output
+         double precision psat,qsat
+
+! JL12 variables for tetens formula
+         double precision,parameter :: Pref_solid_liquid=611.14d0
+         double precision,parameter :: Trefvaporization=35.86D0
+         double precision,parameter :: Trefsublimation=7.66d0
+         double precision,parameter :: Tmin=8.d0
+         double precision,parameter :: r3vaporization=17.269d0
+         double precision,parameter :: r3sublimation=21.875d0
+
+! checked vs. old watersat data 14/05/2012 by JL.
+
+         if (T.gt.T_h2O_ice_liq) then 
+            psat = Pref_solid_liquid*Exp(r3vaporization*(T-T_h2O_ice_liq)/(T-Trefvaporization)) ! liquid / vapour
+         else if (T.lt.Tmin) then
+            print*, "careful, T<Tmin in psat water"
+            psat = Pref_solid_liquid*Exp(r3sublimation*(Tmin-T_h2O_ice_liq)/(Tmin-Trefsublimation)) ! min psat  
+         else                 
+            psat = Pref_solid_liquid*Exp(r3sublimation*(T-T_h2O_ice_liq)/(T-Trefsublimation)) ! solid / vapour
+         endif
+         if(psat.gt.p) then
+            qsat=1.d0
+         else
+            qsat=epsi*psat/(p-(1.d0-epsi)*psat)
+         endif
+         return
+      end subroutine Psat_waterDP
+
+
+
+
+!==================================================================
+      subroutine Lcpdqsat_waterDP(T,p,psat,qsat,dqsat,dlnpsat)
+
+         implicit none
+
+!==================================================================
+!     Purpose
+!     -------
+!     Compute dqsat=L/cp*d (q_sat)/d T and dlnpsat=L/cp d(ln Psat)/d T
+!     for a given temperature (K)! 
+!     Based on the Tetens formula from L.Li physical parametrization manual
+!
+!     Authors
+!     -------
+!     Jeremy Leconte (2012)
+!
+!==================================================================
+
+!        input
+         double precision T, p, psat, qsat
+  
+!        output
+         double precision dqsat,dlnpsat
+
+! JL12 variables for tetens formula
+         double precision,parameter :: Pref_solid_liquid=611.14d0
+         double precision,parameter :: Trefvaporization=35.86d0
+         double precision,parameter :: Tmin=8.d0
+         double precision,parameter :: Trefsublimation=7.66d0
+         double precision,parameter :: r3vaporization=17.269d0
+         double precision,parameter :: r3sublimation=21.875d0
+
+         double precision :: dummy
+
+         if (psat.gt.p) then
+            dqsat=0.d0
+            return
+         endif
+
+         if (T.gt.T_h2O_ice_liq) then
+            dummy = r3vaporization*(T_h2O_ice_liq-Trefvaporization)/(T-Trefvaporization)**2  ! liquid / vapour
+         else if (T.lt.Tmin) then
+            print*, "careful, T<Tmin in Lcp psat water"
+            dummy = r3sublimation*(T_h2O_ice_liq-Trefsublimation)/(Tmin-Trefsublimation)**2  ! solid / vapour
+         else               
+            dummy = r3sublimation*(T_h2O_ice_liq-Trefsublimation)/(T-Trefsublimation)**2  ! solid / vapour
+         endif
+
+         dqsat=RLVTT/RCPD*qsat*(p/(p-(1.d0-epsi)*psat))*dummy
+         dlnpsat=RLVTT/RCPD*dummy
+         return
+      end subroutine Lcpdqsat_waterDP
+
 
 end module watercommon_h