[ncl-talk] Error with "time_axis_labels(cd_date, resT1, restick)"
Rick Brownrigg
brownrig at ucar.edu
Wed Jul 10 17:11:42 MDT 2019
Hi all,
I'm confused by this thread -- where does cd_string enter? I don't see it
referenced in the script.
If the error messages are saying time_axis_labels is undefined, then as
Dave has pointed out, the docs say to load:
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/time_axis_labels.ncl"
See:
http://ncl.ucar.edu/Document/Functions/User_contributed/time_axis_labels.shtml
Rick
On Wed, Jul 10, 2019 at 3:07 PM Zilore Mumba via ncl-talk <ncl-talk at ucar.edu>
wrote:
> Thanks Barry, I did not realise that the contribute script I used was from
> WRF.
> However even with the right script I still get the same error. Perhaps you
> could help me, look at the script and see where am doing it wrong,
> especially where I use cd_calendar and date&time.
> The script is below. I know it is cumbersome because you do not have my
> data. But I will appreciate any help
> ;-----------------------------------------------------------------
> ; Adapted from DKRZ - NCL Example: NCL_xy_multiple_lines_axis.ncl
> ;-----------------------------------------------------------------
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_string.ncl"
>
> begin
>
> filein =
> "/home/zmumba/DATA/tasmax_AFR-44_CNRM-CERFACS-CNRM-CM5_historical_r1i1p1_CLMcom-CCLM4-8-17_v1_day_19960101-20001231.nc"
> ; file name
>
> f = addfile(filein,"r") ; open file
>
> tmax = f->tasmax(:,{-15},{25}) ; read tmax, all timesteps,
> 15S/25E
>
> t = f->time ; for use in x-axis label
> cd_date = cd_calendar(t, -2)
>
> ; open a workstation
> wks_type = "X11" ; graphics output type
> wks_type at wkWidth = 1024 ; graphics output width in
> pixels
> wks_type at wkHeight = 1024 ; graphics output height in
> pixels
>
> wks = gsn_open_wks(wks_type, "plot_xy_tie_series")
>
> ; resources for variable tmax y-axis
> res = True
> res at xyLineThicknessF = 3.0 ; set line
> thickness
> res at xyLineColor = "red" ; set line color
> res at xyDashPatterns = 0 ; solid line
> res at tmXBLabelFontHeightF = 0.010
> res at tiYAxisString = "Daily Max. 2m" + " " + tmax at standard_name
> +" "+ "at 15S/25E" ; axis
> string
>
> restick = True ; resource object
> restick at ttmFormat = "%d %c"
> restick at ttmAxis = "YL"
> restick at ttmMajorStride = 20
> restick at ttmValues = cd_date
>
> time_axis_labels(date&time,res,restick)
>
> ; create the plot
> plot = gsn_csm_xy(wks, t, tmax, res)
>
> end
>
> On Tue, Jul 9, 2019 at 7:33 PM Barry Lynn <barry.h.lynn at gmail.com> wrote:
>
>> Try this one.
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_string.ncl"
>>
>> On Tue, Jul 9, 2019 at 7:18 PM Dave Allured - NOAA Affiliate via ncl-talk
>> <ncl-talk at ucar.edu> wrote:
>>
>>> That is not the correct load command for this function.
>>> Time_axis_labels is loaded separately. Please go to the documentation page
>>> for time_axis_labels, find the correct load command at the top of that
>>> page, and add it to your script.
>>>
>>>
>>> On Tue, Jul 9, 2019 at 12:27 AM Zilore Mumba <zmumba at gmail.com> wrote:
>>>
>>>> Thanks Dave for the hint. i did load the contributed package. I still
>>>> get the same error. The package I loaded is
>>>> "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
>>>> When I check this script I do not see anywhere where there is a
>>>> function like "time_axis_labels".
>>>> I will continue to probe and see where am going wrong
>>>> Incidentally, example6 on
>>>> https://www.ncl.ucar.edu/Document/Functions/User_contributed/time_axis_labels.shtml,
>>>> gives
>>>>
>>>> *time_axis_labels*( data&time, res, restick )
>>>>
>>>> I assumed there was a mistake on data&time and corrected it to
>>>> *time_axis_labels*( date&time, res, restick )
>>>>
>>>> Thanks
>>>>
>>>> Zilore
>>>>
>>>>
>>>> On Mon, Jul 8, 2019 at 10:00 PM Dave Allured - NOAA Affiliate <
>>>> dave.allured at noaa.gov> wrote:
>>>>
>>>>> This is one of the contributed functions that is not loaded
>>>>> automatically. You are missing a "load" command. NCL gave you a clue:
>>>>> "possibly an undefined procedure". Please look at the top of the function
>>>>> documentation page.
>>>>>
>>>>>
>>>>> On Mon, Jul 8, 2019 at 9:14 AM Zilore Mumba via ncl-talk <
>>>>> ncl-talk at ucar.edu> wrote:
>>>>>
>>>>>> I am trying to plot a time series. I have converted the dates in the
>>>>>> file to the format YYYYMMDD with the function
>>>>>> t = f->time ; to get time in file
>>>>>> cd_date = cd_calendar(t, -2) ; print(cd_date) prints correctly,
>>>>>> eg 19960101,...
>>>>>> Then to prepare the format for the x-axis label, I have:
>>>>>>
>>>>>> resT1 = true
>>>>>>
>>>>>> some resources
>>>>>> ..............
>>>>>>
>>>>>> restick = True
>>>>>> restick at ttmFormat = "%d %c"
>>>>>> restick at ttmAxis = "YL"
>>>>>> restick at ttmMajorStride = 20
>>>>>> *time_axis_labels*( data&time, res, restick )
>>>>>>
>>>>>> I am getting the error:
>>>>>> fatal:syntax error: line 57 in file two-line-plot.ncl before or near \n
>>>>>> time_axis_labels( data&time, resT1, restick )
>>>>>> -------------------------------------------------^
>>>>>>
>>>>>> fatal:syntax error: possibly an undefined procedure
>>>>>> fatal:syntax error: line 72 in file two-line-plot.ncl before or near \n
>>>>>> time_axis_labels(cd_date,resT1,restick)
>>>>>> ------------------------------------------^
>>>>>>
>>>>>> fatal:syntax error: possibly an undefined procedure
>>>>>> fatal:Syntax Error in block, block not executed
>>>>>> fatal:error at line 76 in file two-line-plot.ncl
>>>>>>
>>>>>> Is the
>>>>>>
>>>>>> *time_axis_labels function still valid, or is there an alternative function?*
>>>>>>
>>>>>>
>>>>>> *Help will be appreciated.*
>>>>>>
>>>>>> *Zilore*
>>>>>>
>>>>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>>
>> --
>> Barry H. Lynn, Ph.D
>> Senior Associate Scientist, Lecturer,
>> The Institute of the Earth Science,
>> The Hebrew University of Jerusalem,
>> Givat Ram, Jerusalem 91904, Israel
>> Tel: 972 547 231 170
>> Fax: (972)-25662581
>>
>> C.E.O, Weather It Is, LTD
>> Weather and Climate Focus
>> http://weather-it-is.com
>> Jerusalem, Israel
>> Local: 02 930 9525
>> Cell: 054 7 231 170
>> Int-IS: x972 2 930 9525
>>
>> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190710/3d973524/attachment.html>
More information about the ncl-talk
mailing list