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.
Jon, > Just thought of another question – is it possible to find > out what the min-max range of “natural” coordinates is, or would I have to > perform the transformation (theoretically for every timestep) and search the > returned 3D array? For ocean models, to compute the min-max Z range for the dataset, you would need to read the max value over all time of the ocean surface variable, and then find the deepest depth of the depth (bathymetry) variable. So you need to read a 3D surface field (Y,X,T) and the depth field (X,Y), but you wouldn't need to do any of the vertical tranformations. > I suppose that I could do the transformation for one > timestep only and assume the range won’t be very different for other > timesteps since the free surface (or surface pressure) isn’t likely to vary > too much? Yes, you could probably get away with that. > > > > Rich – yes, I probably could use the IDV routines although that means > another dependency for the sake of one routine, plus I need to look at how > to instantiate the IDV objects correctly in all cases. Just requires a bit > more thought. Okay, well, at least use the NJ routines to return the vertical coordinate info (let it do the transformations for you). Nearest neighbor is way too crude for vertical interpolation. Linear will be fine for most visualization, and easy to implement. For each x,y, just find the levels below and above the requested z value and linearly interpolate between them. -- Dr. Richard P. Signell (508) 457-2229 USGS, 384 Woods Hole Rd. Woods Hole, MA 02543-1598
netcdf-java
archives: