Changeset 411 for trunk/LMDZ.MARS/libf/phymars/aeropacity.F

Timestamp:

Nov 22, 2011, 5:34:42 PM (13 years ago)

Author:

tnavarro

Message:

changed scavenging in improvedclouds.F, updated ice radius in callsedim.F and commented outputs in suaer.F90 & aeropacity.F

File:

• trunk/LMDZ.MARS/libf/phymars/aeropacity.F (modified) (2 diffs)

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

@@ -428 +428 @@
           tauscaling(ig) = tauref(ig) * 
      &                     pplev(ig,1) / 700.E0 / taudusttmp(ig)
+c          tauscaling(ig) = 1.e-4
       ENDDO
 
@@ -435 +436 @@
           DO ig=1,ngrid
             aerosol(ig,l,iaerdust(iaer)) = max(1E-20,
-     &                   tauref(ig) *
-     &                   pplev(ig,1) / 700.E0 *
-     &                   aerosol(ig,l,iaerdust(iaer)) / 
-     &                   taudusttmp(ig)
-     &                                        )
+     &                aerosol(ig,l,iaerdust(iaer))* tauscaling(ig))
           ENDDO
         ENDDO
       ENDDO
-
+      
+c output for debug
+        IF (ngrid.NE.1) THEN
+             CALL WRITEDIAGFI(ngridmx,'taudusttmp','virtual tau dust',
+     &      '#',2,taudusttmp)
+             CALL WRITEDIAGFI(ngridmx,'tausca','tauscaling',
+     &      '#',2,tauscaling)
+        ELSE
+             CALL WRITEDIAGFI(ngridmx,'taudusttmp','virtual tau dust',
+     &      '#',0,taudusttmp)
+             CALL WRITEDIAGFI(ngridmx,'tausca','tauscaling',
+     &      '#',0,tauscaling)
+        ENDIF
 c -----------------------------------------------------------------
 c Column integrated visible optical depth in each point