[ncl-talk] Coordinate Subscripting

[Adam Phillips]

Hi Jalal,
I just wrote this short test to see if I could reproduce your error, as I
suspect you have repeating times in your time coordinate variable:
ncl 0> a = (/3.4,2.1,1.,3./)
ncl 1> a!0 = "time"
ncl 2> a&time = (/19980101,19980102,19980102,19980103/)
ncl 3> print(a({19980103}))
fatal:NclOneDValGetClosestIndex: Non-monotonic coordinate value being used,
can't complete coordinate subscript
fatal:Could not obtain coordinate indexes, unable to perform subscript
fatal:["Execute.c":8635]:Execute: Error occurred at or near line 3

That says to me that you need to look very closely at your time coordinate
variable as I recommended before:
prc_reorder&time = cd_calendar(time, -2)
print(prc_reorder&time)
I suspect you will find at least one instance of a repeating time. If you
do have multiple times on the same day, you can try this to create times in
the format of YYYYMMDDHH:
prc_reorder&time = cd_calendar(time, -3)
plot = gsn_csm_contour_map(wks,prc_reorder({1998042200},:,:),res)
Adam





On Fri, Jan 18, 2019 at 6:16 PM Md. Jalal Uddin <dmjalal90 at gmail.com> wrote:

> Hi Adam,
>
> Thanks for your kind help. I tried to follow your instructions; however, I
> am getting the following errors.
>
> fatal:NclOneDValGetClosestIndex: Non-monotonic coordinate value being
> used, can't complete coordinate subscript
> fatal:Could not obtain coordinate indexes, unable to perform subscript
> fatal:["Execute.c":8637]:Execute: Error occurred at or near line 54 in
> file subs_time.ncl
>
> 54 line:
>  plot = gsn_csm_contour_map(wks,prc_reorder({19980501},:,:),res)
>
> When I used yours suggestion like plot =
> gsn_csm_contour_map(wks,prc_reorder({19980422},:,:),res), I got the same
> problem.
>
> Any suggestions, please?
>
> Jalal
>
> On Sat, Jan 19, 2019 at 3:19 AM Adam Phillips <asphilli at ucar.edu> wrote:
>
>> Hi Jalal,
>> Within your script you are calling cd_calendar, and then subsequently
>> calling day_of_year.
>> time at units = "days since 1998-1-1 12:00:00"
>> utc_date = cd_calendar(time, -5)
>> dyear = day_of_year(utc_date(:,0),utc_date(:,1),utc_date(:,2))
>>
>> day_of_year will return numbers ranging from 1-365, and thus the values
>> can repeat if you have multiple years. This is not the type of array that
>> should be used for coordinate subscripting. I would recommend calling
>> cd_calendar, but retrieving the values in YYYYMMDD format, then assign this
>> array as a coordinate variable to your prc_reorder array:
>> delete(prc_reorder&time)   ; may not be needed
>> prc_reorder&time = cd_calendar(time, -2)
>>
>> then you should be able to use the following coordinate subscripting to
>> plot:
>> plot = gsn_csm_contour_map(wks,prc_reorder({19980422},:,:),res)
>>
>> Hope that helps. If you have any further questions let ncl-talk know.
>> Adam
>>
>>
>>
>> On Wed, Jan 16, 2019 at 11:11 PM Md. Jalal Uddin <dmjalal90 at gmail.com>
>> wrote:
>>
>>> Hi all,
>>>
>>> I used cd_calendar to convert utc_date and day_of_year to find the day
>>> of year. Now, I want to subscript the time coordinate and plot. Could
>>> anyone help me? Please find the script in the attached.
>>>
>>> Variable: time
>>> Type: double
>>> Total Size: 22776 bytes
>>>             2847 values
>>> Number of Dimensions: 1
>>> Dimensions and sizes: [time | 2847]
>>> Coordinates:
>>>             time: [ 120..6939]
>>> Number Of Attributes: 5
>>>   standard_name : time
>>>   long_name : Time
>>>   units : days since 1998-1-1 12:00:00
>>>   calendar : standard
>>>   axis : T
>>>
>>> Best regards,
>>> Jalal
>>>
>>> --
>>> *Md. Jalal Uddin*
>>> MSc in Applied Meteorology (English Language)
>>> Nanjing University of Information, Science and Technology, China
>>> B.Sc. in Disaster Management (Hons.)
>>> Patuakhali Science and Technology University, Bangladesh.
>>> Cell: +8613260859092, +8801792052662
>>> Web: www.dmjalal90.weebly.com
>>> Facebook: jalal.hossen.39
>>> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
>>> Twitter: dmjalal90
>>> Skype: dmjalal90
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>
>>
>> --
>> Adam Phillips
>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
>> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>>
>> <http://www.cgd.ucar.edu/staff/asphilli>
>>
>
>
> --
> *Md. Jalal Uddin*
> MSc in Applied Meteorology (English Language)
> Nanjing University of Information, Science and Technology, China
> B.Sc. in Disaster Management (Hons.)
> Patuakhali Science and Technology University, Bangladesh.
> Cell: +8613260859092, +8801792052662
> Web: www.dmjalal90.weebly.com
> Facebook: jalal.hossen.39
> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
> Twitter: dmjalal90
> Skype: dmjalal90
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190122/21ca0bec/attachment.html>