UPP  001
 All Data Structures Files Functions Pages
ALLOCATE_ALL.f
Go to the documentation of this file.
1 
2 !
31  SUBROUTINE allocate_all()
32 !
33  use vrbls4d
34  use vrbls3d
35  use vrbls2d
36  use soil
37  use masks
38  use upp_math, only: ddvdx, ddudy, uuavg
39 !
40  !use params_mod
41  use ctlblk_mod
42 !- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
43  implicit none
44 !
45  include 'mpif.h'
46 !
47  integer ierr,jsx,jex
48  integer i,j,l,k
49 ! Allocate arrays
50  allocate(u(im+1,jsta_2l:jend_2u,lm))
51  allocate(v(im,jsta_2l:jvend_2u,lm))
52  allocate(t(im,jsta_2l:jend_2u,lm))
53 ! CHUANG ADD POTENTIAL TEMP BECAUSE WRF OUTPUT THETA
54 ! allocate(th(im,jsta_2l:jend_2u,lm))
55  allocate(q(im,jsta_2l:jend_2u,lm))
56 ! allocate(w(im,jsta_2l:jend_2u,lp1))
57  allocate(uh(im,jsta_2l:jend_2u,lm))
58  allocate(vh(im,jsta_2l:jend_2u,lm))
59  allocate(wh(im,jsta_2l:jend_2u,lm))
60  allocate(pmid(im,jsta_2l:jend_2u,lm))
61  allocate(pmidv(im,jsta_2l:jend_2u,lm))
62  allocate(pint(im,jsta_2l:jend_2u,lp1))
63  allocate(alpint(im,jsta_2l:jend_2u,lp1))
64  allocate(zmid(im,jsta_2l:jend_2u,lm))
65  allocate(zint(im,jsta_2l:jend_2u,lp1))
66 ! allocate(rainw(im,jsta_2l:jend_2u,lm))
67  allocate(q2(im,jsta_2l:jend_2u,lm))
68  allocate(omga(im,jsta_2l:jend_2u,lm))
69  allocate(dpres(im,jsta_2l:jend_2u,lm))
70  allocate(t_adj(im,jsta_2l:jend_2u,lm))
71  allocate(ttnd(im,jsta_2l:jend_2u,lm))
72  allocate(rswtt(im,jsta_2l:jend_2u,lm))
73  allocate(rlwtt(im,jsta_2l:jend_2u,lm))
74  allocate(exch_h(im,jsta_2l:jend_2u,lm))
75  allocate(train(im,jsta_2l:jend_2u,lm))
76  allocate(tcucn(im,jsta_2l:jend_2u,lm))
77  allocate(el_pbl(im,jsta_2l:jend_2u,lm))
78 
79 !Initialization
80 !$omp parallel do private(i,j,l)
81  do l=1,lm
82  do j=jsta_2l,jend_2u
83  do i=1,im+1
84  u(i,j,l)=0.
85  enddo
86  enddo
87  enddo
88 !$omp parallel do private(i,j,l)
89  do l=1,lm
90  do j=jsta_2l,jvend_2u
91  do i=1,im
92  v(i,j,l)=0.
93  enddo
94  enddo
95  enddo
96 !$omp parallel do private(i,j,l)
97  do l=1,lm
98  do j=jsta_2l,jend_2u
99  do i=1,im
100  t(i,j,l)=spval
101  q(i,j,l)=spval
102  uh(i,j,l)=spval
103  vh(i,j,l)=spval
104  wh(i,j,l)=spval
105  pmid(i,j,l)=spval
106  pmidv(i,j,l)=spval
107  zmid(i,j,l)=spval
108  q2(i,j,l)=spval
109  omga(i,j,l)=spval
110  dpres(i,j,l)=spval
111  t_adj(i,j,l)=spval
112  ttnd(i,j,l)=spval
113  rswtt(i,j,l)=spval
114  rlwtt(i,j,l)=spval
115  exch_h(i,j,l)=spval
116  train(i,j,l)=spval
117  tcucn(i,j,l)=spval
118  el_pbl(i,j,l)=spval
119  enddo
120  enddo
121  enddo
122 !$omp parallel do private(i,j,l)
123  do l=1,lp1
124  do j=jsta_2l,jend_2u
125  do i=1,im
126  pint(i,j,l)=spval
127  alpint(i,j,l)=spval
128  zint(i,j,l)=spval
129  enddo
130  enddo
131  enddo
132 
133 ! MP FIELD
134  allocate(cwm(im,jsta_2l:jend_2u,lm))
135  allocate(f_ice(im,jsta_2l:jend_2u,lm))
136  allocate(f_rain(im,jsta_2l:jend_2u,lm))
137  allocate(f_rimef(im,jsta_2l:jend_2u,lm))
138  allocate(qqw(im,jsta_2l:jend_2u,lm))
139  allocate(qrimef(im,jsta_2l:jend_2u,lm))
140  allocate(qqi(im,jsta_2l:jend_2u,lm))
141  allocate(qqr(im,jsta_2l:jend_2u,lm))
142  allocate(qqs(im,jsta_2l:jend_2u,lm))
143  allocate(qqg(im,jsta_2l:jend_2u,lm))
144  allocate(qqnw(im,jsta_2l:jend_2u,lm))
145  allocate(qqni(im,jsta_2l:jend_2u,lm))
146  allocate(qqnr(im,jsta_2l:jend_2u,lm))
147  allocate(qqnwfa(im,jsta_2l:jend_2u,lm))
148  allocate(qqnifa(im,jsta_2l:jend_2u,lm))
149  allocate(taod5503d(im,jsta_2l:jend_2u,lm))
150  allocate(aextc55(im,jsta_2l:jend_2u,lm))
151  allocate(extcof55(im,jsta_2l:jend_2u,lm))
152  allocate(qc_bl(im,jsta_2l:jend_2u,lm))
153  allocate(cfr(im,jsta_2l:jend_2u,lm))
154  allocate(cfr_raw(im,jsta_2l:jend_2u,lm))
155  allocate(dbz(im,jsta_2l:jend_2u,lm))
156  allocate(dbzr(im,jsta_2l:jend_2u,lm))
157  allocate(dbzi(im,jsta_2l:jend_2u,lm))
158  allocate(dbzc(im,jsta_2l:jend_2u,lm))
159  allocate(mcvg(im,jsta_2l:jend_2u,lm))
160  allocate(nlice(im,jsta_2l:jend_2u,lm))
161 !Initialization
162 !$omp parallel do private(i,j,l)
163  do l=1,lm
164  do j=jsta_2l,jend_2u
165  do i=1,im
166  cwm(i,j,l)=spval
167  f_ice(i,j,l)=spval
168  f_rain(i,j,l)=spval
169  f_rimef(i,j,l)=spval
170  qqw(i,j,l)=spval
171  qrimef(i,j,l)=spval
172  qqi(i,j,l)=spval
173  qqr(i,j,l)=spval
174  qqs(i,j,l)=spval
175  qqg(i,j,l)=spval
176  qqnw(i,j,l)=spval
177  qqni(i,j,l)=spval
178  qqnr(i,j,l)=spval
179  qqnwfa(i,j,l)=spval
180  qqnifa(i,j,l)=spval
181  taod5503d(i,j,l)=spval
182  aextc55(i,j,l)=0.
183  extcof55(i,j,l)=0.
184  qc_bl(i,j,l)=spval
185  cfr(i,j,l)=spval
186  cfr_raw(i,j,l)=spval
187  dbz(i,j,l)=spval
188  dbzr(i,j,l)=spval
189  dbzi(i,j,l)=spval
190  dbzc(i,j,l)=spval
191  mcvg(i,j,l)=spval
192  nlice(i,j,l)=spval
193  enddo
194  enddo
195  enddo
196 ! Wm Lewis: added
197  allocate(nrain(im,jsta_2l:jend_2u,lm))
198  allocate(radius_cloud(im,jsta_2l:jend_2u,lm))
199  allocate(radius_ice(im,jsta_2l:jend_2u,lm))
200  allocate(radius_snow(im,jsta_2l:jend_2u,lm))
201 ! KRS: HWRF Addition for thompson reflectivity
202 ! or non-ferrier physics. wrf-derived
203  allocate(refl_10cm(im,jsta_2l:jend_2u,lm))
204 !GFS FIELD
205  allocate(o3(im,jsta_2l:jend_2u,lm))
206  allocate(o(im,jsta_2l:jend_2u,lm))
207  allocate(o2(im,jsta_2l:jend_2u,lm))
208  allocate(tcucns(im,jsta_2l:jend_2u,lm))
209 !Initialization
210 !$omp parallel do private(i,j,l)
211  do l=1,lm
212  do j=jsta_2l,jend_2u
213  do i=1,im
214  nrain(i,j,l)=spval
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
219  o3(i,j,l)=spval
220  o(i,j,l)=spval
221  o2(i,j,l)=spval
222  tcucns(i,j,l)=spval
223  enddo
224  enddo
225  enddo
226 ! Add GFS d3d fields
227  if (me == 0) print *,' d3d_on=',d3d_on
228  if (d3d_on) then
229  allocate(vdifftt(im,jsta_2l:jend_2u,lm))
230 ! allocate(tcucns(im,jsta_2l:jend_2u,lm))
231  allocate(vdiffmois(im,jsta_2l:jend_2u,lm))
232  allocate(dconvmois(im,jsta_2l:jend_2u,lm))
233  allocate(sconvmois(im,jsta_2l:jend_2u,lm))
234  allocate(nradtt(im,jsta_2l:jend_2u,lm))
235  allocate(o3vdiff(im,jsta_2l:jend_2u,lm))
236  allocate(o3prod(im,jsta_2l:jend_2u,lm))
237  allocate(o3tndy(im,jsta_2l:jend_2u,lm))
238  allocate(mwpv(im,jsta_2l:jend_2u,lm))
239  allocate(unknown(im,jsta_2l:jend_2u,lm))
240  allocate(vdiffzacce(im,jsta_2l:jend_2u,lm))
241  allocate(zgdrag(im,jsta_2l:jend_2u,lm))
242  allocate(cnvctummixing(im,jsta_2l:jend_2u,lm))
243  allocate(vdiffmacce(im,jsta_2l:jend_2u,lm))
244  allocate(mgdrag(im,jsta_2l:jend_2u,lm))
245  allocate(cnvctvmmixing(im,jsta_2l:jend_2u,lm))
246  allocate(ncnvctcfrac(im,jsta_2l:jend_2u,lm))
247  allocate(cnvctumflx(im,jsta_2l:jend_2u,lm))
248  allocate(cnvctdmflx(im,jsta_2l:jend_2u,lm))
249  allocate(cnvctdetmflx(im,jsta_2l:jend_2u,lm))
250  allocate(cnvctzgdrag(im,jsta_2l:jend_2u,lm))
251  allocate(cnvctmgdrag(im,jsta_2l:jend_2u,lm))
252 !Initialization
253 !$omp parallel do private(i,j,l)
254  do l=1,lm
255  do j=jsta_2l,jend_2u
256  do i=1,im
257  vdifftt(i,j,l)=spval
258  vdiffmois(i,j,l)=spval
259  dconvmois(i,j,l)=spval
260  sconvmois(i,j,l)=spval
261  nradtt(i,j,l)=spval
262  o3vdiff(i,j,l)=spval
263  o3prod(i,j,l)=spval
264  o3tndy(i,j,l)=spval
265  mwpv(i,j,l)=spval
266  unknown(i,j,l)=spval
267  vdiffzacce(i,j,l)=spval
268  zgdrag(i,j,l)=spval
269  cnvctummixing(i,j,l)=spval
270  vdiffmacce(i,j,l)=spval
271  mgdrag(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
279  enddo
280  enddo
281  enddo
282  endif
283 !
284  allocate(htm(im,jsta_2l:jend_2u,lm))
285  allocate(vtm(im,jsta_2l:jend_2u,lm))
286 ! add GFIP ICING
287  allocate(icing_gfip(im,jsta_2l:jend_2u,lm))
288  allocate(icing_gfis(im,jsta_2l:jend_2u,lm))
289 !
290 ! add GTG turbulence
291  allocate(catedr(im,jsta_2l:jend_2u,lm))
292  allocate(mwt(im,jsta_2l:jend_2u,lm))
293  allocate(gtg(im,jsta_2l:jend_2u,lm))
294 !Initialization
295 !$omp parallel do private(i,j,l)
296  do l=1,lm
297  do j=jsta_2l,jend_2u
298  do i=1,im
299  htm(i,j,l)=spval
300  vtm(i,j,l)=spval
301  icing_gfip(i,j,l)=spval
302  icing_gfis(i,j,l)=spval
303  catedr(i,j,l)=spval
304  mwt(i,j,l)=spval
305  gtg(i,j,l)=spval
306  enddo
307  enddo
308  enddo
309 !
310 ! FROM SOIL
311 !
312  allocate(smc(im,jsta_2l:jend_2u,nsoil))
313  allocate(stc(im,jsta_2l:jend_2u,nsoil))
314  allocate(sh2o(im,jsta_2l:jend_2u,nsoil))
315  allocate(sldpth(nsoil))
316  allocate(rtdpth(nsoil))
317  allocate(sllevel(nsoil))
318 !Initialization
319 !$omp parallel do private(i,j,l)
320  do l=1,nsoil
321  do j=jsta_2l,jend_2u
322  do i=1,im
323  smc(i,j,l)=spval
324  stc(i,j,l)=spval
325  sh2o(i,j,l)=spval
326  enddo
327  enddo
328  enddo
329 !$omp parallel do private(i)
330  do i=1,nsoil
331  sldpth(i)=spval
332  rtdpth(i)=spval
333  sllevel(i)=spval
334  enddo
335 !
336 ! FROM VRBLS2D
337 !
338 ! SRD
339  allocate(wspd10max(im,jsta_2l:jend_2u))
340  allocate(w_up_max(im,jsta_2l:jend_2u))
341  allocate(w_dn_max(im,jsta_2l:jend_2u))
342  allocate(w_mean(im,jsta_2l:jend_2u))
343  allocate(refd_max(im,jsta_2l:jend_2u))
344  allocate(prate_max(im,jsta_2l:jend_2u))
345  allocate(fprate_max(im,jsta_2l:jend_2u))
346  allocate(up_heli_max(im,jsta_2l:jend_2u))
347  allocate(up_heli_max16(im,jsta_2l:jend_2u))
348  allocate(up_heli_min(im,jsta_2l:jend_2u))
349  allocate(up_heli_min16(im,jsta_2l:jend_2u))
350  allocate(up_heli_max02(im,jsta_2l:jend_2u))
351  allocate(up_heli_min02(im,jsta_2l:jend_2u))
352  allocate(up_heli_max03(im,jsta_2l:jend_2u))
353  allocate(up_heli_min03(im,jsta_2l:jend_2u))
354 !Initialization
355 !$omp parallel do private(i,j)
356  do j=jsta_2l,jend_2u
357  do i=1,im
358  wspd10max(i,j)=spval
359  w_up_max(i,j)=spval
360  w_dn_max(i,j)=spval
361  w_mean(i,j)=spval
362  refd_max(i,j)=spval
363  prate_max(i,j)=spval
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
373  enddo
374  enddo
375  allocate(rel_vort_max(im,jsta_2l:jend_2u))
376  allocate(rel_vort_max01(im,jsta_2l:jend_2u))
377  allocate(rel_vort_maxhy1(im,jsta_2l:jend_2u))
378  allocate(wspd10umax(im,jsta_2l:jend_2u))
379  allocate(wspd10vmax(im,jsta_2l:jend_2u))
380  allocate(refdm10c_max(im,jsta_2l:jend_2u))
381  allocate(hail_max2d(im,jsta_2l:jend_2u))
382  allocate(hail_maxk1(im,jsta_2l:jend_2u))
383  allocate(hail_maxhailcast(im,jsta_2l:jend_2u))
384  allocate(grpl_max(im,jsta_2l:jend_2u))
385  allocate(up_heli(im,jsta_2l:jend_2u))
386  allocate(up_heli16(im,jsta_2l:jend_2u))
387  allocate(ltg1_max(im,jsta_2l:jend_2u))
388  allocate(ltg2_max(im,jsta_2l:jend_2u))
389  allocate(ltg3_max(im,jsta_2l:jend_2u))
390  allocate(nci_ltg(im,jsta_2l:jend_2u))
391  allocate(nca_ltg(im,jsta_2l:jend_2u))
392  allocate(nci_wq(im,jsta_2l:jend_2u))
393  allocate(nca_wq(im,jsta_2l:jend_2u))
394  allocate(nci_refd(im,jsta_2l:jend_2u))
395  allocate(nca_refd(im,jsta_2l:jend_2u))
396 !Initialization
397 !$omp parallel do private(i,j)
398  do j=jsta_2l,jend_2u
399  do i=1,im
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
409  grpl_max(i,j)=spval
410  up_heli(i,j)=spval
411  up_heli16(i,j)=spval
412  ltg1_max(i,j)=spval
413  ltg2_max(i,j)=spval
414  ltg3_max(i,j)=spval
415  nci_ltg(i,j)=spval
416  nca_ltg(i,j)=spval
417  nci_wq(i,j)=spval
418  nca_wq(i,j)=spval
419  nci_refd(i,j)=spval
420  nca_refd(i,j)=spval
421  enddo
422  enddo
423 ! SRD
424 ! CRA
425  allocate(ref_10cm(im,jsta_2l:jend_2u,lm))
426 !Initialization
427 !$omp parallel do private(i,j,l)
428  do l=1,lm
429  do j=jsta_2l,jend_2u
430  do i=1,im
431  ref_10cm(i,j,l)=spval
432  enddo
433  enddo
434  enddo
435  allocate(refc_10cm(im,jsta_2l:jend_2u))
436  allocate(ref1km_10cm(im,jsta_2l:jend_2u))
437  allocate(ref4km_10cm(im,jsta_2l:jend_2u))
438 !Initialization
439 !$omp parallel do private(i,j)
440  do j=jsta_2l,jend_2u
441  do i=1,im
442  refc_10cm(i,j)=spval
443  ref1km_10cm(i,j)=spval
444  ref4km_10cm(i,j)=spval
445  enddo
446  enddo
447 ! CRA
448  allocate(u10(im,jsta_2l:jend_2u))
449  allocate(v10(im,jsta_2l:jend_2u))
450  allocate(tshltr(im,jsta_2l:jend_2u))
451  allocate(qshltr(im,jsta_2l:jend_2u))
452  allocate(mrshltr(im,jsta_2l:jend_2u))
453  allocate(smstav(im,jsta_2l:jend_2u))
454  allocate(ssroff(im,jsta_2l:jend_2u))
455  allocate(bgroff(im,jsta_2l:jend_2u))
456  allocate(vegfrc(im,jsta_2l:jend_2u))
457  allocate(shdmin(im,jsta_2l:jend_2u))
458  allocate(shdmax(im,jsta_2l:jend_2u))
459  allocate(lai(im,jsta_2l:jend_2u))
460  allocate(acsnow(im,jsta_2l:jend_2u))
461  allocate(acgraup(im,jsta_2l:jend_2u))
462  allocate(acfrain(im,jsta_2l:jend_2u))
463  allocate(acsnom(im,jsta_2l:jend_2u))
464  allocate(cmc(im,jsta_2l:jend_2u))
465  allocate(sst(im,jsta_2l:jend_2u))
466  allocate(qz0(im,jsta_2l:jend_2u))
467  allocate(thz0(im,jsta_2l:jend_2u))
468  allocate(uz0(im,jsta_2l:jend_2u))
469  allocate(vz0(im,jsta_2l:jend_2u))
470  allocate(qs(im,jsta_2l:jend_2u))
471  allocate(ths(im,jsta_2l:jend_2u))
472  allocate(sno(im,jsta_2l:jend_2u))
473  allocate(snonc(im,jsta_2l:jend_2u))
474  allocate(ti(im,jsta_2l:jend_2u))
475 !Initialization
476 !$omp parallel do private(i,j)
477  do j=jsta_2l,jend_2u
478  do i=1,im
479  u10(i,j)=spval
480  v10(i,j)=spval
481  tshltr(i,j)=spval
482  qshltr(i,j)=spval
483  mrshltr(i,j)=spval
484  smstav(i,j)=spval
485  ssroff(i,j)=spval
486  bgroff(i,j)=spval
487  vegfrc(i,j)=spval
488  shdmin(i,j)=spval
489  shdmax(i,j)=spval
490  lai(i,j)=spval
491  acsnow(i,j)=spval
492  acgraup(i,j)=spval
493  acfrain(i,j)=spval
494  acsnom(i,j)=spval
495  cmc(i,j)=spval
496  sst(i,j)=spval
497  qz0(i,j)=spval
498  thz0(i,j)=spval
499  uz0(i,j)=spval
500  vz0(i,j)=spval
501  qs(i,j)=spval
502  ths(i,j)=spval
503  sno(i,j)=spval
504  snonc(i,j)=spval
505  ti(i,j)=spval
506  enddo
507  enddo
508 ! Time-averaged fileds
509  allocate(u10mean(im,jsta_2l:jend_2u))
510  allocate(v10mean(im,jsta_2l:jend_2u))
511  allocate(spduv10mean(im,jsta_2l:jend_2u))
512  allocate(swradmean(im,jsta_2l:jend_2u))
513  allocate(swnormmean(im,jsta_2l:jend_2u))
514 !Initialization
515 !$omp parallel do private(i,j)
516  do j=jsta_2l,jend_2u
517  do i=1,im
518  u10mean(i,j)=spval
519  v10mean(i,j)=spval
520  spduv10mean(i,j)=spval
521  swradmean(i,j)=spval
522  swnormmean(i,j)=spval
523  enddo
524  enddo
525 !NAMstart
526  allocate(snoavg(im,jsta_2l:jend_2u))
527  allocate(psfcavg(im,jsta_2l:jend_2u))
528  allocate(t10m(im,jsta_2l:jend_2u))
529  allocate(t10avg(im,jsta_2l:jend_2u))
530  allocate(akmsavg(im,jsta_2l:jend_2u))
531  allocate(akhsavg(im,jsta_2l:jend_2u))
532  allocate(u10max(im,jsta_2l:jend_2u))
533  allocate(v10max(im,jsta_2l:jend_2u))
534  allocate(u10h(im,jsta_2l:jend_2u))
535  allocate(v10h(im,jsta_2l:jend_2u))
536 !Initialization
537 !$omp parallel do private(i,j)
538  do j=jsta_2l,jend_2u
539  do i=1,im
540  snoavg(i,j)=spval
541  psfcavg(i,j)=spval
542  t10m(i,j)=spval
543  t10avg(i,j)=spval
544  akmsavg(i,j)=spval
545  akhsavg(i,j)=spval
546  u10max(i,j)=spval
547  v10max(i,j)=spval
548  u10h(i,j)=spval
549  v10h(i,j)=spval
550  enddo
551  enddo
552 !NAMend
553  allocate(akms(im,jsta_2l:jend_2u))
554  allocate(akhs(im,jsta_2l:jend_2u))
555  allocate(cuprec(im,jsta_2l:jend_2u))
556  allocate(acprec(im,jsta_2l:jend_2u))
557  allocate(ancprc(im,jsta_2l:jend_2u))
558  allocate(cuppt(im,jsta_2l:jend_2u))
559 !Initialization
560 !$omp parallel do private(i,j)
561  do j=jsta_2l,jend_2u
562  do i=1,im
563  akms(i,j)=spval
564  akhs(i,j)=spval
565  cuprec(i,j)=spval
566  acprec(i,j)=spval
567  ancprc(i,j)=spval
568  cuppt(i,j)=spval
569  enddo
570  enddo
571 ! GSDstart
572  allocate(rainc_bucket(im,jsta_2l:jend_2u))
573  allocate(rainc_bucket1(im,jsta_2l:jend_2u))
574  allocate(rainnc_bucket(im,jsta_2l:jend_2u))
575  allocate(rainnc_bucket1(im,jsta_2l:jend_2u))
576  allocate(pcp_bucket(im,jsta_2l:jend_2u))
577  allocate(pcp_bucket1(im,jsta_2l:jend_2u))
578  allocate(snow_bucket(im,jsta_2l:jend_2u))
579  allocate(snow_bucket1(im,jsta_2l:jend_2u))
580  allocate(graup_bucket(im,jsta_2l:jend_2u))
581  allocate(graup_bucket1(im,jsta_2l:jend_2u))
582  allocate(qrmax(im,jsta_2l:jend_2u))
583  allocate(tmax(im,jsta_2l:jend_2u))
584  allocate(snownc(im,jsta_2l:jend_2u))
585  allocate(graupelnc(im,jsta_2l:jend_2u))
586  allocate(tsnow(im,jsta_2l:jend_2u))
587  allocate(qvg(im,jsta_2l:jend_2u))
588  allocate(qv2m(im,jsta_2l:jend_2u))
589  allocate(qvl1(im,jsta_2l:jend_2u))
590  allocate(snfden(im,jsta_2l:jend_2u))
591  allocate(sndepac(im,jsta_2l:jend_2u))
592  allocate(int_smoke(im,jsta_2l:jend_2u))
593  allocate(mean_frp(im,jsta_2l:jend_2u))
594  allocate(int_aod(im,jsta_2l:jend_2u))
595 !Initialization
596 !$omp parallel do private(i,j)
597  do j=jsta_2l,jend_2u
598  do i=1,im
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
609  qrmax(i,j)=spval
610  tmax(i,j)=spval
611  snownc(i,j)=spval
612  graupelnc(i,j)=spval
613  tsnow(i,j)=spval
614  qvg(i,j)=spval
615  qv2m(i,j)=spval
616  qvl1(i,j)=spval
617  snfden(i,j)=spval
618  sndepac(i,j)=spval
619  int_smoke(i,j)=spval
620  mean_frp(i,j)=spval
621  int_aod(i,j)=spval
622  enddo
623  enddo
624  allocate(smoke(im,jsta_2l:jend_2u,lm,nbin_sm))
625 !$omp parallel do private(i,j,l,k)
626  do k=1,nbin_sm
627  do l=1,lm
628  do j=jsta_2l,jend_2u
629  do i=1,im
630  smoke(i,j,l,k)=spval
631  enddo
632  enddo
633  enddo
634  enddo
635 ! GSDend
636  allocate(rswin(im,jsta_2l:jend_2u))
637  allocate(swddni(im,jsta_2l:jend_2u))
638  allocate(swddif(im,jsta_2l:jend_2u))
639  allocate(swdnbc(im,jsta_2l:jend_2u))
640  allocate(swddnic(im,jsta_2l:jend_2u))
641  allocate(swddifc(im,jsta_2l:jend_2u))
642  allocate(swupbc(im,jsta_2l:jend_2u))
643  allocate(swupt(im,jsta_2l:jend_2u))
644  allocate(taod5502d(im,jsta_2l:jend_2u))
645  allocate(aerasy2d(im,jsta_2l:jend_2u))
646  allocate(aerssa2d(im,jsta_2l:jend_2u))
647  allocate(lwp(im,jsta_2l:jend_2u))
648  allocate(iwp(im,jsta_2l:jend_2u))
649  allocate(rlwin(im,jsta_2l:jend_2u))
650  allocate(lwdnbc(im,jsta_2l:jend_2u))
651  allocate(lwupbc(im,jsta_2l:jend_2u))
652  allocate(rlwtoa(im,jsta_2l:jend_2u))
653  allocate(rswtoa(im,jsta_2l:jend_2u))
654 !Initialization
655 !$omp parallel do private(i,j)
656  do j=jsta_2l,jend_2u
657  do i=1,im
658  rswin(i,j)=spval
659  swddni(i,j)=spval
660  swddif(i,j)=spval
661  swdnbc(i,j)=spval
662  swddnic(i,j)=spval
663  swddifc(i,j)=spval
664  swupbc(i,j)=spval
665  swupt(i,j)=spval
666  taod5502d(i,j)=spval
667  aerasy2d(i,j)=spval
668  aerssa2d(i,j)=spval
669  lwp(i,j)=spval
670  iwp(i,j)=spval
671  rlwin(i,j)=spval
672  lwdnbc(i,j)=spval
673  lwupbc(i,j)=spval
674  rlwtoa(i,j)=spval
675  rswtoa(i,j)=spval
676  enddo
677  enddo
678  allocate(tg(im,jsta_2l:jend_2u))
679  allocate(sfcshx(im,jsta_2l:jend_2u))
680  allocate(sfclhx(im,jsta_2l:jend_2u))
681  allocate(fis(im,jsta_2l:jend_2u))
682  allocate(t500(im,jsta_2l:jend_2u))
683  allocate(t700(im,jsta_2l:jend_2u))
684  allocate(z500(im,jsta_2l:jend_2u))
685  allocate(z700(im,jsta_2l:jend_2u))
686  allocate(teql(im,jsta_2l:jend_2u))
687  allocate(ieql(im,jsta_2l:jend_2u))
688  allocate(cfracl(im,jsta_2l:jend_2u))
689  allocate(cfracm(im,jsta_2l:jend_2u))
690  allocate(cfrach(im,jsta_2l:jend_2u))
691  allocate(acfrst(im,jsta_2l:jend_2u))
692  allocate(acfrcv(im,jsta_2l:jend_2u))
693  allocate(hbot(im,jsta_2l:jend_2u))
694  allocate(htop(im,jsta_2l:jend_2u))
695  allocate(aswin(im,jsta_2l:jend_2u))
696  allocate(alwin(im,jsta_2l:jend_2u))
697  allocate(aswout(im,jsta_2l:jend_2u))
698  allocate(alwout(im,jsta_2l:jend_2u))
699  allocate(aswtoa(im,jsta_2l:jend_2u))
700  allocate(alwtoa(im,jsta_2l:jend_2u))
701 !Initialization
702 !$omp parallel do private(i,j)
703  do j=jsta_2l,jend_2u
704  do i=1,im
705  tg(i,j)=spval
706  sfcshx(i,j)=spval
707  sfclhx(i,j)=spval
708  fis(i,j)=spval
709  t500(i,j)=spval
710  t700(i,j)=spval
711  z700(i,j)=spval
712  teql(i,j)=spval
713  ieql(i,j)=0
714  cfracl(i,j)=spval
715  cfracm(i,j)=spval
716  cfrach(i,j)=spval
717  acfrst(i,j)=spval
718  acfrcv(i,j)=spval
719  hbot(i,j)=spval
720  htop(i,j)=spval
721  aswin(i,j)=spval
722  alwin(i,j)=spval
723  aswout(i,j)=spval
724  alwout(i,j)=spval
725  aswtoa(i,j)=spval
726  alwtoa(i,j)=spval
727  enddo
728  enddo
729  allocate(czen(im,jsta_2l:jend_2u))
730  allocate(czmean(im,jsta_2l:jend_2u))
731  allocate(sigt4(im,jsta_2l:jend_2u))
732  allocate(rswout(im,jsta_2l:jend_2u))
733  allocate(radot(im,jsta_2l:jend_2u))
734  allocate(ncfrst(im,jsta_2l:jend_2u)) ! real
735  allocate(ncfrcv(im,jsta_2l:jend_2u)) ! real
736  allocate(smstot(im,jsta_2l:jend_2u))
737  allocate(pctsno(im,jsta_2l:jend_2u))
738  allocate(pshltr(im,jsta_2l:jend_2u))
739  allocate(th10(im,jsta_2l:jend_2u))
740  allocate(q10(im,jsta_2l:jend_2u))
741  allocate(sr(im,jsta_2l:jend_2u))
742  allocate(prec(im,jsta_2l:jend_2u))
743  allocate(subshx(im,jsta_2l:jend_2u))
744  allocate(snopcx(im,jsta_2l:jend_2u))
745  allocate(sfcuvx(im,jsta_2l:jend_2u))
746  allocate(sfcevp(im,jsta_2l:jend_2u))
747  allocate(potevp(im,jsta_2l:jend_2u))
748  allocate(z0(im,jsta_2l:jend_2u))
749  allocate(ustar(im,jsta_2l:jend_2u))
750  allocate(pblh(im,jsta_2l:jend_2u))
751  allocate(pblhgust(im,jsta_2l:jend_2u))
752  allocate(mixht(im,jsta_2l:jend_2u))
753  allocate(twbs(im,jsta_2l:jend_2u))
754  allocate(qwbs(im,jsta_2l:jend_2u))
755 !Initialization
756 !$omp parallel do private(i,j)
757  do j=jsta_2l,jend_2u
758  do i=1,im
759  czen(i,j)=spval
760  czmean(i,j)=spval
761  sigt4(i,j)=spval
762  rswout(i,j)=spval
763  radot(i,j)=spval
764  ncfrst(i,j)=spval
765  ncfrcv(i,j)=spval
766  smstot(i,j)=spval
767  pctsno(i,j)=spval
768  pshltr(i,j)=spval
769  th10(i,j)=spval
770  q10(i,j)=spval
771  sr(i,j)=spval
772  prec(i,j)=spval
773  subshx(i,j)=spval
774  snopcx(i,j)=spval
775  sfcuvx(i,j)=spval
776  sfcevp(i,j)=spval
777  potevp(i,j)=spval
778  z0(i,j)=spval
779  ustar(i,j)=spval
780  pblh(i,j)=spval
781  pblhgust(i,j)=spval
782  mixht(i,j)=spval
783  twbs(i,j)=spval
784  qwbs(i,j)=spval
785  enddo
786  enddo
787  allocate(sfcexc(im,jsta_2l:jend_2u))
788  allocate(grnflx(im,jsta_2l:jend_2u))
789  allocate(soiltb(im,jsta_2l:jend_2u))
790  allocate(z1000(im,jsta_2l:jend_2u))
791  allocate(slp(im,jsta_2l:jend_2u))
792  allocate(pslp(im,jsta_2l:jend_2u))
793  allocate(f(im,jsta_2l:jend_2u))
794  allocate(albedo(im,jsta_2l:jend_2u))
795  allocate(albase(im,jsta_2l:jend_2u))
796  allocate(cldfra(im,jsta_2l:jend_2u))
797  allocate(cprate(im,jsta_2l:jend_2u))
798  allocate(cnvcfr(im,jsta_2l:jend_2u))
799  allocate(ivgtyp(im,jsta_2l:jend_2u))
800  allocate(isltyp(im,jsta_2l:jend_2u))
801  allocate(hbotd(im,jsta_2l:jend_2u))
802  allocate(htopd(im,jsta_2l:jend_2u))
803  allocate(hbots(im,jsta_2l:jend_2u))
804  allocate(htops(im,jsta_2l:jend_2u))
805  allocate(cldefi(im,jsta_2l:jend_2u))
806  allocate(islope(im,jsta_2l:jend_2u))
807  allocate(si(im,jsta_2l:jend_2u))
808  allocate(lspa(im,jsta_2l:jend_2u))
809  allocate(rswinc(im,jsta_2l:jend_2u))
810  allocate(vis(im,jsta_2l:jend_2u))
811  allocate(pd(im,jsta_2l:jend_2u))
812  allocate(mxsnal(im,jsta_2l:jend_2u))
813  allocate(epsr(im,jsta_2l:jend_2u))
814 !Initialization
815 !$omp parallel do private(i,j)
816  do j=jsta_2l,jend_2u
817  do i=1,im
818  sfcexc(i,j)=spval
819  grnflx(i,j)=spval
820  soiltb(i,j)=spval
821  z1000(i,j)=spval
822  slp(i,j)=spval
823  pslp(i,j)=spval
824  f(i,j)=spval
825  albedo(i,j)=spval
826  albase(i,j)=spval
827  cldfra(i,j)=spval
828  cprate(i,j)=spval
829  cnvcfr(i,j)=spval
830  ivgtyp(i,j)=0
831  isltyp(i,j)=0
832  hbotd(i,j)=spval
833  htopd(i,j)=spval
834  hbots(i,j)=spval
835  htops(i,j)=spval
836  cldefi(i,j)=spval
837  islope(i,j)=0
838  si(i,j)=spval
839  lspa(i,j)=spval
840  rswinc(i,j)=spval
841  vis(i,j)=spval
842  pd(i,j)=spval
843  mxsnal(i,j)=spval
844  epsr(i,j)=spval
845  enddo
846  enddo
847 ! add GFS fields
848  allocate(sfcux(im,jsta_2l:jend_2u))
849  allocate(sfcvx(im,jsta_2l:jend_2u))
850  allocate(sfcuxi(im,jsta_2l:jend_2u))
851  allocate(sfcvxi(im,jsta_2l:jend_2u))
852  allocate(avgalbedo(im,jsta_2l:jend_2u))
853  allocate(avgcprate(im,jsta_2l:jend_2u))
854  allocate(avgprec(im,jsta_2l:jend_2u))
855  allocate(avgprec_cont(im,jsta_2l:jend_2u))
856  allocate(avgcprate_cont(im,jsta_2l:jend_2u))
857  allocate(ptop(im,jsta_2l:jend_2u))
858  allocate(pbot(im,jsta_2l:jend_2u))
859  allocate(avgcfrach(im,jsta_2l:jend_2u))
860  allocate(avgcfracm(im,jsta_2l:jend_2u))
861  allocate(avgcfracl(im,jsta_2l:jend_2u))
862  allocate(avgtcdc(im,jsta_2l:jend_2u))
863  allocate(auvbin(im,jsta_2l:jend_2u))
864  allocate(auvbinc(im,jsta_2l:jend_2u))
865  allocate(ptopl(im,jsta_2l:jend_2u))
866  allocate(pbotl(im,jsta_2l:jend_2u))
867  allocate(ttopl(im,jsta_2l:jend_2u))
868  allocate(ptopm(im,jsta_2l:jend_2u))
869  allocate(pbotm(im,jsta_2l:jend_2u))
870  allocate(ttopm(im,jsta_2l:jend_2u))
871  allocate(ptoph(im,jsta_2l:jend_2u))
872  allocate(pboth(im,jsta_2l:jend_2u))
873  allocate(ttoph(im,jsta_2l:jend_2u))
874  allocate(sfcugs(im,jsta_2l:jend_2u))
875  allocate(sfcvgs(im,jsta_2l:jend_2u))
876  allocate(pblcfr(im,jsta_2l:jend_2u))
877  allocate(cldwork(im,jsta_2l:jend_2u))
878  allocate(gtaux(im,jsta_2l:jend_2u))
879  allocate(gtauy(im,jsta_2l:jend_2u))
880  allocate(cd10(im,jsta_2l:jend_2u))
881  allocate(ch10(im,jsta_2l:jend_2u))
882  allocate(mdltaux(im,jsta_2l:jend_2u))
883  allocate(mdltauy(im,jsta_2l:jend_2u))
884  allocate(runoff(im,jsta_2l:jend_2u))
885 !Initialization
886 !$omp parallel do private(i,j)
887  do j=jsta_2l,jend_2u
888  do i=1,im
889  sfcux(i,j)=spval
890  sfcvx(i,j)=spval
891  sfcuxi(i,j)=spval
892  sfcvxi(i,j)=spval
893  avgalbedo(i,j)=spval
894  avgcprate(i,j)=spval
895  avgprec(i,j)=spval
896  avgprec_cont(i,j)=spval
897  avgcprate_cont(i,j)=spval
898  ptop(i,j)=spval
899  pbot(i,j)=spval
900  avgcfrach(i,j)=spval
901  avgcfracm(i,j)=spval
902  avgcfracl(i,j)=spval
903  avgtcdc(i,j)=spval
904  auvbin(i,j)=spval
905  auvbinc(i,j)=spval
906  ptopl(i,j)=spval
907  pbotl(i,j)=spval
908  ttopl(i,j)=spval
909  ptopm(i,j)=spval
910  pbotm(i,j)=spval
911  ttopm(i,j)=spval
912  ptoph(i,j)=spval
913  pboth(i,j)=spval
914  ttoph(i,j)=spval
915  sfcugs(i,j)=spval
916  sfcvgs(i,j)=spval
917  pblcfr(i,j)=spval
918  cldwork(i,j)=spval
919  gtaux(i,j)=spval
920  gtauy(i,j)=spval
921  cd10(i,j)=spval
922  ch10(i,j)=spval
923  mdltaux(i,j)=spval
924  mdltauy(i,j)=spval
925  runoff(i,j)=spval
926  enddo
927  enddo
928  allocate(maxtshltr(im,jsta_2l:jend_2u))
929  allocate(mintshltr(im,jsta_2l:jend_2u))
930  allocate(maxrhshltr(im,jsta_2l:jend_2u))
931  allocate(minrhshltr(im,jsta_2l:jend_2u))
932  allocate(maxqshltr(im,jsta_2l:jend_2u))
933  allocate(minqshltr(im,jsta_2l:jend_2u))
934  allocate(dzice(im,jsta_2l:jend_2u))
935  allocate(alwinc(im,jsta_2l:jend_2u))
936  allocate(alwoutc(im,jsta_2l:jend_2u))
937  allocate(alwtoac(im,jsta_2l:jend_2u))
938  allocate(aswinc(im,jsta_2l:jend_2u))
939  allocate(aswoutc(im,jsta_2l:jend_2u))
940  allocate(aswtoac(im,jsta_2l:jend_2u))
941  allocate(aswintoa(im,jsta_2l:jend_2u))
942  allocate(smcwlt(im,jsta_2l:jend_2u))
943  allocate(suntime(im,jsta_2l:jend_2u))
944  allocate(fieldcapa(im,jsta_2l:jend_2u))
945  allocate(avisbeamswin(im,jsta_2l:jend_2u))
946  allocate(avisdiffswin(im,jsta_2l:jend_2u))
947  allocate(airbeamswin(im,jsta_2l:jend_2u))
948  allocate(airdiffswin(im,jsta_2l:jend_2u))
949  allocate(snowfall(im,jsta_2l:jend_2u))
950  allocate(acond(im,jsta_2l:jend_2u))
951  allocate(edir(im,jsta_2l:jend_2u))
952  allocate(ecan(im,jsta_2l:jend_2u))
953  allocate(etrans(im,jsta_2l:jend_2u))
954  allocate(esnow(im,jsta_2l:jend_2u))
955  allocate(avgedir(im,jsta_2l:jend_2u))
956  allocate(avgecan(im,jsta_2l:jend_2u))
957  allocate(avgetrans(im,jsta_2l:jend_2u))
958  allocate(avgesnow(im,jsta_2l:jend_2u))
959  allocate(avgpotevp(im,jsta_2l:jend_2u))
960  allocate(aod550(im,jsta_2l:jend_2u))
961  allocate(du_aod550(im,jsta_2l:jend_2u))
962  allocate(ss_aod550(im,jsta_2l:jend_2u))
963  allocate(su_aod550(im,jsta_2l:jend_2u))
964  allocate(oc_aod550(im,jsta_2l:jend_2u))
965  allocate(bc_aod550(im,jsta_2l:jend_2u))
966  allocate(landfrac(im,jsta_2l:jend_2u))
967  allocate(paha(im,jsta_2l:jend_2u))
968  allocate(pahi(im,jsta_2l:jend_2u))
969  allocate(tecan(im,jsta_2l:jend_2u))
970  allocate(tetran(im,jsta_2l:jend_2u))
971  allocate(tedir(im,jsta_2l:jend_2u))
972  allocate(twa(im,jsta_2l:jend_2u))
973  allocate(fdnsst(im,jsta_2l:jend_2u))
974  allocate(pwat(im,jsta_2l:jend_2u))
975 !Initialization
976 !$omp parallel do private(i,j)
977  do j=jsta_2l,jend_2u
978  do i=1,im
979  maxtshltr(i,j)=spval
980  mintshltr(i,j)=spval
981  maxrhshltr(i,j)=spval
982  minrhshltr(i,j)=spval
983  maxqshltr(i,j)=spval
984  minqshltr(i,j)=spval
985  dzice(i,j)=spval
986  alwinc(i,j)=spval
987  alwoutc(i,j)=spval
988  alwtoac(i,j)=spval
989  aswinc(i,j)=spval
990  aswoutc(i,j)=spval
991  aswtoac(i,j)=spval
992  aswintoa(i,j)=spval
993  smcwlt(i,j)=spval
994  suntime(i,j)=spval
995  fieldcapa(i,j)=spval
996  avisbeamswin(i,j)=spval
997  avisdiffswin(i,j)=spval
998  airbeamswin(i,j)=spval
999  airdiffswin(i,j)=spval
1000  snowfall(i,j)=spval
1001  acond(i,j)=spval
1002  edir(i,j)=spval
1003  ecan(i,j)=spval
1004  etrans(i,j)=spval
1005  esnow(i,j)=spval
1006  avgedir(i,j)=spval
1007  avgecan(i,j)=spval
1008  avgetrans(i,j)=spval
1009  avgesnow(i,j)=spval
1010  avgpotevp(i,j)=spval
1011  aod550(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
1017  landfrac(i,j)=spval
1018  paha(i,j)=spval
1019  pahi(i,j)=spval
1020  tecan(i,j)=spval
1021  tetran(i,j)=spval
1022  tedir(i,j)=spval
1023  twa(i,j)=spval
1024  fdnsst(i,j)=spval
1025  pwat(i,j)=spval
1026  enddo
1027  enddo
1028 !
1029 ! FROM MASKS
1030 !
1031  allocate(hbm2(im,jsta_2l:jend_2u))
1032  allocate(sm(im,jsta_2l:jend_2u))
1033  allocate(sice(im,jsta_2l:jend_2u))
1034  allocate(lmh(im,jsta_2l:jend_2u)) ! real
1035  allocate(lmv(im,jsta_2l:jend_2u)) ! real
1036  allocate(gdlat(im,jsta_2l:jend_2u))
1037  allocate(gdlon(im,jsta_2l:jend_2u))
1038  allocate(dx(im,jsta_2l:jend_2u))
1039  allocate(dy(im,jsta_2l:jend_2u))
1040 !Initialization
1041 !$omp parallel do private(i,j)
1042  do j=jsta_2l,jend_2u
1043  do i=1,im
1044  hbm2(i,j)=spval
1045  sm(i,j)=spval
1046  sice(i,j)=spval
1047  lmh(i,j)=spval
1048  lmv(i,j)=spval
1049  gdlat(i,j)=spval
1050  gdlon(i,j)=spval
1051  dx(i,j)=spval
1052  dy(i,j)=spval
1053  enddo
1054  enddo
1055 
1056  if (me == 0) print *,' gocart_on=',gocart_on
1057  if (gocart_on) then
1058 !
1059 ! Add GOCART fields
1060 ! vrbls4d
1061  allocate(dust(im,jsta_2l:jend_2u,lm,nbin_du))
1062  allocate(salt(im,jsta_2l:jend_2u,lm,nbin_ss))
1063  allocate(soot(im,jsta_2l:jend_2u,lm,nbin_bc))
1064  allocate(waso(im,jsta_2l:jend_2u,lm,nbin_oc))
1065  allocate(suso(im,jsta_2l:jend_2u,lm,nbin_su))
1066  allocate(pp25(im,jsta_2l:jend_2u,lm,nbin_su))
1067  allocate(pp10(im,jsta_2l:jend_2u,lm,nbin_su))
1068 !Initialization
1069 !$omp parallel do private(i,j,l,k)
1070  do k=1,nbin_du
1071  do l=1,lm
1072  do j=jsta_2l,jend_2u
1073  do i=1,im
1074  dust(i,j,l,k)=spval
1075  enddo
1076  enddo
1077  enddo
1078  enddo
1079 !$omp parallel do private(i,j,l,k)
1080  do k=1,nbin_ss
1081  do l=1,lm
1082  do j=jsta_2l,jend_2u
1083  do i=1,im
1084  salt(i,j,l,k)=spval
1085  enddo
1086  enddo
1087  enddo
1088  enddo
1089 !$omp parallel do private(i,j,l,k)
1090  do k=1,nbin_bc
1091  do l=1,lm
1092  do j=jsta_2l,jend_2u
1093  do i=1,im
1094  soot(i,j,l,k)=spval
1095  enddo
1096  enddo
1097  enddo
1098  enddo
1099 !$omp parallel do private(i,j,l,k)
1100  do k=1,nbin_oc
1101  do l=1,lm
1102  do j=jsta_2l,jend_2u
1103  do i=1,im
1104  waso(i,j,l,k)=spval
1105  enddo
1106  enddo
1107  enddo
1108  enddo
1109 !$omp parallel do private(i,j,l,k)
1110  do k=1,nbin_su
1111  do l=1,lm
1112  do j=jsta_2l,jend_2u
1113  do i=1,im
1114  suso(i,j,l,k)=spval
1115  pp25(i,j,l,k)=spval
1116  pp10(i,j,l,k)=spval
1117  enddo
1118  enddo
1119  enddo
1120  enddo
1121 ! vrbls3d
1122  allocate(ext(im,jsta_2l:jend_2u,lm))
1123  allocate(asy(im,jsta_2l:jend_2u,lm))
1124  allocate(ssa(im,jsta_2l:jend_2u,lm))
1125  allocate(sca(im,jsta_2l:jend_2u,lm))
1126 !Initialization
1127 !$omp parallel do private(i,j)
1128  do l=1,lm
1129  do j=jsta_2l,jend_2u
1130  do i=1,im
1131  ext(i,j,l)=spval
1132  asy(i,j,l)=spval
1133  ssa(i,j,l)=spval
1134  sca(i,j,l)=spval
1135  enddo
1136  enddo
1137  enddo
1138  allocate(duem(im,jsta_2l:jend_2u,nbin_du))
1139  allocate(dusd(im,jsta_2l:jend_2u,nbin_du))
1140  allocate(dudp(im,jsta_2l:jend_2u,nbin_du))
1141  allocate(duwt(im,jsta_2l:jend_2u,nbin_du))
1142  allocate(dusv(im,jsta_2l:jend_2u,nbin_du))
1143  allocate(suem(im,jsta_2l:jend_2u,nbin_su))
1144  allocate(susd(im,jsta_2l:jend_2u,nbin_su))
1145  allocate(sudp(im,jsta_2l:jend_2u,nbin_su))
1146  allocate(suwt(im,jsta_2l:jend_2u,nbin_su))
1147  allocate(ocem(im,jsta_2l:jend_2u,nbin_oc))
1148  allocate(ocsd(im,jsta_2l:jend_2u,nbin_oc))
1149  allocate(ocdp(im,jsta_2l:jend_2u,nbin_oc))
1150  allocate(ocwt(im,jsta_2l:jend_2u,nbin_oc))
1151  allocate(ocsv(im,jsta_2l:jend_2u,nbin_oc))
1152  allocate(bcem(im,jsta_2l:jend_2u,nbin_bc))
1153  allocate(bcsd(im,jsta_2l:jend_2u,nbin_bc))
1154  allocate(bcdp(im,jsta_2l:jend_2u,nbin_bc))
1155  allocate(bcwt(im,jsta_2l:jend_2u,nbin_bc))
1156  allocate(bcsv(im,jsta_2l:jend_2u,nbin_bc))
1157  allocate(ssem(im,jsta_2l:jend_2u,nbin_ss))
1158  allocate(sssd(im,jsta_2l:jend_2u,nbin_ss))
1159  allocate(ssdp(im,jsta_2l:jend_2u,nbin_ss))
1160  allocate(sswt(im,jsta_2l:jend_2u,nbin_ss))
1161  allocate(sssv(im,jsta_2l:jend_2u,nbin_ss))
1162 !Initialization
1163 !$omp parallel do private(i,j,l)
1164  do l=1,nbin_du
1165  do j=jsta_2l,jend_2u
1166  do i=1,im
1167  duem(i,j,l)=spval
1168  dusd(i,j,l)=spval
1169  dudp(i,j,l)=spval
1170  duwt(i,j,l)=spval
1171  dusv(i,j,l)=spval
1172  enddo
1173  enddo
1174  enddo
1175 
1176  do l=1,nbin_su
1177  do j=jsta_2l,jend_2u
1178  do i=1,im
1179  suem(i,j,l)=spval
1180  susd(i,j,l)=spval
1181  sudp(i,j,l)=spval
1182  suwt(i,j,l)=spval
1183  enddo
1184  enddo
1185  enddo
1186 
1187  do l=1,nbin_oc
1188  do j=jsta_2l,jend_2u
1189  do i=1,im
1190  ocem(i,j,l)=spval
1191  ocsd(i,j,l)=spval
1192  ocdp(i,j,l)=spval
1193  ocwt(i,j,l)=spval
1194  ocsv(i,j,l)=spval
1195  enddo
1196  enddo
1197  enddo
1198 
1199  do l=1,nbin_bc
1200  do j=jsta_2l,jend_2u
1201  do i=1,im
1202  bcem(i,j,l)=spval
1203  bcsd(i,j,l)=spval
1204  bcdp(i,j,l)=spval
1205  bcwt(i,j,l)=spval
1206  bcsv(i,j,l)=spval
1207  enddo
1208  enddo
1209  enddo
1210 
1211  do l=1,nbin_ss
1212  do j=jsta_2l,jend_2u
1213  do i=1,im
1214  ssem(i,j,l)=spval
1215  sssd(i,j,l)=spval
1216  ssdp(i,j,l)=spval
1217  sswt(i,j,l)=spval
1218  sssv(i,j,l)=spval
1219  enddo
1220  enddo
1221  enddo
1222  allocate(rhomid(im,jsta_2l:jend_2u,lm))
1223 !Initialization
1224 !$omp parallel do private(i,j,l)
1225  do l=1,lm
1226  do j=jsta_2l,jend_2u
1227  do i=1,im
1228  rhomid(i,j,l)=spval
1229  enddo
1230  enddo
1231  enddo
1232 ! vrbls2d
1233  allocate(dusmass(im,jsta_2l:jend_2u))
1234  allocate(ducmass(im,jsta_2l:jend_2u))
1235  allocate(dusmass25(im,jsta_2l:jend_2u))
1236  allocate(ducmass25(im,jsta_2l:jend_2u))
1237  allocate(susmass(im,jsta_2l:jend_2u))
1238  allocate(sucmass(im,jsta_2l:jend_2u))
1239  allocate(susmass25(im,jsta_2l:jend_2u))
1240  allocate(sucmass25(im,jsta_2l:jend_2u))
1241  allocate(ocsmass(im,jsta_2l:jend_2u))
1242  allocate(occmass(im,jsta_2l:jend_2u))
1243  allocate(ocsmass25(im,jsta_2l:jend_2u))
1244  allocate(occmass25(im,jsta_2l:jend_2u))
1245  allocate(bcsmass(im,jsta_2l:jend_2u))
1246  allocate(bccmass(im,jsta_2l:jend_2u))
1247  allocate(bcsmass25(im,jsta_2l:jend_2u))
1248  allocate(bccmass25(im,jsta_2l:jend_2u))
1249  allocate(sssmass(im,jsta_2l:jend_2u))
1250  allocate(sscmass(im,jsta_2l:jend_2u))
1251  allocate(sssmass25(im,jsta_2l:jend_2u))
1252  allocate(sscmass25(im,jsta_2l:jend_2u))
1253  allocate(dustcb(im,jsta_2l:jend_2u))
1254  allocate(occb(im,jsta_2l:jend_2u))
1255  allocate(bccb(im,jsta_2l:jend_2u))
1256  allocate(sulfcb(im,jsta_2l:jend_2u))
1257  allocate(pp25cb(im,jsta_2l:jend_2u))
1258  allocate(pp10cb(im,jsta_2l:jend_2u))
1259  allocate(sscb(im,jsta_2l:jend_2u))
1260  allocate(dustallcb(im,jsta_2l:jend_2u))
1261  allocate(ssallcb(im,jsta_2l:jend_2u))
1262  allocate(dustpm(im,jsta_2l:jend_2u))
1263  allocate(dustpm10(im,jsta_2l:jend_2u))
1264  allocate(sspm(im,jsta_2l:jend_2u))
1265  allocate(maod(im,jsta_2l:jend_2u))
1266 !Initialization
1267 !$omp parallel do private(i,j)
1268  do j=jsta_2l,jend_2u
1269  do i=1,im
1270  dusmass(i,j)=spval
1271  ducmass(i,j)=spval
1272  dusmass25(i,j)=spval
1273  ducmass25(i,j)=spval
1274  susmass(i,j)=spval
1275  sucmass(i,j)=spval
1276  susmass25(i,j)=spval
1277  sucmass25(i,j)=spval
1278  ocsmass(i,j)=spval
1279  occmass(i,j)=spval
1280  ocsmass25(i,j)=spval
1281  occmass25(i,j)=spval
1282  bcsmass(i,j)=spval
1283  bccmass(i,j)=spval
1284  bcsmass25(i,j)=spval
1285  bccmass25(i,j)=spval
1286  sssmass(i,j)=spval
1287  sscmass(i,j)=spval
1288  sssmass25(i,j)=spval
1289  sscmass25(i,j)=spval
1290  dustcb(i,j)=spval
1291  occb(i,j)=spval
1292  bccb(i,j)=spval
1293  sulfcb(i,j)=spval
1294  pp25cb(i,j)=spval
1295  pp10cb(i,j)=spval
1296  sscb(i,j)=spval
1297  dustallcb(i,j)=spval
1298  ssallcb(i,j)=spval
1299  dustpm(i,j)=spval
1300  dustpm10(i,j)=spval
1301  sspm(i,j)=spval
1302  maod(i,j)=spval
1303  enddo
1304  enddo
1305  endif
1306 ! HWRF RRTMG output
1307  allocate(acswupt(im,jsta_2l:jend_2u))
1308  allocate(swdnt(im,jsta_2l:jend_2u))
1309  allocate(acswdnt(im,jsta_2l:jend_2u))
1310 !Initialization
1311 !$omp parallel do private(i,j)
1312  do j=jsta_2l,jend_2u
1313  do i=1,im
1314  acswupt(i,j)=spval
1315  swdnt(i,j)=spval
1316  acswdnt(i,j)=spval
1317  enddo
1318  enddo
1319 
1320 ! UPP_MATH MODULE DIFFERENTIAL EQUATIONS
1321  allocate(ddvdx(im,jsta_2l:jend_2u))
1322  allocate(ddudy(im,jsta_2l:jend_2u))
1323  allocate(uuavg(im,jsta_2l:jend_2u))
1324 !Initialization
1325 !$omp parallel do private(i,j)
1326  do j=jsta_2l,jend_2u
1327  do i=1,im
1328  ddvdx(i,j)=spval
1329  ddudy(i,j)=spval
1330  uuavg(i,j)=spval
1331  enddo
1332  enddo
1333 
1334 ! AQF
1335  if (me == 0) print *,'aqfcmaq_on= ', aqfcmaq_on
1336  if (aqfcmaq_on) then
1337 
1338  allocate(ozcon(im,jsta_2l:jend_2u,lm))
1339  allocate(pmtf(im,jsta_2l:jend_2u,lm))
1340 
1341 !Initialization
1342 !$omp parallel do private(i,j,l)
1343  do l=1,lm
1344  do j=jsta_2l,jend_2u
1345  do i=1,im
1346  ozcon(i,j,l)=0.
1347  pmtf(i,j,l)=0.
1348  enddo
1349  enddo
1350  enddo
1351 
1352  endif
1353 !
1354  end
masks
Definition: MASKS_mod.f:1
soil
Definition: SOIL_mod.f:1
upp_math
dvdxdudy() computes dudy, dvdx, uwnd
Definition: UPP_MATH.f:17