[ncl-talk] reading station data

[Dennis Shea]

Please look carefully at the code I sent.
The directory and the file name are separate.

ncl hazra_time.ncl > out.hazra_test_daily

Note the the directory name and the filename are separate.
If your file name has changed structure *you* must accomodate.



On Thu, Dec 20, 2018 at 12:40 PM Debasish Hazra <debasish.hazra5 at gmail.com>
wrote:

> Thanks Dennis. But I am getting warning messages when I run your script :
> warning:toint: A bad value was passed (string); input strings must contain
> numeric digits, replacing with missing value
> warning:toint: A bad value was passed (string); input strings must contain
> numeric digits, replacing with missing value
> warning:tofloat: A bad value was passed to (string) tofloat, input strings
> must contain numeric digits, replacing with missing value
> warning:tofloat: A bad value was passed to (string) tofloat, input strings
> must contain numeric digits, replacing with missing value
>
> and finally Error :
> fatal:The result of the conditional expression yields a missing value. NCL
> can not determine branch, see ismissing function
> fatal:["Execute.c":8575]:Execute: Error occurred at or near line 4740 in
> file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
> fatal:["Execute.c":8575]:Execute: Error occurred at or near line 30 in
> file test_daily2.ncl which is :
>
> yyyymmdd = yyyyddd_to_yyyymmdd(yyyyddd)
>
> Any help.
> Thanks
> Debasish
>
> On Thu, Dec 20, 2018 at 10:48 AM Dennis Shea <shea at ucar.edu> wrote:
>
>> Attached is a a method to get appropriate times.
>>
>> You can write code that finds indices that occur within the desired time
>> range.
>>
>> Say:good luck
>>
>>    1 hour: + /- 0.5 hours
>>    3 hour: +/- 1.5 hours
>> etc
>>
>>     ii = ind(...)
>>     if (.not.ismissing(ii(0))) then
>>         AOT??? = avg(aot(ii))
>>     end if
>>
>> good luck
>>
>> On Wed, Dec 19, 2018 at 3:22 PM Debasish Hazra <debasish.hazra5 at gmail.com>
>> wrote:
>>
>>> Hi,
>>> I am trying to use NCL to read station data observation which spans over
>>> a month and at irregular intervals. Idea is to produce observation averages
>>> around fixed time intervals in a day (like 3,6,9,12,15,18 and 21z) using 1
>>> hour windows. However, I am getting fill values in result, although valid
>>> observations are available at those times. Both input and ncl code is
>>> attached. Any help is appreciated.
>>> Thank you
>>> Debasish.
>>> _______________________________________________
>>> 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/20181220/c6254d58/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: out.hazra_test_daily
Type: application/octet-stream
Size: 116919 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181220/c6254d58/attachment.obj>