suaer15.F90

Go to the documentation of this file.

SUBROUTINE suaer15

!**** *SUAER15*   - INITIALIZE COMMON YOMAER15

!     PURPOSE.
!     --------
!           INITIALIZE YOMAER15, THE COMMON THAT CONTAINS THE
!           RADIATIVE CHARACTERISTICS OF THE AEROSOLS

!**   INTERFACE.
!     ----------
!              -----        -----

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

!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMON YOMAER15

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

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

!     REFERENCE.
!     ----------

!     04-06: F. Bouyssel. Meteo-France
!     05-09: A. Alias - PAESOD (black_carbon) is added to PAELAN (oraganic)
!                     - the sulfate is put in a separate type 
!                     (see Hu Ron Ming work) - P.Marquet

!=======================================================================
!-- The (old) five aerosol types were respectively:

!  1/ continental average (+desert)       2/ maritime
!  3/ urban                               4/ volcanic active
!  5/ stratospheric background

!=======================================================================

!-- The (new) six aerosol types are respectively:

!  1/ continental average                 2/ maritime
!  3/ desert                              4/ volcanic active
!  5/ stratospheric background            6/ sulfate

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

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

USE yomrad15 , ONLY : lnewaer15
USE yomaer15 , ONLY : raer15    ,taua15   ,rpiza15  ,rcga15

IMPLICIT NONE

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

REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook('SUAER15',0,zhook_handle)
IF (.NOT.lnewaer15) THEN
  
  raer15=reshape((/&
   & .038520_jprb, .037196_jprb, .040532_jprb, .054934_jprb, .038520_jprb ,&
   & .12613_jprb , .18313_jprb , .10357_jprb , .064106_jprb, .126130_jprb ,&
   & .012579_jprb, .013649_jprb, .018652_jprb, .025181_jprb, .012579_jprb ,&
   & .011890_jprb, .016142_jprb, .021105_jprb, .028908_jprb, .011890_jprb ,&
   & .013792_jprb, .026810_jprb, .052203_jprb, .066338_jprb, .013792_jprb ,&
   & .012579_jprb, .013649_jprb, .018652_jprb, .025181_jprb, .012579_jprb /)&
   & ,shape=(/5,6/))  

  taua15=reshape((/&
   & .730719_jprb,.730719_jprb,&
   & .912819_jprb,.912819_jprb,&
   & .725059_jprb,.725059_jprb,&
   & .745405_jprb,.745405_jprb,&
   & .682188_jprb,.682188_jprb,&
   & 1.35059_jprb,.725059_jprb /)&
   & ,shape=(/2,6/))  

  rpiza15=reshape((/&
   & .872212_jprb,.872212_jprb,&
   & .982545_jprb,.982545_jprb,&
   & .623143_jprb,.623143_jprb,&
   & .944887_jprb,.944887_jprb,&
   & .997975_jprb,.997975_jprb,&
   & .999999_jprb,.999999_jprb /)&
   & ,shape=(/2,6/))  

  rcga15=reshape((/&
   & .647596_jprb,.647596_jprb,&
   & .739002_jprb,.739002_jprb,&
   & .580845_jprb,.580845_jprb,&
   & .662657_jprb,.662657_jprb,&
   & .624246_jprb,.624246_jprb,&
   & .680845_jprb,.680845_jprb /)&
   & ,shape=(/2,6/))  

ELSE

  raer15( :, 1)= (/&
   & .036271_jprb, .030153_jprb, .017343_jprb, .015002_jprb, .008806_jprb /)  
  raer15( :, 2)= (/&
   & .026561_jprb, .032657_jprb, .017977_jprb, .014210_jprb, .016775_jprb /)  
  raer15( :, 3)= (/&
   & .014897_jprb, .016359_jprb, .019789_jprb, .030777_jprb, .013341_jprb /)  
  raer15( :, 4)= (/&
   & .011890_jprb, .016142_jprb, .021105_jprb, .028908_jprb, .011890_jprb /)  
  raer15( :, 5)= (/&
   & .013792_jprb, .026810_jprb, .052203_jprb, .066338_jprb, .013792_jprb /)  
  raer15( :, 6)= (/&
   & .012579_jprb, .013649_jprb, .018652_jprb, .025181_jprb, .012579_jprb /)

  taua15(1, :)= (/&
   & 1.69446_jprb, 1.11855_jprb, 1.09212_jprb, 1.03858_jprb, 1.12044_jprb, 1.35059_jprb  /)  
  taua15(2, :)= (/&
   & 0.40174_jprb, 0.89383_jprb, 0.89546_jprb, 0.51143_jprb, 0.32646_jprb, .725059_jprb /)  
 
  rpiza15(1, :)= (/&
   & .9148907_jprb, .9956173_jprb, .7504584_jprb, .9401905_jprb, .9999999_jprb, .999999_jprb/)  
  rpiza15(2, :)= (/&
   & .8814597_jprb, .9920407_jprb, .9239428_jprb, .9515548_jprb, .9938563_jprb, .999999_jprb/)  
 
  rcga15(1, :)= (/&
   & 0.729019_jprb, 0.803129_jprb, 0.784592_jprb, .7008249_jprb, .7270548_jprb, .680845_jprb/)  
  rcga15(2, :)= (/&
   & 0.663224_jprb, 0.793746_jprb, 0.696315_jprb, .6608509_jprb, .6318786_jprb, .680845_jprb/)  

ENDIF

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

IF (lhook) CALL dr_hook('SUAER15',1,zhook_handle)
END SUBROUTINE suaer15