格点插值

自由行

现在有一列数据是土壤湿度，我想在第一列加上对应的经纬度，但是又需要插值为0.5*0.5的分辨率，原数据是1*·1分辨率我的nc数据比较大，所以不能通过grads打开，就先通过ncl转为了txt，但是转出来以后是一列数据该怎么插值呈0.5*0.5喃，我在原始数据里面就进行了经纬度插值，不知道对不对哦
integer i,j,k,m,h
character*80 filename(1140)
parameter(m=160,n=320)
dimension  x(m,n)
open(30,file='F:\work\ch4.2\yugu\5\filenames.txt')
do j=1,1140
   read(30,*)filename(j)
   print*,trim(filename(j))   
enddo
do k=1,1140
open(1,file='F:\work\ch4.2\yugu\5\'//filename(k))
open(2,file='F:\work\ch4.2\yugu\6\'//filename(k))              !提取后的数据       修改处~~~~~~~~~~~~~~
read(1,*) ((x(i,j),j=1,n),i=1,m)
!print*,x(1,1)
close(1)     
do h=1,160
   do i=h,m,160
    do j=1,n
write(2,"(2f10.4,2x,f20.2)")(h-1)*0.5-89.14,(j-1)*0.5+0,x(h,j)
  end do
end do
   !f(k<696)print*,'剩余',696-k,'个文件文件'
   !if(k==696)print*,'数据处理完成'
enddo
end do
close(2)
end

lqouc

请先尝试自己编程，遇到具体问题再来提问。