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

Timestamp:

Sep 23, 2009, 4:47:40 PM (15 years ago)

Author:

Laurent Fairhead

Message:

• En deconnectant les aérosols (ok_ade=ok_aie=n) on a les mêmes

résultats avant et après les modifs.

• preindustrial readin fields are used to compute natural aerosol fields

to allow for clean double calls to radiation

• full forcing diagnostics (NAT, ANT, ZERO, Cloud forcing, CS,AS) are

activated with lev_histmth 4, If lev_histmth is not 4, the call to the
radiation is minimized, for efficiency, but ade and aie are computed and
applied (however for species wise forcing one would need to do
difference runs) (still quite a bit new forcing info, requires probably

some more explanation)

• there is a hardcoded key in sw_aeroAR4.F90 which lets you choose to use the zero aerosol, or natural aerosol perturbation acting on the meteorology, but still would put out the full forcing diagnostics.
• aod fields from offline aerosol fields are also output in histmth for

all aerosol tracers read in and available for evaluation

• aeropt contains the ss humidity correction from nicolas&yves

File:

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

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

@@ -63 +63 @@
   REAL, DIMENSION(klon,klev)   :: pdel
   REAL, DIMENSION(klon,klev,naero_spc) :: m_allaer
-  REAL, DIMENSION(klon,naero_tot)      :: fractnat_allaer
+  REAL, DIMENSION(klon,klev,naero_spc) :: m_allaer_pi !RAF  
+!  REAL, DIMENSION(klon,naero_tot)      :: fractnat_allaer !RAF delete??
 
   INTEGER :: k, i
@@ -143 +144 @@
   m_allaer(:,:,id_AIPOMM) = pomins(:,:)       ! AIPOMM
 
+!RAF
+  m_allaer_pi(:,:,1)  = bcsol_pi(:,:)        ! ASBCM pre-ind
+  m_allaer_pi(:,:,2)  = pomsol_pi(:,:)       ! ASPOMM pre-ind
+  m_allaer_pi(:,:,3)  = sulfate_pi(:,:)      ! ASSO4M (= SO4) pre-ind
+  m_allaer_pi(:,:,4)  = 0.                ! CSSO4M pre-ind
+  m_allaer_pi(:,:,5)  = sssupco_pi(:,:)      ! SSSSM pre-ind
+  m_allaer_pi(:,:,6)  = sscoarse_pi(:,:)     ! CSSSM pre-ind
+  m_allaer_pi(:,:,7)  = ssacu_pi(:,:)        ! ASSSM pre-ind
+  m_allaer_pi(:,:,8)  = cidust_pi(:,:)       ! CIDUSTM pre-ind
+  m_allaer_pi(:,:,9)  = bcins_pi(:,:)        ! AIBCM pre-ind
+  m_allaer_pi(:,:,10) = pomins_pi(:,:)       ! AIPOMM pre-ind
+
 !
 ! Calculate the total mass of all soluble aersosols
 !
-  mass_solu_aero(:,:)    = sulfate(:,:)    + bcsol(:,:)    + pomsol(:,:)    + &
-       sscoarse(:,:)    + ssacu(:,:)    + sssupco(:,:) 
-  mass_solu_aero_pi(:,:) = sulfate_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + &
-       sscoarse_pi(:,:) + ssacu_pi(:,:) + sssupco_pi(:,:)
+  mass_solu_aero(:,:)    = sulfate(:,:)    + bcsol(:,:)    + pomsol(:,:) !   + &
+!       sscoarse(:,:)    + ssacu(:,:)    + sssupco(:,:) 
+  mass_solu_aero_pi(:,:) = sulfate_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) ! + &
+!       sscoarse_pi(:,:) + ssacu_pi(:,:) + sssupco_pi(:,:)
 
 !****************************************************************************************
@@ -163 +176 @@
   IF (new_aod) THEN 
 
-     fractnat_allaer(:,:) = 0.
-     CALL aeropt_2bands(                 &
+! RAF delete??     fractnat_allaer(:,:) = 0.
+! RAF fractnat_allaer -> m_allaer_pi
+
+     CALL aeropt_2bands( &
           pdel, m_allaer, pdtphys, rhcl, & 
           tau_aero, piz_aero, cg_aero,   &
-          fractnat_allaer, flag_aerosol, &
+          m_allaer_pi, flag_aerosol, &
           pplay, t_seri, presnivs)