1 | MODULE YOESW |
---|
2 | |
---|
3 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
---|
4 | |
---|
5 | IMPLICIT NONE |
---|
6 | |
---|
7 | SAVE |
---|
8 | |
---|
9 | ! ------------------------------------------------------------------ |
---|
10 | !* ** *YOESW* - COEFFICIENTS FOR SHORTWAVE RADIATION TRANSFER |
---|
11 | ! ------------------------------------------------------------------ |
---|
12 | |
---|
13 | REAL(KIND=JPRB) :: APAD(6,3,7) |
---|
14 | REAL(KIND=JPRB) :: BPAD(6,3,7) |
---|
15 | REAL(KIND=JPRB) :: RRAY(6,6) |
---|
16 | REAL(KIND=JPRB), ALLOCATABLE :: RSUN(:) |
---|
17 | !$OMP THREADPRIVATE(RSUN) |
---|
18 | REAL(KIND=JPRB) :: RPDH1 |
---|
19 | REAL(KIND=JPRB) :: RPDU1 |
---|
20 | REAL(KIND=JPRB) :: RPNH |
---|
21 | REAL(KIND=JPRB) :: RPNU |
---|
22 | REAL(KIND=JPRB) :: RSWCE(6) |
---|
23 | REAL(KIND=JPRB) :: RSWCP(6) |
---|
24 | REAL(KIND=JPRB) :: RTDH2O |
---|
25 | REAL(KIND=JPRB) :: RTDUMG |
---|
26 | REAL(KIND=JPRB) :: RTH2O |
---|
27 | REAL(KIND=JPRB) :: RTUMG |
---|
28 | REAL(KIND=JPRB) :: D(6,3) |
---|
29 | REAL(KIND=JPRB) :: REXPO3(6,2,7) |
---|
30 | INTEGER(KIND=JPIM) :: NEXPO3(6) |
---|
31 | |
---|
32 | REAL(KIND=JPRB) :: RYFWCA(6) |
---|
33 | REAL(KIND=JPRB) :: RYFWCB(6) |
---|
34 | REAL(KIND=JPRB) :: RYFWCC(6) |
---|
35 | REAL(KIND=JPRB) :: RYFWCD(6) |
---|
36 | REAL(KIND=JPRB) :: RYFWCE(6) |
---|
37 | REAL(KIND=JPRB) :: RYFWCF(6) |
---|
38 | |
---|
39 | REAL(KIND=JPRB) :: REBCUA(6) |
---|
40 | REAL(KIND=JPRB) :: REBCUB(6) |
---|
41 | REAL(KIND=JPRB) :: REBCUC(6) |
---|
42 | REAL(KIND=JPRB) :: REBCUD(6) |
---|
43 | REAL(KIND=JPRB) :: REBCUE(6) |
---|
44 | REAL(KIND=JPRB) :: REBCUF(6) |
---|
45 | REAL(KIND=JPRB) :: REBCUG(16) |
---|
46 | REAL(KIND=JPRB) :: REBCUH(16) |
---|
47 | REAL(KIND=JPRB) :: REBCUI(6) |
---|
48 | REAL(KIND=JPRB) :: REBCUJ(6) |
---|
49 | |
---|
50 | REAL(KIND=JPRB) :: RASWCA(6) |
---|
51 | REAL(KIND=JPRB) :: RASWCB(6) |
---|
52 | REAL(KIND=JPRB) :: RASWCC(6) |
---|
53 | REAL(KIND=JPRB) :: RASWCD(6) |
---|
54 | REAL(KIND=JPRB) :: RASWCE(6) |
---|
55 | REAL(KIND=JPRB) :: RASWCF(6) |
---|
56 | |
---|
57 | REAL(KIND=JPRB) :: RFUETA(16,3),RFUETB(16,4), RFUETC(16,4) |
---|
58 | REAL(KIND=JPRB) :: RFULIO(16,3) |
---|
59 | REAL(KIND=JPRB) :: RHSAVI(16,3) |
---|
60 | REAL(KIND=JPRB) :: RLILIA(16,5),RLILIB(16,4) |
---|
61 | |
---|
62 | REAL(KIND=JPRB) :: RFLAA0(6) |
---|
63 | REAL(KIND=JPRB) :: RFLAA1(6) |
---|
64 | REAL(KIND=JPRB) :: RFLBB0(6) |
---|
65 | REAL(KIND=JPRB) :: RFLBB1(6) |
---|
66 | REAL(KIND=JPRB) :: RFLBB2(6) |
---|
67 | REAL(KIND=JPRB) :: RFLBB3(6) |
---|
68 | REAL(KIND=JPRB) :: RFLCC0(6) |
---|
69 | REAL(KIND=JPRB) :: RFLCC1(6) |
---|
70 | REAL(KIND=JPRB) :: RFLCC2(6) |
---|
71 | REAL(KIND=JPRB) :: RFLCC3(6) |
---|
72 | |
---|
73 | REAL(KIND=JPRB) :: RFUAA0(6) |
---|
74 | REAL(KIND=JPRB) :: RFUAA1(6) |
---|
75 | REAL(KIND=JPRB) :: RFUBB0(6) |
---|
76 | REAL(KIND=JPRB) :: RFUBB1(6) |
---|
77 | REAL(KIND=JPRB) :: RFUBB2(6) |
---|
78 | REAL(KIND=JPRB) :: RFUBB3(6) |
---|
79 | REAL(KIND=JPRB) :: RFUCC0(6) |
---|
80 | REAL(KIND=JPRB) :: RFUCC1(6) |
---|
81 | REAL(KIND=JPRB) :: RFUCC2(6) |
---|
82 | REAL(KIND=JPRB) :: RFUCC3(6) |
---|
83 | REAL(KIND=JPRB) :: RFLDD0(6) |
---|
84 | REAL(KIND=JPRB) :: RFLDD1(6) |
---|
85 | REAL(KIND=JPRB) :: RFLDD2(6) |
---|
86 | REAL(KIND=JPRB) :: RFLDD3(6) |
---|
87 | |
---|
88 | REAL(KIND=JPRB) :: RSUSHE(6) |
---|
89 | REAL(KIND=JPRB) :: RSUSHF(6) |
---|
90 | REAL(KIND=JPRB) :: RSUSHH(6) |
---|
91 | REAL(KIND=JPRB) :: RSUSHK(6) |
---|
92 | REAL(KIND=JPRB) :: RSUSHA(6) |
---|
93 | REAL(KIND=JPRB) :: RSUSHG(6) |
---|
94 | REAL(KIND=JPRB) :: RSUSHFA(4) |
---|
95 | REAL(KIND=JPRB) :: RSUSHC |
---|
96 | REAL(KIND=JPRB) :: RSUSHD |
---|
97 | |
---|
98 | REAL(KIND=JPRB) :: REFFIA |
---|
99 | REAL(KIND=JPRB) :: REFFIB |
---|
100 | REAL(KIND=JPRB) :: RTIW |
---|
101 | REAL(KIND=JPRB) :: RRIW |
---|
102 | REAL(KIND=JPRB) :: RROMA(6) |
---|
103 | REAL(KIND=JPRB) :: RROMB(6) |
---|
104 | REAL(KIND=JPRB) :: RRASY(6) |
---|
105 | |
---|
106 | REAL(KIND=JPRB) :: RHSRA(6) |
---|
107 | REAL(KIND=JPRB) :: RHSRB(6) |
---|
108 | REAL(KIND=JPRB) :: RHSRC(6) |
---|
109 | REAL(KIND=JPRB) :: RHSRD(6) |
---|
110 | REAL(KIND=JPRB) :: RHSRE(6) |
---|
111 | REAL(KIND=JPRB) :: RHSRF(6) |
---|
112 | REAL(KIND=JPRB) :: RHSRTA |
---|
113 | REAL(KIND=JPRB) :: RHSRTB |
---|
114 | |
---|
115 | REAL(KIND=JPRB) :: RTAUA(6,6) |
---|
116 | REAL(KIND=JPRB) :: RPIZA(6,6) |
---|
117 | REAL(KIND=JPRB) :: RCGA(6,6) |
---|
118 | REAL(KIND=JPRB) :: RAER(6,6) |
---|
119 | |
---|
120 | REAL(KIND=JPRB) :: RTWEIGHT(8) |
---|
121 | REAL(KIND=JPRB) :: RWEIGS(6) |
---|
122 | REAL(KIND=JPRB) :: RWEIGV(6) |
---|
123 | |
---|
124 | REAL(KIND=JPRB) :: RWEIGHT(6,8) |
---|
125 | INTEGER(KIND=JPIM) :: NMPSRTM(14), NTYPS |
---|
126 | |
---|
127 | REAL(KIND=JPRB) :: RADJUST |
---|
128 | |
---|
129 | ! * E.C.M.W.F. PHYSICS PACKAGE * |
---|
130 | |
---|
131 | ! J.-J. MORCRETTE E.C.M.W.F. 89/07/14 |
---|
132 | |
---|
133 | ! NAME TYPE PURPOSE |
---|
134 | ! ---- : ---- : --------------------------------------------------- |
---|
135 | ! APAD : REAL PADE APPROXIMANTS NUMERATOR |
---|
136 | ! BPAD : REAL PADE APPROXIMANTS DENOMINATOR |
---|
137 | ! D : REAL TRANSMISSION LIMIT FOR INFINITE ABSORBER AMOUNT |
---|
138 | ! RRAY : REAL RAYLEIGH SCATTERING COEFFICIENTS |
---|
139 | ! RSUN : REAL SOLAR FRACTION IN SPECTRAL INTERVALS |
---|
140 | ! RPDH1 : 1 + EXPONENT PRESSURE DEPENDENCE H2O |
---|
141 | ! RPDU1 : 1 + EXPONENT PRESSURE DEPENDENCE UNIFORMLY MIXED GASES |
---|
142 | ! RPNH : REFERENCE PRESSURE FACTOR FOR H2O |
---|
143 | ! RPNU : REFERENCE PRESSURE FACTOR FOR UNIFORMLY MIXED GASES |
---|
144 | ! RSWCE : E-TYPE, H2O CONTINUUM ABSORPTION COEFFICIENT |
---|
145 | ! RSWCP : P-TYPE, H2O CONTINUUM ABSORPTION COEFFICIENT |
---|
146 | ! RTDH2O: EXPONENT TEMPERATURE DEPENDENCE H2O |
---|
147 | ! RTDUMG: EXPONENT TEMPERATURE DEPENDENCE UNIFORMLY MIXED GASES |
---|
148 | ! RTH2O : REFERENCE TEMPERATURE H2O |
---|
149 | ! RTUMG : REFERENCE TEMPERATURE UNIFORMLY MIXED GASES |
---|
150 | ! ----------------------------------------------------------------- |
---|
151 | |
---|
152 | ! * E.C.M.W.F. PHYSICS PACKAGE * |
---|
153 | |
---|
154 | ! J.-J. MORCRETTE E.C.M.W.F. 89/07/14 |
---|
155 | |
---|
156 | ! NAME TYPE PURPOSE |
---|
157 | ! ---- : ---- : --------------------------------------------------- |
---|
158 | !* FOUQUART (1987) WATER CLOUD OPTICAL PROPERTIES |
---|
159 | |
---|
160 | ! RYFWCA : REAL : C1 IN OPTICAL THICKNESS FORMULA |
---|
161 | ! RYFWCB : REAL : C2 IN OPTICAL THICKNESS FORMULA |
---|
162 | ! RYFWCC : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
---|
163 | ! RYFWCD : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
---|
164 | ! RYFWCE : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
---|
165 | ! RYFWCF : REAL : ASSYMETRY FACTOR |
---|
166 | |
---|
167 | !* SLINGO (1989) WATER CLOUD OPTICAL PROPERTIES |
---|
168 | |
---|
169 | ! RASWCA : REAL : C1 IN OPTICAL THICKNESS FORMULA |
---|
170 | ! RASWCB : REAL : C2 IN OPTICAL THICKNESS FORMULA |
---|
171 | ! RASWCC : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
---|
172 | ! RASWCD : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
---|
173 | ! RASWCE : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
---|
174 | ! RASWCF : REAL : ASSYMETRY FACTOR |
---|
175 | |
---|
176 | !* SAVIJARVI (1998) WATER CLOUD OPTICAL PROPERTIES (RRTM) |
---|
177 | |
---|
178 | ! RHSAVI : REAL : MASS ABSORPTION COEFFICIENTS (POLYNOMIAL DEVELOPM) |
---|
179 | |
---|
180 | !* LINDNER,LI (2000) WATER CLOUD OPTICAL PROPERTIES (RRTM) |
---|
181 | |
---|
182 | ! RLILIA : REAL : MASS ABSORPTION COEFFICIENTS (POLYNOMIAL DEVELOPM) |
---|
183 | ! RLILIB : REAL : 1-SSA COEFFICIENTS (POLYNOMIAL DEVELOPM) |
---|
184 | |
---|
185 | !* ICE CLOUD OPTICAL PROPERTIES DERIVED FROM EBERT-CURRY (1992) |
---|
186 | |
---|
187 | ! REBCUA : REAL : C1 IN OPTICAL THICKNESS FORMULA |
---|
188 | ! REBCUB : REAL : C2 IN OPTICAL THICKNESS FORMULA |
---|
189 | ! REBCUC : REAL : 1-C3 IN SINGLE SCATTERING ALBEDO FORMULA |
---|
190 | ! REBCUD : REAL : C4 IN SINGLE SCATTERING ALBEDO FORMULA |
---|
191 | ! REBCUE : REAL : C5 IN ASSYMETRY FACTOR FORMULA |
---|
192 | ! REBCUF : REAL : C6 IN ASSYMETRY FACTOR FORMULA |
---|
193 | ! REBCUG : REAL : C7 IN MASS ABSORPTION COEFFICIENT FORMULA |
---|
194 | ! REBCUH : REAL : C8 IN MASS ABSORPTION COEFFICIENT FORMULA |
---|
195 | ! REBCUI : REAL : C7 IN MASS ABSORPTION COEFFICIENT SPECTRAL FORMULA |
---|
196 | ! REBCUJ : REAL : C8 IN MASS ABSORPTION COEFFICIENT SPECTRAL FORMULA |
---|
197 | |
---|
198 | !* ICE CLOUD OPTICAL PROPERTIES DERIVED FROM SUN-SHINE (1995) |
---|
199 | |
---|
200 | ! RSHSUE : REAL : E IN SINGLE SCATTERING ALBEDO FORMULA |
---|
201 | ! RSHSUF : REAL : F IN SINGLE SCATTERING ALBEDO FORMULA |
---|
202 | ! RSHSUH : REAL : H IN ASSYMETRY FACTOR FORMULA |
---|
203 | ! RSHSUK : REAL : K IN ASSYMETRY FACTOR FORMULA |
---|
204 | ! RSHSUA : REAL : ALPHA IN SSA CORRECTION FACTOR FORMULA |
---|
205 | ! RSHSUG : REAL : GAMMA IN ASSYMETRY CORRECTION FACTOR FORMULA |
---|
206 | ! RSHSUFA: REAL : COEFFICIENTS IN TEMPERATURE CORRECTION FACTOR |
---|
207 | |
---|
208 | ! REFFIA : REAL : C9 IN EFFECTIVE RADIUS FORMULA |
---|
209 | ! REFFIB : REAL : C10 IN EFFECTIVE RADIUS FORMULA |
---|
210 | |
---|
211 | !* ICE CLOUD OPTICAL PROPERTIES DERIVED FROM FU-LIOU (1993) |
---|
212 | |
---|
213 | ! RFULIO : REAL : COEFFICIENTS IN EXPRESSION FOR LW EXTINCTION COEFF. |
---|
214 | ! RFLAA : REAL : COEFFICIENTS IN EXPRESSION FOR SW EXTINCTION COEFF. |
---|
215 | ! RFLBB : REAL : COEFFICIENTS IN EXPRESSION FOR SW SINGLE SCATT.ALB. |
---|
216 | ! RFLCC : REAL : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR |
---|
217 | ! RFLDD : REAL : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR |
---|
218 | |
---|
219 | !* ICE CLOUD OPTICAL PROPERTIES DERIVED FROM FU (1996) & FU ET AL. (1998) |
---|
220 | |
---|
221 | ! RFUETA : REAL : COEFFICIENTS IN EXPRESSION FOR LW EXTINCTION COEFF. |
---|
222 | ! RFUAA : REAL : COEFFICIENTS IN EXPRESSION FOR SW EXTINCTION COEFF. |
---|
223 | ! RFUBB : REAL : COEFFICIENTS IN EXPRESSION FOR SW SINGLE SCATT.ALB. |
---|
224 | ! RFUCC : REAL : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR |
---|
225 | |
---|
226 | !* TRANSITION BETWEEN LIQUID AND SOLID WATER |
---|
227 | |
---|
228 | ! RTIW : REAL : TEMPERATURE THRESHOLD |
---|
229 | ! RRIW : REAL : TRANSITION RANGE |
---|
230 | |
---|
231 | !* RAIN OPTICAL PROPERTIES FROM SAVIJARVI (1996) |
---|
232 | |
---|
233 | ! RROMA : REAL : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO |
---|
234 | ! RROMB : REAL : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO |
---|
235 | ! RRASY : REAL : COEFFICIENTS FOR ASSYMETRY FACTOR |
---|
236 | ! RHSRA : REAL : COEFFICIENTS FOR OPTICAL THICKNESS |
---|
237 | ! RHSRB : REAL : COEFFICIENTS FOR OPTICAL THICKNESS |
---|
238 | ! RHSRC : REAL : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO |
---|
239 | ! RHSRD : REAL : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO |
---|
240 | ! RHSRE : REAL : COEFFICIENTS FOR ASSYMETRY FACTOR |
---|
241 | ! RHSRF : REAL : COEFFICIENTS FOR ASSYMETRY FACTOR |
---|
242 | ! RHSRTA : REAL : COEFFICIENTS FOR OPTICAL THICKNESS |
---|
243 | ! RHSRTB : REAL : COEFFICIENTS FOR OPTICAL THICKNESS |
---|
244 | ! ----------------------------------------------------------------- |
---|
245 | |
---|
246 | ! * E.C.M.W.F. PHYSICS PACKAGE * |
---|
247 | |
---|
248 | ! J.-J. MORCRETTE E.C.M.W.F. 89/07/14 |
---|
249 | |
---|
250 | ! NAME TYPE PURPOSE |
---|
251 | ! ---- : ---- : ------- |
---|
252 | ! RTAUA : REAL S.W. NORMALIZED OPTICAL THICKNESS AT 0.55 MICRON |
---|
253 | ! RPIZA : REAL S.W. SINGLE SCATTERING ALBEDO |
---|
254 | ! RCGA : REAL S.W. ASSYMETRY FACTOR |
---|
255 | ! RAER : REAL L.W. ABSORPTION COEFFICIENTS |
---|
256 | ! ----------------------------------------------------------------- |
---|
257 | |
---|
258 | ! * E.C.M.W.F. PHYSICS PACKAGE * |
---|
259 | |
---|
260 | ! J.-J. MORCRETTE E.C.M.W.F. 89/07/14 |
---|
261 | |
---|
262 | ! NAME TYPE PURPOSE |
---|
263 | ! ---- : ---- : ------- |
---|
264 | !RTWEIGHT: REAL S.W. INTEGRATED WEIGHT |
---|
265 | ! NMPSRTM: INTEGER : Indices for mapping SW[1:6] albedo into SRTM[1:14] |
---|
266 | ! ----------------------------------------------------------------- |
---|
267 | !$OMP THREADPRIVATE(apad,bpad,d,nexpo3,nmpsrtm,ntyps,radjust,raer,raswca) |
---|
268 | !$OMP THREADPRIVATE(raswcb,raswcc,raswcd,raswce,raswcf,rcga,rebcua,rebcub) |
---|
269 | !$OMP THREADPRIVATE(rebcuc,rebcud,rebcue,rebcuf,rebcug,rebcuh,rebcui,rebcuj) |
---|
270 | !$OMP THREADPRIVATE(reffia,reffib,rexpo3,rflaa0,rflaa1,rflbb0,rflbb1,rflbb2) |
---|
271 | !$OMP THREADPRIVATE(rflbb3,rflcc0,rflcc1,rflcc2,rflcc3,rfldd0,rfldd1,rfldd2) |
---|
272 | !$OMP THREADPRIVATE(rfldd3,rfuaa0,rfuaa1,rfubb0,rfubb1,rfubb2,rfubb3,rfucc0) |
---|
273 | !$OMP THREADPRIVATE(rfucc1,rfucc2,rfucc3,rfueta,rfuetb,rfuetc,rfulio,rhsavi) |
---|
274 | !$OMP THREADPRIVATE(rhsra,rhsrb,rhsrc,rhsrd,rhsre,rhsrf,rhsrta,rhsrtb,rlilia) |
---|
275 | !$OMP THREADPRIVATE(rlilib,rpdh1,rpdu1,rpiza,rpnh,rpnu,rrasy,rray,rriw,rroma) |
---|
276 | !$OMP THREADPRIVATE(rromb,rsusha,rsushc,rsushd,rsushe,rsushf,rsushfa,rsushg) |
---|
277 | !$OMP THREADPRIVATE(rsushh,rsushk,rswce,rswcp,rtaua,rtdh2o,rtdumg,rth2o,rtiw) |
---|
278 | !$OMP THREADPRIVATE(rtumg,rtweight,rweight,rweigs,rweigv,ryfwca,ryfwcb,ryfwcc,ryfwcd,ryfwce,ryfwcf) |
---|
279 | END MODULE YOESW |
---|