[ncl-talk] Creating Date String

Tue Nov 1 22:11:37 MDT 2016

It looks like you're still passing the old time array version.

time_var=cd_string(time_array,format)

printvarsummary doesn't print the actual variable data, the number you see
there is how the fillvalue is represented. That all looks good. Edit that
above line and you should be golden.

Alan.

On 1 Nov 2016 11:32 p.m., "Barry Lynn" <barry.h.lynn at gmail.com> wrote:

> Dear Alan:
>
> Thank you for helping me push this along.
>
> My variables are already integers
>
> new_time_units = "hours since 1800-01-01 00:00"
> new_time_array = cd_inv_calendar( year(1,1,1),  month(1,1,1), day(1,1,1),
> hour(1,1,1), minute(1,1,1), 0, new_time_units,0)
> print("new_time_array = " + new_time_array)
> format = "" ;; defaults to "%H%M UTC %d %c %Y"
> time_var=cd_string(time_array,format)
> print("time_var = " + time_var)
>
> However, I still get the unit error problem.
>
> Note below, please: new time_array from a print("new_time_array" +
> new_time_array") produces
>
> (0) new_time_array = 1900752.166666667
>
> but (below), the printVarSummary produces a fill value equivalent to a NaN.
>
> Here are the print outs.
>
> (0) year(1,1,1) = 2016
> (0) month(1,1,1) = 11
> (0) day(1,1,1) = 2
> (0) hour(1,1,1) = 0
> (0) new_time_array = 1900752.166666667
>
> Variable: new_time_array
> Type: double
> Total Size: 8 bytes
>             1 values
> Number of Dimensions: 1
> Dimensions and sizes: [1]
> Coordinates:
> Number Of Attributes: 3
>   _FillValue : 9.969209968386869e+36
>   calendar : standard
>   units : hours since 1800-01-01 00:00
> (0) cd_string: The variable containing time values must have a 'units'
> attribute.
> fatal:Illegal right-hand side type for assignment
> fatal:["Execute.c":8575]:Execute: Error occurred at or near line 130 in
> file ./obs_light_12km.ncl
>
>
> On Tue, Nov 1, 2016 at 10:53 PM, Alan Brammer <abrammer at albany.edu> wrote:
>
>>
>> The time variable passed to cd_string() should be in the format of “hours
>> since 1800-01-01 00:00:00” or similar.   The format passed to the argument
>> is for the output string not the input variable.
>>
>>
>> I don’t really understand how you got to the first line of your sample
>> code but trying to build from that is some code below:
>>
>>  ; time_array(0)=hour(1,1,1)
>>  ; time_array(1)=minute(1,1,1)
>>  ; time_array(2)=day(1,1,1)
>>  ; time_array(3)=month(1,1,1)
>>  ; time_array(4)=year(1,1,1)
>>  ;;   time_array(5)=units    ;; are the above strings or numbers?  If
>> they’re numbers this line won’t work.
>>
>> ;; first make a ut_coord complaint time variable from the arrays of time
>> stuff you already have.
>>
>> *time_array = cd_inv_calendar( toint(year(1,1,1)), toint (month(1,1,1)),
>> toint( day(1,1,1)), toint (hour(1,1,1)), toint (minute(1,1,1)), 0, “hours
>> since 1800-01-01 00:00:00”,0)*
>> *printVarSummary(time_array)*
>> ;;;    Should be a single number of double type with a units attribute.
>>
>> *format = “" ;; defaults to "%H%M UTC %d %c %Y"*
>> *time_var=cd_string(time_array,format)*
>>
>>
>>
>> p.s. beware of weird quotation marks if copy/pasting from email.
>>
>> Good luck,
>> Alan
>>
>>
>> ##############################
>> Alan Brammer,
>> Post-Doc Researcher
>>
>> Department of Atmospheric and Environmental Sciences,
>> University at Albany, State University of New York, Albany, NY, 12222
>> abrammer at albany.edu
>> ##############################
>>
>> On 1 Nov 2016, at 06:33, Barry Lynn <barry.h.lynn at gmail.com> wrote:
>>
>> Good Morning:
>>
>> I went to this page
>>
>> https://www.ncl.ucar.edu/Document/Functions/User_contributed
>> /cd_string.shtml
>>
>> and tried to do this:
>>
>> For all the following examples, assume that *time* is an array of time
>> values.
>>
>> *Example 1*
>>
>> Using the default format string of "" will result in "HHMM UTC D Mon
>> YYYY" (e.g., 1800 UTC 4 Jul 1776):
>>
>>   format = ""
>>   stime  = cd_string(time,format)
>>
>> The goal is to plot the string "time_var" as shown below, as a single
>> variable.
>>
>> But I get an error; please see below.
>>
>> (Note: I thought (as a try) to stick this information in the "unitlabel"
>> parameter below, but I would really prefer to plot the date on the right
>> hand side at
>> the top of the map).
>>
>> Here's my code.
>>
>> ;HHMM UTC D Mon YYYY
>>   time_array(0)=hour(1,1,1)
>>   time_array(1)=minute(1,1,1)
>>   time_array(2)=day(1,1,1)
>>   time_array(3)=month(1,1,1)
>>   time_array(4)=year(1,1,1)
>>   time_array(5)=units
>> ; format = "HHMM UTC D Mon YYYY"
>>   format = ""
>>   *time_var=cd_string(time_array,format)*
>>   print("time_var = " + time_var)
>>   opts_ic_cg at FieldTitle           = "Lightning Rate"
>> ; opts_ic_cg at UnitLabel            = "/10-min"
>>  * opts_ic_cg at UnitLabel            = time_var*
>>
>> (0) Working on time: 2016-10-30_00:00:00
>> (0) cd_string: The variable containing time values must have a 'units'
>> attribute.
>> fatal:Illegal right-hand side type for assignment
>>
>> After searching and finding:
>>
>> http://www.ncl.ucar.edu/Document/Functions/User_contributed/
>> time_axis_labels.shtml
>>
>> *time*
>>
>> An array containing the values of time. This variable should have an
>> attribute named "units". In other words, this variable holds the values of
>> the coordinate axis 'time' along with the units.
>>
>>
>> I am not sure what this error means or how to address the error.
>>
>> Thank you.
>>
>>
>> --
>> Barry H. Lynn, Ph.D
>> Senior 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
>> US 914 432 3108
>> _______________________________________________
>> 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 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
> US 914 432 3108
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161102/b43c3bc1/attachment.html 