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.
On Tue, 2006-03-14 at 09:21, mpborkow@xxxxxxxx wrote:
I need to be able to compute the gradient between two points for which I will only know the lat/lon values. I'm looking for the easiest way to do this. Is there a way to spit out the value for a parameter by providing asingle lat/lon?
Matt, Yes, you could use the GDPOINT program to output a value at a specific GPOINT. You can also use the ddx(), ddy(), and/or grad() grid diagnostics at an intervening lat/lon.
Or alternatively, is there a simple way of obtaining the actual gridpoints for a location if I only know the latitude/longitude coordinates so that I can easily reference the grid points in a Fortran program?
The GEMPAK library provides a function GR_PLOC that accepts the character string for GPOINT (which may be grid point locations, lat/lon, station id etc), and provides back the projected grid point location as well as lat/lon values (you could add a write(*,*) statement tothe $NAWIPS/comet/gdpoint/gdpoint.f routine to output the rgy and rgy location following the GR_PLOC call if you wanted to use that program as a starting point). Steve Chiswell Unidata User Support
Matt Borkowski Graduate Student North Carolina State University Department of Marine, Earth, and Atmospheric Science
gembud
archives: