提取GrADS插值结果（oacres/gr2stn/grads2ascii）

liutaoz21

楼主您好，我用您的方法处理出类似的数据，不过我用的是1°*1°的格点，最后是501行701列的数据，我想通过fortran转换成数据在一列的情况，即“行数 列数 数据”，用如下程序：
     program main
      PARAMETER(n=501,m=701)
      DIMENSION r(n,m)

!        读入grd格点数据

        open(1,file='D:\wx\gd\dat\gddats\jmgd\2012070915jm-rain.txt')  
!    read(1)((r(i,j),i=1,n),j=1,m)
    do i=1,n
          do j=1,m
       read(1)r(i,j)
        enddo
        enddo
    open(10,file='D:\wx\gd\dat\gddats\shiyan\2012070915jmgd.txt')

!     输出grd格点数据
     
         do j=1,m
          do i=1,n
       write(10,*)i,j,r(i,j)
        enddo
        enddo
           close(10)
           close(15)
end program
运行时出错，显示“forrtl: severe<256>: unformatted I/0 to unit open for formatted transfers,unit 1,file D:\wx\gd\dat\gddats\2012070915jm-rain.txt...”
我实在搞不明白哪里出了问题，求楼主指教。原来生成的格点数据在附件里。不胜感激！

tintin

我想请问一下，用matlab的二维插值函数能够将格点插值为站点么？

mofangbao

liutaoz21 发表于 2012-8-28 21:52

                               
登录/注册后可看大图

楼主您好，我用您的方法处理出类似的数据，不过我用的是1°*1°的格点，最后是501行701列的数据，我想通过f ...

去读的read(1)改为read(1,*)，这是文本文件，是有格式的，不能用无格式的方法来读，好好看下文件读写的章节

mofangbao

tintin 发表于 2012-8-30 08:41

                               
登录/注册后可看大图

我想请问一下，用matlab的二维插值函数能够将格点插值为站点么？

没怎么用过matlab

tintin

那要是每个站点对应了多个时次的数据插值，这个还能适用么？有没有批量的算法

liutaoz21

mofangbao 发表于 2012-8-30 08:44

                               
登录/注册后可看大图

去读的read(1)改为read(1,*)，这是文本文件，是有格式的，不能用无格式的方法来读，好好看下文件读写的章 ...

谢谢楼主，问题已解决，还是自己基本功不扎实问题

mofangbao

tintin 发表于 2012-8-30 09:33

                               
登录/注册后可看大图

那要是每个站点对应了多个时次的数据插值，这个还能适用么？有没有批量的算法

自己举一反三一下吧

tintin

mofangbao 发表于 2012-8-30 12:36

                               
登录/注册后可看大图

自己举一反三一下吧

多谢O(∩_∩)O~

tintin

竟然是小穷神

yuanyuan99

正需要 谢谢