用grads对28年的逐6小时nc格式资料进行读取并输出日平均的二进制格式

杨柳松

这是代码：
'reinit'
year=1985
while(year<=2012)
month=1
while(month<=12)

    if(month<10)
        'sdfopen I:\JRA25\tmpprs\tmpprs_'year'0'month'.nc'
        'set x 1 288'
        'set y 1 145'
        'set fwrite I:\JRA25\tmpprs_'year'0'month'.dat'
        'set gxout fwrite'
      else
        'sdfopen I:\JRA25\tmpprs\tmpprs_'year''month'.nc'
        'set x 1 288'
        'set y 1 145'
        'set fwrite I:\JRA25\tmpprs_'year''month'.dat'
        'set gxout fwrite'
    endif
    if(month=1|month=3|month=5|month=7|month=8|month=10|month=12)
         data=31
       else
         if(month=4|month=6|month=9|month=11)
           data=30
             else
               if(month=2&(year=year=1988|year=1992|year=1996|year=2000|year=2004|year=2008|year=2012))
                 data=29
                else
                  if(month=2)
                    data=28
                  endif
               endif
         endif
     endif


m=1
while(m<=data)
t1=(m-1)*4+1
t2=m*4
   n=1
   while(n<=17)
  'set z 'n''
  'd ave(tmpprs,t='t1',t='t2')'

  n=n+1
  endwhile
m=m+1
endwhile
'reinit'
month=month+1
endwhile
year=year+1
endwhile

杨柳松

有个语句错误了 if(month=2&(year=1988|year=1992|year=1996|year=2000|year=2004|year=2008|year=2012))

sjysjxyf66

向楼主学习！