1 |
|
|
SUBROUTINE 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 |
|
|
|
45 |
|
|
USE PARKIND1 ,ONLY : JPIM ,JPRB |
46 |
|
|
USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
47 |
|
|
|
48 |
|
|
USE YOESRTAER , ONLY : RSRTAUA ,RSRPIZA ,RSRASYA |
49 |
|
|
|
50 |
|
|
! ---------------------------------------------------------------- |
51 |
|
|
|
52 |
|
|
IMPLICIT NONE |
53 |
|
|
|
54 |
|
|
! DUMMY INTEGER REAL |
55 |
|
|
REAL(KIND=JPRB) :: ZTAUA14(14,6),ZPIZA14(14,6),ZCGA14(14,6) |
56 |
|
|
|
57 |
|
|
INTEGER(KIND=JPIM) :: JAER, JNU |
58 |
|
|
REAL(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 |
|
|
|
110 |
|
|
IF (LHOOK) CALL DR_HOOK('SUSRTAER',0,ZHOOK_HANDLE) |
111 |
|
|
ZTAUA14( 1, :)= (/& |
112 |
|
|
& 0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /) |
113 |
|
|
ZTAUA14( 2, :)= (/& |
114 |
|
|
& 0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /) |
115 |
|
|
ZTAUA14( 3, :)= (/& |
116 |
|
|
& 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) |
117 |
|
|
ZTAUA14( 4, :)= (/& |
118 |
|
|
& 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) |
119 |
|
|
ZTAUA14( 5, :)= (/& |
120 |
|
|
& 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) |
121 |
|
|
ZTAUA14( 6, :)= (/& |
122 |
|
|
& 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) |
123 |
|
|
ZTAUA14( 7, :)= (/& |
124 |
|
|
& 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) |
125 |
|
|
ZTAUA14( 8, :)= (/& |
126 |
|
|
& 0.52838_JPRB , 0.93285_JPRB , 0.93449_JPRB , 0.53078_JPRB , 0.67148_JPRB , 0.46608_JPRB /) |
127 |
|
|
ZTAUA14( 9, :)= (/& |
128 |
|
|
& 0.52838_JPRB , 0.93285_JPRB , 0.93449_JPRB , 0.53078_JPRB , 0.67148_JPRB , 0.46608_JPRB /) |
129 |
|
|
ZTAUA14(10, :)= (/& |
130 |
|
|
& 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) |
131 |
|
|
ZTAUA14(11, :)= (/& |
132 |
|
|
& 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) |
133 |
|
|
ZTAUA14(12, :)= (/& |
134 |
|
|
& 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) |
135 |
|
|
ZTAUA14(13, :)= (/& |
136 |
|
|
& 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) |
137 |
|
|
ZTAUA14(14, :)= (/& |
138 |
|
|
& 0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /) |
139 |
|
|
|
140 |
|
|
ZPIZA14( 1, :)= (/& |
141 |
|
|
& .5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/) |
142 |
|
|
ZPIZA14( 2, :)= (/& |
143 |
|
|
& .5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/) |
144 |
|
|
ZPIZA14( 3, :)= (/& |
145 |
|
|
& .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) |
146 |
|
|
ZPIZA14( 4, :)= (/& |
147 |
|
|
& .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) |
148 |
|
|
ZPIZA14( 5, :)= (/& |
149 |
|
|
& .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) |
150 |
|
|
ZPIZA14( 6, :)= (/& |
151 |
|
|
& .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) |
152 |
|
|
ZPIZA14( 7, :)= (/& |
153 |
|
|
& .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) |
154 |
|
|
ZPIZA14( 8, :)= (/& |
155 |
|
|
& .8970131_JPRB, .9984940_JPRB, .9245594_JPRB, .7768385_JPRB, .9532763_JPRB, .9999999_JPRB/) |
156 |
|
|
ZPIZA14( 9, :)= (/& |
157 |
|
|
& .8970131_JPRB, .9984940_JPRB, .9245594_JPRB, .7768385_JPRB, .9532763_JPRB, .9999999_JPRB/) |
158 |
|
|
ZPIZA14(10, :)= (/& |
159 |
|
|
& .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) |
160 |
|
|
ZPIZA14(11, :)= (/& |
161 |
|
|
& .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) |
162 |
|
|
ZPIZA14(12, :)= (/& |
163 |
|
|
& .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) |
164 |
|
|
ZPIZA14(13, :)= (/& |
165 |
|
|
& .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) |
166 |
|
|
ZPIZA14(14, :)= (/& |
167 |
|
|
& .5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/) |
168 |
|
|
|
169 |
|
|
ZCGA14( 1, :)= (/& |
170 |
|
|
& 0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/) |
171 |
|
|
ZCGA14( 2, :)= (/& |
172 |
|
|
& 0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/) |
173 |
|
|
ZCGA14( 3, :)= (/& |
174 |
|
|
& 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) |
175 |
|
|
ZCGA14( 4, :)= (/& |
176 |
|
|
& 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) |
177 |
|
|
ZCGA14( 5, :)= (/& |
178 |
|
|
& 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) |
179 |
|
|
ZCGA14( 6, :)= (/& |
180 |
|
|
& 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) |
181 |
|
|
ZCGA14( 7, :)= (/& |
182 |
|
|
& 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) |
183 |
|
|
ZCGA14( 8, :)= (/& |
184 |
|
|
& 0.668431_JPRB, 0.788530_JPRB, 0.698682_JPRB, 0.657422_JPRB, .6735182_JPRB, .6519706_JPRB/) |
185 |
|
|
ZCGA14( 9, :)= (/& |
186 |
|
|
& 0.668431_JPRB, 0.788530_JPRB, 0.698682_JPRB, 0.657422_JPRB, .6735182_JPRB, .6519706_JPRB/) |
187 |
|
|
ZCGA14(10, :)= (/& |
188 |
|
|
& 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) |
189 |
|
|
ZCGA14(11, :)= (/& |
190 |
|
|
& 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) |
191 |
|
|
ZCGA14(12, :)= (/& |
192 |
|
|
& 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) |
193 |
|
|
ZCGA14(13, :)= (/& |
194 |
|
|
& 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) |
195 |
|
|
ZCGA14(14, :)= (/& |
196 |
|
|
& 0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/) |
197 |
|
|
|
198 |
|
|
! ---------------------------------------------------------------- |
199 |
|
|
|
200 |
|
|
DO 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 |
206 |
|
|
ENDDO |
207 |
|
|
|
208 |
|
|
! ---------------------------------------------------------------- |
209 |
|
|
|
210 |
|
|
IF (LHOOK) CALL DR_HOOK('SUSRTAER',1,ZHOOK_HANDLE) |
211 |
|
|
END SUBROUTINE SUSRTAER |
212 |
|
|
|