CDMS2 axis.getBounds() returns wrong leading and trailing bounds for monotonic decreasing axis

[JimBiardCics]

I have a data set with a polar stereographic grid. The y axis has 448 elements and is monotonic decreasing from 5837500.0 m to -5337500.0 m. When the getBounds method is called it produces a leading edge bound of (90.0, 90.0) and a trailing edge bound of (-90.0, -90.0). One of the files from the data set is attached.

sample.zip

[doutriaux1]

@JimBiardCics thanks for reporting. For future ref please note that we splits UV-CDAT into various repos, for example cdms2 bug should be reported to the http://github.com/UV-CDAT/cdms repo

[JimBiardCics]

I did some further research and found that the first and last bounds are wrong even if the axis is monotonic increasing. The values on the first and last bounds are the same (with appropriate sign reversal).

[dnadeau4]

@JimBiardCics can you try again. I get Latitude from 31 to 34.

>>> s.getLatitude()[:]
latitude
masked_array(data =
 [[ 31.10267175  31.19941382  31.29580182 ...,  31.67773544  31.58280258
   31.48750008]
 [ 31.25010968  31.3472803   31.44409653 ...,  31.82773957  31.73237994
   31.63665027]
 [ 31.39756424  31.49516539  31.59241187 ...,  31.97777244  31.88198399
   31.78582511]
 ..., 
 [ 34.34589323  34.45253044  34.5588102  ...,  34.98026512  34.87546
   34.77027946]
 [ 34.19869696  34.30486281  34.41067119 ...,  34.83024101  34.72590705
   34.62119777]
 [ 34.05145898  34.15715555  34.26249465 ...,  34.68018796  34.57632298
   34.4720828 ]],
             mask =
 False,
       fill_value = 1e+20)

[dnadeau4]

I also fixed getBounds for Latitude if not in degrees_north units.

[JimBiardCics]

@dnadeau4 I still see the problem with cdms2 version 2.8.1 (that's what's available on the uvcdat conda channel). Why don't we close this and move to the migrated issue that @chaosphere2112 made.

[jasonb5]

Closing this issue, another issue was created in the cdms repo which has been closed.