Changeset 1880 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Dec 22, 2017, 3:05:23 PM (7 years ago)

Author:

mvals

Message:

Mars GCM
-In the sub-grid scale cloud scheme zt is replaced is by ztclf and zq by zqclf ('clf' for 'cloud fraction') in order to avoid any confusion for the further schemes, which need initialization.
-zteff and pqeff are initialized in the first part of the CLFvarying scheme, section 0 of the code, instead of been initialized in section 1 (tendencies) with the sub-timesteps.
-The cloud fraction cannot be lower than the settled value mincloud.(MV)

File:

• trunk/LMDZ.MARS/libf/phymars/watercloud_mod.F (modified) (5 diffs)

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

@@ -109 +109 @@
 
 ! Representation of sub-grid water ice clouds A. Pottier 2013
-!      REAL :: zt(ngrid, nlay)
-      REAL :: zq(ngrid, nlay,nq)
+      REAL :: ztclf(ngrid, nlay)
+      REAL :: zqclf(ngrid, nlay,nq)
       REAL :: zdelt  
       REAL :: norm
       REAL :: ponder
       REAL :: tcond(ngrid,nlay)
-      REAL ::  zqvap(ngrid,nlay)
+      REAL :: zqvap(ngrid,nlay)
       REAL :: zqice(ngrid,nlay)
       REAL ::  spant ! delta T for the temperature distribution
@@ -167 +167 @@
 c   0.  Representation of sub-grid water ice clouds
 c------------------
+c-----Initialization
+      zteff(1:ngrid,1:nlay) = 0
+      pqeff(1:ngrid,1:nlay,1:nq) = 0
+      DO l=1,nlay
+        DO ig=1,ngrid
+             zteff(ig,l)=pt(ig,l)
+        END DO
+      END DO
+      DO l=1,nlay
+        DO ig=1,ngrid
+          DO iq=1,nq
+             pqeff(ig,l,iq)=pq(ig,l,iq)
+          ENDDO
+        ENDDO
+      ENDDO
 c-----Tendencies
       DO l=1,nlay
-       DO ig=1,ngrid
-          zt(ig,l)=pt(ig,l)+ pdt(ig,l)*ptimestep
-       ENDDO
+        DO ig=1,ngrid
+          ztclf(ig,l)=pt(ig,l)+ pdt(ig,l)*ptimestep
+        ENDDO
       ENDDO
       DO l=1,nlay
         DO ig=1,ngrid
           DO iq=1,nq
-             zq(ig,l,iq)=pq(ig,l,iq)+pdq(ig,l,iq)*ptimestep 
+             zqclf(ig,l,iq)=pq(ig,l,iq)+pdq(ig,l,iq)*ptimestep 
           ENDDO
         ENDDO
@@ -189 +204 @@
              flagcloud=.false.
          END IF
-         CALL watersat(ngrid*nlay,zt,pplay,zqsat)
-         zqvap=zq(:,:,igcm_h2o_vap)
-         zqice=zq(:,:,igcm_h2o_ice)
+         !CALL watersat(ngrid*nlay,ztclf,pplay,zqsat) !MV17: we dont need zqsat in the CLFvarying scheme
+         zqvap=zqclf(:,:,igcm_h2o_vap)
+         zqice=zqclf(:,:,igcm_h2o_ice)
          CALL tcondwater(ngrid*nlay,pplay,zqvap+zqice,tcond)
          DO l=1,nlay
            DO ig=1,ngrid
-              zdelt=spant !MAX(spant*zt(ig,l),1.e-12), now totally in K. Fixed width
-              IF (tcond(ig,l) .ge. (zt(ig,l)+zdelt)) THEN
-                 zteff(ig,l)=zt(ig,l)
+              zdelt=spant !MAX(spant*ztclf(ig,l),1.e-12), now totally in K. Fixed width
+              IF (tcond(ig,l) .ge. (ztclf(ig,l)+zdelt)) THEN
+                 zteff(ig,l)=ztclf(ig,l)
                  cloudfrac(ig,l)=1.
-              ELSE IF (tcond(ig,l) .le. (zt(ig,l)-zdelt)) THEN
-                 zteff(ig,l)=zt(ig,l)-zdelt
+              ELSE IF (tcond(ig,l) .le. (ztclf(ig,l)-zdelt)) THEN
+                 zteff(ig,l)=ztclf(ig,l)-zdelt
                  cloudfrac(ig,l)=mincloud
               ELSE
-                 cloudfrac(ig,l)=(tcond(ig,l)-zt(ig,l)+zdelt)/
+                 cloudfrac(ig,l)=(tcond(ig,l)-ztclf(ig,l)+zdelt)/
      &                           (2.0*zdelt)
-                 zteff(ig,l)=(tcond(ig,l)+zt(ig,l)-zdelt)/2.
+                 zteff(ig,l)=(tcond(ig,l)+ztclf(ig,l)-zdelt)/2.
               END IF
               zteff(ig,l)=zteff(ig,l)-pdt(ig,l)*ptimestep
-              IF (cloudfrac(ig,l).le.0) THEN
+              IF (cloudfrac(ig,l).le.mincloud) THEN !MV17: replaced .le.0 by .le.mincloud
                  cloudfrac(ig,l)=mincloud
               ELSE IF (cloudfrac(ig,l).gt.1) THEN
@@ -227 +242 @@
             totcloudfrac(ig)=MAX(totcloudfrac(ig)/norm,1.e-12) ! min value if NaNs
          ENDDO
-c-----No sub-grid cloud representation (CLFvarying=false)
-      ELSE
-         DO l=1,nlay
-            DO ig=1,ngrid
-               zteff(ig,l)=pt(ig,l)
-            END DO
-         END DO
+c-----Effective tracers quantities in the cloud fraction
+         IF (microphys) THEN
+            pqeff(:,:,igcm_ccn_mass)=pq(:,:,igcm_ccn_mass)/
+     &                              cloudfrac(:,:)
+            pqeff(:,:,igcm_ccn_number)=pq(:,:,igcm_ccn_number)/
+     &                              cloudfrac(:,:)
+         END IF ! end if (microphys)
+         pqeff(:,:,igcm_h2o_ice)=pq(:,:,igcm_h2o_ice)/
+     &                           cloudfrac(:,:)
       END IF ! end if (CLFvarying)
-
+      !! MV17: CLFvarying outputs
+      CALL WRITEDIAGFI(ngrid,'pqeffice','pqeffice',
+     &             'kg/kg',3,pqeff(:,:,igcm_h2o_ice))
+      CALL WRITEDIAGFI(ngrid,'zteff','zteff',
+     &             'K',3,zteff(:,:))
+      CALL WRITEDIAGFI(ngrid,'tcond','tcond',
+     &             'K',3,tcond(:,:))
+      CALL WRITEDIAGFI(ngrid,'cloudfrac','cloudfrac',
+     &             'K',3,cloudfrac(:,:))
 c------------------------------------------------------------------
-c Time subsampling for microphysics 
+c Time subsampling for microphysics
+c------------------------------------------------------------------ 
       rhocloud(1:ngrid,1:nlay) = rho_dust
       DO microstep=1,imicro 
@@ -284 +310 @@
      &        + pdq(ig,l,igcm_h2o_vap)
           ENDDO
-        ENDDO
-c------ Effective tracers quantities in the cloud fraction
-        IF (CLFvarying) THEN     
-           pqeff(:,:,:)=pq(:,:,:) ! prevent from buggs (A. Pottier)
-           pqeff(:,:,igcm_ccn_mass) =pq(:,:,igcm_ccn_mass)/
-     &                                cloudfrac(:,:)
-           pqeff(:,:,igcm_ccn_number)= 
-     &     pq(:,:,igcm_ccn_number)/cloudfrac(:,:)
-           pqeff(:,:,igcm_h2o_ice)= pq(:,:,igcm_h2o_ice)/
-     &                               cloudfrac(:,:)
-        ELSE
-           pqeff(:,:,:)=pq(:,:,:)
-           pqeff(:,:,igcm_ccn_mass)= pq(:,:,igcm_ccn_mass)
-           pqeff(:,:,igcm_ccn_number)= pq(:,:,igcm_ccn_number)
-           pqeff(:,:,igcm_h2o_ice)= pq(:,:,igcm_h2o_ice)
-        END IF      
+        ENDDO      
         
 c-------------------------------------------------------------------