Changeset 1337 for LMDZ4/trunk/libf/phylmd/aeropt_5wv.F90

Timestamp:

Apr 2, 2010, 1:31:05 PM (14 years ago)

Author:

Laurent Fairhead

Message:

Additions to aerosol outputs for CMIP5 exercise

---

Additions aux sorties aérosols pour l'exercice CMIP5

Michael, Anne

File:

• LMDZ4/trunk/libf/phylmd/aeropt_5wv.F90 (modified) (4 diffs)

LMDZ4/trunk/libf/phylmd/aeropt_5wv.F90

@@ -11 +11 @@
   USE DIMPHY
   USE aero_mod
+  USE phys_local_var_mod, only: od550aer,od865aer,ec550aer,od550lt1aer
 
   !
@@ -141 +142 @@
   REAL :: tau3d(KLON,KLEV), piz3d(KLON,KLEV), cg3d(KLON,KLEV)
   REAL :: abs3d(KLON,KLEV)     ! epaisseur optique d'absorption
-
+  REAL :: dh(KLON,KLEV)
   
   REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble)   ! ext. coeff. Soluble comp. units *** m2/g 
@@ -624 +625 @@
 !      IF (pplay(i,k).EQ.0) stop  'stop aeropt_5wv p '
       zrho=pplay(i,k)/t_seri(i,k)/RD                  ! kg/m3
+      dh(i,k)=pdel(i,k)/(gravit*zrho)
 !CDIR UNROLL=naero_spc
       mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9
@@ -848 +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
+   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)+ &
+        0.03*tausum(:,2,id_CSSSM)+0.4*tausum(:,2,id_CIDUSTM)
+
+
+
   DEALLOCATE(aerosol_name)