LMDZ
susrtaer.F90
Go to the documentation of this file.
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 
yoesrtaer::rsrpiza
real(kind=jprb), dimension(14, 6) rsrpiza
Definition: yoesrtaer.F90:18
susrtaer
subroutine susrtaer
Definition: susrtaer.F90:2
yoesrtaer::rsrtaua
real(kind=jprb), dimension(14, 6) rsrtaua
Definition: yoesrtaer.F90:17
parkind1::jprb
integer, parameter jprb
Definition: parkind1.F90:31
yoesrtaer::rsrasya
real(kind=jprb), dimension(14, 6) rsrasya
Definition: yoesrtaer.F90:19
yomhook::lhook
logical lhook
Definition: yomhook.F90:12
yomhook::dr_hook
subroutine dr_hook(CDNAME, KSWITCH, PKEY)
Definition: yomhook.F90:17
parkind1::jpim
integer, parameter jpim
Definition: parkind1.F90:13