31 SUBROUTINE allocate_all()
38 use upp_math, only: ddvdx, ddudy, uuavg
50 allocate(u(ista_2l:iend_2u+1,jsta_2l:jend_2u,lm))
51 allocate(v(ista_2l:iend_2u,jsta_2l:jvend_2u,lm))
52 allocate(t(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
55 allocate(q(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
57 allocate(uh(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
58 allocate(vh(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
59 allocate(wh(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
60 allocate(pmid(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
61 allocate(pmidv(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
62 allocate(pint(ista_2l:iend_2u,jsta_2l:jend_2u,lp1))
63 allocate(alpint(ista_2l:iend_2u,jsta_2l:jend_2u,lp1))
64 allocate(zmid(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
65 allocate(zint(ista_2l:iend_2u,jsta_2l:jend_2u,lp1))
67 allocate(q2(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
68 allocate(omga(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
69 allocate(dpres(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
70 allocate(t_adj(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
71 allocate(ttnd(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
72 allocate(rswtt(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
73 allocate(rlwtt(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
74 allocate(exch_h(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
75 allocate(train(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
76 allocate(tcucn(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
77 allocate(el_pbl(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
83 do i=ista_2l,iend_2u+1
134 allocate(cwm(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
135 allocate(f_ice(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
136 allocate(f_rain(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
137 allocate(f_rimef(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
138 allocate(qqw(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
139 allocate(qrimef(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
140 allocate(qqi(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
141 allocate(qqr(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
142 allocate(qqs(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
143 allocate(qqg(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
144 allocate(qqnw(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
145 allocate(qqni(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
146 allocate(qqnr(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
147 allocate(qqnwfa(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
148 allocate(qqnifa(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
149 allocate(taod5503d(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
150 allocate(aextc55(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
151 allocate(extcof55(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
152 allocate(qc_bl(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
153 allocate(cfr(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
154 allocate(cfr_raw(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
155 allocate(dbz(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
156 allocate(dbzr(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
157 allocate(dbzi(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
158 allocate(dbzc(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
159 allocate(mcvg(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
160 allocate(nlice(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
181 taod5503d(i,j,l)=spval
197 allocate(nrain(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
198 allocate(radius_cloud(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
199 allocate(radius_ice(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
200 allocate(radius_snow(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
203 allocate(refl_10cm(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
205 allocate(o3(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
206 allocate(o(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
207 allocate(o2(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
208 allocate(tcucns(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
215 radius_cloud(i,j,l)=spval
216 radius_ice(i,j,l)=spval
217 radius_snow(i,j,l)=spval
218 refl_10cm(i,j,l)=spval
227 if (me == 0) print *,
' d3d_on=',d3d_on
229 allocate(vdifftt(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
231 allocate(vdiffmois(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
232 allocate(dconvmois(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
233 allocate(sconvmois(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
234 allocate(nradtt(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
235 allocate(o3vdiff(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
236 allocate(o3prod(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
237 allocate(o3tndy(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
238 allocate(mwpv(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
239 allocate(unknown(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
240 allocate(vdiffzacce(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
241 allocate(zgdrag(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
242 allocate(cnvctummixing(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
243 allocate(vdiffmacce(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
244 allocate(mgdrag(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
245 allocate(cnvctvmmixing(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
246 allocate(ncnvctcfrac(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
247 allocate(cnvctumflx(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
248 allocate(cnvctdmflx(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
249 allocate(cnvctdetmflx(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
250 allocate(cnvctzgdrag(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
251 allocate(cnvctmgdrag(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
258 vdiffmois(i,j,l)=spval
259 dconvmois(i,j,l)=spval
260 sconvmois(i,j,l)=spval
267 vdiffzacce(i,j,l)=spval
269 cnvctummixing(i,j,l)=spval
270 vdiffmacce(i,j,l)=spval
272 cnvctvmmixing(i,j,l)=spval
273 ncnvctcfrac(i,j,l)=spval
274 cnvctumflx(i,j,l)=spval
275 cnvctdmflx(i,j,l)=spval
276 cnvctdetmflx(i,j,l)=spval
277 cnvctzgdrag(i,j,l)=spval
278 cnvctmgdrag(i,j,l)=spval
284 allocate(htm(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
285 allocate(vtm(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
287 allocate(icing_gfip(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
288 allocate(icing_gfis(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
291 allocate(catedr(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
292 allocate(mwt(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
293 allocate(gtg(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
301 icing_gfip(i,j,l)=spval
302 icing_gfis(i,j,l)=spval
312 allocate(smc(ista_2l:iend_2u,jsta_2l:jend_2u,nsoil))
313 allocate(stc(ista_2l:iend_2u,jsta_2l:jend_2u,nsoil))
314 allocate(sh2o(ista_2l:iend_2u,jsta_2l:jend_2u,nsoil))
315 allocate(sldpth(nsoil))
316 allocate(rtdpth(nsoil))
317 allocate(sllevel(nsoil))
339 allocate(wspd10max(ista_2l:iend_2u,jsta_2l:jend_2u))
340 allocate(w_up_max(ista_2l:iend_2u,jsta_2l:jend_2u))
341 allocate(w_dn_max(ista_2l:iend_2u,jsta_2l:jend_2u))
342 allocate(w_mean(ista_2l:iend_2u,jsta_2l:jend_2u))
343 allocate(refd_max(ista_2l:iend_2u,jsta_2l:jend_2u))
344 allocate(prate_max(ista_2l:iend_2u,jsta_2l:jend_2u))
345 allocate(fprate_max(ista_2l:iend_2u,jsta_2l:jend_2u))
346 allocate(up_heli_max(ista_2l:iend_2u,jsta_2l:jend_2u))
347 allocate(up_heli_max16(ista_2l:iend_2u,jsta_2l:jend_2u))
348 allocate(up_heli_min(ista_2l:iend_2u,jsta_2l:jend_2u))
349 allocate(up_heli_min16(ista_2l:iend_2u,jsta_2l:jend_2u))
350 allocate(up_heli_max02(ista_2l:iend_2u,jsta_2l:jend_2u))
351 allocate(up_heli_min02(ista_2l:iend_2u,jsta_2l:jend_2u))
352 allocate(up_heli_max03(ista_2l:iend_2u,jsta_2l:jend_2u))
353 allocate(up_heli_min03(ista_2l:iend_2u,jsta_2l:jend_2u))
364 fprate_max(i,j)=spval
365 up_heli_max(i,j)=spval
366 up_heli_max16(i,j)=spval
367 up_heli_min(i,j)=spval
368 up_heli_min16(i,j)=spval
369 up_heli_max02(i,j)=spval
370 up_heli_min02(i,j)=spval
371 up_heli_max03(i,j)=spval
372 up_heli_min03(i,j)=spval
375 allocate(rel_vort_max(ista_2l:iend_2u,jsta_2l:jend_2u))
376 allocate(rel_vort_max01(ista_2l:iend_2u,jsta_2l:jend_2u))
377 allocate(rel_vort_maxhy1(ista_2l:iend_2u,jsta_2l:jend_2u))
378 allocate(wspd10umax(ista_2l:iend_2u,jsta_2l:jend_2u))
379 allocate(wspd10vmax(ista_2l:iend_2u,jsta_2l:jend_2u))
380 allocate(refdm10c_max(ista_2l:iend_2u,jsta_2l:jend_2u))
381 allocate(hail_max2d(ista_2l:iend_2u,jsta_2l:jend_2u))
382 allocate(hail_maxk1(ista_2l:iend_2u,jsta_2l:jend_2u))
383 allocate(hail_maxhailcast(ista_2l:iend_2u,jsta_2l:jend_2u))
384 allocate(grpl_max(ista_2l:iend_2u,jsta_2l:jend_2u))
385 allocate(up_heli(ista_2l:iend_2u,jsta_2l:jend_2u))
386 allocate(up_heli16(ista_2l:iend_2u,jsta_2l:jend_2u))
387 allocate(ltg1_max(ista_2l:iend_2u,jsta_2l:jend_2u))
388 allocate(ltg2_max(ista_2l:iend_2u,jsta_2l:jend_2u))
389 allocate(ltg3_max(ista_2l:iend_2u,jsta_2l:jend_2u))
390 allocate(nci_ltg(ista_2l:iend_2u,jsta_2l:jend_2u))
391 allocate(nca_ltg(ista_2l:iend_2u,jsta_2l:jend_2u))
392 allocate(nci_wq(ista_2l:iend_2u,jsta_2l:jend_2u))
393 allocate(nca_wq(ista_2l:iend_2u,jsta_2l:jend_2u))
394 allocate(nci_refd(ista_2l:iend_2u,jsta_2l:jend_2u))
395 allocate(nca_refd(ista_2l:iend_2u,jsta_2l:jend_2u))
400 rel_vort_max(i,j)=spval
401 rel_vort_max01(i,j)=spval
402 rel_vort_maxhy1(i,j)=spval
403 wspd10umax(i,j)=spval
404 wspd10vmax(i,j)=spval
405 refdm10c_max(i,j)=spval
406 hail_max2d(i,j)=spval
407 hail_maxk1(i,j)=spval
408 hail_maxhailcast(i,j)=spval
425 allocate(ref_10cm(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
431 ref_10cm(i,j,l)=spval
435 allocate(refc_10cm(ista_2l:iend_2u,jsta_2l:jend_2u))
436 allocate(ref1km_10cm(ista_2l:iend_2u,jsta_2l:jend_2u))
437 allocate(ref4km_10cm(ista_2l:iend_2u,jsta_2l:jend_2u))
443 ref1km_10cm(i,j)=spval
444 ref4km_10cm(i,j)=spval
448 allocate(u10(ista_2l:iend_2u,jsta_2l:jend_2u))
449 allocate(v10(ista_2l:iend_2u,jsta_2l:jend_2u))
450 allocate(tshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
451 allocate(qshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
452 allocate(mrshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
453 allocate(smstav(ista_2l:iend_2u,jsta_2l:jend_2u))
454 allocate(ssroff(ista_2l:iend_2u,jsta_2l:jend_2u))
455 allocate(bgroff(ista_2l:iend_2u,jsta_2l:jend_2u))
456 allocate(vegfrc(ista_2l:iend_2u,jsta_2l:jend_2u))
457 allocate(shdmin(ista_2l:iend_2u,jsta_2l:jend_2u))
458 allocate(shdmax(ista_2l:iend_2u,jsta_2l:jend_2u))
459 allocate(lai(ista_2l:iend_2u,jsta_2l:jend_2u))
460 allocate(acsnow(ista_2l:iend_2u,jsta_2l:jend_2u))
461 allocate(acgraup(ista_2l:iend_2u,jsta_2l:jend_2u))
462 allocate(acfrain(ista_2l:iend_2u,jsta_2l:jend_2u))
463 allocate(acsnom(ista_2l:iend_2u,jsta_2l:jend_2u))
464 allocate(cmc(ista_2l:iend_2u,jsta_2l:jend_2u))
465 allocate(sst(ista_2l:iend_2u,jsta_2l:jend_2u))
466 allocate(qz0(ista_2l:iend_2u,jsta_2l:jend_2u))
467 allocate(thz0(ista_2l:iend_2u,jsta_2l:jend_2u))
468 allocate(uz0(ista_2l:iend_2u,jsta_2l:jend_2u))
469 allocate(vz0(ista_2l:iend_2u,jsta_2l:jend_2u))
470 allocate(qs(ista_2l:iend_2u,jsta_2l:jend_2u))
471 allocate(ths(ista_2l:iend_2u,jsta_2l:jend_2u))
472 allocate(sno(ista_2l:iend_2u,jsta_2l:jend_2u))
473 allocate(snonc(ista_2l:iend_2u,jsta_2l:jend_2u))
474 allocate(ti(ista_2l:iend_2u,jsta_2l:jend_2u))
509 allocate(u10mean(ista_2l:iend_2u,jsta_2l:jend_2u))
510 allocate(v10mean(ista_2l:iend_2u,jsta_2l:jend_2u))
511 allocate(spduv10mean(ista_2l:iend_2u,jsta_2l:jend_2u))
512 allocate(swradmean(ista_2l:iend_2u,jsta_2l:jend_2u))
513 allocate(swnormmean(ista_2l:iend_2u,jsta_2l:jend_2u))
520 spduv10mean(i,j)=spval
522 swnormmean(i,j)=spval
526 allocate(snoavg(ista_2l:iend_2u,jsta_2l:jend_2u))
527 allocate(psfcavg(ista_2l:iend_2u,jsta_2l:jend_2u))
528 allocate(t10m(ista_2l:iend_2u,jsta_2l:jend_2u))
529 allocate(t10avg(ista_2l:iend_2u,jsta_2l:jend_2u))
530 allocate(akmsavg(ista_2l:iend_2u,jsta_2l:jend_2u))
531 allocate(akhsavg(ista_2l:iend_2u,jsta_2l:jend_2u))
532 allocate(u10max(ista_2l:iend_2u,jsta_2l:jend_2u))
533 allocate(v10max(ista_2l:iend_2u,jsta_2l:jend_2u))
534 allocate(u10h(ista_2l:iend_2u,jsta_2l:jend_2u))
535 allocate(v10h(ista_2l:iend_2u,jsta_2l:jend_2u))
553 allocate(akms(ista_2l:iend_2u,jsta_2l:jend_2u))
554 allocate(akhs(ista_2l:iend_2u,jsta_2l:jend_2u))
555 allocate(cuprec(ista_2l:iend_2u,jsta_2l:jend_2u))
556 allocate(acprec(ista_2l:iend_2u,jsta_2l:jend_2u))
557 allocate(ancprc(ista_2l:iend_2u,jsta_2l:jend_2u))
558 allocate(cuppt(ista_2l:iend_2u,jsta_2l:jend_2u))
572 allocate(rainc_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
573 allocate(rainc_bucket1(ista_2l:iend_2u,jsta_2l:jend_2u))
574 allocate(rainnc_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
575 allocate(rainnc_bucket1(ista_2l:iend_2u,jsta_2l:jend_2u))
576 allocate(pcp_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
577 allocate(pcp_bucket1(ista_2l:iend_2u,jsta_2l:jend_2u))
578 allocate(snow_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
579 allocate(snow_bucket1(ista_2l:iend_2u,jsta_2l:jend_2u))
580 allocate(graup_bucket(ista_2l:iend_2u,jsta_2l:jend_2u))
581 allocate(graup_bucket1(ista_2l:iend_2u,jsta_2l:jend_2u))
582 allocate(qrmax(ista_2l:iend_2u,jsta_2l:jend_2u))
583 allocate(tmax(ista_2l:iend_2u,jsta_2l:jend_2u))
584 allocate(snownc(ista_2l:iend_2u,jsta_2l:jend_2u))
585 allocate(graupelnc(ista_2l:iend_2u,jsta_2l:jend_2u))
586 allocate(tsnow(ista_2l:iend_2u,jsta_2l:jend_2u))
587 allocate(qvg(ista_2l:iend_2u,jsta_2l:jend_2u))
588 allocate(qv2m(ista_2l:iend_2u,jsta_2l:jend_2u))
589 allocate(qvl1(ista_2l:iend_2u,jsta_2l:jend_2u))
590 allocate(snfden(ista_2l:iend_2u,jsta_2l:jend_2u))
591 allocate(sndepac(ista_2l:iend_2u,jsta_2l:jend_2u))
592 allocate(int_smoke(ista_2l:iend_2u,jsta_2l:jend_2u))
593 allocate(mean_frp(ista_2l:iend_2u,jsta_2l:jend_2u))
594 allocate(int_aod(ista_2l:iend_2u,jsta_2l:jend_2u))
599 rainc_bucket(i,j)=spval
600 rainc_bucket1(i,j)=spval
601 rainnc_bucket(i,j)=spval
602 rainnc_bucket1(i,j)=spval
603 pcp_bucket(i,j)=spval
604 pcp_bucket1(i,j)=spval
605 snow_bucket(i,j)=spval
606 snow_bucket1(i,j)=spval
607 graup_bucket(i,j)=spval
608 graup_bucket1(i,j)=spval
624 allocate(smoke(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_sm))
636 allocate(rswin(ista_2l:iend_2u,jsta_2l:jend_2u))
637 allocate(swddni(ista_2l:iend_2u,jsta_2l:jend_2u))
638 allocate(swddif(ista_2l:iend_2u,jsta_2l:jend_2u))
639 allocate(swdnbc(ista_2l:iend_2u,jsta_2l:jend_2u))
640 allocate(swddnic(ista_2l:iend_2u,jsta_2l:jend_2u))
641 allocate(swddifc(ista_2l:iend_2u,jsta_2l:jend_2u))
642 allocate(swupbc(ista_2l:iend_2u,jsta_2l:jend_2u))
643 allocate(swupt(ista_2l:iend_2u,jsta_2l:jend_2u))
644 allocate(taod5502d(ista_2l:iend_2u,jsta_2l:jend_2u))
645 allocate(aerasy2d(ista_2l:iend_2u,jsta_2l:jend_2u))
646 allocate(aerssa2d(ista_2l:iend_2u,jsta_2l:jend_2u))
647 allocate(lwp(ista_2l:iend_2u,jsta_2l:jend_2u))
648 allocate(iwp(ista_2l:iend_2u,jsta_2l:jend_2u))
649 allocate(rlwin(ista_2l:iend_2u,jsta_2l:jend_2u))
650 allocate(lwdnbc(ista_2l:iend_2u,jsta_2l:jend_2u))
651 allocate(lwupbc(ista_2l:iend_2u,jsta_2l:jend_2u))
652 allocate(rlwtoa(ista_2l:iend_2u,jsta_2l:jend_2u))
653 allocate(rswtoa(ista_2l:iend_2u,jsta_2l:jend_2u))
678 allocate(tg(ista_2l:iend_2u,jsta_2l:jend_2u))
679 allocate(sfcshx(ista_2l:iend_2u,jsta_2l:jend_2u))
680 allocate(sfclhx(ista_2l:iend_2u,jsta_2l:jend_2u))
681 allocate(fis(ista_2l:iend_2u,jsta_2l:jend_2u))
682 allocate(t500(ista_2l:iend_2u,jsta_2l:jend_2u))
683 allocate(t700(ista_2l:iend_2u,jsta_2l:jend_2u))
684 allocate(z500(ista_2l:iend_2u,jsta_2l:jend_2u))
685 allocate(z700(ista_2l:iend_2u,jsta_2l:jend_2u))
686 allocate(teql(ista_2l:iend_2u,jsta_2l:jend_2u))
687 allocate(ieql(ista_2l:iend_2u,jsta_2l:jend_2u))
688 allocate(cfracl(ista_2l:iend_2u,jsta_2l:jend_2u))
689 allocate(cfracm(ista_2l:iend_2u,jsta_2l:jend_2u))
690 allocate(cfrach(ista_2l:iend_2u,jsta_2l:jend_2u))
691 allocate(acfrst(ista_2l:iend_2u,jsta_2l:jend_2u))
692 allocate(acfrcv(ista_2l:iend_2u,jsta_2l:jend_2u))
693 allocate(hbot(ista_2l:iend_2u,jsta_2l:jend_2u))
694 allocate(htop(ista_2l:iend_2u,jsta_2l:jend_2u))
695 allocate(aswin(ista_2l:iend_2u,jsta_2l:jend_2u))
696 allocate(alwin(ista_2l:iend_2u,jsta_2l:jend_2u))
697 allocate(aswout(ista_2l:iend_2u,jsta_2l:jend_2u))
698 allocate(alwout(ista_2l:iend_2u,jsta_2l:jend_2u))
699 allocate(aswtoa(ista_2l:iend_2u,jsta_2l:jend_2u))
700 allocate(alwtoa(ista_2l:iend_2u,jsta_2l:jend_2u))
729 allocate(czen(ista_2l:iend_2u,jsta_2l:jend_2u))
730 allocate(czmean(ista_2l:iend_2u,jsta_2l:jend_2u))
731 allocate(sigt4(ista_2l:iend_2u,jsta_2l:jend_2u))
732 allocate(rswout(ista_2l:iend_2u,jsta_2l:jend_2u))
733 allocate(radot(ista_2l:iend_2u,jsta_2l:jend_2u))
734 allocate(ncfrst(ista_2l:iend_2u,jsta_2l:jend_2u))
735 allocate(ncfrcv(ista_2l:iend_2u,jsta_2l:jend_2u))
736 allocate(smstot(ista_2l:iend_2u,jsta_2l:jend_2u))
737 allocate(pctsno(ista_2l:iend_2u,jsta_2l:jend_2u))
738 allocate(pshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
739 allocate(th10(ista_2l:iend_2u,jsta_2l:jend_2u))
740 allocate(q10(ista_2l:iend_2u,jsta_2l:jend_2u))
741 allocate(sr(ista_2l:iend_2u,jsta_2l:jend_2u))
742 allocate(prec(ista_2l:iend_2u,jsta_2l:jend_2u))
743 allocate(subshx(ista_2l:iend_2u,jsta_2l:jend_2u))
744 allocate(snopcx(ista_2l:iend_2u,jsta_2l:jend_2u))
745 allocate(sfcuvx(ista_2l:iend_2u,jsta_2l:jend_2u))
746 allocate(sfcevp(ista_2l:iend_2u,jsta_2l:jend_2u))
747 allocate(potevp(ista_2l:iend_2u,jsta_2l:jend_2u))
748 allocate(z0(ista_2l:iend_2u,jsta_2l:jend_2u))
749 allocate(ustar(ista_2l:iend_2u,jsta_2l:jend_2u))
750 allocate(pblh(ista_2l:iend_2u,jsta_2l:jend_2u))
751 allocate(pblhgust(ista_2l:iend_2u,jsta_2l:jend_2u))
752 allocate(mixht(ista_2l:iend_2u,jsta_2l:jend_2u))
753 allocate(twbs(ista_2l:iend_2u,jsta_2l:jend_2u))
754 allocate(qwbs(ista_2l:iend_2u,jsta_2l:jend_2u))
787 allocate(sfcexc(ista_2l:iend_2u,jsta_2l:jend_2u))
788 allocate(grnflx(ista_2l:iend_2u,jsta_2l:jend_2u))
789 allocate(soiltb(ista_2l:iend_2u,jsta_2l:jend_2u))
790 allocate(z1000(ista_2l:iend_2u,jsta_2l:jend_2u))
791 allocate(slp(ista_2l:iend_2u,jsta_2l:jend_2u))
792 allocate(pslp(ista_2l:iend_2u,jsta_2l:jend_2u))
793 allocate(f(ista_2l:iend_2u,jsta_2l:jend_2u))
794 allocate(albedo(ista_2l:iend_2u,jsta_2l:jend_2u))
795 allocate(albase(ista_2l:iend_2u,jsta_2l:jend_2u))
796 allocate(cldfra(ista_2l:iend_2u,jsta_2l:jend_2u))
797 allocate(cprate(ista_2l:iend_2u,jsta_2l:jend_2u))
798 allocate(cnvcfr(ista_2l:iend_2u,jsta_2l:jend_2u))
799 allocate(ivgtyp(ista_2l:iend_2u,jsta_2l:jend_2u))
800 allocate(isltyp(ista_2l:iend_2u,jsta_2l:jend_2u))
801 allocate(hbotd(ista_2l:iend_2u,jsta_2l:jend_2u))
802 allocate(htopd(ista_2l:iend_2u,jsta_2l:jend_2u))
803 allocate(hbots(ista_2l:iend_2u,jsta_2l:jend_2u))
804 allocate(htops(ista_2l:iend_2u,jsta_2l:jend_2u))
805 allocate(cldefi(ista_2l:iend_2u,jsta_2l:jend_2u))
806 allocate(islope(ista_2l:iend_2u,jsta_2l:jend_2u))
807 allocate(si(ista_2l:iend_2u,jsta_2l:jend_2u))
808 allocate(lspa(ista_2l:iend_2u,jsta_2l:jend_2u))
809 allocate(rswinc(ista_2l:iend_2u,jsta_2l:jend_2u))
810 allocate(vis(ista_2l:iend_2u,jsta_2l:jend_2u))
811 allocate(pd(ista_2l:iend_2u,jsta_2l:jend_2u))
812 allocate(mxsnal(ista_2l:iend_2u,jsta_2l:jend_2u))
813 allocate(epsr(ista_2l:iend_2u,jsta_2l:jend_2u))
848 allocate(sfcux(ista_2l:iend_2u,jsta_2l:jend_2u))
849 allocate(sfcvx(ista_2l:iend_2u,jsta_2l:jend_2u))
850 allocate(sfcuxi(ista_2l:iend_2u,jsta_2l:jend_2u))
851 allocate(sfcvxi(ista_2l:iend_2u,jsta_2l:jend_2u))
852 allocate(avgalbedo(ista_2l:iend_2u,jsta_2l:jend_2u))
853 allocate(avgcprate(ista_2l:iend_2u,jsta_2l:jend_2u))
854 allocate(avgprec(ista_2l:iend_2u,jsta_2l:jend_2u))
855 allocate(avgprec_cont(ista_2l:iend_2u,jsta_2l:jend_2u))
856 allocate(avgcprate_cont(ista_2l:iend_2u,jsta_2l:jend_2u))
857 allocate(ptop(ista_2l:iend_2u,jsta_2l:jend_2u))
858 allocate(pbot(ista_2l:iend_2u,jsta_2l:jend_2u))
859 allocate(avgcfrach(ista_2l:iend_2u,jsta_2l:jend_2u))
860 allocate(avgcfracm(ista_2l:iend_2u,jsta_2l:jend_2u))
861 allocate(avgcfracl(ista_2l:iend_2u,jsta_2l:jend_2u))
862 allocate(avgtcdc(ista_2l:iend_2u,jsta_2l:jend_2u))
863 allocate(auvbin(ista_2l:iend_2u,jsta_2l:jend_2u))
864 allocate(auvbinc(ista_2l:iend_2u,jsta_2l:jend_2u))
865 allocate(ptopl(ista_2l:iend_2u,jsta_2l:jend_2u))
866 allocate(pbotl(ista_2l:iend_2u,jsta_2l:jend_2u))
867 allocate(ttopl(ista_2l:iend_2u,jsta_2l:jend_2u))
868 allocate(ptopm(ista_2l:iend_2u,jsta_2l:jend_2u))
869 allocate(pbotm(ista_2l:iend_2u,jsta_2l:jend_2u))
870 allocate(ttopm(ista_2l:iend_2u,jsta_2l:jend_2u))
871 allocate(ptoph(ista_2l:iend_2u,jsta_2l:jend_2u))
872 allocate(pboth(ista_2l:iend_2u,jsta_2l:jend_2u))
873 allocate(ttoph(ista_2l:iend_2u,jsta_2l:jend_2u))
874 allocate(sfcugs(ista_2l:iend_2u,jsta_2l:jend_2u))
875 allocate(sfcvgs(ista_2l:iend_2u,jsta_2l:jend_2u))
876 allocate(pblcfr(ista_2l:iend_2u,jsta_2l:jend_2u))
877 allocate(cldwork(ista_2l:iend_2u,jsta_2l:jend_2u))
878 allocate(gtaux(ista_2l:iend_2u,jsta_2l:jend_2u))
879 allocate(gtauy(ista_2l:iend_2u,jsta_2l:jend_2u))
880 allocate(cd10(ista_2l:iend_2u,jsta_2l:jend_2u))
881 allocate(ch10(ista_2l:iend_2u,jsta_2l:jend_2u))
882 allocate(mdltaux(ista_2l:iend_2u,jsta_2l:jend_2u))
883 allocate(mdltauy(ista_2l:iend_2u,jsta_2l:jend_2u))
884 allocate(runoff(ista_2l:iend_2u,jsta_2l:jend_2u))
896 avgprec_cont(i,j)=spval
897 avgcprate_cont(i,j)=spval
928 allocate(maxtshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
929 allocate(mintshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
930 allocate(maxrhshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
931 allocate(minrhshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
932 allocate(maxqshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
933 allocate(minqshltr(ista_2l:iend_2u,jsta_2l:jend_2u))
934 allocate(dzice(ista_2l:iend_2u,jsta_2l:jend_2u))
935 allocate(alwinc(ista_2l:iend_2u,jsta_2l:jend_2u))
936 allocate(alwoutc(ista_2l:iend_2u,jsta_2l:jend_2u))
937 allocate(alwtoac(ista_2l:iend_2u,jsta_2l:jend_2u))
938 allocate(aswinc(ista_2l:iend_2u,jsta_2l:jend_2u))
939 allocate(aswoutc(ista_2l:iend_2u,jsta_2l:jend_2u))
940 allocate(aswtoac(ista_2l:iend_2u,jsta_2l:jend_2u))
941 allocate(aswintoa(ista_2l:iend_2u,jsta_2l:jend_2u))
942 allocate(smcwlt(ista_2l:iend_2u,jsta_2l:jend_2u))
943 allocate(suntime(ista_2l:iend_2u,jsta_2l:jend_2u))
944 allocate(fieldcapa(ista_2l:iend_2u,jsta_2l:jend_2u))
945 allocate(avisbeamswin(ista_2l:iend_2u,jsta_2l:jend_2u))
946 allocate(avisdiffswin(ista_2l:iend_2u,jsta_2l:jend_2u))
947 allocate(airbeamswin(ista_2l:iend_2u,jsta_2l:jend_2u))
948 allocate(airdiffswin(ista_2l:iend_2u,jsta_2l:jend_2u))
949 allocate(snowfall(ista_2l:iend_2u,jsta_2l:jend_2u))
950 allocate(acond(ista_2l:iend_2u,jsta_2l:jend_2u))
951 allocate(edir(ista_2l:iend_2u,jsta_2l:jend_2u))
952 allocate(ecan(ista_2l:iend_2u,jsta_2l:jend_2u))
953 allocate(etrans(ista_2l:iend_2u,jsta_2l:jend_2u))
954 allocate(esnow(ista_2l:iend_2u,jsta_2l:jend_2u))
955 allocate(avgedir(ista_2l:iend_2u,jsta_2l:jend_2u))
956 allocate(avgecan(ista_2l:iend_2u,jsta_2l:jend_2u))
957 allocate(avgetrans(ista_2l:iend_2u,jsta_2l:jend_2u))
958 allocate(avgesnow(ista_2l:iend_2u,jsta_2l:jend_2u))
959 allocate(avgpotevp(ista_2l:iend_2u,jsta_2l:jend_2u))
960 allocate(aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
961 allocate(du_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
962 allocate(ss_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
963 allocate(su_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
964 allocate(oc_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
965 allocate(bc_aod550(ista_2l:iend_2u,jsta_2l:jend_2u))
966 allocate(landfrac(ista_2l:iend_2u,jsta_2l:jend_2u))
967 allocate(paha(ista_2l:iend_2u,jsta_2l:jend_2u))
968 allocate(pahi(ista_2l:iend_2u,jsta_2l:jend_2u))
969 allocate(tecan(ista_2l:iend_2u,jsta_2l:jend_2u))
970 allocate(tetran(ista_2l:iend_2u,jsta_2l:jend_2u))
971 allocate(tedir(ista_2l:iend_2u,jsta_2l:jend_2u))
972 allocate(twa(ista_2l:iend_2u,jsta_2l:jend_2u))
973 allocate(fdnsst(ista_2l:iend_2u,jsta_2l:jend_2u))
974 allocate(pwat(ista_2l:iend_2u,jsta_2l:jend_2u))
981 maxrhshltr(i,j)=spval
982 minrhshltr(i,j)=spval
996 avisbeamswin(i,j)=spval
997 avisdiffswin(i,j)=spval
998 airbeamswin(i,j)=spval
999 airdiffswin(i,j)=spval
1008 avgetrans(i,j)=spval
1010 avgpotevp(i,j)=spval
1012 du_aod550(i,j)=spval
1013 ss_aod550(i,j)=spval
1014 su_aod550(i,j)=spval
1015 oc_aod550(i,j)=spval
1016 bc_aod550(i,j)=spval
1031 allocate(hbm2(ista_2l:iend_2u,jsta_2l:jend_2u))
1032 allocate(sm(ista_2l:iend_2u,jsta_2l:jend_2u))
1033 allocate(sice(ista_2l:iend_2u,jsta_2l:jend_2u))
1034 allocate(lmh(ista_2l:iend_2u,jsta_2l:jend_2u))
1035 allocate(lmv(ista_2l:iend_2u,jsta_2l:jend_2u))
1036 allocate(gdlat(ista_2l:iend_2u,jsta_2l:jend_2u))
1037 allocate(gdlon(ista_2l:iend_2u,jsta_2l:jend_2u))
1038 allocate(dx(ista_2l:iend_2u,jsta_2l:jend_2u))
1039 allocate(dy(ista_2l:iend_2u,jsta_2l:jend_2u))
1042 do j=jsta_2l,jend_2u
1043 do i=ista_2l,iend_2u
1056 if (me == 0) print *,
' gocart_on=',gocart_on
1061 allocate(dust(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_du))
1062 allocate(salt(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_ss))
1063 allocate(soot(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_bc))
1064 allocate(waso(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_oc))
1065 allocate(suso(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_su))
1066 allocate(pp25(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_su))
1067 allocate(pp10(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_su))
1072 do j=jsta_2l,jend_2u
1073 do i=ista_2l,iend_2u
1082 do j=jsta_2l,jend_2u
1083 do i=ista_2l,iend_2u
1092 do j=jsta_2l,jend_2u
1093 do i=ista_2l,iend_2u
1102 do j=jsta_2l,jend_2u
1103 do i=ista_2l,iend_2u
1112 do j=jsta_2l,jend_2u
1113 do i=ista_2l,iend_2u
1122 allocate(ext(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1123 allocate(asy(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1124 allocate(ssa(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1125 allocate(sca(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1129 do j=jsta_2l,jend_2u
1130 do i=ista_2l,iend_2u
1138 allocate(duem(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_du))
1139 allocate(dusd(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_du))
1140 allocate(dudp(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_du))
1141 allocate(duwt(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_du))
1142 allocate(dusv(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_du))
1143 allocate(suem(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_su))
1144 allocate(susd(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_su))
1145 allocate(sudp(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_su))
1146 allocate(suwt(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_su))
1147 allocate(ocem(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_oc))
1148 allocate(ocsd(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_oc))
1149 allocate(ocdp(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_oc))
1150 allocate(ocwt(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_oc))
1151 allocate(ocsv(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_oc))
1152 allocate(bcem(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_bc))
1153 allocate(bcsd(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_bc))
1154 allocate(bcdp(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_bc))
1155 allocate(bcwt(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_bc))
1156 allocate(bcsv(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_bc))
1157 allocate(ssem(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_ss))
1158 allocate(sssd(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_ss))
1159 allocate(ssdp(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_ss))
1160 allocate(sswt(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_ss))
1161 allocate(sssv(ista_2l:iend_2u,jsta_2l:jend_2u,nbin_ss))
1165 do j=jsta_2l,jend_2u
1166 do i=ista_2l,iend_2u
1177 do j=jsta_2l,jend_2u
1178 do i=ista_2l,iend_2u
1188 do j=jsta_2l,jend_2u
1189 do i=ista_2l,iend_2u
1200 do j=jsta_2l,jend_2u
1201 do i=ista_2l,iend_2u
1212 do j=jsta_2l,jend_2u
1213 do i=ista_2l,iend_2u
1222 allocate(rhomid(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1226 do j=jsta_2l,jend_2u
1227 do i=ista_2l,iend_2u
1233 allocate(dusmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1234 allocate(ducmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1235 allocate(dusmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1236 allocate(ducmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1237 allocate(susmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1238 allocate(sucmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1239 allocate(susmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1240 allocate(sucmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1241 allocate(ocsmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1242 allocate(occmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1243 allocate(ocsmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1244 allocate(occmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1245 allocate(bcsmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1246 allocate(bccmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1247 allocate(bcsmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1248 allocate(bccmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1249 allocate(sssmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1250 allocate(sscmass(ista_2l:iend_2u,jsta_2l:jend_2u))
1251 allocate(sssmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1252 allocate(sscmass25(ista_2l:iend_2u,jsta_2l:jend_2u))
1253 allocate(dustcb(ista_2l:iend_2u,jsta_2l:jend_2u))
1254 allocate(occb(ista_2l:iend_2u,jsta_2l:jend_2u))
1255 allocate(bccb(ista_2l:iend_2u,jsta_2l:jend_2u))
1256 allocate(sulfcb(ista_2l:iend_2u,jsta_2l:jend_2u))
1257 allocate(pp25cb(ista_2l:iend_2u,jsta_2l:jend_2u))
1258 allocate(pp10cb(ista_2l:iend_2u,jsta_2l:jend_2u))
1259 allocate(sscb(ista_2l:iend_2u,jsta_2l:jend_2u))
1260 allocate(dustallcb(ista_2l:iend_2u,jsta_2l:jend_2u))
1261 allocate(ssallcb(ista_2l:iend_2u,jsta_2l:jend_2u))
1262 allocate(dustpm(ista_2l:iend_2u,jsta_2l:jend_2u))
1263 allocate(dustpm10(ista_2l:iend_2u,jsta_2l:jend_2u))
1264 allocate(sspm(ista_2l:iend_2u,jsta_2l:jend_2u))
1265 allocate(maod(ista_2l:iend_2u,jsta_2l:jend_2u))
1268 do j=jsta_2l,jend_2u
1269 do i=ista_2l,iend_2u
1272 dusmass25(i,j)=spval
1273 ducmass25(i,j)=spval
1276 susmass25(i,j)=spval
1277 sucmass25(i,j)=spval
1280 ocsmass25(i,j)=spval
1281 occmass25(i,j)=spval
1284 bcsmass25(i,j)=spval
1285 bccmass25(i,j)=spval
1288 sssmass25(i,j)=spval
1289 sscmass25(i,j)=spval
1297 dustallcb(i,j)=spval
1307 allocate(acswupt(ista_2l:iend_2u,jsta_2l:jend_2u))
1308 allocate(swdnt(ista_2l:iend_2u,jsta_2l:jend_2u))
1309 allocate(acswdnt(ista_2l:iend_2u,jsta_2l:jend_2u))
1312 do j=jsta_2l,jend_2u
1313 do i=ista_2l,iend_2u
1321 allocate(ddvdx(ista_2l:iend_2u,jsta_2l:jend_2u))
1322 allocate(ddudy(ista_2l:iend_2u,jsta_2l:jend_2u))
1323 allocate(uuavg(ista_2l:iend_2u,jsta_2l:jend_2u))
1326 do j=jsta_2l,jend_2u
1327 do i=ista_2l,iend_2u
1335 if (me == 0) print *,
'aqfcmaq_on= ', aqfcmaq_on
1336 if (aqfcmaq_on)
then
1338 allocate(ozcon(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1339 allocate(pmtf(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
1344 do j=jsta_2l,jend_2u
1345 do i=ista_2l,iend_2u
dvdxdudy() computes dudy, dvdx, uwnd