Re: missing value issue

Alexander

The function 'readAsciiTable' doesn't assing a '_FillValue to to it's output, (i.e. variables simdata and obsdata in your script) so you need to assign the attribute your self, i.e.

> obsdata@_FillValue = -999.

For future reference, a useful test for this kind of problem is the 'printMinMax()' function, since it won't return any recognised missing values, so if

printMinMax(obsdata, True) returns '-999.' you know that the _FillValue attribute needs to be set.

Will

**********************************************************
Will Hobbs, Ph.D. William.R.Hobbs@jpl.nasa.gov
Jet Propulsion Laboratory
4800 Oak Grove Dr. office: 300-324a
M/S 300-323 phone: (818) 354-0466
Pasadena, CA 91109 fax: (818) 354-0966
**********************************************************

From: Alexander Semenov <asemenov@alaska.edu<mailto:asemenov@alaska.edu>>
Date: Fri, 13 Apr 2012 10:29:12 -0800
To: <ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>>
Subject: missing value issue

Hi,

I have a problem with missing values. I'm comparing model output and observations (both txt files). Some observations have missing values (-999.000) and apparently ncl thinks of them as of real values and calculates bias, RMSE....

I can't figure out how I can do it in ncl for the particular code (!!): to get it to understand that if the observations have a missing value = -999 to skip this value and go to the next line.

--
regards
*******************************************************
Alexander Semenov
PhD Student - Research Assistant
International Arctic Research Center
Department of atmospheric sciences
University of Alaska Fairbanks
930 Koyukuk dr, 408c3
Fairbanks, AK, USA, 99775
work phone +19074742672
_______________________________________________ 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 Fri Apr 13 12:44:42 2012