关于如何画风廓线的提问

schLltt

！！！！！！！！！！！！！！！！！！！风速风向由.txt转换成二进制的fortran程序
program ex01
integer i
integer,parameter:: nt=10
real:: wnd1(nt),wnd2(nt),dir1(nt),dir2(nt)

open(101,file='D:\fengkuoxian\wind.txt',status='old')
do i=1,nt
   read(101,*) wnd1(i),wnd2(i)
enddo
close(101)
open(102,file='D:\fengkuoxian\dire.txt',status='old')
do i=1,nt
   read(102,*) dir1(i),dir2(i)
enddo
close(102)

open(103,file='D:\fengkuoxian\wind.dat',form='unformatted')   
  do i=1,nt
    write(103) wnd1(nt),wnd2(nt)
  enddo
close(103)

open(104,file='D:\fengkuoxian\dire.dat',form='unformatted')   
  do i=1,nt
    write(104) dir1(nt),dir2(nt)
  enddo
close(104)
end

！！！！！！！！！！！！！！！！！！！！！ctl文件
dset D:\fengkuoxian\dire.dat
title direction
undef 9999
xdef 25 linear 0 1
ydef 10 levels 75 140 290 540 695 890 1005 1145 1330 1590
zdef 0 levels 0
tdef 25 linear 00Z01jan2010 1hr
vars 1
dir  10 99 direction  
ENDVARS

dset  D:\fengkuoxian\wind.dat
title wind speed
undef 9999
xdef 25 linear 0 1
ydef 10 levels 75 140 290 540 695 890 1005 1145 1330 1590
zdef 0 levels 0
tdef 25 linear 00Z01jan2010 1hr
vars 1
wnd 10 99  wind
ENDVARS

！！！！！！！！！！！！！！！！！！

画图的时候d wnd;dir 出现如图所示问题，请神们帮我看看啊

mofangbao

你该好好学习下GrADS官方文档~还有可以论坛搜索学习一下Fortran与GrADS的一些常见使用方式，不能按自己以为的方式去写东西

schLltt

mofangbao 发表于 2015-3-4 12:59
你该好好学习下GrADS官方文档~还有可以论坛搜索学习一下Fortran与GrADS的一些常见使用方式，不能按自己以为 ...

看文档只能理解到这样了。。能不能再给点建议呢

river

没懂，你的资料你不说谁知道是什么样格式的，站点，还是格点？从fortran程序里只能看出来是个一维的，只有时间维度，没有层次和格点的变化。。。但是从你的ctl一看，不仅有时间维，X和Y也都有变化。。。根本是对不上的东西。所以清风说的是没问题的，你缺的是grads的最基本的知识。如果像你说的“看文档只能理解到这样了”，那我只能说你根本就没仔细认真的看。不然绝对不会只有这种程度······

另外再给你一个地址，可以去学习参考一下http://bbs.06climate.com/forum.php?mod=viewthread&tid=30220&highlight=%B7%E7%C0%AA%CF%DF

schLltt

river 发表于 2015-3-4 16:26
没懂，你的资料你不说谁知道是什么样格式的，站点，还是格点？从fortran程序里只能看出来是个一维的，只有 ...

格点数据，两列十行，需要画的是横坐标为时间，纵坐标为高度的风廓线，程序是在将文本文件转换为二进制文件，那就是说ctl写的不对？

schLltt

river 发表于 2015-3-4 16:26
没懂，你的资料你不说谁知道是什么样格式的，站点，还是格点？从fortran程序里只能看出来是个一维的，只有 ...

就是根据那个帖子改的。。