【已解决】Fortran提取出数据后，用grads画图，总出现Z的维数有问题，请求前辈指导

shirly · 发表于 2013-11-28 09:23:20

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 shirly 于 2013-12-1 11:49 编辑

用fortran提取数据，编码如下：
program Pro_dat
implicit none
integer,parameter:: n=4000
integer i,j,k,e,z,day,it
character start*80
REAL T(n),P(n),W(n)

!读数
open(1,file='d:/kuoxian/2009tk/2009yitiqu/全部探空数据s5298320090109.19.txt',status='old')
do i=1,14
read(1,*)start
print*,start
end do
do i=1,n
read(1,'(12x,f6.2,4x,f6.2,7x,f3.0)',iostat=it) T(i),P(i),W(i) !,P(i),W(i)
if(it<0) exit
write(*,'(2(f6.2,4x),6x,f3.0)')T(i),P(i),W(i)
end do
!写入提取出来的数据
open(2,file='d:/kuoxian/2009tk/tiquTPW/TPW.dat',form='unformatted',access='direct',recl=1)
do i=1,n-1
write(2,rec=i) T(i)
end do
end

生成的数据说明文件如下：
dset d:/kuoxian/2009tk/tiquTPW/TPW.dat
title T
undef 9.96921e+36
xdef 1 linear 72.75 0.5
ydef 1 linear 16.75 0.5
zdef 2942 linear 1 1
tdef 1 linear 19Z09Jan2009 1mn
vars 1
T 0 99 t
endvars

运行的gs文件如下：
'reinit'
'open d:/kuoxian/2009tk/tiquTPw/TPW.ctl'
'set x 1'
'set y 1'
'set t 1'
'set z 1 100'
'd T'
'printim d:/kuoxian/gs/tu/T.png white'
'reinit'

运行出错信息：

出错信息