Re: Coordinate subscripting problems

Hi Seth,

Thanks for providing such a good analysis of your problem.

There was a bug fixed in version a034 that dealt with single value
coordinate subscripting. see:

   http://www.ncl.ucar.edu/new_release.shtml#ImportantBugFixes

If you tell me what system you have, I can give you a new ncl
executable to try to see if this makes a difference. If it does, then
I can provide you with the full software, or you can download it
yourself from www.earthsystemgrid.org.

--Mary

On Thu, 12 Apr 2007 mcginnis_at_ucar.edu wrote:

> I am encountering some very strange behavior using variable subscripting on
> some NCEP-II data. If anyone can give me some insight into what's going on,
> I'd appreciate it.
>
> I've attached a script that demonstrates the problem I'm talking and the
> output from that script, but here's a short summary:
>
> I'm running NCL 4.2.0.a033 on a linux machine. I have NetCDF file with data
> from the NCEP-II reanalysis (air temperature, in my example, but it happens
> with other files as well).
>
> I'm reading in data with a line like the following:
>
> temp = short2flt(fin->air(time|:, level|0, {lat|LAT}, {lon|LON}))
>
> I'm selecting a single point, so NCL selects the closest value. If I slowly
> change LAT from a value near one point in the latitude coordinate variable
> array (a grid center) to another, the data values change when it's halfway
> between, just as I'd expect.
>
> Here's the weird thing: the data values *also* change when I'm changing from
> one side of a grid center to the other, and, moreover, I get a third value if
> my latitude is exactly equal to the grid center latitude. This doesn't
> happen with longitude.
>
> Any ideas what's going on?
>
> Thanks,
>
> --Seth McGinnis
>
>
> ----------------------------------------------------------------
>
>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Apr 12 2007 - 16:30:49 MDT