Changeset 1467

Timestamp:

Aug 4, 2015, 12:05:18 PM (9 years ago)

Author:

tnavarro

Message:

new option supersat to allow for supersaturation of water

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/callkeys.h (modified) (2 diffs)
• libf/phymars/conf_phys.F (modified) (1 diff)
• libf/phymars/physiq.F (modified) (1 diff)
• libf/phymars/watercloud.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -2190 +2190 @@
 - Moved "iniprint.h" back to dyn3d (and added link to it in dyn1d), as it
   is preferable to not share modules/commons between physics and dynamics.
+
+== 04/08/2015 == TN
+- New option supersat to allow for the supersaturation of water vapor.
+  Default value is True. The goal is to disable supersaturation even 
+  with microphysics for comparison purposes.

trunk/LMDZ.MARS/libf/phymars/callkeys.h

@@ -12 +12 @@
      &   ,callg2d,linear,rayleigh,tracer,active,doubleq,submicron       &
      &   ,lifting,freedust,callddevil,scavenging,sedimentation          &
-     &   ,activice,water,tifeedback,microphys,caps,photochem,calltherm  &
-     &   ,callrichsl,callslope,tituscap,callyamada4
+     &   ,activice,water,tifeedback,microphys,supersat,caps,photochem   &
+     &   ,calltherm,callrichsl,callslope,tituscap,callyamada4
      
       COMMON/callkeys_i/iradia,iaervar,iddist,ilwd,ilwb,ilwn,ncouche    &
@@ -57 +57 @@
       logical active,doubleq,submicron,lifting,callddevil,scavenging
       logical sedimentation
-      logical water,activice,tifeedback,microphys,caps
+      logical water,activice,tifeedback,microphys,supersat,caps
       logical photochem
       integer nltemodel

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

@@ -477 +477 @@
          write(*,*)" microphys = ",microphys
 
+! supersat
+         write(*,*)"Allow super-saturation of water vapor?"
+         supersat=.true. ! default value
+         call getin("supersat",supersat)
+         write(*,*)"supersat = ",supersat
+
 ! microphysical parameter contact       
          write(*,*) "water contact parameter ?"

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

@@ -2162 +2162 @@
      &                       'Mean reff',
      &                       'm',2,rave)
+            call WRITEDIAGFI(ngrid,'saturation',
+     &           'h2o vap saturation ratio','dimless',3,satu)
             if (scavenging) then
               CALL WRITEDIAGFI(ngrid,"Nccntot",

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

@@ -94 +94 @@
       REAL subpdt(ngrid,nlay)         ! cf. pdtcloud
 
+      ! no supersaturation when option supersat is false
+      REAL zt(ngrid,nlay)       ! local value of temperature
+      REAL zqsat(ngrid,nlay)    ! saturation
 
       INTEGER iq,ig,l
@@ -391 +394 @@
         ENDDO
       ENDIF
-        
-        
+
       IF(microphys) THEN
+
+       ! In case one does not want to allow supersatured water when using microphysics.
+       ! Not done by default.
+       IF(.not.supersat) THEN     
+        zt  = pt + (pdt+pdtcloud)*ptimestep
+        call watersat(ngrid*nlay,zt,pplay,zqsat)
+        DO l=1, nlay
+         DO ig=1,ngrid
+          IF (pq(ig,l,igcm_h2o_vap)
+     &      + (pdq(ig,l,igcm_h2o_vap) + pdqcloud(ig,l,igcm_h2o_vap))
+     &      * ptimestep .ge. zqsat(ig,l)) THEN
+             pdqcloud(ig,l,igcm_h2o_vap) = 
+     &         (zqsat(ig,l) - pq(ig,l,igcm_h2o_vap))/ptimestep
+     &        - pdq(ig,l,igcm_h2o_vap)
+             pdqcloud(ig,l,igcm_h2o_ice) = 
+     &         -pdqcloud(ig,l,igcm_h2o_vap)
+             ! no need to correct ccn_number, updaterad can handle this properly.
+          ENDIF
+         ENDDO
+        ENDDO        
+       ENDIF
        
        DO l=1, nlay