Re: Missing Value Error

Chris,

The WRF scripts are not set up to handle data that is all missing, which appears to the case here.

I've CC'ed wrfhelp here so they can issue a fix. For the wrfhelp folks, the offending code is:

  mx = max(data)
  mn = min(data)

  if(mn.ne.mx.and.opts.and.isatt(opts,"ContourParameters")) then

Meanwhile, one possible work-around is to set your data to some valid
value, if it is all missing. Then, you will just get a "constant field" type of plot:

   if(all(ismissing(data))) then
     print("Warning: data is all missing, setting it to 1.0")
     data = 0.0
   end if

--Mary

On Mar 26, 2012, at 11:19 AM, Chris Klich wrote:

> Hello, I am currently trying to make a contour plot of several of my WRF output files. I am interpolating model levels onto "nice" heights (1km, 2km, etc.). However, when I try to plot domain 2 at 1km, I get the following error:
>
> fatal:The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing function
> fatal:Execute: Error occurred at or near line 1896 in file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
>
> fatal:Execute: Error occurred at or near line 2671 in file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
>
> fatal:Execute: Error occurred at or near line 81 in file co_interp.ncl
>
> This only happens with domain 2, domains 1 and 3 plot fine. Over this area there is terrain higher than 1km, but this is the case in domains 1 and 3 also, and it just plots these areas as white boxes. Why is it not doing this for domain 2? I thought if it had a missing value it would plot it as a white box and move on, but something seems to be stopping this.
>
> Does anyone have any suggestions that may help getting this to run properly? I attached the program if it will be helpful.
>
> Thanks,
> Chris
> <co_interp.ncl>_______________________________________________
> 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
Received on Mon Mar 26 14:48:42 2012