用fortran提取站点数据出不来图

鸾依雨觞 · 发表于 2015-4-16 10:57:57

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

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

x

本人想提取站点数据，编辑了程序，程序可以运行，stnmap后有map文件，但是用open语句打开ctl后图上显示没有把数据写进去，这是什么情况啊，请教高手啊
Fortran程序：
real rlat(31),rlon(31),tim
real avet(31)
integer nlev,nflag,i
character*8 stid(31)
open(1,file='E:\1\shuju\stid.txt',status='old')
open(2,file='E:\1\shuju\avet.txt',status='old')
open(18,file='E:\1\shuju\a.dat',status='unknown',form='binary')
open(19,file='E:\1\shuju\a.txt',status='unknown',form='formatted')
read(1,*)(stid(i),rlon(i),rlat(i),i=1,31)
read(2,*)(avet(i),i=1,31)
tim=0.0
nlev=1
nflag=1
do i=1,31
write(18)stid(i),rlat(i),rlon(i),tim,nlev,nflag,avet(i)
write(19,*)stid(i),rlat(i),rlon(i),tim,nlev,nflag,avet(i)
end do
nlev=0
write(18)stid(1),rlat(1),rlon(1),tim,nlev,nflag
write(19,*)stid(1),rlat(1),rlon(1),tim,nlev,nflag
close(1)
close(2)
close(18)
close(19)
end
ctl文件：
DSET e:\1\shuju\a.dat
DTYPE station
STNMAP e:\1\shuju\a.map
UNDEF -9999.0
TITLE shuju
TDEF 1 LINEAR JUN1951 1dy
VARS 1
avet 0 99 avet
ENDVARS

grads出图情况如图片所示

liuchao0120 · 发表于 2015-4-16 15:19:09

你的图呢？能贴出来你的stid.txt和avet.txt吗？

鸾依雨觞 · 发表于 2015-4-17 11:19:03

liuchao0120 发表于 2015-4-16 15:19
你的图呢？能贴出来你的stid.txt和avet.txt吗？

好的。
图上没有站点数据：

图上没有站点数据，站点主要是云南地区的

stid：
56688 26.1000 103.1667
56772 25.2333 102.5000
56775 25.5833 102.4333
56778 25.0167 102.6833
56781 25.5500 103.2667
56785 25.3333 103.0333
56863 24.9333 102.4833
56871 24.6500 102.6167
56872 24.9500 102.6167
56880 24.9167 103.1667
56881 24.7333 103.2667
56882 24.8833 102.8000
56483 28.6000 103.9500
56489 28.2333 103.6333
56497 28.0667 104.2500
56582 27.7500 103.9000
56585 27.1833 103.5333
56586 27.3500 103.7167
56594 27.6333 104.0500
56595 27.4333 104.9667
56596 27.8500 105.0500
56673 26.9167 102.9167
56684 26.4167 103.2833
56697 26.2167 104.0833
56782 25.4167 103.5500
56783 25.5000 103.8000
56786 23.5833 103.8333
56788 25.0333 103.6667
56790 25.6667 104.2500
56883 24.8333 103.9833
56891 24.9833 104.3167
avet:
197
158
157
156
148
141
153
156
125
164
163
150
179
164
177
143
120
117
170
115
135
206
127
135
137
150
145
150
143
142
151

liuchao0120 · 发表于 2015-4-20 10:42:54

最后一个没有写进去吧，你看你最后一个a.txt文件

鸾依雨觞 · 发表于 2015-4-22 21:22:47

liuchao0120 发表于 2015-4-20 10:42
最后一个没有写进去吧，你看你最后一个a.txt文件

写进去了啊
56688    103.1667    26.10000    0.0000000E+00          1          1
197.0000
56772    102.5000    25.23330    0.0000000E+00          1          1
158.0000
56775    102.4333    25.58330    0.0000000E+00          1          1
157.0000
56778    102.6833    25.01670    0.0000000E+00          1          1
156.0000
56781    103.2667    25.55000    0.0000000E+00          1          1
148.0000
56785    103.0333    25.33330    0.0000000E+00          1          1
141.0000
56863    102.4833    24.93330    0.0000000E+00          1          1
153.0000
56871    102.6167    24.65000    0.0000000E+00          1          1
156.0000
56872    102.6167    24.95000    0.0000000E+00          1          1
125.0000
56880    103.1667    24.91670    0.0000000E+00          1          1
164.0000
56881    103.2667    24.73330    0.0000000E+00          1          1
163.0000
56882    102.8000    24.88330    0.0000000E+00          1          1
150.0000
56483    103.9500    28.60000    0.0000000E+00          1          1
179.0000
56489    103.6333    28.23330    0.0000000E+00          1          1
164.0000
56497    104.2500    28.06670    0.0000000E+00          1          1
177.0000
56582    103.9000    27.75000    0.0000000E+00          1          1
143.0000
56585    103.5333    27.18330    0.0000000E+00          1          1
120.0000
56586    103.7167    27.35000    0.0000000E+00          1          1
117.0000
56594    104.0500    27.63330    0.0000000E+00          1          1
170.0000
56595    104.9667    27.43330    0.0000000E+00          1          1
115.0000
56596    105.0500    27.85000    0.0000000E+00          1          1
135.0000
56673    102.9167    26.91670    0.0000000E+00          1          1
206.0000
56684    103.2833    26.41670    0.0000000E+00          1          1
127.0000
56697    104.0833    26.21670    0.0000000E+00          1          1
135.0000
56782    103.5500    25.41670    0.0000000E+00          1          1
137.0000
56783    103.8000    25.50000    0.0000000E+00          1          1
150.0000
56786    103.8333    23.58330    0.0000000E+00          1          1
145.0000
56788    103.6667    25.03330    0.0000000E+00          1          1
150.0000
56790    104.2500    25.66670    0.0000000E+00          1          1
143.0000
56883    103.9833    24.83330    0.0000000E+00          1          1
142.0000
56891    104.3167    24.98330    0.0000000E+00          1          1
151.0000
56688    103.1667    26.10000    0.0000000E+00          0          1

hzgloria · 发表于 2015-8-21 15:57:09

问题解决了吗？

zhuzhejl · 发表于 2015-9-24 18:32:24

想问一下各个站点平均值怎么求循环顺序老是很混乱谢谢

下是一种 · 发表于 2015-9-25 09:51:29

正在学习。。。

surelee · 发表于 2016-4-22 19:33:27

生成.map文件了么？

kaijieyue · 发表于 2016-5-17 17:30:19

没有写入数据

		自动登录	找回密码
密码			立即注册

[求助] 用fortran提取站点数据出不来图

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