source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/suaerl.F90 @ 3331

Last change on this file since 3331 was 3331, checked in by acozic, 6 years ago

Add modification for isotopes

  • Property svn:executable set to *
File size: 3.6 KB
Line 
1SUBROUTINE SUAERL
2
3!**** *SUAERL*   - INITIALIZE COMMON YOEAER
4
5!     PURPOSE.
6!     --------
7!           INITIALIZE YOEAER, THE COMMON THAT CONTAINS THE
8!           RADIATIVE CHARACTERISTICS OF THE AEROSOLS
9
10!**   INTERFACE.
11!     ----------
12!              -----        -----
13
14!        EXPLICIT ARGUMENTS :
15!        --------------------
16!        NONE
17
18!        IMPLICIT ARGUMENTS :
19!        --------------------
20!        COMMON YOEAER
21
22!     METHOD.
23!     -------
24!        SEE DOCUMENTATION
25
26!     EXTERNALS.
27!     ----------
28
29!     REFERENCE.
30!     ----------
31!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE "IFS MODEL"
32
33!     AUTHOR.
34!     -------
35!        JEAN-JACQUES MORCRETTE *ECMWF*
36
37!     MODIFICATIONS.
38!     --------------
39!        ORIGINAL : 88-02-15
40!        99-05-25   JJMorcrette   Revised Aerosol LW properties (OPAC)
41!        M.Hamrud      01-Oct-2003 CY28 Cleaning
42
43!     ------------------------------------------------------------------
44
45USE PARKIND1  ,ONLY : JPIM     ,JPRB
46USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
47
48USE YOESW    , ONLY : RAER
49
50!      ----------------------------------------------------------------
51
52!*       1.    SHORTWAVE COEFFICIENTS
53!              ----------------------
54
55!      Initialised in SUAERSn depending on NSW number of SW spectral
56!      intervals
57
58!      ----------------------------------------------------------------
59
60!*       2.    LONGWAVE COEFFICIENTS
61!              ---------------------
62
63!=======================================================================
64!-- The (old) five aerosol types were respectively:
65
66!  1/ continental average (+desert)       2/ maritime
67!  3/ urban                               4/ volcanic active
68!  5/ stratospheric background
69
70!      RAER  = RESHAPE((/
71!     &          .038520, .037196, .040532, .054934, .038520
72!     &        , .126130, .18313 , .10357 , .064106, .126130
73!     &        , .012579, .013649, .018652, .025181, .012579
74!     &        , .011890, .016142, .021105, .028908, .011890
75!     &        , .013792, .026810, .052203, .066338, .013792 /)
76!     & ,SHAPE=(/5,5/))
77
78!=======================================================================
79
80!-- The six aerosol types are respectively:
81
82!  1/ continental average                 2/ maritime
83!  3/ desert                              4/ urban
84!  5/ volcanic active                     6/ stratospheric background
85
86! The quantities given are:
87! TAU : ratio of average optical thickness in interval to that at 0.55
88!       micron
89! PIZA: average single scattering albedo
90! CGA : average asymmetry factor
91
92! computed from Hess and Koepke (con, mar, des, urb)
93!          from Bonnel et al.   (vol, str)
94
95!-- data are entered for the 6 spectral intervals of the LW scheme (line)
96!   and the different types of aerosols (column)
97
98IMPLICIT NONE
99REAL(KIND=JPRB) :: ZHOOK_HANDLE
100IF (LHOOK) CALL DR_HOOK('SUAERL',0,ZHOOK_HANDLE)
101RAER( :, 1)= (/&
102 & .036271_JPRB, .030153_JPRB, .017343_JPRB, .015002_JPRB, .008806_JPRB, .006865_JPRB /) 
103RAER( :, 2)= (/&
104 & .026561_JPRB, .032657_JPRB, .017977_JPRB, .014210_JPRB, .016775_JPRB, .022123_JPRB /) 
105RAER( :, 3)= (/&
106 & .014897_JPRB, .016359_JPRB, .019789_JPRB, .030777_JPRB, .013341_JPRB, .014321_JPRB /) 
107RAER( :, 4)= (/&
108 & .001863_JPRB, .002816_JPRB, .002355_JPRB, .002557_JPRB, .001774_JPRB, .001780_JPRB /) 
109RAER( :, 5)= (/&
110 & .011890_JPRB, .016142_JPRB, .021105_JPRB, .028908_JPRB, .011890_JPRB, .011890_JPRB /) 
111RAER( :, 6)= (/&
112 & .013792_JPRB, .026810_JPRB, .052203_JPRB, .066338_JPRB, .013792_JPRB, .013792_JPRB /) 
113
114!      ----------------------------------------------------------------
115
116IF (LHOOK) CALL DR_HOOK('SUAERL',1,ZHOOK_HANDLE)
117END SUBROUTINE SUAERL
Note: See TracBrowser for help on using the repository browser.