[ncl-talk] If statements and Fill_Values

[Rashed Mahmood]

I think there is some logical issue(s) in your if statements:
Issues:
1)  if(CL(n).eq.CL at _FillValue)then ; This would not work because the
"@_FillValue" attribute cannot be used inside if statements, I think. You
can define a missing value variable before hand, e.g:

missV = CL at _FillValue
print(missV)
 if(CL(n).eq.missV)then
...

2) if(all(ismissing(CL)))then : For this statement since you are asking "if
ALL Values  of CL are missing which might not be the case since some values
would be valid.

Does your variable CL has _FillValue attribute assigned? What is
printVarSummary(CL) suggest?

*The simplest statement would be( if CL has "_FillValue" attribute assigned
to it)*:

*if(ismissing(CL(n)))then*
  continue
else
 ...
end if

hope that helps.
Rashed






On Fri, Jan 10, 2020 at 1:21 AM Beáta Szabó-Takács via ncl-talk <
ncl-talk at ucar.edu> wrote:

> Dear Karin and Barry,
>
> I checked different ways to skip the _FillValues. Here is my brief script:
>
>
>
> begin
>
>
>
>   input = addfile("NetCDF_2017_07_18.nc","r")
>
>   CL=input->cLayer1
>                                                   ;First cloud layer, its
> _FillValue is -9900
>
>   CL1 = new((/5400/), typeof(CL), 1e20)
>
>   CL2 = new((/5400/), typeof(CL), 1e20)
>
>   CL3 = new((/5400/), typeof(CL), 1e20)
>
>   CL4 = new((/5400/), typeof(CL), 1e20)
>
>
>
>    do n=1,5399
>
>
>
>              if(CL(n).eq.CL at _FillValue)then
>
>                   continue
>
>              else
>
>                  CL1(n) = CL(n)
>
>               end if
>
>
>
>             if(all(ismissing(CL)))then
>
>                    continue
>
>            else
>
>                   CL2(n) = CL(n)
>
>            end if
>
>
>
>             if(CL(n).ne.CL at _FillValue)then
>
>                   CL3(n) = CL(n)
>
>             end if
>
>
>
>            if(.not.ismissing(CL(n)))then
>
>                 CL4(n) = CL(n)
>
>            end if
>
>
>
>    end do
>
>
>
>    print(CL1(1:500))
>
>    print(CL2(1:500))
>
>    print(CL3(1:500))
>
>    print(CL4(1:500))
>
>
>
> end
>
>
>
> Unfortunately, none of the if statements skip the _FillValues. I do not
> understand why.
>
> Beata
>
>
>
> *From:* ncl-talk [mailto:ncl-talk-bounces at ucar.edu] *On Behalf Of *Karin
> Meier-Fleischer via ncl-talk
> *Sent:* Friday, January 10, 2020 9:26 AM
> *To:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] If statements and Fill_Values
>
>
>
> Hi,
>
> what about
>
> if(all(ismissing(ratio_hr))) then
>    continue
> end if
>
> -Karin
>
> Am 10.01.20 um 09:17 schrieb Barry Lynn via ncl-talk:
>
> Hi:
>
>
>
> I will give it a try.
>
>
>
> In the meantime, I am wondering why you suggest the ne, rather than the
> equal to with a continue and an else?
>
>
>
> On Fri, Jan 10, 2020 at 9:52 AM Beáta Szabó-Takács <szabo.b at czechglobe.cz>
> wrote:
>
> Dear Barry,
>
> The continue will proceed directly to the next iteration. Try to skip the
> _FillValues by:
>
>
>
> if(max(ratio_hr).ne.ratio_hr at _FillValue.and. min(ratio_hr)
> .ne.ratio_hr at _FillValue)then
>
> statement(s)
>
> end if
>
>
>
> I hope it solves your issue.
>
> Kind regards,
>
> Beata
>
>
>
> *From:* ncl-talk [mailto:ncl-talk-bounces at ucar.edu] *On Behalf Of *Barry
> Lynn via ncl-talk
> *Sent:* Friday, January 10, 2020 7:40 AM
> *To:* ncl-talk <ncl-talk at ucar.edu>
> *Subject:* [ncl-talk] If statements and Fill_Values
>
>
>
> Hi:
>
>
>
> I am trying to skip a loop where all the numbers (min/max) are _FillValues.
>
>
>
> min=9.96921e+36   max=9.96921e+36
>
>
>
> However, I get this error, even though my if statement is checking this.
>
>
>
> fatal:The result of the conditional expression yields a missing value. NCL
> can not determine branch, see ismissing function
>
>
>
>   if (max(ratio_hr)
> .eq.ratio_hr at _FillValue.and.min(ratio_hr).eq.ratio_hr at _FillValue)then
>
>    continue
>
>   else
>
>
>
> Is there someway around this?
>
>
>
> Thank you
>
> --
>
> Barry H. Lynn, Ph.D
>
> Senior Associate Scientist, 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
>
>
>
>
> --
>
> Barry H. Lynn, Ph.D
>
> Senior Associate Scientist, 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
>
>
>
> _______________________________________________
>
> ncl-talk mailing list
>
> ncl-talk at ucar.edu
>
> List instructions, subscriber options, unsubscribe:
>
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> _______________________________________________
> 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/20200110/7ee67b9b/attachment.html>