LMDZ
suaersn.F90
Go to the documentation of this file.
1 SUBROUTINE suaersn (KTSW,KSW)
2 
3 !**** *SUAERS* - 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 ! 96-01-27 JJ Morcrette Various spectral resolutions
41 ! 99-05-25 JJMorcrette Revised aerosol optical properties
42 ! 00-10-25 JJMorcrette 6 spectral intervals
43 
44 ! ------------------------------------------------------------------
45 
46 #include "tsmbkind.h"
47 
48 USE yoesw , ONLY : rtaua ,rpiza ,rcga
49 
50 ! ----------------------------------------------------------------
51 
52 IMPLICIT NONE
53 
54 
55 ! DUMMY INTEGER SCALARS
56 integer_m :: ksw
57 integer_m :: ktsw
58 
59 real_b :: ztaua2(2,6) ,zpiza2(2,6) ,zcga2(2,6)
60 real_b :: ztaua4(4,6) ,zpiza4(4,6) ,zcga4(4,6)
61 real_b :: ztaua6(6,6) ,zpiza6(6,6) ,zcga6(6,6)
62 
63 ! LOCAL INTEGER SCALARS
64 integer_m :: jaer, jnu
65 
66 ! ----------------------------------------------------------------
67 
68 !* 1. SHORTWAVE COEFFICIENTS
69 ! ----------------------
70 !=======================================================================
71 !-- The (old) five aerosol types were respectively:
72 
73 ! 1/ continental average (+desert) 2/ maritime
74 ! 3/ urban 4/ volcanic active
75 ! 5/ stratospheric background
76 
77 !-- old values were not spectrally defined:
78 ! ZTAU2 = .730719, .912819, .725059, .745405, .682188
79 ! ZPIZA2 = .872212, .982545, .623143, .944887, .997975
80 ! ZCGA2 = .647596, .739002, .580845, .662657, .624246
81 !=======================================================================
82 
83 !-- The six aerosol types are respectively:
84 
85 ! 1/ continental average 2/ maritime
86 ! 3/ desert 4/ urban
87 ! 5/ volcanic active 6/ stratospheric background
88 
89 ! The quantities given are:
90 ! TAU : ratio of average optical thickness in interval to that at 0.55
91 ! micron
92 ! PIZA: average single scattering albedo
93 ! CGA : average asymmetry factor
94 
95 ! computed from Hess and Koepke (con, mar, des, urb)
96 ! from Bonnel et al. (vol, str)
97 
98 
99 ! 1.1 TWO SPECTRAL INTERVALS (0.25-0.69-4.00microns)
100 
101 ztaua2(1, :)= (/&
102  &1.69446_jprb , 1.11855_jprb , 1.09212_jprb , 1.72145_jprb , 1.03858_jprb , 1.12044_jprb /)
103 ztaua2(2, :)= (/&
104  &0.40174_jprb , 0.89383_jprb , 0.89546_jprb , 0.40741_jprb , 0.51143_jprb , 0.32646_jprb /)
105 
106 zpiza2(1, :)= (/&
107  &.9148907_jprb, .9956173_jprb, .7504584_jprb, .8131335_jprb, .9401905_jprb, .9999999_jprb/)
108 zpiza2(2, :)= (/&
109  &.8814597_jprb, .9920407_jprb, .9239428_jprb, .7546879_jprb, .9515548_jprb, .9938563_jprb/)
110 
111 zcga2(1, :)= (/&
112  &0.729019_jprb, 0.803129_jprb, 0.784592_jprb, 0.712208_jprb, .7008249_jprb, .7270548_jprb/)
113 zcga2(2, :)= (/&
114  &0.663224_jprb, 0.793746_jprb, 0.696315_jprb, 0.652612_jprb, .6608509_jprb, .6318786_jprb/)
115 
116 
117 ! 1.2 FOUR SPECTRAL INTERVALS (0.25-0.69-1.19-2.38-4.00microns)
118 
119 ztaua4(1, :)= (/&
120  &1.69446_jprb , 1.11855_jprb , 1.09212_jprb , 1.72145_jprb , 1.03858_jprb , 1.12044_jprb /)
121 ztaua4(2, :)= (/&
122  &0.52838_jprb , 0.93285_jprb , 0.93449_jprb , 0.53078_jprb , 0.67148_jprb , 0.46608_jprb /)
123 ztaua4(3, :)= (/&
124  &0.20543_jprb , 0.84642_jprb , 0.84958_jprb , 0.21673_jprb , 0.28270_jprb , 0.10915_jprb /)
125 ztaua4(4, :)= (/&
126  &0.10849_jprb , 0.66699_jprb , 0.65255_jprb , 0.11600_jprb , 0.06529_jprb , 0.04468_jprb /)
127 
128 zpiza4(1, :)= (/&
129  &.9148907_jprb, .9956173_jprb, .7504584_jprb, .8131335_jprb, .9401905_jprb, .9999999_jprb/)
130 zpiza4(2, :)= (/&
131  &.8970131_jprb, .9984940_jprb, .9245594_jprb, .7768385_jprb, .9532763_jprb, .9999999_jprb/)
132 zpiza4(3, :)= (/&
133  &.8287144_jprb, .9949396_jprb, .9279543_jprb, .6765051_jprb, .9467578_jprb, .9955938_jprb/)
134 zpiza4(4, :)= (/&
135  &.5230504_jprb, .7868518_jprb, .8531531_jprb, .4048149_jprb, .8748231_jprb, .2355667_jprb/)
136 
137 zcga4(1, :)= (/&
138  &0.729019_jprb, 0.803129_jprb, 0.784592_jprb, 0.712208_jprb, .7008249_jprb, .7270548_jprb/)
139 zcga4(2, :)= (/&
140  &0.668431_jprb, 0.788530_jprb, 0.698682_jprb, 0.657422_jprb, .6735182_jprb, .6519706_jprb/)
141 zcga4(3, :)= (/&
142  &0.636342_jprb, 0.802467_jprb, 0.691305_jprb, 0.627497_jprb, .6105750_jprb, .4760794_jprb/)
143 zcga4(4, :)= (/&
144  &0.700610_jprb, 0.818871_jprb, 0.702399_jprb, 0.689886_jprb, .4629866_jprb, .1907639_jprb/)
145 
146 
147 ! 1.3 SIX SPECTRAL INTERVALS (0.185-0.25-0.44-0.69-1.19-2.38-4.00microns)
148 
149 ztaua6(1, :)= (/&
150  &1.69446_jprb , 1.11855_jprb , 1.09212_jprb , 1.72145_jprb , 1.03858_jprb , 1.12044_jprb /)
151 ztaua6(2, :)= (/&
152  &1.69446_jprb , 1.11855_jprb , 1.09212_jprb , 1.72145_jprb , 1.03858_jprb , 1.12044_jprb /)
153 ztaua6(3, :)= (/&
154  &1.69446_jprb , 1.11855_jprb , 1.09212_jprb , 1.72145_jprb , 1.03858_jprb , 1.12044_jprb /)
155 ztaua6(4, :)= (/&
156  &0.52838_jprb , 0.93285_jprb , 0.93449_jprb , 0.53078_jprb , 0.67148_jprb , 0.46608_jprb /)
157 ztaua6(5, :)= (/&
158  &0.20543_jprb , 0.84642_jprb , 0.84958_jprb , 0.21673_jprb , 0.28270_jprb , 0.10915_jprb /)
159 ztaua6(6, :)= (/&
160  &0.10849_jprb , 0.66699_jprb , 0.65255_jprb , 0.11600_jprb , 0.06529_jprb , 0.04468_jprb /)
161 
162 zpiza6(1, :)= (/&
163  &.9148907_jprb, .9956173_jprb, .7504584_jprb, .8131335_jprb, .9401905_jprb, .9999999_jprb/)
164 zpiza6(2, :)= (/&
165  &.9148907_jprb, .9956173_jprb, .7504584_jprb, .8131335_jprb, .9401905_jprb, .9999999_jprb/)
166 zpiza6(3, :)= (/&
167  &.9148907_jprb, .9956173_jprb, .7504584_jprb, .8131335_jprb, .9401905_jprb, .9999999_jprb/)
168 zpiza6(4, :)= (/&
169  &.8970131_jprb, .9984940_jprb, .9245594_jprb, .7768385_jprb, .9532763_jprb, .9999999_jprb/)
170 zpiza6(5, :)= (/&
171  &.8287144_jprb, .9949396_jprb, .9279543_jprb, .6765051_jprb, .9467578_jprb, .9955938_jprb/)
172 zpiza6(6, :)= (/&
173  &.5230504_jprb, .7868518_jprb, .8531531_jprb, .4048149_jprb, .8748231_jprb, .2355667_jprb/)
174 
175 zcga6(1, :)= (/&
176  &0.729019_jprb, 0.803129_jprb, 0.784592_jprb, 0.712208_jprb, .7008249_jprb, .7270548_jprb/)
177 zcga6(2, :)= (/&
178  &0.729019_jprb, 0.803129_jprb, 0.784592_jprb, 0.712208_jprb, .7008249_jprb, .7270548_jprb/)
179 zcga6(3, :)= (/&
180  &0.729019_jprb, 0.803129_jprb, 0.784592_jprb, 0.712208_jprb, .7008249_jprb, .7270548_jprb/)
181 zcga6(4, :)= (/&
182  &0.668431_jprb, 0.788530_jprb, 0.698682_jprb, 0.657422_jprb, .6735182_jprb, .6519706_jprb/)
183 zcga6(5, :)= (/&
184  &0.636342_jprb, 0.802467_jprb, 0.691305_jprb, 0.627497_jprb, .6105750_jprb, .4760794_jprb/)
185 zcga6(6, :)= (/&
186  &0.700610_jprb, 0.818871_jprb, 0.702399_jprb, 0.689886_jprb, .4629866_jprb, .1907639_jprb/)
187 
188 
189 ! ----------------------------------------------------------------
190 
191 IF (ksw == 2) THEN
192  DO jnu=1,ksw
193  DO jaer=1,6
194  rtaua(jnu,jaer)=ztaua2(jnu,jaer)
195  rpiza(jnu,jaer)=zpiza2(jnu,jaer)
196  rcga(jnu,jaer) =zcga2(jnu,jaer)
197  ENDDO
198  ENDDO
199 ELSEIF (ksw == 4) THEN
200  DO jnu=1,ksw
201  DO jaer=1,6
202  rtaua(jnu,jaer)=ztaua4(jnu,jaer)
203  rpiza(jnu,jaer)=zpiza4(jnu,jaer)
204  rcga(jnu,jaer) =zcga4(jnu,jaer)
205  ENDDO
206  ENDDO
207 ELSEIF (ksw == 6) THEN
208  DO jnu=1,ksw
209  DO jaer=1,6
210  rtaua(jnu,jaer)=ztaua6(jnu,jaer)
211  rpiza(jnu,jaer)=zpiza6(jnu,jaer)
212  rcga(jnu,jaer) =zcga6(jnu,jaer)
213  ENDDO
214  ENDDO
215 ELSE
216  stop 'SUAERSN: WRONG NUMBER OF SPECTRAL INTERVALS'
217 ENDIF
218 
219 ! ----------------------------------------------------------------
220 
221 RETURN
222 END SUBROUTINE suaersn
yoesw
Definition: yoesw.F90:1
yoesw::rtaua
real(kind=jprb), dimension(6, 6) rtaua
Definition: yoesw.F90:115
suaersn
subroutine suaersn(KTSW, KSW)
Definition: suaersn.F90:2
yoesw::rcga
real(kind=jprb), dimension(6, 6) rcga
Definition: yoesw.F90:117
yoesw::rpiza
real(kind=jprb), dimension(6, 6) rpiza
Definition: yoesw.F90:116