UPP  11.0.0
 All Data Structures Files Functions Pages
CALHEL.f File Reference

Subroutine that computes storm relative helicity. More...

Go to the source code of this file.

Functions/Subroutines

subroutine calhel (DEPTH, UST, VST, HELI, USHR1, VSHR1, USHR6, VSHR6)
 This routine computes estimated storm motion and storm-relative environmental helicity. More...
 

Detailed Description

Subroutine that computes storm relative helicity.

Definition in file CALHEL.f.

Function/Subroutine Documentation

subroutine calhel ( real, dimension(2), intent(in)  DEPTH,
real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u), intent(out)  UST,
real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u), intent(out)  VST,
real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u,2), intent(out)  HELI,
real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u)  USHR1,
real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u)  VSHR1,
real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u)  USHR6,
real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u)  VSHR6 
)

This routine computes estimated storm motion and storm-relative environmental helicity.

(Davies-Jones et al 1990) the algorithm processd as follows.

The storm motion computation no longer employs the Davies and Johns (1993) method which defined storm motion as 30 degress to the right of the 0-6 km mean wind at 75% of the speed for mean speeds less than 15 m/s and 20 degrees to the right for speeds greater than 15 m/s. Instead, we now use the dynamic method (Bunkers et al. 1988) which has been found to do better in cases with 'non-classic' hodographs (such as Northwest-flow events) and do as well or better than the old method in more classic situations.

Parameters
[in]DEPTHDepth in meters over whcih helicity should be computed; allows one to distinguish 0-3 km and 0-1 km values.
[out]USTEstimated U Component (m/s) Of Storm motion.
[out]VSTEstimated V Component (m/s) Of Storm motion.
[out]HELIStorm-relative heliciry (m**2/s**2).
[out]USHR1U Component (m/s) Of 0-1 km shear.
[out]VSHR1V Component (m/s) Of 0-1 km shear.
[out]USHR6U Component (m/s) Of 0-0.5 to 5.5-6.0 km shear.
[out]VSHR6V Component (m/s) Of 0-0.5 to 5.5-6.0 km shear.

Program history log:

Date Programmer Comments
1994-08-22 Michael Baldwin Initial
1997-03-27 Michael Baldwin Speed up code
1998-06-15 T Black Conversion from 1-D to 2-D
2000-01-04 Jim Tuccillo MPI Version
2000-01-10 G Manikin Changed to Bunkers method
2002-05-22 G Manikin Now allow choice of computing helicity over two different (0-1 and 0-3 km) depths
2003-03-25 G Manikin Modified code to compute mean winds using arithmetic averages instead of mass weighting; differences are minor but want to be consistent with the Bunkers method
2004-04-16 M Pyle Minimal modifications but put into NMM WRFPOST code
2005-02-25 H Chuang Add computation for ARW A grid
2005-07-07 Binbin Zhou Add RSM for A grid
2019-10-30 Bo Cui Remove "goto" statement
2021-09-02 Bo Cui Decompose UPP in X direction
2022-10-07 Tracy Hertneky Add left mover for storm motion in SH
Author
Michael Baldwin W/NP2
Date
1994-08-22

Definition at line 43 of file CALHEL.f.

Referenced by miscln().