coordinate indexing problem

I am having trouble extracting a spatial region from a
3 dimensional variable using coordinate indexes.

My data looks like:

Variable: pct_pft
Type: float
Total Size: 17625600 bytes
             4406400 values
Number of Dimensions: 3
Dimensions and sizes: [pft | 17] x [lat | 360] x [lon | 720]
Coordinates:
             pft: [1..17]
             lat: [-89.75..89.75]
             lon: [0.25..359.75]
Number Of Attributes: 3
   lonFlip : longitude coordinate v has been reordered via lonFlip
   units : unitless
   long_name : percent pft

My region and the line that fails:

slat = 60.
elat = 85.
slon = 290.
elon = 345.

doh2 = pct_pft(0,{slat:elat},{slon:elon})

fatal:NclOneDValGetRangeIndex: Non-monotonic coordinate value being
used, can't complete coordinate subscript
fatal:Could not obtain coordinate indexes, unable to perform subscript
fatal:Execute: Error occurred at or near line 169 in file
mk_mksrf.warmArctic.ncl

My error message is very similar to an archive question, (see below).
Unfortunately, it looks like Dennis was unable to replicate the problem,
so there is no solution in the archive:

http://www.ncl.ucar.edu/Support/talk_archives/2008/0799.html

thanks in advance for any help -

nanr
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Nov 25 12:49:22 2009