Changeset 2346 for trunk/LMDZ.MARS/libf/phymars/co2sat.F

Timestamp:

Jun 9, 2020, 2:20:46 PM (4 years ago)

Author:

aslmd

Message:

Cancelled commits 2344 and 2345. Sorry. No bad commit but wrong username to be fixed.

File:

• trunk/LMDZ.MARS/libf/phymars/co2sat.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/co2sat.F

@@ +1 @@
+      SUBROUTINE co2sat(naersize,t,p,psat)
+c       SUBROUTINE co2sat(naersize,t,p,qsat) JA
+      IMPLICIT NONE
+
 c=======================================================================
-c     SUBROUTINE co2sat
-c-----------------------------------------------------------------------
-c     Aim:
-c     ----
-c     Compute saturated steam pressure (from James et al, 1992)
+c
+c
+c  now:  straight psat of CO2 (or qsat of CO2 but need of mmean)
+c
 c=======================================================================
-      subroutine co2sat(naersize, t, psat)
-     
-      implicit none
-c-----------------------------------------------------------------------
-c     VARIABLES
-c-----------------------------------------------------------------------
-c     Inputs:
-c     ------- 
-      integer, intent(in) ::
-     &   naersize ! dimension of tables t and psat
-     
-      real, intent(in) ::
-     &   t(naersize) ! temperature table
 
-c     Output:
-c     -------
-      real, intent(out) ::
-     &   psat(naersize) ! Saturated steam pressure (Pa)
+c   declarations:
+c   -------------
+c   arguments:
+c   ----------
 
-c     Local:
-c     ------
-      integer ::
-     &   i ! loop on naersize
-c=======================================================================
-c===== BEGIN
-c=======================================================================
-      do i = 1, naersize
-        psat(i) = 1.382 * 1e12 * exp(-3182.48/t(i))
-      end do
-c=======================================================================
-c===== END
-c=======================================================================
+c   INPUT
+      integer naersize 
+      real t(naersize) , p(naersize)
+c   OUTPUT
+c      real qsat(naersize) JA
+      real psat(naersize)
+
+c   local:
+c   ------
+      INTEGER i
+      REAL r2,r3,r4 , To, es
+      SAVE r2,r3,r4
+      DATA r2,r3,r4/611.14,21.875,7.66/
+      SAVE To
+      DATA To/273.16/
+          
+      do i=1,naersize
+
+
+c        pression de vapeur saturante (James et al. 1992):
+
+          psat(i)  = 1.382 * 1e12 * exp(-3182.48/t(i)) !; (Pa)
+
+c         OR:
+
+c         qsat(i) = psat/p(i)*44.01/mmean ! Need of updated information on mmean
+c         qsat(i) = max(qsat(i), 1.e-30)
+
+
+      enddo
+c      qsat=psat JA
+          
+
       RETURN
       END