任意垂直剖面的地形画法

黎大页子 · 发表于 2014-4-27 04:06:15

pandaxuwj 发表于 2014-4-23 13:50
顶一个，不过我想问一下楼主，我用wrf输出数据的垂直方向是等压面坐标，怎么画出类似图形？

把HGT*0.001换成psfc*0.01

  ln=lon1
  i=1
  while(ln<=lon2)
lt=lat1+(lat2-lat1)*(ln-lon1)/(lon2-lon1)
'set lon 'ln
'set lat 'lt
'set z 1'
'd psfc*0.01'
res1=sublin(result,2)
top=subwrd(res1,4)

'q w2xy 'ln' 1000'
res2=sublin(result,1)
tx=subwrd(res2,3)
ty=subwrd(res2,6)

if(top<1000)
'q w2xy 'ln' 'top
res3=sublin(result,1)
mt=subwrd(res3,6)

   'set line 1 1 6'
   'draw line 'tx' 'ty' 'tx' 'mt
endif

ln=ln+0.002
i=i+1
  endwhile

pandaxuwj · 发表于 2014-4-27 09:16:36

黎大页子发表于 2014-4-27 04:06
把HGT*0.001换成psfc*0.01

ln=lon1

谢谢楼主的解答

人走茶不凉 · 发表于 2014-4-27 11:20:45

顶，收藏，有空儿的时候好好看看

唐卡卡55 · 发表于 2014-4-27 13:17:41

谢谢分享谢谢分享谢谢分享

707013105 · 发表于 2014-4-27 13:48:39

好像论文要用到，先收藏了

更大华 · 发表于 2014-4-29 02:08:40

感谢分享。中国的科研就是要共享，才能共同提高

手心里的阳光 · 发表于 2014-4-30 02:30:11

看了楼主的帖子运用了下，但是发现了问题，问题提示：unable to locate endwhile statement for the while statement ...不知道错误在哪。。都copy了楼主的gs，尝试将lat和lon的值带入lt中，也不能解决。。希望得到楼主的帮助~~

pandaxuwj · 发表于 2014-4-30 12:44:36

手心里的阳光发表于 2014-4-30 02:30
看了楼主的帖子运用了下，但是发现了问题，问题提示：unable to locate endwhile statement for the while ...

那个问题是说找不到你的endwhile，你查看一下循环结尾，如果没有endwhile就补上，如果本来有的话，可能是endwhile是最后一行没法识别，试试敲个回车，在下面加上一行

黎大页子 · 发表于 2014-4-30 13:08:04

手心里的阳光发表于 2014-4-30 02:30
看了楼主的帖子运用了下，但是发现了问题，问题提示：unable to locate endwhile statement for the while ...

解决了吗？

dream · 发表于 2014-5-2 18:56:33

有个OROG文件是画地形的，那个是个什么东西

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

[图形美化] 任意垂直剖面的地形画法