This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.
Phillip, The BVSQ() function in the grid functions accepts an input value, which typically is THTA (though nothing would prevent you from using THTE) which you chose to pass to the function. The $GEMPAK/source/gemlib/df/dfbvsq.c routine then uses your input grid for specified thta grid as described in the routine: * BVSQ ( THTA ) = [ GRAVTY * LDF (THTA) ] / [ LAV (THTA) * DZ ] * * * * DZ = change in height across the layer * * = -( RDGAS / GRAVTY ) * LAV (THTA) * * * ( LAV (PRES) / 1000 ) ** KAPPA * LDF (PRES) / LAV (PRES) * * in THTA coordinates * Steve Chiswell Unidata User Support > Hello, > > My question deals with the moist Brunt-Vaisala frequency. Does this > algorithm just replace theta with theta-e or does the algorithm use the > formula used in Durran and Klemp 1982? I probably need to use the > latter, but it would be a little messy if I tried it on gfunc or > gdpfun. I'm currently using version 5.9.4 for reference. Thank you for > you help and time. > > Phillip Marzette > > > Ticket Details =================== Ticket ID: ZAN-221187 Department: Support GEMPAK Priority: Normal Status: Closed