求助！nc资料怎么提取4-6月降水资料

尼罗河星星

我的数据如下   希望提取1978~2013年4-6月的降雨资料，用grads怎么编写程序呢，要判断闰年还是平年吗?求大神协助！！


我编的程序是如下：
'reinit'
'set gxout fwrite'
'set fwrite F:\2\test\allpre4_6.grd'
'sdfopen F:\2\test\allpre.nc'

year=1979
snum=6574
while(year<=2013)
  t1=year-1956
    if(t1=0|tt1=4|t1=8|t1=12|t1=16|t1=20|t1=24|t1=28|t1=32|t1=36|t1=40|t1=44|t1=48|t1=52|t1=56|t1=60|t1=64)
    num=366
    else
    num=365
    endif

say year
say num

*资料截取要用格点，不能用经纬度
    iday=iday+snum+1   
    iiday=iday
        ttt=iiday+90
         while(iday<=ttt)
         'set t 'iday''
         'set x 1 128'
         'set y 1 72'
         'd pre'
         iday=iday+1
         endwhile
year=year+1

endwhile
'disable fwrite'
'reinit'
;



但是提示报错为

影形

请问楼主的降水数据是怎么下的？我只下到降水率，还是每年单独下载的

尼罗河星星

影形 发表于 2015-4-9 08:06
请问楼主的降水数据是怎么下的？我只下到降水率，还是每年单独下载的

这个是老师给我的

qianqian_sunny

影形 发表于 2015-4-9 08:06
请问楼主的降水数据是怎么下的？我只下到降水率，还是每年单独下载的

同问同问，NCEP的月平均降水都是降水率，不知道能在哪里下到降水资料？

qianqian_sunny

尼罗河星星 发表于 2015-4-9 13:12
这个是老师给我的

楼主，如果这个资料是公共资料的话能不能给我一份？

影形

qianqian_sunny 发表于 2015-4-16 17:09
楼主，如果这个资料是公共资料的话能不能给我一份？

我不是楼主，我也求资料来着/(ㄒoㄒ)/~~

qianqian_sunny

影形 发表于 2015-4-17 21:48
我不是楼主，我也求资料来着/(ㄒoㄒ)/~~

我最近下资料的时候发现，月平均降水率的单位是（mm/d)，我觉得是每个月的总降水量除以这个月的天数，所以我们应该能得每个月的总降水量，仅仅是我的想法

影形

qianqian_sunny 发表于 2015-4-18 16:02
我最近下资料的时候发现，月平均降水率的单位是（mm/d)，我觉得是每个月的总降水量除以这个月的天数，所 ...

我差NCEP/NCAR数据只有可降水率（是你说的平均降水率么），量级在10^-4,我*3600*24才得到降水量的数量级。。。

qianqian_sunny

影形 发表于 2015-4-19 15:49
我差NCEP/NCAR数据只有可降水率（是你说的平均降水率么），量级在10^-4,我*3600*24才得到降水量的数量级 ...

我下的是CMAP的降水资料，数值0—70mm/day，看着很正常

mengnannan1314

我是提取一月的几天在计算的，看一下楼主的。