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

Timestamp:

Jun 9, 2020, 2:15:04 PM (4 years ago)

Author:

aslmd

Message:

CO2 cloud microphysics

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=======================================================================
+      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     Local:
+c     ------
+      integer ::
+     &   i ! loop on naersize
 c=======================================================================
-c
-c
-c  now:  straight psat of CO2 (or qsat of CO2 but need of mmean)
-c
+c===== BEGIN
 c=======================================================================
-
-c   declarations:
-c   -------------
-c   arguments:
-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
-          
-
+      do i = 1, naersize
+        psat(i) = 1.382 * 1e12 * exp(-3182.48/t(i))
+      end do
+c=======================================================================
+c===== END
+c=======================================================================
       RETURN
       END