Missing value problem?

From: Jiang, Lifen <lfjiang_at_nyahnyahspammersnyahnyah>
Date: Sat Jun 01 2013 - 10:47:51 MDT

Dear NCL members,

I am calculating annual average of my monthly data using NCL "month_to_annual" function. However, it always returned empty file (all grids are missing value).

I don't think the codes (see below) have problems because they worked very well for other data but for this. I guess the issue is related with missing value of the data file but cann't find what it is.

I have uploaded the data file to UCAR's ftp under folder "incoming". I appreciate if anyone can help check what is the problem.

Thank you!



    load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
    f1 = addfile("npp_Lmon_NorESM1-M_historical_r1i1p1_185001-200512.nc", "r")
    t1 = f1->npp ; t1(time,lat,lon)
    t1@_FillValue = 0
    t1@missing_value = t1@_FillValue
    Yearly = month_to_annual(t1(time|:,lat|:,lon|:), 1) ; TMP(156,nlat,mlon)
    system("/bin/rm -f yearly_npp_NorESM1-M.nc") ; remove any pre-exist file
    fout = addfile("yearly_npp_NorESM1-M.nc", "c") ; new netCDF file
    fout@title = "Month to year"
    fout->Yearly = Yearly

ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
Received on Sat Jun 1 10:48:05 2013

This archive was generated by hypermail 2.1.8 : Tue Jun 04 2013 - 09:50:17 MDT