Dynamic.gs calculates the following variables from the geopotential height, temperature, and horizontal wind.
Variable Name Units
-Geostrophic Wind : ug,vg [m/s]
-Ageostrophic Wind : ua,va [m/s]
-Q-Vectors : Q1,Q1 [pa/m2/s]
-Temp Advection : tadv [K/s]
-Vort Advection : vadv [-]
-Frontogenesis : F [K/m/s]x10^9
-Fn Vector : fnx,fny [K/m/s]x10^9
-Deformation : def1,def2 [m]
The variables are then saved into the name, and you can plot them as normal variables. However, since there are a lot of calculations to make in this script, the variables are only calculated for the current dimensions. To calculate these values at all x,y,z,t points, simply set your dimensions to span the entire 4D domain.
Options: -help: pulls up help page
-var : allows the user to point to specified height, temperature, wind variables
Default variables are: hgtprs, tmpprs, ugrdprs, vgrdprs
Example Usage: dynamic -var hgtprs tmpprs ugrdprs vgrdprs
This example will use the hgtprs, tmpprs, ugrdprs, vgrdprs variables to calculate the dynamic variables.
Example Image: NARR Shaded: Fronogenesis, Contoured in Green: Fn vector divergence, Vectors: Fn Vectors for 06UTC Feb 15, 2003
|850-700mb average: Frontogensis, Fn Vectors, Fn Divergence, plotted for 06UTC Feb 15 2003|
I chose this example to compare to the example used in this presentation. Things seem to match up reasonably well. Though I do urge some caution here, I am not (nor have I ever been) a dynamics expert, and some of the calculations were taken from either old GrADS scripts, or approximations of textbook equations. Also, it is possible that the listed units are not 100% correct. These caveats aside, I have been using this script to successfully diagnose areas of frontogenesis, and ageostrophic circulations from model output for quite some time, so I think overall it serves it's purpose well. I post it now, with a newly added help page, and slightly cleaner code.
Download Dynamic.gs Here
Download Frontogenesis Example Here