[ncl-talk] Error with "time_axis_labels(cd_date, resT1, restick)"
Zilore Mumba
zmumba at gmail.com
Wed Jul 10 15:06:35 MDT 2019
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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190710/cf84c772/attachment.html>
More information about the ncl-talk
mailing list