Changeset 1181 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_optic.F90

Timestamp:

Jun 15, 2009, 6:49:35 PM (15 years ago)

Author:

jghattas

Message:

• Ajoute le possibilite de forcer avec les aerosols poussiere et sel de mer.
• Ajoute le module aerosol_mod qui contient des parmetres pour les aerosols
• flag_aerosol a legerement change signification

A Cozic

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_optic.F90 (modified) (11 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_optic.F90

@@ -4 +4 @@
      pplay, paprs, t_seri, rhcl, &
      mass_ins_aero, mass_ins_aero_pi, &
-     tau_aero, piz_aero, cg_aero )
+     tau_aero, piz_aero, cg_aero, &
+     tausum_aero, tau3d_aero )
 
 ! This routine will :
@@ -12 +13 @@
   
   USE dimphy
+  USE aerosol_mod
   IMPLICIT NONE
 
@@ -30 +32 @@
   REAL, DIMENSION(klon,klev), INTENT(OUT)     :: mass_ins_aero    ! Total mass for all indissoluble aerosols
   REAL, DIMENSION(klon,klev), INTENT(OUT)     :: mass_ins_aero_pi !     -"-     preindustrial values
-  REAL, DIMENSION(klon,klev,9,2), INTENT(OUT) :: tau_aero    ! Aerosol optical thickness
-  REAL, DIMENSION(klon,klev,9,2), INTENT(OUT) :: piz_aero    ! Single scattering albedo aerosol
-  REAL, DIMENSION(klon,klev,9,2), INTENT(OUT) :: cg_aero     ! asymmetry parameter aerosol
+  REAL, DIMENSION(klon,klev,naero_grp,nbands), INTENT(OUT) :: tau_aero    ! Aerosol optical thickness
+  REAL, DIMENSION(klon,klev,naero_grp,nbands), INTENT(OUT) :: piz_aero    ! Single scattering albedo aerosol
+  REAL, DIMENSION(klon,klev,naero_grp,nbands), INTENT(OUT) :: cg_aero     ! asymmetry parameter aerosol
+  REAL, DIMENSION(klon,nwave,naero_spc), INTENT(OUT)       :: tausum_aero
+  REAL, DIMENSION(klon,klev,nwave,naero_spc), INTENT(OUT)  :: tau3d_aero
 
 ! Local variables
 !****************************************************************************************
   REAL, DIMENSION(klon)        :: aerindex ! POLDER aerosol index 
-  REAL, DIMENSION(klon,10)     :: fractnat_allaer
   REAL, DIMENSION(klon,klev)   :: sulfate  ! SO4 aerosol concentration [ug/m3]
   REAL, DIMENSION(klon,klev)   :: bcsol    ! BC soluble concentration [ug/m3]
@@ -43 +46 @@
   REAL, DIMENSION(klon,klev)   :: pomsol   ! POM soluble concentration [ug/m3]
   REAL, DIMENSION(klon,klev)   :: pomins   ! POM insoluble concentration [ug/m3]
+  REAL, DIMENSION(klon,klev)   :: cidust    ! DUST aerosol concentration  [ug/m3]
+  REAL, DIMENSION(klon,klev)   :: sscoarse  ! SS Coarse concentration [ug/m3]
+  REAL, DIMENSION(klon,klev)   :: sssupco   ! SS Super Coarse concentration [ug/m3]
+  REAL, DIMENSION(klon,klev)   :: ssacu     ! SS Acumulation concentration [ug/m3]
   REAL, DIMENSION(klon,klev)   :: sulfate_pi
   REAL, DIMENSION(klon,klev)   :: bcsol_pi
@@ -48 +55 @@
   REAL, DIMENSION(klon,klev)   :: pomsol_pi
   REAL, DIMENSION(klon,klev)   :: pomins_pi
+  REAL, DIMENSION(klon,klev)   :: cidust_pi
+  REAL, DIMENSION(klon,klev)   :: sscoarse_pi
+  REAL, DIMENSION(klon,klev)   :: sssupco_pi
+  REAL, DIMENSION(klon,klev)   :: ssacu_pi
   REAL, DIMENSION(klon,klev)   :: pdel
-  REAL, DIMENSION(klon,klev,8) :: m_allaer
+  REAL, DIMENSION(klon,klev,naero_spc) :: m_allaer
+  REAL, DIMENSION(klon,naero_tot)      :: fractnat_allaer
 
   INTEGER :: k, i
@@ -59 +71 @@
 ! Read and interpolate sulfate
   IF ( flag_aerosol .EQ. 1 .OR. &
-       flag_aerosol .EQ. 4 .OR. &
        flag_aerosol .EQ. 6 ) THEN 
 
@@ -69 +80 @@
 ! Read and interpolate bcsol and bcins
   IF ( flag_aerosol .EQ. 2 .OR. &
-       flag_aerosol .EQ. 4 .OR. &
-       flag_aerosol .EQ. 5 ) THEN 
+       flag_aerosol .EQ. 6 ) THEN 
 
      ! Get bc aerosol distribution 
@@ -83 +93 @@
 ! Read and interpolate pomsol and pomins
   IF ( flag_aerosol .EQ. 3 .OR. &
-       flag_aerosol .EQ. 4 .OR. &
-       flag_aerosol .EQ. 5 .OR. &
        flag_aerosol .EQ. 6 ) THEN
 
@@ -95 +103 @@
 
 
+! Read and interpolate csssm, ssssm, assssm
+  IF (flag_aerosol .EQ. 4 .OR. &
+      flag_aerosol .EQ. 6 ) THEN 
+
+      CALL readaerosol_interp(1 ,rjourvrai, debut, pplay, paprs, sscoarse,sscoarse_pi) 
+      CALL readaerosol_interp(2 ,rjourvrai, debut, pplay, paprs, ssacu, ssacu_pi) 
+
+  ELSE
+     sscoarse(:,:) = 0. ; sscoarse_pi(:,:) = 0. 
+     ssacu(:,:) = 0.    ; ssacu_pi(:,:) = 0. 
+  ENDIF
+
+! Read and interpolate cidustm
+  IF (flag_aerosol .EQ. 5 .OR.  &
+      flag_aerosol .EQ. 6 ) THEN 
+
+      CALL readaerosol_interp(6 , rjourvrai, debut, pplay, paprs, cidust, cidust_pi) 
+
+  ELSE
+      cidust(:,:) = 0. ; cidust_pi(:,:) = 0. 
+  ENDIF
+
+!
 ! Store all aerosols in one variable
 !
-! ACo pour couplage aerosol offline 07/04/2009
-! Tableau contenant les masses pour tous les aerosols 
-! les valeurs a zero seront a remplacer par les bons
-! tableaux lorsque les routines de lectures seront
-! ajoutees. 
-  m_allaer(:,:,1) = 0.                ! SSSSM || CSSSM ! Coarse Soluble Sea Salt Mass
-  m_allaer(:,:,2) = 0.                ! ASSSM
+  m_allaer(:,:,1) = sscoarse(:,:)     ! SSSSM || CSSSM
+  m_allaer(:,:,2) = ssacu(:,:)        ! ASSSM
   m_allaer(:,:,3) = bcsol(:,:)        ! ASBCM
   m_allaer(:,:,4) = pomsol(:,:)       ! ASPOMM
   m_allaer(:,:,5) = sulfate(:,:)      ! ASSO4M || CSSO4M   
-  m_allaer(:,:,6) = 0.                ! CIDUSTM        ! Coarse Insoluble DUST Mass
+  m_allaer(:,:,6) = cidust(:,:)       ! CIDUSTM
   m_allaer(:,:,7) = bcins(:,:)        ! AIBCM
   m_allaer(:,:,8) = pomins(:,:)       ! AIPOMM
@@ -114 +140 @@
 ! Calculate the total mass of all indissoluble aersosols
 !
-  mass_ins_aero(:,:)    = sulfate(:,:)    + bcsol(:,:)    + pomsol(:,:) 
-  mass_ins_aero_pi(:,:) = sulfate_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) 
+  mass_ins_aero(:,:)    = sulfate(:,:)    + bcsol(:,:)    + pomsol(:,:)    + &
+       sscoarse(:,:)   + ssacu(:,:)    + cidust(:,:) 
+  mass_ins_aero_pi(:,:) = sulfate_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + &
+       sscoarse_pi(:,:) + ssacu_pi(:,:) + cidust_pi(:,:) 
 
 !****************************************************************************************
@@ -137 +165 @@
      
      ! aeropt_5wv only for validation and diagnostics.
-     ! In this version no diagnostics are set. 
-     ! jg : may be desactivated if no diagnostics added.
      CALL aeropt_5wv( &
           pdel, m_allaer, &
           pdtphys, rhcl, aerindex, & 
-          flag_aerosol, pplay, t_seri)
+          flag_aerosol, pplay, t_seri, &
+          tausum_aero, tau3d_aero)
   ELSE