[ncl-talk] problem with missing values

[Ipshita Majhi]

The fill value is -9999 for x364 , and I had put
set_default_fillvalue("float", -9999.0) at the start of the code.
Thank You
Best Regards
Ipshita


On Friday, December 12, 2014, Dennis Shea <shea at ucar.edu> wrote:

> printVarSummary(x364)
>
> What is the _FillValue ?
>
> ---
>
> et_default_fillvalue("float", -9999.0)
> pweek=new((/47,52/),"float")            ; pweek at _FillValue
>
>
>
> do i=0,47-1
> do j=0,357,7
>
> pweek(i,j/7)=avg(x364(i,j:j+6))
>
>       if x364 has an _FillValue of 9.92226e+36
>       and the result of avg is a missing value,
>       the _FillValue of x354 will transfer to pweek
>      Also, 9.92226e+36 is the 'printed' value, This is not the
>      actual value. Rounding has occurred for printing purposes.
>
>       if x364 at _FillValue = 9.92226e+36
>       and you prefer -9999.0 then set
>       x364 at _FillValue == pweek at _FillValue *before* the loop
>
>
> ;  end if
> end do
> end do
>
> I think the above should fix your issue.
>
> The whole loop below could be done much more efficiently
>
>
>   pweek = where(pweek.gt. 1e25, pweek at _FillValue, pweek)
>
>
> do i=0,47-1
> do j=0,51
> print(i)
> print(j)
> if (pweek(i,j).eq.9.92226ee+36)then
> pweek(i,j)=9999.0
> end if
> end do
> end do
>
>
>
>
> On Fri, Dec 12, 2014 at 3:10 PM, Ipshita Majhi <ipmajhi at alaska.edu> wrote:
>>
>> Thank you
>>
>> Yes it is a typo in the code I fixed it yet it still gives max value as
>> 9.6921e+36 , and if I write to substitute it then it gives 9.92226e+36
>> as the max value.
>>
>> Best Regards
>> Ipshita
>>
>>
>>
>> On Friday, December 12, 2014, Gus Correa <gus at ldeo.columbia.edu> wrote:
>>
>>> Hi Ipshita
>>>
>>> Is this a typo in your email or in the code?
>>>
>>>  > if (pweek(i,j).eq.9.92226ee+36)then
>>>
>>> Note the double "ee".
>>>
>>> Gus Correa
>>>
>>> On 12/12/2014 04:29 PM, Ipshita Majhi wrote:
>>> > Dear NCL ,
>>> >
>>> > I have written a code to write to convert 365 days into 52 weeks avg.
>>> > I am facing error regarding missing value. I have written statements
>>> > which have substituted missing value as -9999.0
>>> > Yet it gives the maximum value as 9.96921e+36 and when I write to
>>> > substitute that with -9999.0, it gives maximum value as 9.969222e+36.
>>> > Here is the snippet of my code. I am very confused about this. It will
>>> > be great if somebody could help me
>>> >
>>> > set_default_fillvalue("float", -9999.0)
>>> >
>>> > pweek=new((/47,52/),"float")
>>> >
>>> >
>>> > do i=0,47-1
>>> > do j=0,357,7
>>> >
>>> > pweek(i,j/7)=avg(x364(i,j:j+6))
>>> >
>>> > ;  end if
>>> > end do
>>> > end do
>>> >
>>> >
>>> > do i=0,47-1
>>> > do j=0,51
>>> > print(i)
>>> > print(j)
>>> > if (pweek(i,j).eq.9.92226ee+36)then
>>> > pweek(i,j)=9999.0
>>> > end if
>>> > end do
>>> > end do
>>> >
>>> > printMinMax(pweek, 0)
>>> >
>>> > Thank You
>>> > Best Regards
>>> > Ipshita
>>> >
>>> >
>>> >
>>> > _______________________________________________
>>> > ncl-talk mailing list
>>> > List instructions, subscriber options, unsubscribe:
>>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>> >
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> 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/20141212/0419e3c8/attachment.html