Changeset 1417 for LMDZ4

Timestamp:

Jul 16, 2010, 2:06:55 PM (14 years ago)

Author:

jghattas

Message:

Corrected bug in calculation of the diagnostic variable ec550aer
"Extinction at 550nm". (Maria Raffaella Vuolo, LSCE)

File:

• LMDZ4/branches/LMDZ4_AR5/libf/phylmd/aeropt_5wv.F90 (modified) (1 diff)

LMDZ4/branches/LMDZ4_AR5/libf/phylmd/aeropt_5wv.F90

@@ -850 +850 @@
 !  ENDDO
 
-   DO i=1, klon
-      od550aer(i)=SUM(tausum(i,2,:))
-      od865aer(i)=SUM(tausum(i,5,:))
-      DO k=1, KLEV
-         ec550aer(i,k)=SUM(tau(i,k,2,:))/dh(i,k)
-      END DO    
-   END DO
+  DO i=1, klon
+     od550aer(i)=0.
+     DO m=1,naero_spc
+        od550aer(i)=od550aer(i)+tausum(i,2,m)
+     END DO
+  END DO
+  DO i=1, klon
+     od865aer(i)=0.
+     DO m=1,naero_spc
+        od865aer(i)=od865aer(i)+tausum(i,5,m)
+     END DO
+  END DO
+  DO i=1, klon
+     DO k=1, KLEV
+        ec550aer(i,k)=0.
+        DO m=1,naero_spc
+           ec550aer(i,k)=ec550aer(i,k)+tau(i,k,2,m)/dh(i,k)
+        END DO
+     END DO
+  END DO
+  
    od550lt1aer(:)=tausum(:,2,id_ASSO4M)+tausum(:,2,id_ASBCM)+tausum(:,2,id_AIBCM)+ &
         tausum(:,2,id_ASPOMM)+tausum(:,2,id_AIPOMM)+tausum(:,2,id_ASSSM)+ &