source: LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90 @ 2709

Last change on this file since 2709 was 2693, checked in by oboucher, 8 years ago

Cosmetic changes
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:keywords set to Author Date Id Revi
File size: 15.1 KB
RevLine 
[2004]1! $Id: readaerosol_optic_rrtm.F90 2693 2016-11-01 09:08:00Z fairhead $
2!
[2146]3SUBROUTINE readaerosol_optic_rrtm(debut, aerosol_couple,  &
[2644]4     new_aod, flag_aerosol, flag_bc_internal_mixture, itap, rjourvrai, &
[2004]5     pdtphys, pplay, paprs, t_seri, rhcl, presnivs, &
[2146]6     tr_seri, mass_solu_aero, mass_solu_aero_pi, &
[2004]7     tau_aero, piz_aero, cg_aero, &
8     tausum_aero, tau3d_aero )
9
[2146]10  ! This routine will :
11  ! 1) recevie the aerosols(already read and interpolated) corresponding to flag_aerosol
12  ! 2) calculate the optical properties for the aerosols
13  !
14
[2004]15  USE dimphy
16  USE aero_mod
[2146]17  USE phys_local_var_mod, only: sconcso4,sconcno3,sconcoa,sconcbc,sconcss,sconcdust, &
18       concso4,concno3,concoa,concbc,concss,concdust,loadso4,loadoa,loadbc,loadss,loaddust, &
19       load_tmp1,load_tmp2,load_tmp3,load_tmp4,load_tmp5,load_tmp6,load_tmp7
20
[2504]21  USE infotrac_phy
[2693]22  USE YOMCST
[2146]23
[2004]24  IMPLICIT NONE
[2693]25
[2004]26  include "clesphys.h"
27
[2146]28  ! Input arguments
29  !****************************************************************************************
[2004]30  LOGICAL, INTENT(IN)                      :: debut
[2146]31  LOGICAL, INTENT(IN)                      :: aerosol_couple
[2004]32  LOGICAL, INTENT(IN)                      :: new_aod
33  INTEGER, INTENT(IN)                      :: flag_aerosol
[2644]34  LOGICAL, INTENT(IN)                      :: flag_bc_internal_mixture
[2004]35  INTEGER, INTENT(IN)                      :: itap
36  REAL, INTENT(IN)                         :: rjourvrai
37  REAL, INTENT(IN)                         :: pdtphys
38  REAL, DIMENSION(klon,klev), INTENT(IN)   :: pplay
39  REAL, DIMENSION(klon,klev+1), INTENT(IN) :: paprs
40  REAL, DIMENSION(klon,klev), INTENT(IN)   :: t_seri
41  REAL, DIMENSION(klon,klev), INTENT(IN)   :: rhcl   ! humidite relative ciel clair
42  REAL, DIMENSION(klev), INTENT(IN)        :: presnivs
[2146]43  REAL, DIMENSION(klon,klev,nbtr), INTENT(IN) :: tr_seri ! concentration tracer
[2004]44
[2146]45  ! Output arguments
46  !****************************************************************************************
[2004]47  REAL, DIMENSION(klon,klev), INTENT(OUT)     :: mass_solu_aero    ! Total mass for all soluble aerosols
48  REAL, DIMENSION(klon,klev), INTENT(OUT)     :: mass_solu_aero_pi !     -"-     preindustrial values
49  REAL, DIMENSION(klon,klev,2,NSW), INTENT(OUT) :: tau_aero    ! Aerosol optical thickness
50  REAL, DIMENSION(klon,klev,2,NSW), INTENT(OUT) :: piz_aero    ! Single scattering albedo aerosol
51  REAL, DIMENSION(klon,klev,2,NSW), INTENT(OUT) :: cg_aero     ! asymmetry parameter aerosol
52  REAL, DIMENSION(klon,nwave,naero_tot), INTENT(OUT)       :: tausum_aero
53  REAL, DIMENSION(klon,klev,nwave,naero_tot), INTENT(OUT)  :: tau3d_aero
54
[2146]55  ! Local variables
56  !****************************************************************************************
57  REAL, DIMENSION(klon)        :: aerindex      ! POLDER aerosol index
58  REAL, DIMENSION(klon,klev)   :: sulfacc       ! SO4 accumulation concentration [ug/m3]
59  REAL, DIMENSION(klon,klev)   :: sulfcoarse    ! SO4 coarse concentration [ug/m3]
60  REAL, DIMENSION(klon,klev)   :: bcsol         ! BC soluble concentration [ug/m3]
61  REAL, DIMENSION(klon,klev)   :: bcins         ! BC insoluble concentration [ug/m3]
62  REAL, DIMENSION(klon,klev)   :: pomsol        ! POM soluble concentration [ug/m3]
63  REAL, DIMENSION(klon,klev)   :: pomins        ! POM insoluble concentration [ug/m3]
64  REAL, DIMENSION(klon,klev)   :: cidust        ! DUST aerosol concentration  [ug/m3]
65  REAL, DIMENSION(klon,klev)   :: sscoarse      ! SS Coarse concentration [ug/m3]
66  REAL, DIMENSION(klon,klev)   :: sssupco       ! SS Super Coarse concentration [ug/m3]
67  REAL, DIMENSION(klon,klev)   :: ssacu         ! SS Acumulation concentration [ug/m3]
68  REAL, DIMENSION(klon,klev)   :: nitracc       ! nitrate accumulation concentration [ug/m3]
69  REAL, DIMENSION(klon,klev)   :: nitrcoarse    ! nitrate coarse concentration [ug/m3]
70  REAL, DIMENSION(klon,klev)   :: nitrinscoarse ! nitrate insoluble coarse concentration [ug/m3]
71  REAL, DIMENSION(klon,klev)   :: sulfacc_pi
72  REAL, DIMENSION(klon,klev)   :: sulfcoarse_pi
[2004]73  REAL, DIMENSION(klon,klev)   :: bcsol_pi
74  REAL, DIMENSION(klon,klev)   :: bcins_pi
75  REAL, DIMENSION(klon,klev)   :: pomsol_pi
76  REAL, DIMENSION(klon,klev)   :: pomins_pi
77  REAL, DIMENSION(klon,klev)   :: cidust_pi
78  REAL, DIMENSION(klon,klev)   :: sscoarse_pi
79  REAL, DIMENSION(klon,klev)   :: sssupco_pi
80  REAL, DIMENSION(klon,klev)   :: ssacu_pi
[2146]81  REAL, DIMENSION(klon,klev)   :: nitracc_pi
82  REAL, DIMENSION(klon,klev)   :: nitrcoarse_pi
83  REAL, DIMENSION(klon,klev)   :: nitrinscoarse_pi
84  REAL, DIMENSION(klon,klev)   :: pdel, zrho
85  REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer
86  REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer_pi !RAF 
[2180]87
[2146]88  integer :: id_ASBCM, id_ASPOMM, id_ASSO4M, id_ASMSAM, id_CSSO4M, id_CSMSAM, id_SSSSM
89  integer :: id_CSSSM, id_ASSSM, id_CIDUSTM, id_AIBCM, id_AIPOMM, id_ASNO3M, id_CSNO3M, id_CINO3M
[2004]90  INTEGER :: k, i
91
[2146]92  !--air density
93  zrho(:,:)=pplay(:,:)/t_seri(:,:)/RD                     !--kg/m3
[2004]94
[2146]95  !****************************************************************************************
96  ! 1) Get aerosol mass
97  !   
98  !****************************************************************************************
99  !
100  !
101  IF (aerosol_couple) THEN   !--we get aerosols from tr_seri array from INCA
102     !
103     !--copy fields from INCA tr_seri
104     !--convert to ug m-3 unit for consistency with offline fields
105     !
[2693]106     DO i=1,nbtr
107        SELECT CASE(trim(solsym(i)))
108           CASE ("ASBCM")
[2146]109              id_ASBCM = i
[2693]110           CASE ("ASPOMM")
[2146]111              id_ASPOMM = i
[2693]112           CASE ("ASSO4M")
[2146]113              id_ASSO4M = i
[2693]114           CASE ("ASMSAM")
[2146]115              id_ASMSAM = i
[2693]116           CASE ("CSSO4M")
[2146]117              id_CSSO4M = i
[2693]118           CASE ("CSMSAM")
[2146]119              id_CSMSAM = i
[2693]120           CASE ("SSSSM")
[2146]121              id_SSSSM = i
[2693]122           CASE ("CSSSM")
[2146]123              id_CSSSM = i
[2693]124           CASE ("ASSSM")
[2146]125              id_ASSSM = i
[2693]126           CASE ("CIDUSTM")
[2146]127              id_CIDUSTM = i
[2693]128           CASE ("AIBCM")
[2146]129              id_AIBCM = i
[2693]130           CASE ("AIPOMM")
[2146]131              id_AIPOMM = i
[2693]132           CASE ("ASNO3M")
[2146]133              id_ASNO3M = i
[2693]134           CASE ("CSNO3M")
[2146]135              id_CSNO3M = i
[2693]136           CASE ("CINO3M")
[2146]137              id_CINO3M = i
[2693]138           END SELECT
139     ENDDO
[2004]140
[2146]141     bcsol(:,:)        =   tr_seri(:,:,id_ASBCM)                         *zrho(:,:)*1.e9  ! ASBCM
142     pomsol(:,:)       =   tr_seri(:,:,id_ASPOMM)                        *zrho(:,:)*1.e9  ! ASPOMM
143     sulfacc(:,:)      =  (tr_seri(:,:,id_ASSO4M)+tr_seri(:,:,id_ASMSAM))*zrho(:,:)*1.e9  ! ASSO4M (=SO4) + ASMSAM (=MSA)
144     sulfcoarse(:,:)   =  (tr_seri(:,:,id_CSSO4M)+tr_seri(:,:,id_CSMSAM))*zrho(:,:)*1.e9  ! CSSO4M (=SO4) + CSMSAM (=MSA)
145     sssupco(:,:)      =   tr_seri(:,:,id_SSSSM)                         *zrho(:,:)*1.e9  ! SSSSM
146     sscoarse(:,:)     =   tr_seri(:,:,id_CSSSM)                         *zrho(:,:)*1.e9  ! CSSSM
147     ssacu(:,:)        =   tr_seri(:,:,id_ASSSM)                         *zrho(:,:)*1.e9  ! ASSSM
148     cidust(:,:)       =   tr_seri(:,:,id_CIDUSTM)                       *zrho(:,:)*1.e9  ! CIDUSTM
149     bcins(:,:)        =   tr_seri(:,:,id_AIBCM)                         *zrho(:,:)*1.e9  ! AIBCM
150     pomins(:,:)       =   tr_seri(:,:,id_AIPOMM)                        *zrho(:,:)*1.e9  ! AIPOMM
151     nitracc(:,:)      =   tr_seri(:,:,id_ASNO3M)                        *zrho(:,:)*1.e9  ! ASNO3M
152     nitrcoarse(:,:)   =   tr_seri(:,:,id_CSNO3M)                        *zrho(:,:)*1.e9  ! CSNO3M
153     nitrinscoarse(:,:)=   tr_seri(:,:,id_CINO3M)                        *zrho(:,:)*1.e9  ! CINO3M
154     !
155     bcsol_pi(:,:)        =   0.0 ! ASBCM pre-ind
156     pomsol_pi(:,:)       =   0.0 ! ASPOMM pre-ind
157     sulfacc_pi(:,:)      =   0.0 ! ASSO4M (=SO4) + ASMSAM (=MSA) pre-ind
158     sulfcoarse_pi(:,:)   =   0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) pre-ind
159     sssupco_pi(:,:)      =   0.0 ! SSSSM pre-ind
160     sscoarse_pi(:,:)     =   0.0 ! CSSSM pre-ind
161     ssacu_pi(:,:)        =   0.0 ! ASSSM pre-ind
162     cidust_pi(:,:)       =   0.0 ! CIDUSTM pre-ind
163     bcins_pi(:,:)        =   0.0 ! AIBCM pre-ind
164     pomins_pi(:,:)       =   0.0 ! AIPOMM pre-ind
165     nitracc_pi(:,:)      =   0.0 ! ASNO3M pre-ind
166     nitrcoarse_pi(:,:)   =   0.0 ! CSNO3M pre-ind
167     nitrinscoarse_pi(:,:)=   0.0 ! CINO3M
168     !
169  ELSE !--not aerosol_couple
170     !
171     ! Read and interpolate sulfate
[2693]172     IF ( flag_aerosol .EQ. 1 .OR. flag_aerosol .EQ. 6 ) THEN
[2004]173
[2146]174        CALL readaerosol_interp(id_ASSO4M_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, sulfacc, sulfacc_pi,loadso4)
175     ELSE
176        sulfacc(:,:) = 0. ; sulfacc_pi(:,:) = 0.
177        loadso4=0.
[2693]178     ENDIF
[2004]179
[2146]180     ! Read and interpolate bcsol and bcins
[2693]181     IF ( flag_aerosol .EQ. 2 .OR. flag_aerosol .EQ. 6 ) THEN
[2004]182
[2146]183        ! Get bc aerosol distribution
184        CALL readaerosol_interp(id_ASBCM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, bcsol, bcsol_pi, load_tmp1 )
185        CALL readaerosol_interp(id_AIBCM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, bcins, bcins_pi, load_tmp2 )
186        loadbc(:)=load_tmp1(:)+load_tmp2(:)
187     ELSE
188        bcsol(:,:) = 0. ; bcsol_pi(:,:) = 0.
189        bcins(:,:) = 0. ; bcins_pi(:,:) = 0.
190        loadbc=0.
[2693]191     ENDIF
[2004]192
[2146]193     ! Read and interpolate pomsol and pomins
[2693]194     IF ( flag_aerosol .EQ. 3 .OR. flag_aerosol .EQ. 6 ) THEN
[2004]195
[2146]196        CALL readaerosol_interp(id_ASPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomsol, pomsol_pi, load_tmp3)
197        CALL readaerosol_interp(id_AIPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomins, pomins_pi, load_tmp4)
198        loadoa(:)=load_tmp3(:)+load_tmp4(:)
199     ELSE
200        pomsol(:,:) = 0. ; pomsol_pi(:,:) = 0.
201        pomins(:,:) = 0. ; pomins_pi(:,:) = 0.
202        loadoa=0.
[2693]203     ENDIF
[2004]204
[2146]205     ! Read and interpolate csssm, ssssm, assssm
[2693]206     IF (flag_aerosol .EQ. 4 .OR. flag_aerosol .EQ. 6 ) THEN
[2004]207
[2148]208        CALL readaerosol_interp(id_SSSSM_phy ,itap, pdtphys,rjourvrai, &
209        debut, pplay, paprs, t_seri, sssupco, sssupco_pi, load_tmp5)
210        CALL readaerosol_interp(id_CSSSM_phy ,itap, pdtphys,rjourvrai, &
211        debut, pplay, paprs, t_seri, sscoarse,sscoarse_pi, load_tmp6)
212        CALL readaerosol_interp(id_ASSSM_phy ,itap, pdtphys, rjourvrai, &
213        debut, pplay, paprs, t_seri, ssacu, ssacu_pi, load_tmp7)
[2146]214        loadss(:)=load_tmp5(:)+load_tmp6(:)+load_tmp7(:)
215     ELSE
216        sscoarse(:,:) = 0. ; sscoarse_pi(:,:) = 0.
217        ssacu(:,:)    = 0. ; ssacu_pi(:,:) = 0.
218        sssupco(:,:)  = 0. ; sssupco_pi = 0.
219        loadss=0.
220     ENDIF
[2004]221
[2146]222     ! Read and interpolate cidustm
[2693]223     IF (flag_aerosol .EQ. 5 .OR. flag_aerosol .EQ. 6 ) THEN
[2004]224
[2146]225        CALL readaerosol_interp(id_CIDUSTM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, cidust, cidust_pi, loaddust)
226
227     ELSE
228        cidust(:,:) = 0. ; cidust_pi(:,:) = 0.
229        loaddust=0.
230     ENDIF
231     !
232     sulfcoarse(:,:)      =   0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA)
233     sulfcoarse_pi(:,:)   =   0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) pre-ind
234     !
235     !--placeholder for offline nitrate   
236     !
237     nitracc(:,:)         =   0.0
238     nitracc_pi(:,:)      =   0.0
239     nitrcoarse(:,:)      =   0.0
240     nitrcoarse_pi(:,:)   =   0.0
241     nitrinscoarse(:,:)   =   0.0
242     nitrinscoarse_pi(:,:)=   0.0
243
244  ENDIF !--not aerosol_couple
245
246  !
247  ! Store all aerosols in one variable
248  !
249  m_allaer(:,:,id_ASBCM_phy)  = bcsol(:,:)        ! ASBCM
250  m_allaer(:,:,id_ASPOMM_phy) = pomsol(:,:)       ! ASPOMM
251  m_allaer(:,:,id_ASSO4M_phy) = sulfacc(:,:)      ! ASSO4M (= SO4)
252  m_allaer(:,:,id_CSSO4M_phy) = sulfcoarse(:,:)   ! CSSO4M
253  m_allaer(:,:,id_SSSSM_phy)  = sssupco(:,:)      ! SSSSM
254  m_allaer(:,:,id_CSSSM_phy)  = sscoarse(:,:)     ! CSSSM
255  m_allaer(:,:,id_ASSSM_phy)  = ssacu(:,:)        ! ASSSM
256  m_allaer(:,:,id_CIDUSTM_phy)= cidust(:,:)       ! CIDUSTM
257  m_allaer(:,:,id_AIBCM_phy)  = bcins(:,:)        ! AIBCM
258  m_allaer(:,:,id_ASNO3M_phy) = nitracc(:,:)      ! ASNO3M
259  m_allaer(:,:,id_CSNO3M_phy) = nitrcoarse(:,:)   ! CSNO3M
260  m_allaer(:,:,id_CINO3M_phy) = nitrinscoarse(:,:)! CINO3M
261  m_allaer(:,:,id_AIPOMM_phy) = pomins(:,:)       ! AIPOMM
262  m_allaer(:,:,id_STRAT_phy)  = 0.0
263
264  !RAF
265  m_allaer_pi(:,:,id_ASBCM_phy)  = bcsol_pi(:,:)        ! ASBCM pre-ind
266  m_allaer_pi(:,:,id_ASPOMM_phy) = pomsol_pi(:,:)       ! ASPOMM pre-ind
267  m_allaer_pi(:,:,id_ASSO4M_phy) = sulfacc_pi(:,:)      ! ASSO4M (= SO4) pre-ind
268  m_allaer_pi(:,:,id_CSSO4M_phy) = sulfcoarse_pi(:,:)   ! CSSO4M pre-ind
269  m_allaer_pi(:,:,id_SSSSM_phy)  = sssupco_pi(:,:)      ! SSSSM pre-ind
270  m_allaer_pi(:,:,id_CSSSM_phy)  = sscoarse_pi(:,:)     ! CSSSM pre-ind
271  m_allaer_pi(:,:,id_ASSSM_phy)  = ssacu_pi(:,:)        ! ASSSM pre-ind
272  m_allaer_pi(:,:,id_CIDUSTM_phy)= cidust_pi(:,:)       ! CIDUSTM pre-ind
273  m_allaer_pi(:,:,id_AIBCM_phy)  = bcins_pi(:,:)        ! AIBCM pre-ind
274  m_allaer_pi(:,:,id_ASNO3M_phy) = nitracc_pi(:,:)      ! ASNO3M pre-ind
275  m_allaer_pi(:,:,id_CSNO3M_phy) = nitrcoarse_pi(:,:)   ! CSNO3M pre-ind
276  m_allaer_pi(:,:,id_CINO3M_phy) = nitrinscoarse_pi(:,:)! CINO3M pre-ind
277  m_allaer_pi(:,:,id_AIPOMM_phy) = pomins_pi(:,:)       ! AIPOMM pre-ind
278  m_allaer_pi(:,:,id_STRAT_phy)  = 0.0
279
280  !
281  ! Calculate the total mass of all soluble aersosols
282  ! to be revisited for AR6
[2323]283  mass_solu_aero(:,:)    = sulfacc(:,:)    + bcsol(:,:)    + pomsol(:,:)   + nitracc(:,:) + ssacu(:,:)
284  mass_solu_aero_pi(:,:) = sulfacc_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + nitracc_pi(:,:) + ssacu_pi(:,:)
[2146]285
286  !****************************************************************************************
287  ! 2) Calculate optical properties for the aerosols
288  !
289  !****************************************************************************************
[2004]290  DO k = 1, klev
291     DO i = 1, klon
292        pdel(i,k) = paprs(i,k) - paprs (i,k+1)
[2693]293     ENDDO
294  ENDDO
[2004]295
[2325]296!--new aerosol properties
[2146]297  ! aeropt_6bands for rrtm
[2644]298  CALL aeropt_6bands_rrtm(          &
299       pdel, m_allaer, rhcl,        &
300       tau_aero, piz_aero, cg_aero, &
301       m_allaer_pi, flag_aerosol,   &
302       flag_bc_internal_mixture, zrho )
[2004]303
[2146]304  ! aeropt_5wv only for validation and diagnostics
[2644]305  CALL aeropt_5wv_rrtm(              &
306       pdel, m_allaer,               &
307       rhcl, aerindex,               &
308       flag_aerosol,                 &
309       flag_bc_internal_mixture,     &
310       pplay, t_seri,                &
[2146]311       tausum_aero, tau3d_aero )
[2004]312
[2146]313  ! Diagnostics calculation for CMIP5 protocol
314  sconcso4(:)  =m_allaer(:,1,id_ASSO4M_phy)*1.e-9
315  sconcno3(:)  =(m_allaer(:,1,id_ASNO3M_phy)+m_allaer(:,1,id_CSNO3M_phy)+m_allaer(:,1,id_CINO3M_phy))*1.e-9
316  sconcoa(:)   =(m_allaer(:,1,id_ASPOMM_phy)+m_allaer(:,1,id_AIPOMM_phy))*1.e-9
317  sconcbc(:)   =(m_allaer(:,1,id_ASBCM_phy)+m_allaer(:,1,id_AIBCM_phy))*1.e-9
318  sconcss(:)   =(m_allaer(:,1,id_ASSSM_phy)+m_allaer(:,1,id_CSSSM_phy)+m_allaer(:,1,id_SSSSM_phy))*1.e-9
319  sconcdust(:) =m_allaer(:,1,id_CIDUSTM_phy)*1.e-9
320  concso4(:,:) =m_allaer(:,:,id_ASSO4M_phy)*1.e-9
321  concno3(:,:) =(m_allaer(:,:,id_ASNO3M_phy)+m_allaer(:,:,id_CSNO3M_phy)+m_allaer(:,:,id_CINO3M_phy))*1.e-9
322  concoa(:,:)  =(m_allaer(:,:,id_ASPOMM_phy)+m_allaer(:,:,id_AIPOMM_phy))*1.e-9
323  concbc(:,:)  =(m_allaer(:,:,id_ASBCM_phy)+m_allaer(:,:,id_AIBCM_phy))*1.e-9
324  concss(:,:)  =(m_allaer(:,:,id_ASSSM_phy)+m_allaer(:,:,id_CSSSM_phy)+m_allaer(:,:,id_SSSSM_phy))*1.e-9
325  concdust(:,:)=m_allaer(:,:,id_CIDUSTM_phy)*1.e-9
326
[2004]327END SUBROUTINE readaerosol_optic_rrtm
Note: See TracBrowser for help on using the repository browser.