56 SUBROUTINE calhel3(LLOW,LUPP,UST,VST,HELI)
59 use vrbls3d, only: zmid, uh, vh, u, v, zint
60 use vrbls2d, only: fis, u10, v10
64 use ctlblk_mod, only: jsta, jend, jsta_m, jend_m, jsta_2l, jend_2u, &
65 lm, im, jm, me, spval, &
66 ista, iend, ista_m, iend_m, ista_2l, iend_2u
67 use gridspec_mod
, only: gridtype
71 real,
PARAMETER :: p150=15000.0,p300=30000.0,s15=15.0
72 real,
PARAMETER :: d3000=3000.0,pi6=0.5235987756,pi9=0.34906585
73 real,
PARAMETER :: d5500=5500.0,d6000=6000.0,d7000=7000.0
74 real,
PARAMETER :: d500=500.0
76 real,
PARAMETER :: d1000=1000.0
77 real,
PARAMETER :: d1500=1500.0
79 REAL,
PARAMETER :: pi = 3.1415927
84 integer,
dimension(ista_2l:iend_2u,jsta_2l:jend_2u),
intent(in) :: llow, lupp
85 REAL,
dimension(ista_2l:iend_2u,jsta_2l:jend_2u),
intent(out) :: ust,vst
86 REAL,
dimension(ista_2l:iend_2u,jsta_2l:jend_2u),
intent(out) :: heli
88 real,
dimension(ista_2l:iend_2u,jsta_2l:jend_2u) :: htsfc, ust6, vst6, ust5, vst5, &
89 ust1, vst1, ushr1, vshr1, &
90 ushr6, vshr6, u1, v1, u2, v2, &
91 hgt1, hgt2, umean, vmean
92 real,
dimension(ista_2l:iend_2u,jsta_2l:jend_2u) :: ushr05,vshr05,elt,elb
105 integer,
dimension(ista_2l:iend_2u,jsta_2l:jend_2u) :: count6, count5, count1, l1, l2
111 INTEGER ive(jm),ivw(jm)
112 integer i,j,iw,ie,js,jn,jvn,jvs,l,n,lv
113 integer istart,istop,jstart,jstop
114 real z2,dzabv,umean5,vmean5,umean1,vmean1,umean6,vmean6, &
115 denom,z1,z3,dz,dz1,dz2,du1,du2,dv1,dv2
158 IF(gridtype ==
'E')
THEN
169 ELSE IF(gridtype ==
'B')
THEN
203 IF(gridtype /=
'A') CALL exch(fis(ista_2l:iend_2u,jsta_2l:jend_2u))
205 IF(gridtype /=
'A') CALL exch(zmid(ista_2l:iend_2u,jsta_2l:jend_2u,l))
215 IF (gridtype==
'B')
THEN
216 htsfc(i,j) = (0.25/g)*(fis(iw,j)+fis(ie,j)+fis(i,jn)+fis(ie,jn))
221 z2 = 0.25*(zmid(iw,j,l)+zmid(ie,j,l)+zmid(i,jn,l)+zmid(ie,jn,l))
223 htsfc(i,j) = (0.25/g)*(fis(iw,j)+fis(ie,j)+fis(i,jn)+fis(i,js))
228 z2 = 0.25*(zmid(iw,j,l)+zmid(ie,j,l)+zmid(i,jn,l)+zmid(i,js,l))
230 dzabv = z2-htsfc(i,j)
233 IF (dzabv <= d6000 .AND. l <= lv)
THEN
234 ust6(i,j) = ust6(i,j) + uh(i,j,l)
235 vst6(i,j) = vst6(i,j) + vh(i,j,l)
236 count6(i,j) = count6(i,j) + 1
239 IF (dzabv < d6000 .AND. dzabv >= d5500 .AND. l <= lv)
THEN
240 ust5(i,j) = ust5(i,j) + uh(i,j,l)
241 vst5(i,j) = vst5(i,j) + vh(i,j,l)
242 count5(i,j) = count5(i,j) + 1
245 IF (dzabv < d500 .AND. l <= lv)
THEN
246 ust1(i,j) = ust1(i,j) + uh(i,j,l)
247 vst1(i,j) = vst1(i,j) + vh(i,j,l)
248 count1(i,j) = count1(i,j) + 1
251 IF (dzabv >= d1000 .AND. dzabv <= d1500 .AND. l <= lv)
THEN
258 IF (dzabv >= d500 .AND. dzabv < d1000 .AND. &
259 l <= lv .AND. l1(i,j) <= l2(i,j))
THEN
275 IF (count5(i,j) == 0)
THEN
281 IF (gridtype==
'B')
THEN
282 z2 = 0.25*(zmid(iw,j,l)+zmid(ie,j,l)+zmid(i,jn,l)+zmid(ie,jn,l))
284 z2 = 0.25*(zmid(iw,j,l)+zmid(ie,j,l)+zmid(i,jn,l)+zmid(i,js,l))
289 IF (dzabv < d7000 .AND. dzabv >= d6000)
THEN
290 ust5(i,j) = ust5(i,j) + uh(i,j,l)
291 vst5(i,j) = vst5(i,j) + vh(i,j,l)
306 IF (count6(i,j) > 0 .AND. count1(i,j) > 0 .AND. count5(i,j) > 0)
THEN
307 umean5 = ust5(i,j) / count5(i,j)
308 vmean5 = vst5(i,j) / count5(i,j)
309 umean1 = ust1(i,j) / count1(i,j)
310 vmean1 = vst1(i,j) / count1(i,j)
311 umean6 = ust6(i,j) / count6(i,j)
312 vmean6 = vst6(i,j) / count6(i,j)
325 ushr6(i,j) = umean5 - umean1
326 vshr6(i,j) = vmean5 - vmean1
328 denom = ushr6(i,j)*ushr6(i,j)+vshr6(i,j)*vshr6(i,j)
329 IF (denom /= 0.0)
THEN
330 ust(i,j) = umean6 + (7.5*vshr6(i,j)/sqrt(denom))
331 vst(i,j) = vmean6 - (7.5*ushr6(i,j)/sqrt(denom))
343 IF(l1(i,j) > 0 .AND. l2(i,j) > 0)
THEN
344 umean(i,j) = u1(i,j) + (d1000 - hgt1(i,j))*(u2(i,j) - &
345 u1(i,j))/(hgt2(i,j) - hgt1(i,j))
346 vmean(i,j) = v1(i,j) + (d1000 - hgt1(i,j))*(v2(i,j) - &
347 v1(i,j))/(hgt2(i,j) - hgt1(i,j))
348 ELSE IF(l1(i,j) > 0 .AND. l2(i,j) == 0)
THEN
351 ELSE IF(l1(i,j) == 0 .AND. l2(i,j) > 0)
THEN
359 IF(l1(i,j) > 0 .OR. l2(i,j) > 0)
THEN
360 ushr05(i,j) = umean1 - u10(i,j)
361 vshr05(i,j) = vmean1 - v10(i,j)
362 ushr1(i,j) = umean(i,j) - u10(i,j)
363 vshr1(i,j) = vmean(i,j) - v10(i,j)
388 if(gridtype /=
'A')
then
389 call exch(zint(1,jsta_2l,l))
390 call exch(zint(1,jsta_2l,l+1))
398 IF (gridtype==
'B')
THEN
399 z2=0.25*(zmid(iw,j,l)+zmid(ie,j,l)+ &
400 zmid(i,jn,l)+zmid(ie,jn,l))
402 z2=0.25*(zmid(iw,j,l)+zmid(ie,j,l)+ &
403 zmid(i,jn,l)+zmid(i,js,l))
406 elt(i,j) = zint(i,j,lupp(i,j))-htsfc(i,j)
407 elb(i,j) = zint(i,j,llow(i,j))-htsfc(i,j)
410 IF(dzabv <= elt(i,j) .AND. dzabv >= elb(i,j) .AND. l <= nint(lmv(i,j)))
THEN
411 IF (gridtype==
'B')
THEN
412 z1 = 0.25*(zmid(iw,j,l+1)+zmid(ie,j,l+1)+ &
413 zmid(i,jn,l+1)+zmid(ie,jn,l+1))
414 z3 = 0.25*(zmid(iw,j,l-1)+zmid(ie,j,l-1)+ &
415 zmid(i,jn,l-1)+zmid(ie,jn,l-1))
416 dz = 0.25*((zint(iw,j,l)+zint(ie,j,l)+ &
417 zint(i,jn,l)+zint(ie,jn,l))- &
418 (zint(iw,j,l+1)+zint(ie,j,l+1)+ &
419 zint(i,jn,l+1)+zint(ie,jn,l+1)))
421 z1 = 0.25*(zmid(iw,j,l+1)+zmid(ie,j,l+1)+ &
422 zmid(i,jn,l+1)+zmid(i,js,l+1))
423 z3 = 0.25*(zmid(iw,j,l-1)+zmid(ie,j,l-1)+ &
424 zmid(i,jn,l-1)+zmid(i,js,l-1))
425 dz = 0.25*((zint(iw,j,l)+zint(ie,j,l)+ &
426 zint(i,js,l)+zint(i,jn,l))- &
427 (zint(iw,j,l+1)+zint(ie,j,l+1)+ &
428 zint(i,js,l+1)+zint(i,jn,l+1)))
432 du1 = uh(i,j,l+1)-uh(i,j,l)
433 du2 = uh(i,j,l)-uh(i,j,l-1)
434 dv1 = vh(i,j,l+1)-vh(i,j,l)
435 dv2 = vh(i,j,l)-vh(i,j,l-1)
436 IF( l >= lupp(i,j) .AND. l <= llow(i,j) )
THEN
437 IF( vh(i,j,l) <spval.and.uh(i,j,l) <spval.and. &
438 vh(i,j,l+1)<spval.and.uh(i,j,l+1)<spval.and. &
439 vh(i,j,l-1)<spval.and.uh(i,j,l-1)<spval.and. &
440 vst(i,j) <spval.and.ust(i,j) <spval) &
441 heli(i,j) = ((vh(i,j,l)-vst(i,j))* &
442 (dz2*(du1/dz1)+dz1*(du2/dz2)) &
443 - (uh(i,j,l)-ust(i,j))* &
444 (dz2*(dv1/dz1)+dz1*(dv2/dz2))) &
445 *dz/(dz1+dz2)+heli(i,j)
447 IF(lupp(i,j) == llow(i,j)) heli(i,j) = 0.
subroutine calhel3(LLOW, LUPP, UST, VST, HELI)
Subroutine that computes storm relative helicity.