如何求区域平均

呼啦呼啦

我用的是一个全球卫星网格点的逐候亮温资料（如果有做亮温资料处理的，请尽管联系我吧~~），主要是用于对气候事件的诊断分析，所以我就先做了该资料的距平， 得到了一个grd，里面存放的是1999-2009年这11年每个月与该月11年平均值的一个距平值，我现在想做一个亮温异常（即距平）的时间序列，便将全球格点求了个面积平均，可是得到的结果都是0，不知道怎么回事。。。。
  可能我说的比较啰嗦，总的来说就是

   1、如何求全球格点的一个平均，以得到一个一维的时间序列图呢？


   2、亮温可以用于什么样的气候事件中呢？

呼啦呼啦

subroutine sumaverage(jp,ii,jj,m,iy)                                 !全球总格点平均的距平
implicit none
integer i,j,ii,jj,m,iy,k,n
real sumjp(m,iy),jp(ii,jj,m,iy)
do n=1,iy
  do k=1,m
         sumjp(k,n)=0.0
  enddo
enddo
do n=1,iy
   do k=1,m
    do i=5,140
     do j=5,68
           sumjp(k,n)=sumjp(k,n)+jp(i,j,k,n)
         enddo
    enddo
     sumjp(k,n)=sumjp(k,n)/8704
   enddo
  enddo
end subroutine


这是我编的一段求全球格点平均的程序，大家给看看 对不对啊？

river

呼啦呼啦 发表于 2014-8-21 18:58
subroutine sumaverage(jp,ii,jj,m,iy)                                 !全球总格点平均的距平
implicit none
integer i,j,ii,jj,m ...

赋初值的位置应该是这样
do n=1,iy
   do k=1,m
sumjp(k,n)
    do i=5,140
     do j=5,68
           sumjp(k,n)=sumjp(k,n)+jp(i,j,k,n)
         enddo
    enddo
     sumjp(k,n)=sumjp(k,n)/8704
   enddo
  enddo

lqouc

你是指这个帖子么？其实还好啦，我最早的发的帖子都是低级的，一步步走过来，正好可以见证进步啊。

呼啦呼啦

lqouc 发表于 2016-1-4 20:17
你是指这个帖子么？其实还好啦，我最早的发的帖子都是低级的，一步步走过来，正好可以见证进步啊。

才发现我也可以改帖子名字，已经搞定啦！

jintao

{:5_196:}{:5_196:}