source: LMDZ5/branches/testing/libf/phymar/suaersn.F90 @ 5431

SUBROUTINE SUAERSN (KTSW,KSW)

!**** *SUAERS*   - INITIALIZE COMMON YOEAER

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

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

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

!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMON YOEAER

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

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

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

!     AUTHOR.
!     -------
!        JEAN-JACQUES MORCRETTE *ECMWF*

!     MODIFICATIONS.
!     --------------
!        ORIGINAL : 88-02-15
!        96-01-27  JJ Morcrette  Various spectral resolutions
!        99-05-25  JJMorcrette   Revised aerosol optical properties
!        00-10-25  JJMorcrette   6 spectral intervals

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

#include "tsmbkind.h"

USE YOESW    , ONLY : RTAUA     ,RPIZA    ,RCGA

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

IMPLICIT NONE


!     DUMMY INTEGER SCALARS
INTEGER_M :: KSW
INTEGER_M :: KTSW

REAL_B :: ZTAUA2(2,6)  ,ZPIZA2(2,6)  ,ZCGA2(2,6)
REAL_B :: ZTAUA4(4,6)  ,ZPIZA4(4,6)  ,ZCGA4(4,6)
REAL_B :: ZTAUA6(6,6)  ,ZPIZA6(6,6)  ,ZCGA6(6,6)

!     LOCAL INTEGER SCALARS
INTEGER_M :: JAER, JNU

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

!*       1.    SHORTWAVE COEFFICIENTS
!              ----------------------
!=======================================================================
!-- The (old) five aerosol types were respectively:

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

!-- old values were not spectrally defined:
! ZTAU2  = .730719, .912819, .725059, .745405, .682188
! ZPIZA2 = .872212, .982545, .623143, .944887, .997975
! ZCGA2  = .647596, .739002, .580845, .662657, .624246
!=======================================================================

!-- The six aerosol types are respectively:

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

! The quantities given are:
! TAU : ratio of average optical thickness in interval to that at 0.55 
!       micron
! PIZA: average single scattering albedo
! CGA : average asymmetry factor

! computed from Hess and Koepke (con, mar, des, urb)
!          from Bonnel et al.   (vol, str)


!        1.1   TWO SPECTRAL INTERVALS (0.25-0.69-4.00microns)

ZTAUA2(1, :)= (/&
 &1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /)
ZTAUA2(2, :)= (/&
 &0.40174_JPRB , 0.89383_JPRB , 0.89546_JPRB , 0.40741_JPRB , 0.51143_JPRB , 0.32646_JPRB /)
 
ZPIZA2(1, :)= (/&
 &.9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/)
ZPIZA2(2, :)= (/&
 &.8814597_JPRB, .9920407_JPRB, .9239428_JPRB, .7546879_JPRB, .9515548_JPRB, .9938563_JPRB/)
 
ZCGA2(1, :)= (/&
 &0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/)
ZCGA2(2, :)= (/&
 &0.663224_JPRB, 0.793746_JPRB, 0.696315_JPRB, 0.652612_JPRB, .6608509_JPRB, .6318786_JPRB/)


!        1.2   FOUR SPECTRAL INTERVALS (0.25-0.69-1.19-2.38-4.00microns)

ZTAUA4(1, :)= (/&
 &1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /)
ZTAUA4(2, :)= (/&
 &0.52838_JPRB , 0.93285_JPRB , 0.93449_JPRB , 0.53078_JPRB , 0.67148_JPRB , 0.46608_JPRB /)
ZTAUA4(3, :)= (/&
 &0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /)
ZTAUA4(4, :)= (/&
 &0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /)
 
ZPIZA4(1, :)= (/&
 &.9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/)
ZPIZA4(2, :)= (/&
 &.8970131_JPRB, .9984940_JPRB, .9245594_JPRB, .7768385_JPRB, .9532763_JPRB, .9999999_JPRB/)
ZPIZA4(3, :)= (/&
 &.8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/)
ZPIZA4(4, :)= (/&
 &.5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/)
 
ZCGA4(1, :)= (/&
 &0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/)
ZCGA4(2, :)= (/&
 &0.668431_JPRB, 0.788530_JPRB, 0.698682_JPRB, 0.657422_JPRB, .6735182_JPRB, .6519706_JPRB/)
ZCGA4(3, :)= (/&
 &0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/)
ZCGA4(4, :)= (/&
 &0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/)


!        1.3   SIX SPECTRAL INTERVALS (0.185-0.25-0.44-0.69-1.19-2.38-4.00microns)

ZTAUA6(1, :)= (/&
 &1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /)
ZTAUA6(2, :)= (/&
 &1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /)
ZTAUA6(3, :)= (/&
 &1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /)
ZTAUA6(4, :)= (/&
 &0.52838_JPRB , 0.93285_JPRB , 0.93449_JPRB , 0.53078_JPRB , 0.67148_JPRB , 0.46608_JPRB /)
ZTAUA6(5, :)= (/&
 &0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /)
ZTAUA6(6, :)= (/&
 &0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /)
 
ZPIZA6(1, :)= (/&
 &.9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/)
ZPIZA6(2, :)= (/&
 &.9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/)
ZPIZA6(3, :)= (/&
 &.9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/)
ZPIZA6(4, :)= (/&
 &.8970131_JPRB, .9984940_JPRB, .9245594_JPRB, .7768385_JPRB, .9532763_JPRB, .9999999_JPRB/)
ZPIZA6(5, :)= (/&
 &.8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/)
ZPIZA6(6, :)= (/&
 &.5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/)
 
ZCGA6(1, :)= (/&
 &0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/)
ZCGA6(2, :)= (/&
 &0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/)
ZCGA6(3, :)= (/&
 &0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/)
ZCGA6(4, :)= (/&
 &0.668431_JPRB, 0.788530_JPRB, 0.698682_JPRB, 0.657422_JPRB, .6735182_JPRB, .6519706_JPRB/)
ZCGA6(5, :)= (/&
 &0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/)
ZCGA6(6, :)= (/&
 &0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/)


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

IF (KSW == 2) THEN
  DO JNU=1,KSW
    DO JAER=1,6
      RTAUA(JNU,JAER)=ZTAUA2(JNU,JAER)
      RPIZA(JNU,JAER)=ZPIZA2(JNU,JAER)
      RCGA(JNU,JAER) =ZCGA2 (JNU,JAER)
    ENDDO
  ENDDO
ELSEIF (KSW == 4) THEN
  DO JNU=1,KSW
    DO JAER=1,6
      RTAUA(JNU,JAER)=ZTAUA4(JNU,JAER)
      RPIZA(JNU,JAER)=ZPIZA4(JNU,JAER)
      RCGA(JNU,JAER) =ZCGA4 (JNU,JAER)
    ENDDO
  ENDDO
ELSEIF (KSW == 6) THEN
  DO JNU=1,KSW
    DO JAER=1,6
      RTAUA(JNU,JAER)=ZTAUA6(JNU,JAER)
      RPIZA(JNU,JAER)=ZPIZA6(JNU,JAER)
      RCGA(JNU,JAER) =ZCGA6 (JNU,JAER)
    ENDDO
  ENDDO
ELSE
  STOP 'SUAERSN: WRONG NUMBER OF SPECTRAL INTERVALS'
ENDIF

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

RETURN
END SUBROUTINE SUAERSN