处理cmap的候降水资料

飞王子

求助：cmap逐候的降水资料，在它的网站上有从1979至2014年的，以及气候态的1979-2000的候平均的气候态资料，可是我其他的资料使用的气候态是1981-2010年，所以想提取这30年气候态的候平均降水资料。遇到的问题：因为有闰年的存在，导致如下的问题


ga-> set t 148
Time values set: 1981:1:5:0 1981:1:5:0
ga-> set t 221
Time values set: 1982:1:5:0 1982:1:5:0
ga-> set t 294
Time values set: 1983:1:5:0 1983:1:5:0
ga-> set t 367
Time values set: 1984:1:5:0 1984:1:5:0
ga-> set t 440
Time values set: 1985:1:4:0 1985:1:4:0
最直接的体现是
ga-> set t 2338
Time values set: 2010:12:29:0 2010:12:29:0
ga-> set t 148
Time values set: 1981:1:5:0 1981:1:5:0*******************多出了一候

也就是说73*30=2338-148，而不是end-start+1=73*30

希望家园里的朋友们多多指点O(∩_∩)O~

飞王子

问题解决了，在此也说一下吧，人家网上有资料的使用说明就是在网页的下面，引自“http://www.esrl.noaa.gov/psd/data/gridded/data.cmap.html”
Reading in GrADS: If you read the data in GrADS, it will assume the values are 5 days apart the entire time period including leap years. This isn't correct so you may get the wrong dates. To handle this in GrADS, you can use a ddf file rather than using sdfopen and set the calendar options to 365_day_calendar. Or, you can just use the "t" cvalues as each full year has 73 values starting Jan 1.

飞王子

特别要注意一下ddf文件（红色字体）
dset F:\Data\cmap\precip.pentad.mean.include.nc
title CPC Merged Analysis of Precipitation Pentad data (includes NCEP Reanalysis I)1979-2014
undef -9.96921e+36
dtype netcdf
options 365_day_calendar
xdef lon 144 linear 1.25 2.5
ydef lat 72 linear -88.75 2.5
tdef time 2591 linear 00Z01JAN1979 7200mn
vars 1
precip=>precip  0  t,y,x  Average Pentad Rate of Precipitation
endvars
另外，因为是降水的，所以我把set z 删掉了，这样就能用xdfopen打开ddf文件了

皮特

这么好的贴没人顶。。。

皮特

楼主计算侯降水距平百分率，出现了3e+5的结果，是为什么？

皮特

侯平均气候态的1979-2000资料有下载吗，只发现月平均气候态资料

niuda

皮特 发表于 2015-10-5 12:03
楼主计算侯降水距平百分率，出现了3e+5的结果，是为什么？

请教 降水距平百分率 怎么算？

飞王子

皮特 发表于 2015-10-5 12:03
楼主计算侯降水距平百分率，出现了3e+5的结果，是为什么？

额，是不是算错了

飞王子

皮特 发表于 2015-10-5 12:04
侯平均气候态的1979-2000资料有下载吗，只发现月平均气候态资料

1981-2010区间的应该是有的

聆念LN

请问楼主ddf文件中的红色字体是什么意思 没有看懂啊。。