[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