NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
Kevin and IDV users List-As I told Kevin offline, I'm going to add some functionality to the layerDifference and layerAverage functions that will allow the user to specify the units of the levels. The default value is to use the units of the levels in the data. I will add this after the next 3.0 update to make sure it doesn't create any unforseen problems.
So, Kevin will be able to change his formula below to be:thrm(F1,user_Top_of_Layer[default=500],user_Bottom_of_Layer [default=1000],user_Level_Units[default=hPa])
However, one thing to note about this is that the IDV has to download all the levels and then select out the levels that the user requests. For large grids, this can take some time. The user can use the Levels tab in the Data Subset panel to select the range from the top level to bottom level to reduce the amount of data downloaded, but that still can be quite a lot of data.
Another way to do this is to create a formula that takes just the two levels. For thermal wind, the formula is:
THRM ( S ) = [ u (GEO(S)) (level1) - u (GEO(S)) (level2),
v (GEO(S)) (level1) - v (GEO(S)) (level2) ]
basically, subtract the components of the geostrophic wind of the base
level from the top level. This can be rewritten as:
sub(geo(zTop), geo(zBase))and the user would select two parameters - z at the top level and z at the base level. While it would be more selecting, it would reduce the amount of data read and speed up the processing time.
Just some food for thought. Don On 2/3/12 9:52 AM, Kevin R. Tyle wrote:
Hi all, I set up a formula in IDV to plot thermal wind vectors, a la: <derived description="thermal wind in a layer" formula="thrm(F1,user_Top_of_Layer[default=500],user_Bottom_of_Layer [default=1000])" id="thermalWind" isdefault="false" isenduser="true"> Whenever I do this on a GRIB2 format grid file, e.g. realtime NAM or GFS from the motherlode or atm-nsf opendap servers, I end up with all missing data when I load in the formula, specify my levels, and choose the 3D geopotential height from the GRIB2 file. If I do the same on the same model data that has been converted to GEMPAK, and also with the sample model data used in the IDV Workshop (which I believe is GRIB1), the thermal wind vectors plot fine. I suspect there is some issue specifically related to GRIB2 format here ... Thanks for looking into this! --Kevin _____________________________________________ Kevin Tyle, Systems Administrator Dept. of Atmospheric & Environmental Sciences University at Albany Earth Science 235, 1400 Washington Avenue Albany, NY 12222 Email: ktyle@xxxxxxxxxx Phone: 518-442-4578 _____________________________________________ _______________________________________________ idvusers mailing list idvusers@xxxxxxxxxxxxxxxx For list information, to unsubscribe, visit: http://www.unidata.ucar.edu/mailing_lists/
-- Don Murray NOAA/ESRL/PSD and CIRES 303-497-3596 http://www.esrl.noaa.gov/psd/people/don.murray/
idvusers archives: