Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4_AR5/libf/phylmd/aeropt_5wv.F90

    r1337 r1417  
    850850!  ENDDO
    851851
    852    DO i=1, klon
    853       od550aer(i)=SUM(tausum(i,2,:))
    854       od865aer(i)=SUM(tausum(i,5,:))
    855       DO k=1, KLEV
    856          ec550aer(i,k)=SUM(tau(i,k,2,:))/dh(i,k)
    857       END DO   
    858    END DO
     852  DO i=1, klon
     853     od550aer(i)=0.
     854     DO m=1,naero_spc
     855        od550aer(i)=od550aer(i)+tausum(i,2,m)
     856     END DO
     857  END DO
     858  DO i=1, klon
     859     od865aer(i)=0.
     860     DO m=1,naero_spc
     861        od865aer(i)=od865aer(i)+tausum(i,5,m)
     862     END DO
     863  END DO
     864  DO i=1, klon
     865     DO k=1, KLEV
     866        ec550aer(i,k)=0.
     867        DO m=1,naero_spc
     868           ec550aer(i,k)=ec550aer(i,k)+tau(i,k,2,m)/dh(i,k)
     869        END DO
     870     END DO
     871  END DO
     872 
    859873   od550lt1aer(:)=tausum(:,2,id_ASSO4M)+tausum(:,2,id_ASBCM)+tausum(:,2,id_AIBCM)+ &
    860874        tausum(:,2,id_ASPOMM)+tausum(:,2,id_AIPOMM)+tausum(:,2,id_ASSSM)+ &
Note: See TracChangeset for help on using the changeset viewer.