[ncl-talk] query
Dennis Shea
shea at ucar.edu
Fri Oct 6 10:11:21 MDT 2017
recline automatically strips out all _FillValue.
Sent from my iPhone
> On Oct 6, 2017, at 1:50 AM, Barry Lynn <barry.h.lynn at gmail.com> wrote:
>
> Hello:
>
> I'd like to use regline, and I don't always know how many points might be in the data. For example, if I want to look at 10 points, but one of them is missing, then when I calculate the slope of the line it will include a missing value.
>
> I am wondering how I might skip a missing value, and only calculate the line for present 9 points.
>
> Thanks,
>
> Barry
>
>> On Fri, Oct 6, 2017 at 7:08 AM, saurabh rathore <rohitsrb2020 at gmail.com> wrote:
>> Thank you Kevin for your help in modifying my script.
>>
>> regards, saurabh
>>
>>> On Fri, Oct 6, 2017 at 10:05 AM, Kevin Hallock <hallock at ucar.edu> wrote:
>>> Hi Saurabh,
>>>
>>> I have attached a modified version of your script, “depth_trend_line_mod.ncl”. One major change is that it now checks the output of the regline_stats() function to make sure the result is not a missing value.
>>>
>>> The issue I found was that during the first iteration of your loop (i=0), the value returned by regline_stats() contained a missing value and did not have the “b95” attribute.
>>>
>>> Also, I chose to explicitly declare temp (“temp = new(1, double)”) before the loop, instead of implicitly declaring it inside the loop (“temp = regline_stats(time,ts(i,:))”). When regline_stats() experiences an error, the value returned is a float instead of a double, which causes issues in later iterations of the loop (with an error such as “fatal:["NclVar.c":1390]:Assignment type mismatch, right hand side can't be coerced to type of left hand side”).
>>>
>>> I hope this helps!
>>> Kevin
>>>
>>>> On Oct 4, 2017, at 1:15 AM, saurabh rathore <rohitsrb2020 at gmail.com> wrote:
>>>>
>>>> Dear NCL users,
>>>>
>>>> I am trying to save the regression coefficient value from a 2d variable which is in depth and time and also its 95% confidence interval i.e. b95 in regliine_stat function. I want to save the upper and lower bound value of b95 but somehow the function is not working. Please help me out for this. I am attaching the data file and the code that I am using.
>>>>
>>>> regards, saurabh
>>>>
>>>> --
>>>>
>>>>
>>>> REGARDS
>>>>
>>>> Saurabh Rathore
>>>> Research Scholar (PhD.)
>>>> Centre For Oceans, Rivers, Atmosphere & Land Science Technology
>>>> Indian Institute Of Technology, Kharagpur
>>>> contact :- 91- 8345984434
>>>> <depth_trend_line.ncl><heat_gd1.nc>_______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>>
>>
>> --
>>
>>
>> REGARDS
>>
>> Saurabh Rathore
>> Research Scholar (PhD.)
>> Centre For Oceans, Rivers, Atmosphere & Land Science Technology
>> Indian Institute Of Technology, Kharagpur
>> contact :- 91- 8345984434
>>
>> _______________________________________________
>> 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
> _______________________________________________
> 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/20171006/8372d07e/attachment.html>
More information about the ncl-talk
mailing list