tlnp图的问题

shzljg

TLNP问题我有探空资料，想画TLNP图，利用fortran 读，然后写成二进制文件，然后写相应二进制文件的ctl文件，最后写gs文件，把plotskew.gs拷贝到gs文件后面，运行gs.但是画出来的图不对，不知道是什么原因。
以下为数据。
51431   81.33   43.95  663   66
     928    9999      17      14       0       0
     925      69      18      14       0       0
     850     143      20      12      80       3
     700     307       7       2     315       5
     500     575     -10     -20     295       8
     400     744     -19     -34      25       3
     300     950     -37     -47     255      17
     250    1074     -46     -54     225      25
     200    1222     -48     -57     220      17
     150    1412     -49     -58     235      20
     100    1672     -57     -66     230       9
fortran :program tankong
real head(5)
integer dat(11,6)
open(1,file='e:\storm\51431.txt')
open(2,file='e:\storm\51431.dat',form='binary',status='new')
open(3,file='e:\storm\51431_check.txt',status='new')
read(1,*) (head(i),i=1,5)
read(1,*) ((dat(i,j),j=1,6),i=1,11)
write(2) (dat(i,3),i=1,11)
write(3,"(i4)") (dat(i,3),i=1,11)
write(2) (dat(i,4),i=1,11)
write(3,"(i4)") (dat(i,4),i=1,11)
write(2) (dat(i,6),i=1,11)
  write(3,"(i4)") (dat(i,6),i=1,11)
write(2) (dat(i,5),i=1,11)
write(3,"(i4)") (dat(i,5),i=1,11)
  close(1)
  close(2)
  close(3)
end
ctl 文件：DSET e:\storm\51431.dat
UNDEF -999.9
XDEF 1 LINEAR 81.33 0.1
YDEF 1 LINEAR 43.95 0.1
ZDEF 11 LEVELS 928 925 850 700 500 400 300 250 200 150 100
TDEF 1 LINEAR 08Z15jul2007 1hr
VARS 4
temp 11 99 temp
dew 11 99 dew
speed 11 99  wind direct
dir 11 99  wind velocity
ENDVARS
gs: 'reinit'
'open e:\storm\51431.ctl'
'set x 1'
'set y 1'
'set t 1'
'set z 1 11'
rc=plotskew(temp,dew,speed,dir)

topmad

笨笨原来做过相关的教程 可以看看他的帖子 他把代码啥的都贴出来了

shzljg

能给详细的网页地址吗？没找到

Miss肿么办

同求画Tlnp图的方法，有数据又不想手绘扫描。。。

kellyding_100

tlogp图右上角方框里的时间为制作该图的时间(电脑默认)，非探空时间 我才知道

hdswy

楼主后来怎么解决的呢，同问

Rplay

请问这位大神，你的问题后来解决来吗，如何解决的？小弟遇到了相似的问题，不知道哪里出错，还需请教大神。

流觞曲水

你这是micaps第五类数据格式，可以直接导入到micap里画，图也更漂亮：

tlogp图右上角方框里的时间为制作该图的时间(电脑默认)，非探空时间

kenever

流觞曲水 发表于 2013-5-19 07:26
你这是micaps第五类数据格式，可以直接导入到micap里画，图也更漂亮：

tlogp图右上角方框里的时间为制作 ...

micaps为什么显示没法直接读取第五类数据 要有高空填图数据才可以用 可是我又没有高空填图数据怎么办

ckyzydys

topmad 发表于 2012-6-26 08:14
笨笨原来做过相关的教程 可以看看他的帖子 他把代码啥的都贴出来了

请问老师，您提到的笨笨的帖子叫什么呀，谢谢