[ncl-talk] Converting time to datestring

[Barry Lynn]

Hi:

You have to define year, month, etc from your data.

Can you do that?


On Thu, 4 Feb 2021 at 20:39 Zilore Mumba <zmumba at gmail.com> wrote:

> Hello Lynn,
> Sorry I have had no success with the code you sent me. On the line
> new_time_array = cd_inv_calendar( year(it,1,1),  month(it,1,1),
> day(it,1,1), hour(it,1,1), minute(it,1,1), 0, new_time_units,0)
> I get the error: fatal:Undefined identifier: (year) is undefined, can't
> continue
>
> On Thu, Feb 4, 2021 at 7:55 AM Barry Lynn <barry.h.lynn at gmail.com> wrote:
>
>> Hi:
>>
>> Try adjusting the third line here:
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_string.ncl"
>>
>> begin code...
>>
>> new_time_units = "hours since 1800-01-01 00:00"
>> new_time_units = new_time_units
>> new_time_array = cd_inv_calendar( year(it,1,1),  month(it,1,1),
>> day(it,1,1), hour(it,1,1), minute(it,1,1), 0, new_time_units,0)
>> print("new_time_array = " + new_time_array)
>> printVarSummary(new_time_array)
>> format = "" ;; defaults to "%H%M IST %d %c %Y"
>>  format = "%H%M IST %d %c %Y"
>>  new_time_array = new_time_array + 2.01
>> ;format = "%H%M IDT %d %c %Y"
>> ;new_time_array = new_time_array + 3.01
>> time_var=cd_string(new_time_array,format)
>> print("time_var = " + time_var)
>>
>> On Thu, Feb 4, 2021 at 7:17 AM Zilore Mumba via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> Hello,
>>> I hope I can get quick help to set me off as it is taking me days to
>>> make sense of what I have to do.
>>> I have historical (1971-200) and model simulation (2021-2050) data. I
>>> have adapted an ncl script from the DKRZ website for plotting a time series.
>>>
>>> With print(times) I get the time range 7701.5 to 18658.5.
>>> When I plot the time series at a point I get the x-axis range 6000 60
>>> 20000
>>> ncdump -h on the file indicates that time:units = "days since 1949-12-01
>>> 00:00:00"
>>> my data has format yyyymmdd (yearmothday). I want the x-axis to show
>>> dates.
>>>
>>> I have been looking in the ncl documentation how to do this. This may
>>> take me days or weeks.
>>> I would appreciate any help in this.
>>> The script I have is below
>>>
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
>>> begin
>>>
>>>    diri = "./"                             ; directory where to find the
>>> data
>>>    fili =
>>> "pr_AFR-44_CNRM-CERFACS-CNRM-CM5_historical_r1i1p1_SMHI-RCA4_v1_day_19710101_20001231Reg.nc"
>>>             ; file name
>>>
>>>    f   = addfile(diri+fili,"r")             ; open file
>>>
>>>    t1 = f->pr(:,{-15},{30})                 ; read precip time series at
>>> 15S/30E
>>>
>>>    times  = f->time                         ; get times in file
>>>    cd_date = cd_calendar(times, -2)
>>>
>>> ; open a workstation
>>> ;----------------------------------------------------------------
>>>    wks_type = "X11"                         ; graphics output type
>>>    wks_type at wkWidth   =  1024               ; graphics output width in
>>> pixels
>>>    wks_type at wkHeight  =  1048               ; graphics output height in
>>> pixels
>>>
>>>    wks = gsn_open_wks(wks_type, "plot_xy_multiple_lines_axis")
>>>
>>> ; resources for variable t1, y-axis
>>> ;----------------------------------------------------------
>>>    resT1                      =  True
>>>    resT1 at ttmValues            =  cd_date
>>>    resT1 at xyLineThicknessF     =  3.0                        ; set line
>>> thickness
>>>    resT1 at xyLineColor          =  "red"                      ; set line
>>> color
>>>    resT1 at xyDashPatterns       =  0                          ; solid line
>>>    resT1 at tmXBLabelFontHeightF = 0.010
>>>    resT1 at tiYAxisString        =  "Daily " + " " + t1 at standard_name +"
>>>  "+ "at 15S/30E"              ; y-axis string
>>>
>>> ; create the plot
>>> ;-----------------------------------------------
>>>    plot = gsn_csm_xy(wks, times&time, t1, resT1)
>>>
>>> end
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at mailman.ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>>
>> --
>> Guest Editor of Special Issue “Extreme Hydroclimatic Events and
>> Prediction”
>> LINK:
>>
>> https://www.mdpi.com/journal/water/special_issues/Extreme_Hydroclimatic_Events
>>
>> 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
>>
>> --
Guest Editor of Special Issue “Extreme Hydroclimatic Events and Prediction”
LINK:
https://www.mdpi.com/journal/water/special_issues/Extreme_Hydroclimatic_Events

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: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210204/dc38fa21/attachment.html>