cressman插值错误，求大神指点

Dai

程序算法部分如下，输入的数据我都检查了应该没错误，但是插值后数据、出图总是不对，算法我觉得也没错，但应该就是算法部分的问题，请大神指点啊！！！！！！！！！
!!!!!!!!!!计算距离
do i=1,numsdlon
  do ii=1,numsdlat
  smaxt=0.0
  total_weight=0.0
  ilevs=0
   do j=1,total
    slon1 = slon(j)*pi/180 !站点经度
    glon1 = glon(i)*pi/180 !格点经度
    slat1 = slat(j)*pi/180 !站点纬度
    glat1 = glat(ii)*pi/180 !格点纬度
    dist=cos(slat1)*cos(glat1)*cos(slon1-glon1)+sin(slat1)*sin(glat1)
    dist=a*acos(dist)                  
if(dist<=radius.and.maxt2(j)/=xmiss)then
!!!!!!!!!!!!!!计算权重!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    ilevs=ilevs+1
          weight=(radius**2-dist**2)/(radius**2+dist**2)
          smaxt=smaxt+(weight*weight*maxt2(j))
          total_weight=total_weight+weight            
    endif
   enddo
  print*,ilevs
   if(ilevs>=1)then
   gmaxt(i,ii)=smaxt/total_weight
   else
   gmaxt(i,ii)=xmiss
   endif
  enddo
enddo

lqouc

算法是自己写的？
如果不是自己写的那估计是ctl或者数据io格式的问题。

Dai

lqouc 发表于 2016-3-1 15:56
算法是自己写的？
如果不是自己写的那估计是ctl或者数据io格式的问题。

算法是找了资料、公式改的，所以算法是对的吗

Dai

lqouc 发表于 2016-3-1 15:56
算法是自己写的？
如果不是自己写的那估计是ctl或者数据io格式的问题。

而且我把数据存到txt看了，感觉图没画错，出来的数据就是这样的，所以觉得是不是算法不对

Dai

lqouc 发表于 2016-3-1 15:56
算法是自己写的？
如果不是自己写的那估计是ctl或者数据io格式的问题。

确实是输出的经纬度反了 现在好了 谢谢