Changeset 1484 for trunk/LMDZ.GENERIC/libf/phystd/condense_co2.F90

Timestamp:

Oct 15, 2015, 2:35:35 PM (9 years ago)

Author:

mturbet

Message:

Harmonization of surface tracer tendencies calculation

File:

• trunk/LMDZ.GENERIC/libf/phystd/condense_co2.F90 (modified) (7 diffs)

trunk/LMDZ.GENERIC/libf/phystd/condense_co2.F90

@@ -2 +2 @@
           pcapcal,pplay,pplev,ptsrf,pt,                  &
           pphi,pdt,pdu,pdv,pdtsrf,pu,pv,pq,pdq,          &
-          piceco2,albedo,pemisurf,                       &
+          piceco2,pdqsurfc,albedo,pemisurf,              &
           albedo_bareground,albedo_co2_ice_SPECTV,       &
           pdtc,pdtsrfc,pdpsrf,pduc,pdvc,                 &
@@ -89 +89 @@
       REAL,INTENT(OUT) :: pdvc(ngrid,nlayer)
       REAL,INTENT(OUT) :: pdqc(ngrid,nlayer,nq)
+      REAL,INTENT(OUT) :: pdqsurfc(ngrid)
 
 !-----------------------------------------------------------------------
@@ -102 +103 @@
       REAL ztnuc (ngrid,nlayer)
       REAL ztcondsol(ngrid) 
-      REAL zdiceco2(ngrid)
       REAL zcondicea(ngrid,nlayer), zcondices(ngrid)
       REAL zfallice(ngrid), Mfallice(ngrid) 
@@ -225 +225 @@
          pdpsrf(ig) = 0.
          condsub(ig) = .false.
-         zdiceco2(ig) = 0.
+         pdqsurfc(ig) = 0.
       ENDDO
 
@@ -435 +435 @@
 !     Changing CO2 ice amount and pressure
 
-            zdiceco2(ig) = zcondices(ig) + zfallice(ig)
-            piceco2(ig)  = piceco2(ig) + zdiceco2(ig)*ptimestep
-            pdpsrf(ig)   = -zdiceco2(ig)*g
+            pdqsurfc(ig) = zcondices(ig) + zfallice(ig)
+            piceco2(ig)  = piceco2(ig) + pdqsurfc(ig)*ptimestep
+            pdpsrf(ig)   = -pdqsurfc(ig)*g
 
             IF(ABS(pdpsrf(ig)*ptimestep).GT.pplev(ig,1)) THEN
@@ -461 +461 @@
          if(.not.piceco2(ig).ge.0.) THEN
             if(piceco2(ig).le.-1.e-8) print*,   &
-                'WARNING in condense_co2cloud: piceco2(',ig,')=', piceco2(ig)
+              'WARNING : in condense_co2cloud: piceco2(',ig,')=', piceco2(ig)
             piceco2(ig)=0.
          endif
-         if (piceco2(ig) .gt. 1.) then  ! CO2 Albedo condition changed to 1 mm coverage. Change by MT2015.
+         if (piceco2(ig) .gt. 1.) then  ! CO2 Albedo condition changed to ~1 mm coverage. Change by MT2015.
             DO nw=1,L_NSPECTV
                albedo(ig,nw) = albedo_co2_ice_SPECTV(nw)
@@ -476 +476 @@
             pemisurf(ig)  = emissiv
          end if
+         
+         piceco2(ig)  = piceco2(ig) - pdqsurfc(ig)*ptimestep ! This line was added so that tendencies are added outside the routine. MT2015.
+         if(.not.piceco2(ig).ge.0.) THEN
+            if(piceco2(ig).le.-1.e-8) print*,   &
+              'WARNING 2 : in condense_co2cloud: piceco2(',ig,')=', piceco2(ig)
+            piceco2(ig)=0.
+         endif
+         
       end do