[ncl-talk] problem with missing values

[Ipshita Majhi]

Thank you , my problem has been resolved by:-

pweek = where(pweek.gt. 1e25, pweek at _FillValue, pweek)

I am very grateful for the help.

Best Regards
Ipshita

On Fri, Dec 12, 2014 at 4:27 PM, Ipshita Majhi <ipmajhi at alaska.edu> wrote:
>
> I will try the above mentioned loop that Dennis suggested and will let you
> all know the outcome.
> Thank you all for the support
> Best Regards
> Ipshita
>
> On Fri, Dec 12, 2014 at 4:20 PM, 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/40799d43/attachment.html