求grd转txt编程方法

资本虾

以前会用fortran将txt转换成grd格式的数据文件用于GRADS作图，现有grd格式的全球月平均气温资料想转换了ASC11码的，可用txt打开的数据文件，求个程序指导,最好是用fortran的 多谢！

小蝌蚪

我编的一个小程序，比楼上的简单一点，不知你是否能看明白啊？这和txt转grd差不多的。
program main
implicit none
integer i,j,n
parameter (n=168)
real ss(1,n)
open(100,file='d:\pjtest\hgt35155.grd',form='binary')
open(200,file='d:\pjtest\hgt35155.txt')
do j=1,n
       read(100) ss(1,j)
        write(200,*) ss(1,j)
     enddo
close (100)
close (200)
end

mofangbao

理解grads的数据存储方式你就会了

river

是一样过程啊，用fortran把grd数据读出来，再写成txt的就行了啊。如果你的资料是grads可以识别的，并且有配套的ctl什么的话，使用兰溪的一个脚本就可以写成txt的了

herenso

跟把txt转成grd相比，grd转txt在fortran里不就是个逆过程么~grd的排放就跟你用fortran写出来的是一样的吧。

topmad

一正一反而已  看看grads里面的数据是如何存放的就知道鸟

资本虾

river 发表于 2013-3-22 08:13

                               
登录/注册后可看大图

是一样过程啊，用fortran把grd数据读出来，再写成txt的就行了啊。如果你的资料是grads可以识别的，并且有配 ...

我的资料是用nc转成的grd，就想给它转成txt的，兰溪的脚本在哪？多谢！

言深深

给你一段代码，看懂了再改！
parameter(m=73,n=144,o=778)
!mn空间格点数目，o时间长度
open(1,file='kaka.dat',access='direct',form='binary',recl=n*m*o*4)
read(1,rec=1) (((x(i,j,k),j=1,n),i=1,m),k=1,o)
close(1)     

open(2,file='lev1000.txt')
do i=1,m
do j=1,n
write(2,"(<o>f10.2)")(x(i,j,k),k=1,o)
enddo
enddo
close(2)

river

资本虾 发表于 2013-3-22 09:23

                               
登录/注册后可看大图

我的资料是用nc转成的grd，就想给它转成txt的，兰溪的脚本在哪？多谢！

在这呢，自己看看吧。http://bbs.06climate.com/forum.php?mod=viewthread&tid=11515

kongfeng0824

言深深 发表于 2013-3-22 09:25

                               
登录/注册后可看大图

给你一段代码，看懂了再改！
parameter(m=73,n=144,o=778)
!mn空间格点数目，o时间长度

深深版主大才！