|
UPP v11.0.0
|
Subroutine that computes the updraft helicity. More...
Go to the source code of this file.
Functions/Subroutines | |
| subroutine | calupdhel (updhel) |
Subroutine that computes the updraft helicity.
| [out] | UPDHEL | Updraft helicity (m^2/s^2). |
| Date | Programmer | Comments |
|---|---|---|
| 2007-10-22 | M Pyle | Initial |
| 2007-10-22 | M Pyle | based on SPC Algorithm courtesy of David Bright |
| 2011-01-11 | M Pyle | converted to F90 for unified post |
| 2011-04-05 | H Chuang | added B grid option |
| 2020-11-06 | J Meng | Use UPP_MATH Module |
| 2021-10-31 | J Meng | 2D DECOMPOSITION |
Definition in file CALUPDHEL.f.
| subroutine calupdhel | ( | real, dimension(ista_2l:iend integer :: l, j, i integer, dimension(jm) :: ihe,ihw ! integer dxval,dyval,cenlat,cenlon,truelat1,truelat2 ! integer latstart,lonstart,latlast,lonlast ! !*************************************************************************** ! start calupdhel here. ! ! write(6,*) 'min/max wh(:,:,20):: ', minval(wh(:,:,20)), & ! maxval(wh(:,:,20)) do l=1,lm call exch(uh(ista_2l,jsta_2l,l)) end do if (gridtype == 'b')then do l=1,lm call exch(vh(ista_2l,jsta_2l,l)) end do end if !$omp parallel do private(i,j) do j=jsta_2l,jend_2u do i=ista_2l,iend_2u updhel(i,j) = d00 enddo enddo do j=jsta_2l,jend_2u ihw(j) = -mod(j,2) ihe(j) = ihw(j)+1 enddo ! integrate (w * relative vorticity * dz) over the 2 km to ! 5 km agl depth. ! initial try without horizontal averaging !$omp parallel do private(i,j) do j=jsta_m,jend_m do i=ista_m,iend_m htsfc(i,j) = zint(i,j,nint(lmh(i,j))+1) | updhel | ) |
Definition at line 17 of file CALUPDHEL.f.