txt转成dat求助【已解决】

希尔博拉斯

小弟最近在做一个经度，高度的剖面图，经度为13个格点，高度17层，原文件为txt格式的，现在要把它转成dat画图程序如下：
parameter(m=17,n=13)
dimension  x(m,n)
integer i,j
open(11,file='d:\fprogram\egvt1.txt')
read(11,*) ((x(i,j),j=1,n),i=1,m)
close(11)  

open(12,file='d:\fprogram\egvt.dat',access='direct',form='binary',recl=n*4)
do i=1,m
do j=1,n
write(12,rec=i) x(i,j)
enddo
enddo
close(12)

可是这样的dat我配一个ctl后出不了图，但是如果我把ctl配成13个经度，17个纬度的话，是可以出图的，求解问题出在哪？

lqouc

不是高度么？怎么出来纬度了。

希尔博拉斯

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

                               
登录/注册后可看大图

不是高度么？怎么出来纬度了。

本来应该是17个高度，1个纬度，可是这样画不出图。但是我如果把在ctl中把这17个高度变成17个纬度，一个高度，那么可以画出图。。。

lqouc

把ctl贴出来。
不过先说说你的Fortran啊，你都用binary了就不要写rec了，如果要写rec那就不要对水平空间进行显式循环。grads的数据如何排列是有严格规则的。

小傅之追梦

CTL写成经度-高度，你画图的时候设置set z 1 17   ????

小傅之追梦

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

                               
登录/注册后可看大图

CTL写成经度-高度，你画图的时候设置set z 1 17   ????

默认是第一层的

希尔博拉斯

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

                               
登录/注册后可看大图

把ctl贴出来。
不过先说说你的Fortran啊，你都用binary了就不要写rec了，如果要写rec那就不要对水平空间进 ...

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 1 t,z,y,x mean daily air temperature
endvars
这是ctl

希尔博拉斯

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

                               
登录/注册后可看大图

把ctl贴出来。
不过先说说你的Fortran啊，你都用binary了就不要写rec了，如果要写rec那就不要对水平空间进 ...

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 17 linear 0 1
zdef 1 levels 1000
tdef 1 linear 00Z01JAN1948 1yr
vars 1
air 1 t,z,y,x mean daily air temperature
endvars
如果我把ctl写成这样子，是可以出图的。。

希尔博拉斯

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

                               
登录/注册后可看大图

CTL写成经度-高度，你画图的时候设置set z 1 17   ????

对啊。。设置
set x 1 13
set y 1
set z 1 17
set t 1
可是就是不行啊。。

lqouc

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

                               
登录/注册后可看大图

dset d:/fprogram/egvt.dat
title qingzang
unpack scale_factor add_offset 1 32066

吧air后面的1改成17