source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/susrtaer.F90 @ 5308

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

Add modification for isotopes

  • Property svn:executable set to *
File size: 8.0 KB
Line 
1SUBROUTINE SUSRTAER
2
3!**** *SUSRTAER*   - INITIALIZE COMMON YOESRTAER
4
5!     PURPOSE.
6!     --------
7!           INITIALIZE YOESRTAER, 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 YOESRTAER
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 : 03-03-06
40!        03-03-06  JJMorcrette   14 spectral intervals for RRTM_SW
41!        M.Hamrud      01-Oct-2003 CY28 Cleaning
42
43!     ------------------------------------------------------------------
44
45USE PARKIND1  ,ONLY : JPIM     ,JPRB
46USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
47
48USE YOESRTAER   , ONLY : RSRTAUA     ,RSRPIZA    ,RSRASYA
49
50!      ----------------------------------------------------------------
51
52IMPLICIT NONE
53
54!     DUMMY INTEGER REAL
55REAL(KIND=JPRB) :: ZTAUA14(14,6),ZPIZA14(14,6),ZCGA14(14,6)
56
57INTEGER(KIND=JPIM) :: JAER, JNU
58REAL(KIND=JPRB) :: ZHOOK_HANDLE
59
60!      ----------------------------------------------------------------
61
62!*       1.    SHORTWAVE COEFFICIENTS
63!              ----------------------
64!=======================================================================
65!-- The (old) five aerosol types were respectively:
66
67!  1/ continental average (+desert)       2/ maritime
68!  3/ urban                               4/ volcanic active
69!  5/ stratospheric background
70
71!-- old values were not spectrally defined:
72! ZTAU2  = .730719, .912819, .725059, .745405, .682188
73! ZPIZA2 = .872212, .982545, .623143, .944887, .997975
74! ZCGA2  = .647596, .739002, .580845, .662657, .624246
75!=======================================================================
76
77!-- The six aerosol types are respectively:
78
79!  1/ continental average                 2/ maritime
80!  3/ desert                              4/ urban
81!  5/ volcanic active                     6/ stratospheric background
82
83! The quantities given are:
84! TAU : ratio of average optical thickness in interval to that at 0.55
85!       micron
86! PIZA: average single scattering albedo
87! CGA : average asymmetry factor
88
89! computed from Hess and Koepke (con, mar, des, urb)
90!          from Bonnel et al.   (vol, str)
91
92!       1.1    FOURTEEN SPECTRAL INTERVALS FOR RRTM_SW
93
94! SW : 14 spectral intervals
95!  3.846 -  3.077
96!  3.077 -  2.500
97!  2.500 -  2.150
98!  2.150 -  1.942
99!  1.942 -  1.626
100!  1.626 -  1.299
101!  1.299 -  1.242
102!  1.242 -  0.7782
103!  0.7782-  0.6250
104!  0.6250-  0.4415
105!  0.4415-  0.3448
106!  0.3448-  0.2632
107!  0.2632-  0.2000
108! 12.195 -  3.846
109
110IF (LHOOK) CALL DR_HOOK('SUSRTAER',0,ZHOOK_HANDLE)
111ZTAUA14( 1, :)= (/&
112 & 0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /) 
113ZTAUA14( 2, :)= (/&
114 & 0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /) 
115ZTAUA14( 3, :)= (/&
116 & 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) 
117ZTAUA14( 4, :)= (/&
118 & 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) 
119ZTAUA14( 5, :)= (/&
120 & 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) 
121ZTAUA14( 6, :)= (/&
122 & 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) 
123ZTAUA14( 7, :)= (/&
124 & 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) 
125ZTAUA14( 8, :)= (/&
126 & 0.52838_JPRB , 0.93285_JPRB , 0.93449_JPRB , 0.53078_JPRB , 0.67148_JPRB , 0.46608_JPRB /) 
127ZTAUA14( 9, :)= (/&
128 & 0.52838_JPRB , 0.93285_JPRB , 0.93449_JPRB , 0.53078_JPRB , 0.67148_JPRB , 0.46608_JPRB /) 
129ZTAUA14(10, :)= (/&
130 & 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) 
131ZTAUA14(11, :)= (/&
132 & 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) 
133ZTAUA14(12, :)= (/&
134 & 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) 
135ZTAUA14(13, :)= (/&
136 & 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) 
137ZTAUA14(14, :)= (/&
138 & 0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /) 
139 
140ZPIZA14( 1, :)= (/&
141 & .5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/) 
142ZPIZA14( 2, :)= (/&
143 & .5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/) 
144ZPIZA14( 3, :)= (/&
145 & .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) 
146ZPIZA14( 4, :)= (/&
147 & .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) 
148ZPIZA14( 5, :)= (/&
149 & .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) 
150ZPIZA14( 6, :)= (/&
151 & .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) 
152ZPIZA14( 7, :)= (/&
153 & .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) 
154ZPIZA14( 8, :)= (/&
155 & .8970131_JPRB, .9984940_JPRB, .9245594_JPRB, .7768385_JPRB, .9532763_JPRB, .9999999_JPRB/) 
156ZPIZA14( 9, :)= (/&
157 & .8970131_JPRB, .9984940_JPRB, .9245594_JPRB, .7768385_JPRB, .9532763_JPRB, .9999999_JPRB/) 
158ZPIZA14(10, :)= (/&
159 & .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) 
160ZPIZA14(11, :)= (/&
161 & .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) 
162ZPIZA14(12, :)= (/&
163 & .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) 
164ZPIZA14(13, :)= (/&
165 & .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) 
166ZPIZA14(14, :)= (/&
167 & .5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/) 
168
169ZCGA14( 1, :)= (/&
170 & 0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/) 
171ZCGA14( 2, :)= (/&
172 & 0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/) 
173ZCGA14( 3, :)= (/&
174 & 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) 
175ZCGA14( 4, :)= (/&
176 & 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) 
177ZCGA14( 5, :)= (/&
178 & 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) 
179ZCGA14( 6, :)= (/&
180 & 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) 
181ZCGA14( 7, :)= (/&
182 & 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) 
183ZCGA14( 8, :)= (/&
184 & 0.668431_JPRB, 0.788530_JPRB, 0.698682_JPRB, 0.657422_JPRB, .6735182_JPRB, .6519706_JPRB/) 
185ZCGA14( 9, :)= (/&
186 & 0.668431_JPRB, 0.788530_JPRB, 0.698682_JPRB, 0.657422_JPRB, .6735182_JPRB, .6519706_JPRB/) 
187ZCGA14(10, :)= (/&
188 & 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) 
189ZCGA14(11, :)= (/&
190 & 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) 
191ZCGA14(12, :)= (/&
192 & 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) 
193ZCGA14(13, :)= (/&
194 & 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) 
195ZCGA14(14, :)= (/&
196 & 0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/) 
197
198!      ----------------------------------------------------------------
199
200DO JNU=1,14
201  DO JAER=1,6
202    RSRTAUA(JNU,JAER)=ZTAUA14(JNU,JAER)
203    RSRPIZA(JNU,JAER)=ZPIZA14(JNU,JAER)
204    RSRASYA(JNU,JAER)=ZCGA14 (JNU,JAER)
205  ENDDO
206ENDDO
207
208!      ----------------------------------------------------------------
209
210IF (LHOOK) CALL DR_HOOK('SUSRTAER',1,ZHOOK_HANDLE)
211END SUBROUTINE SUSRTAER
212
Note: See TracBrowser for help on using the repository browser.