1 | ! |
---|
2 | ! $Id: phys_local_var_mod.F90 1347 2010-04-13 15:12:56Z musat $ |
---|
3 | ! |
---|
4 | MODULE phys_local_var_mod |
---|
5 | |
---|
6 | ! Variables locales pour effectuer les appels en serie |
---|
7 | !====================================================================== |
---|
8 | ! |
---|
9 | ! |
---|
10 | !====================================================================== |
---|
11 | ! Declaration des variables |
---|
12 | |
---|
13 | REAL, SAVE, ALLOCATABLE :: t_seri(:,:), q_seri(:,:) |
---|
14 | !$OMP THREADPRIVATE(t_seri, q_seri) |
---|
15 | REAL, SAVE, ALLOCATABLE :: ql_seri(:,:),qs_seri(:,:) |
---|
16 | !$OMP THREADPRIVATE(ql_seri,qs_seri) |
---|
17 | REAL, SAVE, ALLOCATABLE :: u_seri(:,:), v_seri(:,:) |
---|
18 | !$OMP THREADPRIVATE(u_seri, v_seri) |
---|
19 | |
---|
20 | REAL, SAVE, ALLOCATABLE :: tr_seri(:,:,:) |
---|
21 | !$OMP THREADPRIVATE(tr_seri) |
---|
22 | REAL, SAVE, ALLOCATABLE :: d_t_dyn(:,:), d_q_dyn(:,:) |
---|
23 | !$OMP THREADPRIVATE(d_t_dyn, d_q_dyn) |
---|
24 | REAL, SAVE, ALLOCATABLE :: d_u_dyn(:,:), d_v_dyn(:,:) |
---|
25 | !$OMP THREADPRIVATE(d_u_dyn, d_v_dyn) |
---|
26 | REAL, SAVE, ALLOCATABLE :: d_t_con(:,:),d_q_con(:,:) |
---|
27 | !$OMP THREADPRIVATE(d_t_con,d_q_con) |
---|
28 | REAL, SAVE, ALLOCATABLE :: d_u_con(:,:),d_v_con(:,:) |
---|
29 | !$OMP THREADPRIVATE(d_u_con,d_v_con) |
---|
30 | REAL, SAVE, ALLOCATABLE :: d_t_wake(:,:),d_q_wake(:,:) |
---|
31 | !$OMP THREADPRIVATE( d_t_wake,d_q_wake) |
---|
32 | REAL, SAVE, ALLOCATABLE :: d_t_lsc(:,:),d_q_lsc(:,:),d_ql_lsc(:,:) |
---|
33 | !$OMP THREADPRIVATE(d_t_lsc,d_q_lsc,d_ql_lsc) |
---|
34 | REAL, SAVE, ALLOCATABLE :: d_t_ajsb(:,:), d_q_ajsb(:,:) |
---|
35 | !$OMP THREADPRIVATE(d_t_ajsb, d_q_ajsb) |
---|
36 | REAL, SAVE, ALLOCATABLE :: d_t_ajs(:,:), d_q_ajs(:,:) |
---|
37 | !$OMP THREADPRIVATE(d_t_ajs, d_q_ajs) |
---|
38 | REAL, SAVE, ALLOCATABLE :: d_u_ajs(:,:), d_v_ajs(:,:) |
---|
39 | !$OMP THREADPRIVATE(d_u_ajs, d_v_ajs) |
---|
40 | REAL, SAVE, ALLOCATABLE :: d_t_eva(:,:),d_q_eva(:,:) |
---|
41 | !$OMP THREADPRIVATE(d_t_eva,d_q_eva) |
---|
42 | !tendances dues a oro et lif |
---|
43 | REAL, SAVE, ALLOCATABLE :: d_t_oli(:,:) |
---|
44 | !$OMP THREADPRIVATE(d_t_oli) |
---|
45 | REAL, SAVE, ALLOCATABLE :: d_u_oli(:,:), d_v_oli(:,:) |
---|
46 | !$OMP THREADPRIVATE(d_u_oli, d_v_oli) |
---|
47 | REAL, SAVE, ALLOCATABLE :: d_t_vdf(:,:), d_q_vdf(:,:) |
---|
48 | !$OMP THREADPRIVATE( d_t_vdf, d_q_vdf) |
---|
49 | REAL, SAVE, ALLOCATABLE :: d_u_vdf(:,:), d_v_vdf(:,:) |
---|
50 | !$OMP THREADPRIVATE(d_u_vdf, d_v_vdf) |
---|
51 | REAL, SAVE, ALLOCATABLE :: d_t_oro(:,:) |
---|
52 | !$OMP THREADPRIVATE(d_t_oro) |
---|
53 | REAL, SAVE, ALLOCATABLE :: d_u_oro(:,:), d_v_oro(:,:) |
---|
54 | !$OMP THREADPRIVATE(d_u_oro, d_v_oro) |
---|
55 | REAL, SAVE, ALLOCATABLE :: d_t_lif(:,:) |
---|
56 | !$OMP THREADPRIVATE(d_t_lif) |
---|
57 | REAL, SAVE, ALLOCATABLE :: d_u_lif(:,:), d_v_lif(:,:) |
---|
58 | !$OMP THREADPRIVATE(d_u_lif, d_v_lif) |
---|
59 | ! Tendances Ondes de G non oro (runs strato). |
---|
60 | REAL, SAVE, ALLOCATABLE :: d_u_hin(:,:) |
---|
61 | !$OMP THREADPRIVATE(d_u_hin) |
---|
62 | REAL, SAVE, ALLOCATABLE :: d_v_hin(:,:) |
---|
63 | !$OMP THREADPRIVATE(d_v_hin) |
---|
64 | REAL, SAVE, ALLOCATABLE :: d_t_hin(:,:) |
---|
65 | !$OMP THREADPRIVATE(d_t_hin) |
---|
66 | |
---|
67 | ! tendance du a la conersion Ec -> E thermique |
---|
68 | REAL, SAVE, ALLOCATABLE :: d_t_ec(:,:) |
---|
69 | !$OMP THREADPRIVATE(d_t_ec) |
---|
70 | REAL, SAVE, ALLOCATABLE :: d_ts(:,:), d_tr(:,:,:) |
---|
71 | !$OMP THREADPRIVATE(d_ts, d_tr) |
---|
72 | |
---|
73 | ! diagnostique pour le rayonnement |
---|
74 | REAL, SAVE, ALLOCATABLE :: topswad_aero(:), solswad_aero(:) ! diag |
---|
75 | !$OMP THREADPRIVATE(topswad_aero,solswad_aero) |
---|
76 | REAL, SAVE, ALLOCATABLE :: topswai_aero(:), solswai_aero(:) ! diag |
---|
77 | !$OMP THREADPRIVATE(topswai_aero,solswai_aero) |
---|
78 | REAL, SAVE, ALLOCATABLE :: topswad0_aero(:), solswad0_aero(:) ! diag |
---|
79 | !$OMP THREADPRIVATE(topswad0_aero,solswad0_aero) |
---|
80 | REAL, SAVE, ALLOCATABLE :: topsw_aero(:,:), solsw_aero(:,:) ! diag |
---|
81 | !$OMP THREADPRIVATE(topsw_aero,solsw_aero) |
---|
82 | REAL, SAVE, ALLOCATABLE :: topsw0_aero(:,:), solsw0_aero(:,:) ! diag |
---|
83 | !$OMP THREADPRIVATE(topsw0_aero,solsw0_aero) |
---|
84 | REAL, SAVE, ALLOCATABLE :: topswcf_aero(:,:), solswcf_aero(:,:) ! diag |
---|
85 | !$OMP THREADPRIVATE(topswcf_aero,solswcf_aero) |
---|
86 | REAL, SAVE, ALLOCATABLE :: tausum_aero(:,:,:) |
---|
87 | !$OMP THREADPRIVATE(tausum_aero) |
---|
88 | REAL, SAVE, ALLOCATABLE :: tau3d_aero(:,:,:,:) |
---|
89 | !$OMP THREADPRIVATE(tau3d_aero) |
---|
90 | REAL, SAVE, ALLOCATABLE :: scdnc(:,:) |
---|
91 | !$OMP THREADPRIVATE(scdnc) |
---|
92 | REAL, SAVE, ALLOCATABLE :: cldncl(:) |
---|
93 | !$OMP THREADPRIVATE(cldncl) |
---|
94 | REAL, SAVE, ALLOCATABLE :: reffclwtop(:) |
---|
95 | !$OMP THREADPRIVATE(reffclwtop) |
---|
96 | REAL, SAVE, ALLOCATABLE :: lcc(:) |
---|
97 | !$OMP THREADPRIVATE(lcc) |
---|
98 | REAL, SAVE, ALLOCATABLE :: reffclws(:,:) |
---|
99 | !$OMP THREADPRIVATE(reffclws) |
---|
100 | REAL, SAVE, ALLOCATABLE :: reffclwc(:,:) |
---|
101 | !$OMP THREADPRIVATE(reffclwc) |
---|
102 | REAL, SAVE, ALLOCATABLE :: cldnvi(:) |
---|
103 | !$OMP THREADPRIVATE(cldnvi) |
---|
104 | REAL, SAVE, ALLOCATABLE :: lcc3d(:,:) |
---|
105 | !$OMP THREADPRIVATE(lcc3d) |
---|
106 | REAL, SAVE, ALLOCATABLE :: lcc3dcon(:,:) |
---|
107 | !$OMP THREADPRIVATE(lcc3dcon) |
---|
108 | REAL, SAVE, ALLOCATABLE :: lcc3dstra(:,:) |
---|
109 | !$OMP THREADPRIVATE(lcc3dstra) |
---|
110 | REAL, SAVE, ALLOCATABLE :: od550aer(:) |
---|
111 | !$OMP THREADPRIVATE(od550aer) |
---|
112 | REAL, SAVE, ALLOCATABLE :: absvisaer(:) |
---|
113 | !$OMP THREADPRIVATE(absvisaer) |
---|
114 | REAL, SAVE, ALLOCATABLE :: od865aer(:) |
---|
115 | !$OMP THREADPRIVATE(od865aer) |
---|
116 | REAL, SAVE, ALLOCATABLE :: ec550aer(:,:) |
---|
117 | !$OMP THREADPRIVATE(ec550aer) |
---|
118 | REAL, SAVE, ALLOCATABLE :: od550lt1aer(:) |
---|
119 | !$OMP THREADPRIVATE(od550lt1aer) |
---|
120 | REAL, SAVE, ALLOCATABLE :: sconcso4(:) |
---|
121 | !$OMP THREADPRIVATE(sconcso4) |
---|
122 | REAL, SAVE, ALLOCATABLE :: sconcoa(:) |
---|
123 | !$OMP THREADPRIVATE(sconcoa) |
---|
124 | REAL, SAVE, ALLOCATABLE :: sconcbc(:) |
---|
125 | !$OMP THREADPRIVATE(sconcbc) |
---|
126 | REAL, SAVE, ALLOCATABLE :: sconcss(:) |
---|
127 | !$OMP THREADPRIVATE(sconcss) |
---|
128 | REAL, SAVE, ALLOCATABLE :: sconcdust(:) |
---|
129 | !$OMP THREADPRIVATE(sconcdust) |
---|
130 | REAL, SAVE, ALLOCATABLE :: concso4(:,:) |
---|
131 | !$OMP THREADPRIVATE(concso4) |
---|
132 | REAL, SAVE, ALLOCATABLE :: concoa(:,:) |
---|
133 | !$OMP THREADPRIVATE(concoa) |
---|
134 | REAL, SAVE, ALLOCATABLE :: concbc(:,:) |
---|
135 | !$OMP THREADPRIVATE(concbc) |
---|
136 | REAL, SAVE, ALLOCATABLE :: concss(:,:) |
---|
137 | !$OMP THREADPRIVATE(concss) |
---|
138 | REAL, SAVE, ALLOCATABLE :: concdust(:,:) |
---|
139 | !$OMP THREADPRIVATE(concdust) |
---|
140 | REAL, SAVE, ALLOCATABLE :: loadso4(:) |
---|
141 | !$OMP THREADPRIVATE(loadso4) |
---|
142 | REAL, SAVE, ALLOCATABLE :: loadoa(:) |
---|
143 | !$OMP THREADPRIVATE(loadoa) |
---|
144 | REAL, SAVE, ALLOCATABLE :: loadbc(:) |
---|
145 | !$OMP THREADPRIVATE(loadbc) |
---|
146 | REAL, SAVE, ALLOCATABLE :: loadss(:) |
---|
147 | !$OMP THREADPRIVATE(loadss) |
---|
148 | REAL, SAVE, ALLOCATABLE :: loaddust(:) |
---|
149 | !$OMP THREADPRIVATE(loaddust) |
---|
150 | REAL, SAVE, ALLOCATABLE :: load_tmp1(:) |
---|
151 | !$OMP THREADPRIVATE(load_tmp1) |
---|
152 | REAL, SAVE, ALLOCATABLE :: load_tmp2(:) |
---|
153 | !$OMP THREADPRIVATE(load_tmp2) |
---|
154 | REAL, SAVE, ALLOCATABLE :: load_tmp3(:) |
---|
155 | !$OMP THREADPRIVATE(load_tmp3) |
---|
156 | REAL, SAVE, ALLOCATABLE :: load_tmp4(:) |
---|
157 | !$OMP THREADPRIVATE(load_tmp4) |
---|
158 | REAL, SAVE, ALLOCATABLE :: load_tmp5(:) |
---|
159 | !$OMP THREADPRIVATE(load_tmp5) |
---|
160 | REAL, SAVE, ALLOCATABLE :: load_tmp6(:) |
---|
161 | !$OMP THREADPRIVATE(load_tmp6) |
---|
162 | REAL, SAVE, ALLOCATABLE :: load_tmp7(:) |
---|
163 | !$OMP THREADPRIVATE(load_tmp7) |
---|
164 | |
---|
165 | CONTAINS |
---|
166 | |
---|
167 | !====================================================================== |
---|
168 | SUBROUTINE phys_local_var_init |
---|
169 | use dimphy |
---|
170 | use infotrac, ONLY : nbtr |
---|
171 | USE aero_mod |
---|
172 | |
---|
173 | IMPLICIT NONE |
---|
174 | #include "indicesol.h" |
---|
175 | allocate(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev)) |
---|
176 | allocate(u_seri(klon,klev),v_seri(klon,klev)) |
---|
177 | |
---|
178 | allocate(tr_seri(klon,klev,nbtr)) |
---|
179 | allocate(d_t_dyn(klon,klev),d_q_dyn(klon,klev)) |
---|
180 | allocate(d_u_dyn(klon,klev),d_v_dyn(klon,klev)) |
---|
181 | allocate(d_t_con(klon,klev),d_q_con(klon,klev)) |
---|
182 | allocate(d_u_con(klon,klev),d_v_con(klon,klev)) |
---|
183 | allocate(d_t_wake(klon,klev),d_q_wake(klon,klev)) |
---|
184 | allocate(d_t_lsc(klon,klev),d_q_lsc(klon,klev)) |
---|
185 | allocate(d_ql_lsc(klon,klev)) |
---|
186 | allocate(d_t_ajsb(klon,klev),d_q_ajsb(klon,klev)) |
---|
187 | allocate(d_t_ajs(klon,klev),d_q_ajs(klon,klev)) |
---|
188 | allocate(d_u_ajs(klon,klev),d_v_ajs(klon,klev)) |
---|
189 | allocate(d_t_eva(klon,klev),d_q_eva(klon,klev)) |
---|
190 | allocate(d_t_vdf(klon,klev),d_q_vdf(klon,klev)) |
---|
191 | allocate(d_u_vdf(klon,klev),d_v_vdf(klon,klev)) |
---|
192 | allocate(d_t_oli(klon,klev),d_t_oro(klon,klev)) |
---|
193 | allocate(d_u_oli(klon,klev),d_v_oli(klon,klev)) |
---|
194 | allocate(d_u_oro(klon,klev),d_v_oro(klon,klev)) |
---|
195 | allocate(d_t_lif(klon,klev),d_t_ec(klon,klev)) |
---|
196 | allocate(d_u_lif(klon,klev),d_v_lif(klon,klev)) |
---|
197 | allocate(d_ts(klon,klev), d_tr(klon,klev,nbtr)) |
---|
198 | allocate(topswad_aero(klon), solswad_aero(klon)) |
---|
199 | allocate(topswai_aero(klon), solswai_aero(klon)) |
---|
200 | allocate(topswad0_aero(klon), solswad0_aero(klon)) |
---|
201 | allocate(topsw_aero(klon,naero_grp), solsw_aero(klon,naero_grp)) |
---|
202 | allocate(topsw0_aero(klon,naero_grp), solsw0_aero(klon,naero_grp)) |
---|
203 | allocate(topswcf_aero(klon,3), solswcf_aero(klon,3)) |
---|
204 | allocate(d_u_hin(klon,klev),d_v_hin(klon,klev),d_t_hin(klon,klev)) |
---|
205 | allocate(tausum_aero(klon,nwave,naero_spc)) |
---|
206 | allocate(tau3d_aero(klon,klev,nwave,naero_spc)) |
---|
207 | allocate(scdnc(klon, klev)) |
---|
208 | allocate(cldncl(klon)) |
---|
209 | allocate(reffclwtop(klon)) |
---|
210 | allocate(lcc(klon)) |
---|
211 | allocate(reffclws(klon, klev)) |
---|
212 | allocate(reffclwc(klon, klev)) |
---|
213 | allocate(cldnvi(klon)) |
---|
214 | allocate(lcc3d(klon, klev)) |
---|
215 | allocate(lcc3dcon(klon, klev)) |
---|
216 | allocate(lcc3dstra(klon, klev)) |
---|
217 | allocate(od550aer(klon)) |
---|
218 | allocate(od865aer(klon)) |
---|
219 | allocate(absvisaer(klon)) |
---|
220 | allocate(ec550aer(klon,klev)) |
---|
221 | allocate(od550lt1aer(klon)) |
---|
222 | allocate(sconcso4(klon)) |
---|
223 | allocate(sconcoa(klon)) |
---|
224 | allocate(sconcbc(klon)) |
---|
225 | allocate(sconcss(klon)) |
---|
226 | allocate(sconcdust(klon)) |
---|
227 | allocate(concso4(klon,klev)) |
---|
228 | allocate(concoa(klon,klev)) |
---|
229 | allocate(concbc(klon,klev)) |
---|
230 | allocate(concss(klon,klev)) |
---|
231 | allocate(concdust(klon,klev)) |
---|
232 | allocate(loadso4(klon)) |
---|
233 | allocate(loadoa(klon)) |
---|
234 | allocate(loadbc(klon)) |
---|
235 | allocate(loadss(klon)) |
---|
236 | allocate(loaddust(klon)) |
---|
237 | allocate(load_tmp1(klon)) |
---|
238 | allocate(load_tmp2(klon)) |
---|
239 | allocate(load_tmp3(klon)) |
---|
240 | allocate(load_tmp4(klon)) |
---|
241 | allocate(load_tmp5(klon)) |
---|
242 | allocate(load_tmp6(klon)) |
---|
243 | allocate(load_tmp7(klon)) |
---|
244 | |
---|
245 | END SUBROUTINE phys_local_var_init |
---|
246 | |
---|
247 | !====================================================================== |
---|
248 | SUBROUTINE phys_local_var_end |
---|
249 | use dimphy |
---|
250 | IMPLICIT NONE |
---|
251 | #include "indicesol.h" |
---|
252 | deallocate(t_seri,q_seri,ql_seri,qs_seri) |
---|
253 | deallocate(u_seri,v_seri) |
---|
254 | |
---|
255 | deallocate(tr_seri) |
---|
256 | deallocate(d_t_dyn,d_q_dyn) |
---|
257 | deallocate(d_u_dyn,d_v_dyn) |
---|
258 | deallocate(d_t_con,d_q_con) |
---|
259 | deallocate(d_u_con,d_v_con) |
---|
260 | deallocate(d_t_wake,d_q_wake) |
---|
261 | deallocate(d_t_lsc,d_q_lsc) |
---|
262 | deallocate(d_ql_lsc) |
---|
263 | deallocate(d_t_ajsb,d_q_ajsb) |
---|
264 | deallocate(d_t_ajs,d_q_ajs) |
---|
265 | deallocate(d_u_ajs,d_v_ajs) |
---|
266 | deallocate(d_t_eva,d_q_eva) |
---|
267 | deallocate(d_t_vdf,d_q_vdf) |
---|
268 | deallocate(d_u_vdf,d_v_vdf) |
---|
269 | deallocate(d_t_oli,d_t_oro) |
---|
270 | deallocate(d_u_oli,d_v_oli) |
---|
271 | deallocate(d_u_oro,d_v_oro) |
---|
272 | deallocate(d_t_lif,d_t_ec) |
---|
273 | deallocate(d_u_lif,d_v_lif) |
---|
274 | deallocate(d_ts, d_tr) |
---|
275 | deallocate(topswad_aero,solswad_aero) |
---|
276 | deallocate(topswai_aero,solswai_aero) |
---|
277 | deallocate(topswad0_aero,solswad0_aero) |
---|
278 | deallocate(topsw_aero,solsw_aero) |
---|
279 | deallocate(topsw0_aero,solsw0_aero) |
---|
280 | deallocate(topswcf_aero,solswcf_aero) |
---|
281 | deallocate(tausum_aero) |
---|
282 | deallocate(tau3d_aero) |
---|
283 | deallocate(scdnc) |
---|
284 | deallocate(cldncl) |
---|
285 | deallocate(reffclwtop) |
---|
286 | deallocate(lcc) |
---|
287 | deallocate(reffclws) |
---|
288 | deallocate(reffclwc) |
---|
289 | deallocate(cldnvi) |
---|
290 | deallocate(lcc3d) |
---|
291 | deallocate(lcc3dcon) |
---|
292 | deallocate(lcc3dstra) |
---|
293 | deallocate(od550aer) |
---|
294 | deallocate(od865aer) |
---|
295 | deallocate(absvisaer) |
---|
296 | deallocate(ec550aer) |
---|
297 | deallocate(od550lt1aer) |
---|
298 | deallocate(sconcso4) |
---|
299 | deallocate(sconcoa) |
---|
300 | deallocate(sconcbc) |
---|
301 | deallocate(sconcss) |
---|
302 | deallocate(sconcdust) |
---|
303 | deallocate(concso4) |
---|
304 | deallocate(concoa) |
---|
305 | deallocate(concbc) |
---|
306 | deallocate(concss) |
---|
307 | deallocate(concdust) |
---|
308 | deallocate(loadso4) |
---|
309 | deallocate(loadoa) |
---|
310 | deallocate(loadbc) |
---|
311 | deallocate(loadss) |
---|
312 | deallocate(loaddust) |
---|
313 | deallocate(load_tmp1) |
---|
314 | deallocate(load_tmp2) |
---|
315 | deallocate(load_tmp3) |
---|
316 | deallocate(load_tmp4) |
---|
317 | deallocate(load_tmp5) |
---|
318 | deallocate(load_tmp6) |
---|
319 | deallocate(load_tmp7) |
---|
320 | deallocate(d_u_hin,d_v_hin,d_t_hin) |
---|
321 | |
---|
322 | END SUBROUTINE phys_local_var_end |
---|
323 | |
---|
324 | END MODULE phys_local_var_mod |
---|