40 use upp_ifi_mod,
only: set_ifi_dims
46 use upp_math,
only: ddvdx, ddudy, uuavg
58 allocate(u(ista_2l:iend_2u+1,jsta_2l:jend_2u,lm))
59 allocate(v(ista_2l:iend_2u,jsta_2l:jvend_2u,lm))
60 allocate(t(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
63 allocate(q(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
65 allocate(uh(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
66 allocate(vh(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
67 allocate(wh(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
68 allocate(pmid(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
69 allocate(pmidv(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
70 allocate(pint(ista_2l:iend_2u,jsta_2l:jend_2u,lp1))
71 allocate(alpint(ista_2l:iend_2u,jsta_2l:jend_2u,lp1))
72 allocate(zmid(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
73 allocate(zint(ista_2l:iend_2u,jsta_2l:jend_2u,lp1))
75 allocate(q2(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
76 allocate(omga(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
77 allocate(dpres(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
78 allocate(t_adj(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
79 allocate(ttnd(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
80 allocate(rswtt(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
81 allocate(rlwtt(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
82 allocate(exch_h(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
83 allocate(train(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
84 allocate(tcucn(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
85 allocate(el_pbl(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
93 do i=ista_2l,iend_2u+1
100 do j=jsta_2l,jvend_2u
144 allocate(cwm(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
145 allocate(f_ice(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
146 allocate(f_rain(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
147 allocate(f_rimef(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
148 allocate(qqw(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
149 allocate(qrimef(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
150 allocate(qqi(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
151 allocate(qqr(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
152 allocate(qqs(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
153 allocate(qqg(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
154 allocate(qqh(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
155 allocate(qqnw(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
156 allocate(qqni(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
157 allocate(qqnr(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
158 allocate(qqnwfa(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
159 allocate(qqnifa(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
160 allocate(taod5503d(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
161 allocate(aextc55(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
162 allocate(extcof55(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
163 allocate(qc_bl(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
164 allocate(cfr(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
165 allocate(effri(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
166 allocate(effrl(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
167 allocate(effrs(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
168 allocate(cfr_raw(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
169 allocate(dbz(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
170 allocate(dbzr(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
171 allocate(dbzi(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
172 allocate(dbzc(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
173 allocate(mcvg(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
174 allocate(nlice(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
196 taod5503d(i,j,l)=spval
215 allocate(nrain(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
216 allocate(radius_cloud(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
217 allocate(radius_ice(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
218 allocate(radius_snow(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
221 allocate(refl_10cm(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
223 allocate(o3(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
224 allocate(o(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
225 allocate(o2(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
226 allocate(tcucns(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
233 radius_cloud(i,j,l)=spval
234 radius_ice(i,j,l)=spval
235 radius_snow(i,j,l)=spval
236 refl_10cm(i,j,l)=spval
245 if (me == 0) print *,
' d3d_on=',d3d_on
247 allocate(vdifftt(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
249 allocate(vdiffmois(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
250 allocate(dconvmois(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
251 allocate(sconvmois(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
252 allocate(nradtt(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
253 allocate(o3vdiff(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
254 allocate(o3prod(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
255 allocate(o3tndy(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
256 allocate(mwpv(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
257 allocate(unknown(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
258 allocate(vdiffzacce(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
259 allocate(zgdrag(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
260 allocate(cnvctummixing(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
261 allocate(vdiffmacce(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
262 allocate(mgdrag(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
263 allocate(cnvctvmmixing(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
264 allocate(ncnvctcfrac(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
265 allocate(cnvctumflx(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
266 allocate(cnvctdmflx(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
267 allocate(cnvctdetmflx(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
268 allocate(cnvctzgdrag(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
269 allocate(cnvctmgdrag(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
276 vdiffmois(i,j,l)=spval
277 dconvmois(i,j,l)=spval
278 sconvmois(i,j,l)=spval
285 vdiffzacce(i,j,l)=spval
287 cnvctummixing(i,j,l)=spval
288 vdiffmacce(i,j,l)=spval
290 cnvctvmmixing(i,j,l)=spval
291 ncnvctcfrac(i,j,l)=spval
292 cnvctumflx(i,j,l)=spval
293 cnvctdmflx(i,j,l)=spval
294 cnvctdetmflx(i,j,l)=spval
295 cnvctzgdrag(i,j,l)=spval
296 cnvctmgdrag(i,j,l)=spval
302 allocate(htm(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
303 allocate(vtm(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
305 allocate(icing_gfip(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
306 allocate(icing_gfis(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
314 icing_gfip(i,j,l)=spval
315 icing_gfis(i,j,l)=spval
322 allocate(catedr(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
323 allocate(mwt(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
324 allocate(gtg(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
325 allocate(cit(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
342 allocate(smc(ista_2l:iend_2u,jsta_2l:jend_2u,nsoil))
343 allocate(stc(ista_2l:iend_2u,jsta_2l:jend_2u,nsoil))
344 allocate(sh2o(ista_2l:iend_2u,jsta_2l:jend_2u,nsoil))
345 allocate(sldpth(nsoil))
346 allocate(rtdpth(nsoil))
347 allocate(sllevel(nsoil))
368 allocate(cape(ista_2l:iend_2u,jsta_2l:jend_2u))
369 allocate(cin(ista_2l:iend_2u,jsta_2l:jend_2u))
370 allocate(ifi_apcp(ista_2l:iend_2u,jsta_2l:jend_2u))
372 allocate(wspd10max(ista_2l:iend_2u,jsta_2l:jend_2u))
373 allocate(w_up_max(ista_2l:iend_2u,jsta_2l:jend_2u))
374 allocate(w_dn_max(ista_2l:iend_2u,jsta_2l:jend_2u))
375 allocate(w_mean(ista_2l:iend_2u,jsta_2l:jend_2u))
376 allocate(refd_max(ista_2l:iend_2u,jsta_2l:jend_2u))
377 allocate(prate_max(ista_2l:iend_2u,jsta_2l:jend_2u))
378 allocate(fprate_max(ista_2l:iend_2u,jsta_2l:jend_2u))
379 allocate(up_heli_max(ista_2l:iend_2u,jsta_2l:jend_2u))
380 allocate(up_heli_max16(ista_2l:iend_2u,jsta_2l:jend_2u))
381 allocate(up_heli_min(ista_2l:iend_2u,jsta_2l:jend_2u))
382 allocate(up_heli_min16(ista_2l:iend_2u,jsta_2l:jend_2u))
383 allocate(up_heli_max02(ista_2l:iend_2u,jsta_2l:jend_2u))
384 allocate(up_heli_min02(ista_2l:iend_2u,jsta_2l:jend_2u))
385 allocate(up_heli_max03(ista_2l:iend_2u,jsta_2l:jend_2u))
386 allocate(up_heli_min03(ista_2l:iend_2u,jsta_2l:jend_2u))
400 fprate_max(i,j)=spval
401 up_heli_max(i,j)=spval
402 up_heli_max16(i,j)=spval
403 up_heli_min(i,j)=spval
404 up_heli_min16(i,j)=spval
405 up_heli_max02(i,j)=spval
406 up_heli_min02(i,j)=spval
407 up_heli_max03(i,j)=spval
408 up_heli_min03(i,j)=spval
411 allocate(rel_vort_max(ista_2l:iend_2u,jsta_2l:jend_2u))
412 allocate(rel_vort_max01(ista_2l:iend_2u,jsta_2l:jend_2u))
413 allocate(rel_vort_maxhy1(ista_2l:iend_2u,jsta_2l:jend_2u))
414 allocate(wspd10umax(ista_2l:iend_2u,jsta_2l:jend_2u))
415 allocate(wspd10vmax(ista_2l:iend_2u,jsta_2l:jend_2u))
416 allocate(refdm10c_max(ista_2l:iend_2u,jsta_2l:jend_2u))
417 allocate(hail_max2d(ista_2l:iend_2u,jsta_2l:jend_2u))
418 allocate(hail_maxk1(ista_2l:iend_2u,jsta_2l:jend_2u))
419 allocate(hail_maxhailcast(ista_2l:iend_2u,jsta_2l:jend_2u))
420 allocate(grpl_max(ista_2l:iend_2u,jsta_2l:jend_2u))
421 allocate(up_heli(ista_2l:iend_2u,jsta_2l:jend_2u))
422 allocate(up_heli16(ista_2l:iend_2u,jsta_2l:jend_2u))
423 allocate(ltg1_max(ista_2l:iend_2u,jsta_2l:jend_2u))
424 allocate(ltg2_max(ista_2l:iend_2u,jsta_2l:jend_2u))
425 allocate(ltg3_max(ista_2l:iend_2u,jsta_2l:jend_2u))
426 allocate(nci_ltg(ista_2l:iend_2u,jsta_2l:jend_2u))
427 allocate(nca_ltg(ista_2l:iend_2u,jsta_2l:jend_2u))
428 allocate(nci_wq(ista_2l:iend_2u,jsta_2l:jend_2u))
429 allocate(nca_wq(ista_2l:iend_2u,jsta_2l:jend_2u))
430 allocate(nci_refd(ista_2l:iend_2u,jsta_2l:jend_2u))
431 allocate(nca_refd(ista_2l:iend_2u,jsta_2l:jend_2u))
436 rel_vort_max(i,j)=spval
437 rel_vort_max01(i,j)=spval
438 rel_vort_maxhy1(i,j)=spval
439 wspd10umax(i,j)=spval
440 wspd10vmax(i,j)=spval
441 refdm10c_max(i,j)=spval
442 hail_max2d(i,j)=spval
443 hail_maxk1(i,j)=spval
444 hail_maxhailcast(i,j)=spval
461 allocate(ref_10cm(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
467 ref_10cm(i,j,l)=spval
471 allocate(refc_10cm(ista_2l:iend_2u,jsta_2l:jend_2u))
472 allocate(ref1km_10cm(ista_2l:iend_2u,jsta_2l:jend_2u))
473 allocate(ref4km_10cm(ista_2l:iend_2u,jsta_2l:jend_2u))
479 ref1km_10cm(i,j)=spval
480 ref4km_10cm(i,j)=spval
484 allocate(u10(ista_2l:iend_2u,jsta_2l:jend_2u))
485 allocate(v10(ista_2l:iend_2u,jsta_2l:jend_2u))
486 allocate(tshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
487 allocate(qshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
488 allocate(mrshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
489 allocate(smstav(ista_2l:iend_2u,jsta_2l:jend_2u))
490 allocate(ssroff(ista_2l:iend_2u,jsta_2l:jend_2u))
491 allocate(bgroff(ista_2l:iend_2u,jsta_2l:jend_2u))
492 allocate(vegfrc(ista_2l:iend_2u,jsta_2l:jend_2u))
493 allocate(shdmin(ista_2l:iend_2u,jsta_2l:jend_2u))
494 allocate(shdmax(ista_2l:iend_2u,jsta_2l:jend_2u))
495 allocate(lai(ista_2l:iend_2u,jsta_2l:jend_2u))
496 allocate(xlaixy(ista_2l:iend_2u,jsta_2l:jend_2u))
497 allocate(acsnow(ista_2l:iend_2u,jsta_2l:jend_2u))
498 allocate(acgraup(ista_2l:iend_2u,jsta_2l:jend_2u))
499 allocate(acfrain(ista_2l:iend_2u,jsta_2l:jend_2u))
500 allocate(acsnom(ista_2l:iend_2u,jsta_2l:jend_2u))
501 allocate(cmc(ista_2l:iend_2u,jsta_2l:jend_2u))
502 allocate(sst(ista_2l:iend_2u,jsta_2l:jend_2u))
503 allocate(qz0(ista_2l:iend_2u,jsta_2l:jend_2u))
504 allocate(thz0(ista_2l:iend_2u,jsta_2l:jend_2u))
505 allocate(uz0(ista_2l:iend_2u,jsta_2l:jend_2u))
506 allocate(vz0(ista_2l:iend_2u,jsta_2l:jend_2u))
507 allocate(qs(ista_2l:iend_2u,jsta_2l:jend_2u))
508 allocate(ths(ista_2l:iend_2u,jsta_2l:jend_2u))
509 allocate(sno(ista_2l:iend_2u,jsta_2l:jend_2u))
510 allocate(snonc(ista_2l:iend_2u,jsta_2l:jend_2u))
511 allocate(ti(ista_2l:iend_2u,jsta_2l:jend_2u))
547 allocate(u10mean(ista_2l:iend_2u,jsta_2l:jend_2u))
548 allocate(v10mean(ista_2l:iend_2u,jsta_2l:jend_2u))
549 allocate(spduv10mean(ista_2l:iend_2u,jsta_2l:jend_2u))
550 allocate(swradmean(ista_2l:iend_2u,jsta_2l:jend_2u))
551 allocate(swnormmean(ista_2l:iend_2u,jsta_2l:jend_2u))
558 spduv10mean(i,j)=spval
560 swnormmean(i,j)=spval
564 allocate(snoavg(ista_2l:iend_2u,jsta_2l:jend_2u))
565 allocate(psfcavg(ista_2l:iend_2u,jsta_2l:jend_2u))
566 allocate(t10m(ista_2l:iend_2u,jsta_2l:jend_2u))
567 allocate(t10avg(ista_2l:iend_2u,jsta_2l:jend_2u))
568 allocate(akmsavg(ista_2l:iend_2u,jsta_2l:jend_2u))
569 allocate(akhsavg(ista_2l:iend_2u,jsta_2l:jend_2u))
570 allocate(u10max(ista_2l:iend_2u,jsta_2l:jend_2u))
571 allocate(v10max(ista_2l:iend_2u,jsta_2l:jend_2u))
572 allocate(u10h(ista_2l:iend_2u,jsta_2l:jend_2u))
573 allocate(v10h(ista_2l:iend_2u,jsta_2l:jend_2u))
591 allocate(akms(ista_2l:iend_2u,jsta_2l:jend_2u))
592 allocate(akhs(ista_2l:iend_2u,jsta_2l:jend_2u))
593 allocate(cuprec(ista_2l:iend_2u,jsta_2l:jend_2u))
594 allocate(acprec(ista_2l:iend_2u,jsta_2l:jend_2u))
595 allocate(ancprc(ista_2l:iend_2u,jsta_2l:jend_2u))
596 allocate(cuppt(ista_2l:iend_2u,jsta_2l:jend_2u))
610 allocate(rainc_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
611 allocate(rainc_bucket1(ista_2l:iend_2u,jsta_2l:jend_2u))
612 allocate(rainnc_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
613 allocate(rainnc_bucket1(ista_2l:iend_2u,jsta_2l:jend_2u))
614 allocate(pcp_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
615 allocate(pcp_bucket1(ista_2l:iend_2u,jsta_2l:jend_2u))
616 allocate(snow_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
617 allocate(snow_bucket1(ista_2l:iend_2u,jsta_2l:jend_2u))
618 allocate(graup_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
619 allocate(graup_bucket1(ista_2l:iend_2u,jsta_2l:jend_2u))
620 allocate(frzrn_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
621 allocate(snow_acm(ista_2l:iend_2u,jsta_2l:jend_2u))
622 allocate(snow_bkt(ista_2l:iend_2u,jsta_2l:jend_2u))
623 allocate(qrmax(ista_2l:iend_2u,jsta_2l:jend_2u))
624 allocate(tmax(ista_2l:iend_2u,jsta_2l:jend_2u))
625 allocate(snownc(ista_2l:iend_2u,jsta_2l:jend_2u))
626 allocate(graupelnc(ista_2l:iend_2u,jsta_2l:jend_2u))
627 allocate(tsnow(ista_2l:iend_2u,jsta_2l:jend_2u))
628 allocate(qvg(ista_2l:iend_2u,jsta_2l:jend_2u))
629 allocate(qv2m(ista_2l:iend_2u,jsta_2l:jend_2u))
630 allocate(qvl1(ista_2l:iend_2u,jsta_2l:jend_2u))
631 allocate(snfden(ista_2l:iend_2u,jsta_2l:jend_2u))
632 allocate(sndepac(ista_2l:iend_2u,jsta_2l:jend_2u))
633 allocate(mean_frp(ista_2l:iend_2u,jsta_2l:jend_2u))
634 allocate(hwp(ista_2l:iend_2u,jsta_2l:jend_2u))
635 allocate(smoke_ave(ista_2l:iend_2u,jsta_2l:jend_2u))
636 allocate(dust_ave(ista_2l:iend_2u,jsta_2l:jend_2u))
637 allocate(coarsepm_ave(ista_2l:iend_2u,jsta_2l:jend_2u))
642 rainc_bucket(i,j)=spval
643 rainc_bucket1(i,j)=spval
644 rainnc_bucket(i,j)=spval
645 rainnc_bucket1(i,j)=spval
646 pcp_bucket(i,j)=spval
647 pcp_bucket1(i,j)=spval
648 snow_bucket(i,j)=spval
649 snow_bucket1(i,j)=spval
650 graup_bucket(i,j)=spval
651 graup_bucket1(i,j)=spval
666 coarsepm_ave(i,j)=spval
669 allocate(smoke(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_sm))
670 allocate(fv3dust(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_sm))
671 allocate(coarsepm(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_sm))
672 allocate(ebb(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_sm))
679 fv3dust(i,j,l,k)=spval
680 coarsepm(i,j,l,k)=spval
687 allocate(rswin(ista_2l:iend_2u,jsta_2l:jend_2u))
688 allocate(swddni(ista_2l:iend_2u,jsta_2l:jend_2u))
689 allocate(swddif(ista_2l:iend_2u,jsta_2l:jend_2u))
690 allocate(swdnbc(ista_2l:iend_2u,jsta_2l:jend_2u))
691 allocate(swddnic(ista_2l:iend_2u,jsta_2l:jend_2u))
692 allocate(swddifc(ista_2l:iend_2u,jsta_2l:jend_2u))
693 allocate(swupbc(ista_2l:iend_2u,jsta_2l:jend_2u))
694 allocate(swupt(ista_2l:iend_2u,jsta_2l:jend_2u))
695 allocate(taod5502d(ista_2l:iend_2u,jsta_2l:jend_2u))
696 allocate(aerasy2d(ista_2l:iend_2u,jsta_2l:jend_2u))
697 allocate(aerssa2d(ista_2l:iend_2u,jsta_2l:jend_2u))
698 allocate(lwp(ista_2l:iend_2u,jsta_2l:jend_2u))
699 allocate(iwp(ista_2l:iend_2u,jsta_2l:jend_2u))
700 allocate(rlwin(ista_2l:iend_2u,jsta_2l:jend_2u))
701 allocate(lwdnbc(ista_2l:iend_2u,jsta_2l:jend_2u))
702 allocate(lwupbc(ista_2l:iend_2u,jsta_2l:jend_2u))
703 allocate(rlwtoa(ista_2l:iend_2u,jsta_2l:jend_2u))
704 allocate(rswtoa(ista_2l:iend_2u,jsta_2l:jend_2u))
729 allocate(tg(ista_2l:iend_2u,jsta_2l:jend_2u))
730 allocate(sfcshx(ista_2l:iend_2u,jsta_2l:jend_2u))
731 allocate(sfclhx(ista_2l:iend_2u,jsta_2l:jend_2u))
732 allocate(fis(ista_2l:iend_2u,jsta_2l:jend_2u))
733 allocate(t500(ista_2l:iend_2u,jsta_2l:jend_2u))
734 allocate(t700(ista_2l:iend_2u,jsta_2l:jend_2u))
735 allocate(z500(ista_2l:iend_2u,jsta_2l:jend_2u))
736 allocate(z700(ista_2l:iend_2u,jsta_2l:jend_2u))
737 allocate(teql(ista_2l:iend_2u,jsta_2l:jend_2u))
738 allocate(ieql(ista_2l:iend_2u,jsta_2l:jend_2u))
739 allocate(cfracl(ista_2l:iend_2u,jsta_2l:jend_2u))
740 allocate(cfracm(ista_2l:iend_2u,jsta_2l:jend_2u))
741 allocate(cfrach(ista_2l:iend_2u,jsta_2l:jend_2u))
742 allocate(acfrst(ista_2l:iend_2u,jsta_2l:jend_2u))
743 allocate(acfrcv(ista_2l:iend_2u,jsta_2l:jend_2u))
744 allocate(hbot(ista_2l:iend_2u,jsta_2l:jend_2u))
745 allocate(htop(ista_2l:iend_2u,jsta_2l:jend_2u))
746 allocate(aswin(ista_2l:iend_2u,jsta_2l:jend_2u))
747 allocate(alwin(ista_2l:iend_2u,jsta_2l:jend_2u))
748 allocate(aswout(ista_2l:iend_2u,jsta_2l:jend_2u))
749 allocate(alwout(ista_2l:iend_2u,jsta_2l:jend_2u))
750 allocate(aswtoa(ista_2l:iend_2u,jsta_2l:jend_2u))
751 allocate(alwtoa(ista_2l:iend_2u,jsta_2l:jend_2u))
780 allocate(czen(ista_2l:iend_2u,jsta_2l:jend_2u))
781 allocate(czmean(ista_2l:iend_2u,jsta_2l:jend_2u))
782 allocate(sigt4(ista_2l:iend_2u,jsta_2l:jend_2u))
783 allocate(rswout(ista_2l:iend_2u,jsta_2l:jend_2u))
784 allocate(radot(ista_2l:iend_2u,jsta_2l:jend_2u))
785 allocate(ncfrst(ista_2l:iend_2u,jsta_2l:jend_2u))
786 allocate(ncfrcv(ista_2l:iend_2u,jsta_2l:jend_2u))
787 allocate(smstot(ista_2l:iend_2u,jsta_2l:jend_2u))
788 allocate(pctsno(ista_2l:iend_2u,jsta_2l:jend_2u))
789 allocate(pshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
790 allocate(th10(ista_2l:iend_2u,jsta_2l:jend_2u))
791 allocate(q10(ista_2l:iend_2u,jsta_2l:jend_2u))
792 allocate(sr(ista_2l:iend_2u,jsta_2l:jend_2u))
793 allocate(prec(ista_2l:iend_2u,jsta_2l:jend_2u))
794 allocate(subshx(ista_2l:iend_2u,jsta_2l:jend_2u))
795 allocate(snopcx(ista_2l:iend_2u,jsta_2l:jend_2u))
796 allocate(sfcuvx(ista_2l:iend_2u,jsta_2l:jend_2u))
797 allocate(sfcevp(ista_2l:iend_2u,jsta_2l:jend_2u))
798 allocate(potevp(ista_2l:iend_2u,jsta_2l:jend_2u))
799 allocate(z0(ista_2l:iend_2u,jsta_2l:jend_2u))
800 allocate(ustar(ista_2l:iend_2u,jsta_2l:jend_2u))
801 allocate(pblh(ista_2l:iend_2u,jsta_2l:jend_2u))
802 allocate(pblhgust(ista_2l:iend_2u,jsta_2l:jend_2u))
803 allocate(mixht(ista_2l:iend_2u,jsta_2l:jend_2u))
804 allocate(twbs(ista_2l:iend_2u,jsta_2l:jend_2u))
805 allocate(qwbs(ista_2l:iend_2u,jsta_2l:jend_2u))
838 allocate(sfcexc(ista_2l:iend_2u,jsta_2l:jend_2u))
839 allocate(grnflx(ista_2l:iend_2u,jsta_2l:jend_2u))
840 allocate(soiltb(ista_2l:iend_2u,jsta_2l:jend_2u))
841 allocate(z1000(ista_2l:iend_2u,jsta_2l:jend_2u))
842 allocate(slp(ista_2l:iend_2u,jsta_2l:jend_2u))
843 allocate(pslp(ista_2l:iend_2u,jsta_2l:jend_2u))
844 allocate(f(ista_2l:iend_2u,jsta_2l:jend_2u))
845 allocate(albedo(ista_2l:iend_2u,jsta_2l:jend_2u))
846 allocate(albase(ista_2l:iend_2u,jsta_2l:jend_2u))
847 allocate(cldfra(ista_2l:iend_2u,jsta_2l:jend_2u))
848 allocate(cprate(ista_2l:iend_2u,jsta_2l:jend_2u))
849 allocate(cnvcfr(ista_2l:iend_2u,jsta_2l:jend_2u))
850 allocate(ivgtyp(ista_2l:iend_2u,jsta_2l:jend_2u))
851 allocate(isltyp(ista_2l:iend_2u,jsta_2l:jend_2u))
852 allocate(hbotd(ista_2l:iend_2u,jsta_2l:jend_2u))
853 allocate(htopd(ista_2l:iend_2u,jsta_2l:jend_2u))
854 allocate(hbots(ista_2l:iend_2u,jsta_2l:jend_2u))
855 allocate(htops(ista_2l:iend_2u,jsta_2l:jend_2u))
856 allocate(cldefi(ista_2l:iend_2u,jsta_2l:jend_2u))
857 allocate(islope(ista_2l:iend_2u,jsta_2l:jend_2u))
858 allocate(si(ista_2l:iend_2u,jsta_2l:jend_2u))
859 allocate(lspa(ista_2l:iend_2u,jsta_2l:jend_2u))
860 allocate(rswinc(ista_2l:iend_2u,jsta_2l:jend_2u))
861 allocate(vis(ista_2l:iend_2u,jsta_2l:jend_2u))
862 allocate(pd(ista_2l:iend_2u,jsta_2l:jend_2u))
863 allocate(mxsnal(ista_2l:iend_2u,jsta_2l:jend_2u))
864 allocate(epsr(ista_2l:iend_2u,jsta_2l:jend_2u))
899 allocate(sfcux(ista_2l:iend_2u,jsta_2l:jend_2u))
900 allocate(sfcvx(ista_2l:iend_2u,jsta_2l:jend_2u))
901 allocate(sfcuxi(ista_2l:iend_2u,jsta_2l:jend_2u))
902 allocate(sfcvxi(ista_2l:iend_2u,jsta_2l:jend_2u))
903 allocate(avgalbedo(ista_2l:iend_2u,jsta_2l:jend_2u))
904 allocate(avgcprate(ista_2l:iend_2u,jsta_2l:jend_2u))
905 allocate(avgprec(ista_2l:iend_2u,jsta_2l:jend_2u))
906 allocate(avgprec_cont(ista_2l:iend_2u,jsta_2l:jend_2u))
907 allocate(avgcprate_cont(ista_2l:iend_2u,jsta_2l:jend_2u))
908 allocate(ptop(ista_2l:iend_2u,jsta_2l:jend_2u))
909 allocate(pbot(ista_2l:iend_2u,jsta_2l:jend_2u))
910 allocate(avgcfrach(ista_2l:iend_2u,jsta_2l:jend_2u))
911 allocate(avgcfracm(ista_2l:iend_2u,jsta_2l:jend_2u))
912 allocate(avgcfracl(ista_2l:iend_2u,jsta_2l:jend_2u))
913 allocate(avgtcdc(ista_2l:iend_2u,jsta_2l:jend_2u))
914 allocate(auvbin(ista_2l:iend_2u,jsta_2l:jend_2u))
915 allocate(auvbinc(ista_2l:iend_2u,jsta_2l:jend_2u))
916 allocate(ptopl(ista_2l:iend_2u,jsta_2l:jend_2u))
917 allocate(pbotl(ista_2l:iend_2u,jsta_2l:jend_2u))
918 allocate(ttopl(ista_2l:iend_2u,jsta_2l:jend_2u))
919 allocate(ptopm(ista_2l:iend_2u,jsta_2l:jend_2u))
920 allocate(pbotm(ista_2l:iend_2u,jsta_2l:jend_2u))
921 allocate(ttopm(ista_2l:iend_2u,jsta_2l:jend_2u))
922 allocate(ptoph(ista_2l:iend_2u,jsta_2l:jend_2u))
923 allocate(pboth(ista_2l:iend_2u,jsta_2l:jend_2u))
924 allocate(ttoph(ista_2l:iend_2u,jsta_2l:jend_2u))
925 allocate(sfcugs(ista_2l:iend_2u,jsta_2l:jend_2u))
926 allocate(sfcvgs(ista_2l:iend_2u,jsta_2l:jend_2u))
927 allocate(pblcfr(ista_2l:iend_2u,jsta_2l:jend_2u))
928 allocate(cldwork(ista_2l:iend_2u,jsta_2l:jend_2u))
929 allocate(gtaux(ista_2l:iend_2u,jsta_2l:jend_2u))
930 allocate(gtauy(ista_2l:iend_2u,jsta_2l:jend_2u))
931 allocate(cd10(ista_2l:iend_2u,jsta_2l:jend_2u))
932 allocate(ch10(ista_2l:iend_2u,jsta_2l:jend_2u))
933 allocate(mdltaux(ista_2l:iend_2u,jsta_2l:jend_2u))
934 allocate(mdltauy(ista_2l:iend_2u,jsta_2l:jend_2u))
935 allocate(runoff(ista_2l:iend_2u,jsta_2l:jend_2u))
947 avgprec_cont(i,j)=spval
948 avgcprate_cont(i,j)=spval
979 allocate(maxtshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
980 allocate(mintshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
981 allocate(maxrhshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
982 allocate(minrhshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
983 allocate(maxqshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
984 allocate(minqshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
985 allocate(dzice(ista_2l:iend_2u,jsta_2l:jend_2u))
986 allocate(alwinc(ista_2l:iend_2u,jsta_2l:jend_2u))
987 allocate(alwoutc(ista_2l:iend_2u,jsta_2l:jend_2u))
988 allocate(alwtoac(ista_2l:iend_2u,jsta_2l:jend_2u))
989 allocate(aswinc(ista_2l:iend_2u,jsta_2l:jend_2u))
990 allocate(aswoutc(ista_2l:iend_2u,jsta_2l:jend_2u))
991 allocate(aswtoac(ista_2l:iend_2u,jsta_2l:jend_2u))
992 allocate(aswintoa(ista_2l:iend_2u,jsta_2l:jend_2u))
993 allocate(smcwlt(ista_2l:iend_2u,jsta_2l:jend_2u))
994 allocate(suntime(ista_2l:iend_2u,jsta_2l:jend_2u))
995 allocate(fieldcapa(ista_2l:iend_2u,jsta_2l:jend_2u))
996 allocate(avisbeamswin(ista_2l:iend_2u,jsta_2l:jend_2u))
997 allocate(avisdiffswin(ista_2l:iend_2u,jsta_2l:jend_2u))
998 allocate(airbeamswin(ista_2l:iend_2u,jsta_2l:jend_2u))
999 allocate(airdiffswin(ista_2l:iend_2u,jsta_2l:jend_2u))
1000 allocate(snowfall(ista_2l:iend_2u,jsta_2l:jend_2u))
1001 allocate(acond(ista_2l:iend_2u,jsta_2l:jend_2u))
1002 allocate(edir(ista_2l:iend_2u,jsta_2l:jend_2u))
1003 allocate(ecan(ista_2l:iend_2u,jsta_2l:jend_2u))
1004 allocate(etrans(ista_2l:iend_2u,jsta_2l:jend_2u))
1005 allocate(esnow(ista_2l:iend_2u,jsta_2l:jend_2u))
1006 allocate(avgedir(ista_2l:iend_2u,jsta_2l:jend_2u))
1007 allocate(avgecan(ista_2l:iend_2u,jsta_2l:jend_2u))
1008 allocate(avgetrans(ista_2l:iend_2u,jsta_2l:jend_2u))
1009 allocate(avgesnow(ista_2l:iend_2u,jsta_2l:jend_2u))
1010 allocate(avgpotevp(ista_2l:iend_2u,jsta_2l:jend_2u))
1011 allocate(aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
1012 allocate(du_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
1013 allocate(ss_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
1014 allocate(su_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
1015 allocate(oc_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
1016 allocate(bc_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
1017 allocate(landfrac(ista_2l:iend_2u,jsta_2l:jend_2u))
1018 allocate(paha(ista_2l:iend_2u,jsta_2l:jend_2u))
1019 allocate(pahi(ista_2l:iend_2u,jsta_2l:jend_2u))
1020 allocate(tecan(ista_2l:iend_2u,jsta_2l:jend_2u))
1021 allocate(tetran(ista_2l:iend_2u,jsta_2l:jend_2u))
1022 allocate(tedir(ista_2l:iend_2u,jsta_2l:jend_2u))
1023 allocate(twa(ista_2l:iend_2u,jsta_2l:jend_2u))
1024 allocate(fdnsst(ista_2l:iend_2u,jsta_2l:jend_2u))
1025 allocate(pwat(ista_2l:iend_2u,jsta_2l:jend_2u))
1028 do j=jsta_2l,jend_2u
1029 do i=ista_2l,iend_2u
1030 maxtshltr(i,j)=spval
1031 mintshltr(i,j)=spval
1032 maxrhshltr(i,j)=spval
1033 minrhshltr(i,j)=spval
1034 maxqshltr(i,j)=spval
1035 minqshltr(i,j)=spval
1046 fieldcapa(i,j)=spval
1047 avisbeamswin(i,j)=spval
1048 avisdiffswin(i,j)=spval
1049 airbeamswin(i,j)=spval
1050 airdiffswin(i,j)=spval
1059 avgetrans(i,j)=spval
1061 avgpotevp(i,j)=spval
1063 du_aod550(i,j)=spval
1064 ss_aod550(i,j)=spval
1065 su_aod550(i,j)=spval
1066 oc_aod550(i,j)=spval
1067 bc_aod550(i,j)=spval
1082 allocate(hbm2(ista_2l:iend_2u,jsta_2l:jend_2u))
1083 allocate(sm(ista_2l:iend_2u,jsta_2l:jend_2u))
1084 allocate(sice(ista_2l:iend_2u,jsta_2l:jend_2u))
1085 allocate(lmh(ista_2l:iend_2u,jsta_2l:jend_2u))
1086 allocate(lmv(ista_2l:iend_2u,jsta_2l:jend_2u))
1087 allocate(gdlat(ista_2l:iend_2u,jsta_2l:jend_2u))
1088 allocate(gdlon(ista_2l:iend_2u,jsta_2l:jend_2u))
1089 allocate(dx(ista_2l:iend_2u,jsta_2l:jend_2u))
1090 allocate(dy(ista_2l:iend_2u,jsta_2l:jend_2u))
1093 do j=jsta_2l,jend_2u
1094 do i=ista_2l,iend_2u
1109 if (gocart_on .or.gccpp_on .or. nasa_on)
then
1113 allocate(dust(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_du))
1114 allocate(salt(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_ss))
1115 allocate(soot(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_bc))
1116 allocate(waso(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_oc))
1117 allocate(suso(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_su))
1119 allocate(no3(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_no3))
1120 allocate(nh4(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_nh4))
1122 allocate(pp25(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_su))
1123 allocate(pp10(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_su))
1128 do j=jsta_2l,jend_2u
1129 do i=ista_2l,iend_2u
1138 do j=jsta_2l,jend_2u
1139 do i=ista_2l,iend_2u
1148 do j=jsta_2l,jend_2u
1149 do i=ista_2l,iend_2u
1158 do j=jsta_2l,jend_2u
1159 do i=ista_2l,iend_2u
1169 do j=jsta_2l,jend_2u
1170 do i=ista_2l,iend_2u
1180 do j=jsta_2l,jend_2u
1181 do i=ista_2l,iend_2u
1192 do j=jsta_2l,jend_2u
1193 do i=ista_2l,iend_2u
1202 allocate(ext(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1203 allocate(asy(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1204 allocate(ssa(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1205 allocate(sca(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1209 do j=jsta_2l,jend_2u
1210 do i=ista_2l,iend_2u
1218 if ( d2d_chem )
then
1219 allocate(duem(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_du))
1220 allocate(dusd(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_du))
1221 allocate(dudp(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_du))
1222 allocate(duwt(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_du))
1223 allocate(dusv(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_du))
1224 allocate(suem(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_su))
1225 allocate(susd(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_su))
1226 allocate(sudp(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_su))
1227 allocate(suwt(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_su))
1228 allocate(ocem(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_oc))
1229 allocate(ocsd(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_oc))
1230 allocate(ocdp(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_oc))
1231 allocate(ocwt(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_oc))
1232 allocate(ocsv(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_oc))
1233 allocate(bcem(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_bc))
1234 allocate(bcsd(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_bc))
1235 allocate(bcdp(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_bc))
1236 allocate(bcwt(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_bc))
1237 allocate(bcsv(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_bc))
1238 allocate(ssem(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_ss))
1239 allocate(sssd(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_ss))
1240 allocate(ssdp(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_ss))
1241 allocate(sswt(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_ss))
1242 allocate(sssv(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_ss))
1246 do j=jsta_2l,jend_2u
1247 do i=ista_2l,iend_2u
1258 do j=jsta_2l,jend_2u
1259 do i=ista_2l,iend_2u
1269 do j=jsta_2l,jend_2u
1270 do i=ista_2l,iend_2u
1281 do j=jsta_2l,jend_2u
1282 do i=ista_2l,iend_2u
1293 do j=jsta_2l,jend_2u
1294 do i=ista_2l,iend_2u
1304 allocate(rhomid(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1308 do j=jsta_2l,jend_2u
1309 do i=ista_2l,iend_2u
1315 allocate(dusmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1316 allocate(ducmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1317 allocate(dusmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1318 allocate(ducmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1319 allocate(susmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1320 allocate(sucmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1321 allocate(susmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1322 allocate(sucmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1323 allocate(ocsmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1324 allocate(occmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1325 allocate(ocsmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1326 allocate(occmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1327 allocate(bcsmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1328 allocate(bccmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1329 allocate(bcsmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1330 allocate(bccmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1331 allocate(sssmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1332 allocate(sscmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1333 allocate(sssmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1334 allocate(sscmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1335 allocate(dustcb(ista_2l:iend_2u,jsta_2l:jend_2u))
1337 allocate(no3cb(ista_2l:iend_2u,jsta_2l:jend_2u))
1338 allocate(nh4cb(ista_2l:iend_2u,jsta_2l:jend_2u))
1340 allocate(occb(ista_2l:iend_2u,jsta_2l:jend_2u))
1341 allocate(bccb(ista_2l:iend_2u,jsta_2l:jend_2u))
1342 allocate(sulfcb(ista_2l:iend_2u,jsta_2l:jend_2u))
1343 allocate(pp25cb(ista_2l:iend_2u,jsta_2l:jend_2u))
1344 allocate(pp10cb(ista_2l:iend_2u,jsta_2l:jend_2u))
1345 allocate(sscb(ista_2l:iend_2u,jsta_2l:jend_2u))
1346 allocate(dustallcb(ista_2l:iend_2u,jsta_2l:jend_2u))
1347 allocate(ssallcb(ista_2l:iend_2u,jsta_2l:jend_2u))
1348 allocate(dustpm(ista_2l:iend_2u,jsta_2l:jend_2u))
1349 allocate(dustpm10(ista_2l:iend_2u,jsta_2l:jend_2u))
1350 allocate(sspm(ista_2l:iend_2u,jsta_2l:jend_2u))
1351 allocate(maod(ista_2l:iend_2u,jsta_2l:jend_2u))
1354 do j=jsta_2l,jend_2u
1355 do i=ista_2l,iend_2u
1358 dusmass25(i,j)=spval
1359 ducmass25(i,j)=spval
1362 susmass25(i,j)=spval
1363 sucmass25(i,j)=spval
1366 ocsmass25(i,j)=spval
1367 occmass25(i,j)=spval
1370 bcsmass25(i,j)=spval
1371 bccmass25(i,j)=spval
1374 sssmass25(i,j)=spval
1375 sscmass25(i,j)=spval
1387 dustallcb(i,j)=spval
1397 allocate(acswupt(ista_2l:iend_2u,jsta_2l:jend_2u))
1398 allocate(swdnt(ista_2l:iend_2u,jsta_2l:jend_2u))
1399 allocate(acswdnt(ista_2l:iend_2u,jsta_2l:jend_2u))
1402 do j=jsta_2l,jend_2u
1403 do i=ista_2l,iend_2u
1411 allocate(ddvdx(ista_2l:iend_2u,jsta_2l:jend_2u))
1412 allocate(ddudy(ista_2l:iend_2u,jsta_2l:jend_2u))
1413 allocate(uuavg(ista_2l:iend_2u,jsta_2l:jend_2u))
1416 do j=jsta_2l,jend_2u
1417 do i=ista_2l,iend_2u
1425 if (me == 0) print *,
'aqf_on= ', aqf_on
1428 allocate(avgozcon(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1429 allocate(avgpmtf(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1430 allocate(aqm_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
1435 do j=jsta_2l,jend_2u
1436 do i=ista_2l,iend_2u
1437 avgozcon(i,j,l)=spval
1438 avgpmtf(i,j,l)=spval
1445 do j=jsta_2l,jend_2u
1446 do i=ista_2l,iend_2u
1447 aqm_aod550(i,j)=spval