Changeset 703 for trunk/LMDZ.MARS/libf

Timestamp:

Jun 12, 2012, 4:29:19 PM (13 years ago)

Author:

tnavarro

Message:

corrected bug of tracers handling in watercloud that led to a incorrect CO2 depletion. Actually, all non-water cycle tracers were not modified by processes that occur before watercloud between commit 633 and this one.

File:

• trunk/LMDZ.MARS/libf/phymars/watercloud.F (modified) (1 diff)

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

@@ -262 +262 @@
        ENDDO
 c------ Tracers tendencies pdqcloud
-       DO iq=1,nq
-        DO l=1,nlay
+       DO l=1,nlay
          DO ig=1,ngrid
-            pdqcloud(ig,l,iq) = subpdq(ig,l,iq)/real(imicro)
-     &       - pdq(ig,l,iq)
+            pdqcloud(ig,l,igcm_h2o_ice) = 
+     &        subpdq(ig,l,igcm_h2o_ice)/real(imicro) 
+     &       - pdq(ig,l,igcm_h2o_ice)
+            pdqcloud(ig,l,igcm_h2o_vap) = 
+     &        subpdq(ig,l,igcm_h2o_vap)/real(imicro) 
+     &       - pdq(ig,l,igcm_h2o_vap)
+          IF(microphys) THEN
+            pdqcloud(ig,l,igcm_dust_mass) = 
+     &        subpdq(ig,l,igcm_dust_mass)/real(imicro) 
+     &       - pdq(ig,l,igcm_dust_mass)
+            pdqcloud(ig,l,igcm_dust_number) = 
+     &        subpdq(ig,l,igcm_dust_number)/real(imicro) 
+     &       - pdq(ig,l,igcm_dust_number)
+            pdqcloud(ig,l,igcm_ccn_mass) = 
+     &        subpdq(ig,l,igcm_ccn_mass)/real(imicro) 
+     &       - pdq(ig,l,igcm_ccn_mass)
+            pdqcloud(ig,l,igcm_ccn_number) = 
+     &        subpdq(ig,l,igcm_ccn_number)/real(imicro) 
+     &       - pdq(ig,l,igcm_ccn_number)
+           ENDIF
          ENDDO
-        ENDDO
        ENDDO