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!
Lifen
****************
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)
printVarSummary(t1)
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:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
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