su_aerp.F90

Go to the documentation of this file.

SUBROUTINE su_aerp

!**** *SU_AERP*   - INITIALIZE MODULES YOEAERSRC, YOEAERSNK

!     PURPOSE.
!     --------
!           INITIALIZE YOEAERSRC AND YOEAERSNK, THE MODULES THAT CONTAINS 
!           COEFFICIENTS NEEDED TO RUN THE PROGNOSTIC AEROSOLS

!**   INTERFACE.
!     ----------
!        *CALL* *SU_AERP

!        EXPLICIT ARGUMENTS :
!        --------------------
!        NONE

!        IMPLICIT ARGUMENTS :
!        --------------------
!        YOEAERSRC, YOEAERSNK, YOEAERATM

!     METHOD.
!     -------
!        SEE DOCUMENTATION

!     EXTERNALS.
!     ----------

!     REFERENCE.
!     ----------
!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS

!     AUTHOR.
!     -------
!        JEAN-JACQUES MORCRETTE *ECMWF*
!        from O.BOUCHER (LOA, 1998-03) 

!     MODIFICATIONS.
!     --------------
!        ORIGINAL : 2004-05-10

!     ------------------------------------------------------------------

USE parkind1  ,ONLY : jprb
USE yomhook   ,ONLY : lhook,   dr_hook

USE yoeaersrc , ONLY : rssflx

USE yoeaersnk , ONLY : r_r, r_s, ralphar, ralphas, rfraer, rfrgas, &
  &  rrhmax, rrhtab, rrho_ss, rssgrow, rmmd_ss, rmmd_dd, rrho_dd, &
  &  rfrbc , rfrif , rfrom  , rfrso4 , rfrdd  , rfrss , rho_wat, rho_ice, &
  &  rvdpoce, rvdpsic, rvdplnd, rvdplic, rvsedoce, rvsedsic, rvsedlnd, rvsedlic, &
  &  nbrh

USE yoeaeratm , ONLY : rmasse, rmfmin, ndd1, nss1

IMPLICIT NONE

REAL(KIND=JPRB) :: ZHOOK_HANDLE
!     ----------------------------------------------------------------
IF (lhook) CALL dr_hook('SU_AERP',0,zhook_handle)

!-- For the ECMWF model, the following tables when dimensioned to 12 
!   can refer to 12 values of RH
!      (RHTAB, RHHO_SS, RSSGROW)
!   or to 12 types/bins of aerosols with the following mapping: NTYPAER
!   1- 3  sea-salt  0.03 - 0.5 -  5  - 20 microns                                  1
!   4- 6  dust      0.03 - 0.5 - 0.9 - 20 microns                                  2
!   7- 8  POM       hydrophilic, hydrophobic                                       3
!   9-10  BC        hydrophilic, hydrophobic                                       4
!  11     sulfate                                                                                                  5
!  12     fly ash                                                                                                  6
!  13     stratospheric aerosols                                                                   7
!  14     volcanic aerosols                                                                        8
!  15                                                                                                              9
!      (RVDPOCE, RVDSIC, RVDPLND, RVDPLIC)
!      (RVSEDOCE,RVSEDSIC,RVSEDLND,RVSEDLIC)



!*      1.    PARAMETERS RELATED TO SOURCES
!             -----------------------------

!-- OB's original 12 types and 24 different values had the following mapping
! DMS SO2 SO4 H2S DMSO MSA H2O2
! BC(2) OM(2) FlyAsh DU(2) SS(10)

!-- OB's original SS 10 bins     
! bin sizes: 0.03-0.06-0.13-0.25-0.5-1.0-2.0-5.0-10.-15.-20
!RSSFLX = (/ &
!  &  0.20526E-09_JPRB, 0.49292E-09_JPRB, 0.97079E-09_JPRB, 0.31938E-08_JPRB &
!  &, 0.16245E-07_JPRB, 0.86292E-07_JPRB, 0.31326E-06_JPRB, 0.24671E-06_JPRB &
!  &, 0.14109E-06_JPRB, 0.11784E-06_JPRB /)

! maximum possible number of aerosol types
!NMAXTAER=9   already defined in SU_AERW

!N.B. Fluxes of sea salt for each size bin are given in mg m-2 s-1 at wind 
!     speed of 1 m s-1 at 10m height (at 80% RH) in OB's seasalt.F
! RSSFLX also in mg m-2 s-1       
!-- OB's ECMWF 3 bins of sea salt: 0.03, 0.5, 5, 20 microns  
rssflx = (/ 4.85963536e-09_jprb, 4.15358556e-07_jprb, 5.04905813e-07_jprb /)

! OB's original vdep_oce, vdep_sic, vdep_ter, vdep_lic were given over 24 values 

! following 12 values for 10 SS and 2 DU in cm s-1
!RVDPOCE = (/   0.1_JPRB,   1.2_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB &
!          &,   0.1_JPRB,   1.2_JPRB,   1.2_JPRB,   1.2_JPRB,   1.5_JPRB,   1.5_JPRB /)
!
!RVDPSIC = (/   0.1_JPRB,   1.2_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB &
!          &,   0.1_JPRB,   1.2_JPRB,   1.2_JPRB,   1.2_JPRB,   1.5_JPRB,   1.5_JPRB /)
!
!RVDPLND = (/   0.1_JPRB,   1.2_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB &
!          &,   0.1_JPRB,   1.2_JPRB,   1.2_JPRB,   1.2_JPRB,   1.5_JPRB,   1.5_JPRB /)
!
!RVDPLIC = (/   0.1_JPRB,   1.2_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB &
!          &,   0.1_JPRB,   1.2_JPRB,   1.2_JPRB,   1.2_JPRB,   1.5_JPRB,   1.5_JPRB /)



!*      2.    PARAMETERS RELATED TO SINKS
!             ---------------------------

r_r = 0.001_jprb
r_s = 0.001_jprb

rfraer = 0.5_jprb
rfrgas = 1.0_jprb

!*      2.1   SEA SALT
!             -------- 
!-- parameters related to SEA SALT: 12 relates to 12 values of relative humidity

nbrh=12
rrhmax = 95._jprb
rrhtab = (/ 0._jprb, 10._jprb, 20._jprb, 30._jprb, 40._jprb, 50._jprb &
       & , 60._jprb, 70._jprb, 80._jprb, 85._jprb, 90._jprb, 95._jprb /)
rrho_ss = (/ 2160._jprb, 2160._jprb, 2160._jprb, 2160._jprb, 1451.6_jprb &
     & , 1367.9_jprb, 1302.9_jprb, 1243.2_jprb, 1182.7_jprb, 1149.5_jprb &
     & , 1111.6_jprb, 1063.1_jprb /)
rssgrow = (/ 0.503_jprb, 0.503_jprb, 0.503_jprb, 0.503_jprb, 0.724_jprb &
     & , 0.782_jprb, 0.838_jprb, 0.905_jprb, 1.000_jprb, 1.072_jprb &
     & , 1.188_jprb, 1.447_jprb /)

!-- OB's original 10 bins !RMMD_SS = (/ 0.09_JPRB, 0.19_JPRB, 0.38_JPRB, 0.75_JPRB, 1.50_JPRB &
!         & , 3.00_JPRB, 7.00_JPRB, 15.0_JPRB, 25.0_JPRB, 35.0_JPRB /)

!-- OB's ECMWF 3 bins of sea salt
!  bins are 0.03 - 0.5 -  5.0  - 20 microns

rmmd_ss = (/ 0.30_jprb, 3.00_jprb, 10.00_jprb /)
rfrss   = (/  0.7_jprb,  0.7_jprb,   0.7_jprb /)
rho_wat = 1000._jprb
rho_ice = 500._jprb

!-  computed off-line by gems_ss.f  (m s-1)

rvsedoce(1:3) = (/   0.24e-04_jprb,   0.20e-02_jprb,   0.20e-01_jprb /)
rvsedsic(1:3) = (/   0.24e-04_jprb,   0.20e-02_jprb,   0.20e-01_jprb /)
rvsedlnd(1:3) = (/   0.24e-04_jprb,   0.20e-02_jprb,   0.20e-01_jprb /)
rvsedlic(1:3) = (/   0.24e-04_jprb,   0.20e-02_jprb,   0.20e-01_jprb /)

! adapted from LMDZ   (m s-1)

rvdpoce(1:3) = (/ 0.100e-02_jprb, 0.110e-01_jprb,   0.145e-01_jprb /)
rvdpsic(1:3) = (/ 0.100e-02_jprb, 0.110e-01_jprb,   0.145e-01_jprb /)
rvdplnd(1:3) = (/ 0.100e-02_jprb, 0.110e-01_jprb,   0.145e-01_jprb /)
rvdplic(1:3) = (/ 0.100e-02_jprb, 0.110e-01_jprb,   0.145e-01_jprb /)


!*      2.2   DESERT DUST
!             ----------- 
!- parameters related to DESERT DUST  (OB's ECMWF 3 bins)
!  bins are 0.03 - 0.55 -  0.9  - 20 microns

rmmd_dd = (/  0.32_jprb,  0.75_jprb,   9.0_jprb /)
rrho_dd = (/ 2600._jprb, 2600._jprb, 2600._jprb /)
rfrdd   = (/   0.7_jprb,   0.7_jprb,   0.7_jprb /)

!-  computed off-line by gems_dust.f

rvsedoce(4:6) = (/ 0.70e-04_jprb, 0.20e-03_jprb, 0.20e-02_jprb /)
rvsedsic(4:6) = (/ 0.70e-04_jprb, 0.20e-03_jprb, 0.20e-02_jprb /)
rvsedlnd(4:6) = (/ 0.70e-04_jprb, 0.20e-03_jprb, 0.20e-02_jprb /)
rvsedlic(4:6) = (/ 0.70e-04_jprb, 0.20e-03_jprb, 0.20e-02_jprb /)

! adapted from LMDZ   (m s-1)

rvdpoce(4:6)  = (/ 0.100e-02_jprb, 0.500e-02_jprb, 0.110e-01_jprb /)
rvdpsic(4:6)  = (/ 0.100e-02_jprb, 0.500e-02_jprb, 0.110e-01_jprb /)
rvdplnd(4:6)  = (/ 0.100e-02_jprb, 0.500e-02_jprb, 0.110e-01_jprb /)
rvdplic(4:6)  = (/ 0.100e-02_jprb, 0.500e-02_jprb, 0.110e-01_jprb /)


!*      2.3   OTHER AEROSOLS (to be improved later!)
!             --------------
!- parameters related to other aerosol types
!- particulate organic matter POM
rfrom = (/  0.0_jprb,  0.7_jprb /)

rvsedoce(7:8) = (/   0.10e+00_jprb,   0.10e+00_jprb /)
rvsedsic(7:8) = (/   0.10e+00_jprb,   0.10e+00_jprb /)
rvsedlnd(7:8) = (/   0.10e+00_jprb,   0.10e+00_jprb /)
rvsedlic(7:8) = (/   0.10e+00_jprb,   0.10e+00_jprb /)

! adapted from LMDZ   (m s-1)

rvdpoce(7:8) = (/   0.10e-02_jprb,   0.10e-02_jprb /)
rvdpsic(7:8) = (/   0.10e-02_jprb,   0.10e-02_jprb /)
rvdplnd(7:8) = (/   0.10e-02_jprb,   0.10e-02_jprb /)
rvdplic(7:8) = (/   0.10e-02_jprb,   0.10e-02_jprb /)


!- black carbon
rfrbc = (/  0.0_jprb,  0.7_jprb /)

rvsedoce(9:10) = (/   0.10e+00_jprb,   0.10e+00_jprb /)
rvsedsic(9:10) = (/   0.10e+00_jprb,   0.10e+00_jprb /)
rvsedlnd(9:10) = (/   0.10e+00_jprb,   0.10e+00_jprb /)
rvsedlic(9:10) = (/   0.10e+00_jprb,   0.10e+00_jprb /)

! adapted from LMDZ   (m s-1)

rvdpoce(9:10) = (/   0.10e-02_jprb,   0.10e-02_jprb /)
rvdpsic(9:10) = (/   0.10e-02_jprb,   0.10e-02_jprb /)
rvdplnd(9:10) = (/   0.10e-02_jprb,   0.10e-02_jprb /)
rvdplic(9:10) = (/   0.10e-02_jprb,   0.10e-02_jprb /)

!- sulfate
rfrso4 = 0.7_jprb

rvsedoce(11) = 0.05_jprb
rvsedsic(11) = 0.25_jprb
rvsedlnd(11) = 0.25_jprb
rvsedlic(11) = 0.25_jprb

! adapted from LMDZ   (m s-1)

rvdpoce(11) = 0.05e-02_jprb
rvdpsic(11) = 0.25e-02_jprb
rvdplnd(11) = 0.25e-02_jprb
rvdplic(11) = 0.25e-02_jprb

!- fly ash
rfrif  = 0.7_jprb

rvsedoce(12) = 0.20e+00_jprb
rvsedsic(12) = 0.20e+00_jprb
rvsedlnd(12) = 0.20e+00_jprb
rvsedlic(12) = 0.20e+00_jprb

! adapted from LMDZ   (m s-1)

rvdpoce(12) = 0.20e-02_jprb
rvdpsic(12) = 0.20e-02_jprb
rvdplnd(12) = 0.20e-02_jprb
rvdplic(12) = 0.20e-02_jprb




!- NB: 15 values for all possible types of ECMWF aerosols
ralphar = (/ &
 & 0.001_jprb, 0.001_jprb, 0.001_jprb, 0.001_jprb, 0.001_jprb, 0.001_jprb, 0.001_jprb, 0.001_jprb, &
 & 0.001_jprb, 0.001_jprb, 0.001_jprb, 0.001_jprb, 0.001_jprb, 0.001_jprb, 0.001_jprb /)
ralphas = (/ &
 &  0.01_jprb,  0.01_jprb,  0.01_jprb,  0.01_jprb,  0.01_jprb,  0.01_jprb,  0.01_jprb,  0.01_jprb, &
 &  0.01_jprb,  0.01_jprb,  0.01_jprb,  0.01_jprb,  0.01_jprb,  0.01_jprb,  0.01_jprb /)

                                                                                                                                                                 
!*      3.    PARAMETERS RELATED TO TRANSPORT WITHIN THE FREE ATMOSPHERE
!             ----------------------------------------------------------

ndd1=4
nss1=1

rmfmin = 1.e-10_jprb

rmasse = (/ &
  &  6.02e+23_jprb, 6.02e+23_jprb, 6.02e+23_jprb, 6.02e+23_jprb, 6.02e+23_jprb, 6.02e+23_jprb &
  &, 6.02e+23_jprb, 6.02e+23_jprb, 6.02e+23_jprb, 6.02e+23_jprb, 6.02e+23_jprb, 6.02e+23_jprb &
  &, 6.02e+23_jprb, 6.02e+23_jprb, 6.02e+23_jprb /)


!     ----------------------------------------------------------------
IF (lhook) CALL dr_hook('SU_AERP',1,zhook_handle)
END SUBROUTINE su_aerp