Changeset 2347 for trunk/LMDZ.MARS/libf/phymars/tcondco2.F90

Timestamp:

Jun 9, 2020, 2:52:48 PM (4 years ago)

Author:

aslmd

Message:

CO2 cloud microphysics

File:

• trunk/LMDZ.MARS/libf/phymars/tcondco2.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/tcondco2.F90

@@ -8 +8 @@
 ! Condensation temperature for co2 ice; based on 
 ! the saturation in co2sat.F JA17
-!--------------------------------------------------i
+!--------------------------------------------------
 
 integer, intent(in) :: ngrid,nlay
@@ -14 +14 @@
 double precision, intent(out), dimension(ngrid,nlay):: tcond ! CO2 condensation temperature     (atm)
 double precision:: A,B,pco2
+real :: qco2
 integer:: ig,l
 
 A=dlog(1.382d12)
 B=-3182.48
-
+qco2=0.
 DO l=1,nlay
    DO ig=1,ngrid
+!      qco2 = min(1., q(ig,l))!added by CM not sure, to be verified before
+                             !  cleaning
+!      pco2 = qco2 * (mmean(ig,l)/44.01) * p(ig,l)
       pco2 = q(ig,l) * (mmean(ig,l)/44.01) * p(ig,l)
       tcond(ig,l)=B/(dlog(pco2)-A)