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

Timestamp:

Nov 14, 2013, 12:00:25 PM (11 years ago)

Author:

tnavarro

Message:

bug with freedust mode

File:

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

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

@@ -266 +266 @@
       ENDIF
 
+      IF (dustmult .gt. 0) tauref(:) = tauref(:)*dustmult
+
 c -----------------------------------------------------------------
 c Computing the opacity in each layer
@@ -470 +472 @@
 
       IF (freedust) THEN
-        ! tauref has been initialized to 0 before.
+       ! tauref has been initialized to 0 before.
+       DO iaer=1,naerdust
         DO l=1,nlayer
-          DO ig=1,ngrid
-            tauref(ig) = tauref(ig) + 
+         DO ig=1,ngrid
+          tauref(ig) = tauref(ig) + 
      &                aerosol(ig,l,iaerdust(iaer))
           ENDDO
         ENDDO
+       ENDDO
       ENDIF