txt转成dat求助【已解决】

小傅之追梦

麻烦贴出CTL，并且把画图时的报错截图

小傅之追梦

你的CTL改成
dset d:/fprogram/egvt.dat
title qingzang
unpack scale_factor add_offset 1 32066
undef -9.99e+33
xdef 13 linear 75 2.5
ydef 1 linear 0 1
zdef 17 levels 1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10
tdef 1 linear 00Z01JAN1948 1yr
vars 1
air 17 t,z,y,x mean daily air temperature
endvars
请尝试一下

驍放假了

帮顶。。。。。。

lqouc

希尔博拉斯 发表于 2014-3-22 21:19

                               
登录/注册后可看大图

对啊。。设置
set x 1 13
set y 1

fortran的输出你改了？

希尔博拉斯

lqouc 发表于 2014-3-22 21:31

                               
登录/注册后可看大图

fortran的输出你改了？

求问fortran该怎么改。。

希尔博拉斯

小傅之追梦 发表于 2014-3-22 21:24

                               
登录/注册后可看大图

你的CTL改成
dset d:/fprogram/egvt.dat
title qingzang

现在能出图了，可是图显示有问题，似乎是fortran的问题。

lqouc

希尔博拉斯 发表于 2014-3-22 21:34

                               
登录/注册后可看大图

求问fortran该怎么改。。

把rec=i和recl=4去掉，然后ctl里面air后面写成17 99

小傅之追梦

希尔博拉斯 发表于 2014-3-22 21:34

                               
登录/注册后可看大图

现在能出图了，可是图显示有问题，似乎是fortran的问题。

我改动了你的Fortran输出，现在正常
parameter(m=17,n=13)
dimension  x(m,n)
integer i,j
open(11,file='egvt1.txt')
read(11,*) ((x(i,j),j=1,n),i=1,m)
close(11)  

open(12,file='egvt.grd',form='binary')
write(12) ((x(i,j),j=1,n),i=1,m)
close(12)
end

river

先去看一下grads要求的数据格式，这个循环明显不对，循环由外向内的顺序应该是  t z y x

小傅之追梦

出图如下~~~