做风场和海温序列的回归遇到的问题，求大神指导。。。

距离

我是用65年平均的海温再做区域平均得到一个65年的时间序列，然后把这个时间序列和全球的风场做回归：
1：这是我的海温指数的fortran程序：
do j=1,iy
    do i=1,ix
       do m=1,65
            p=0
             sst1(i,j,m)=0.0
              do t=9,11
                       if(sst(i,j,t)/=32767.00)then
               sst1(i,j,m)=sst1(i,j,m)+sst(i,j,t+(m-1)*12)
                           p=p+1
                           else
                           sst1(i,j,m)=sst1(i,j,m)+32767.00
                           p=p+1
              end if
                        end do
          sst2(i,j,m)=sst1(i,j,m)/p
                end do
        end do
end do
do m=1,65
   sst21(m)=0.0
      p=0
      do j=47,55
       do i=52,60
            if(sst2(i,j,m)/=32767.00)then
           sst21(m)=sst21(m)+sst2(i,j,m)
            p=p+1
                   else
                end if
         end do
  end do
sst2a(m)=sst21(m)/p
end do

irec7=0
do m=1,65
irec7=irec7+1
write(7,rec=irec7)sst2a(m)
end do


我多次检查过了，算出来的区域平均值是对的
2：然后再把这个时间序列和u，v分别计算回归系数bu,bv。以下是我的回归程序：
!!!!!!!!!!!以下是求秋季的回归场，先取出uv的年平均值
do j=1,iy
    do i=1,ix
           do m=1,65
          yu1(i,j,m)=0.0
                     p=0
           do t=9,11
            yu1(i,j,m)=yu1(i,j,m)+u(i,j,t+(m-1)*12)
            p=p+1
                   end do
         end do
     end do
end do
yu=yu1/3


do j=1,iy
    do i=1,ix
           do m=1,65
          yv1(i,j,m)=0.0
                     p=0
           do t=9,11
            yv1(i,j,m)=yv1(i,j,m)+v(i,j,t+(m-1)*12)
            p=p+1
                   end do
         end do
     end do
end do
yv=yv1/3

ys=0.0
！！！！！！！！！！！！开始回归系数的计算bu.....。。。。。。
do j=1,iy
    do i=1,ix
           do m=1,65
        ys(i,j)=ys(i,j)+yu(i,j,m)*sstx(m)
                end do
         end do
end do

yua=0.0
do j=1,iy
   do i=1,ix
       do m=1,65
       yua(i,j)=yua(i,j)+yu(i,j,m)
           end do
        end do
end do

sstx1a=0.0
do m=1,65
sstx1a=sstx1a+sstx(m)
end do

sstx2=0.0
do m=1,65
sstx2=sstx2+sstx(m)**2
end do
do j=1,iy
   do i=1,ix
   bu(i,j)=(ys(i,j)-((1/65)*yua(i,j)*sstx1a))/(sstx2-(1/65)*(sstx1a**2))
   end do
end do

！！！！！！！！！！！！！计算bv..。。。。。。。。。
ysv=0.0
do j=1,iy
    do i=1,ix
           do m=1,im
        ysv(i,j)=ysv(i,j)+yv(i,j,m)*sstx(m)
                end do
         end do
end do
yva=0.0
do j=1,iy
   do i=1,ix
       do m=1,im
       yva(i,j)=yva(i,j)+yv(i,j,m)
           end do
        end do
end do
do j=1,iy
   do i=1,ix
   bv(i,j)=(ysv(i,j)-(1/65)*yva(i,j)*sstx1a)/(sstx2-(1/65)*(sstx1a**2))
   end do
end do
write(6,rec=1)((bu(i,j),i=1,ix),j=1,iy)
write(7,rec=1)((bv(i,j),i=1,ix),j=1,iy)
stop
end

3,计算完毕，然后画出来的bu,bv矢量图全是纬向风。。。。。。没有涡旋什么的。。无语了。。。

river

上图，上CTL

lqouc

楼主是刚接触气象么？
你说的uv是海面风么？
不管是什么风，你觉得年平均尺度上涡旋会很明显么？就拿你说的纬向风，显然一整年的纬向风风向是基本不变的，但是如果是季风区，那么经向风会有一个年循环，这样平均值后显然就没有了。
而且不明白你做回归是什么意图，区域平均的海温序列和该区域每个点的风做回归你想表达什么物理意义呢？

距离

就是做了年平均后，就去掉了海温本身的年变化啊，研究的主要是年际变化，所以要取年平均的海温和风场啊，区域平均的海温序列和周围区域的风做回归，就是要找出正相关的区域对应的风场是辐合还是辐散的。。。。然后再看春夏秋冬的年际变化。。。

最后得到的相关场和风场的回归应该要对应起来；；；
下边是我的相关图和回归图：

距离

我的第一张图设置的i=144，后边的几张图设置的i=145.这样设置的原因是在grads里边q ctlinfo的时候看到的数据是144个点，但是q dim 的时候看到x的变化时1到145.。。。。。不明白这是怎么回事，所以做了以上改动。

小蝌蚪

你好，我有个问题不明白，我是用指数回归风场。做法跟楼主差不多。但是结果的显著性检验怎么做呢？一个要素的回归检验我是懂的，但是风场有uv两个要素，是分别做回归后在合成画的，那么它的显著性怎么求呢？求指导~谢谢~@lqouc@距离

流浪的心

风场显著性检验

赵尼玛

同小蝌蚪一样，想问一下风场的显著性检验如何做？  

guoqu510

同样的问题，求教。

Sarah

风场的回归到底是怎么检验的呢？是用全风速来做相关性检验吗？