IDL垂直剖面图纵坐标为气压的间距问题

莎朗嘿哟YY

我用IDL画了一个随气压变化的垂直剖面图，如下：
我现在的问题是纵坐标的问题，一般越往上100hpa之间的距离比低层的要大，就是下图的样子
求问在IDL里面怎么实现呢？？？谢谢！

chenhongfengye

我这个是对数变化的，我觉得你可以对数分之一变化就可以吧

莎朗嘿哟YY

chenhongfengye 发表于 2014-5-25 14:57
我这个是对数变化的，我觉得你可以对数分之一变化就可以吧

能告诉我具体的命令是什么？刚接触IDL，有的命令不是太清楚。。谢谢！

chenhongfengye

   axis,yaxis=0,/ylog,yticks=15,yticklen=1,$
     ytickv=[10,20,40,60,80,100,200,400,600,800,1000,2000,4000,6000,8000,10000],$

莎朗嘿哟YY

chenhongfengye 发表于 2014-5-25 15:12
axis,yaxis=0,/ylog,yticks=15,yticklen=1,$
     ytickv=[10,20,40,60,80,100,200,400,600,800,1000,2 ...

恩，好的，谢谢，我试试

莎朗嘿哟YY

chenhongfengye 发表于 2014-5-25 15:12
axis,yaxis=0,/ylog,yticks=15,yticklen=1,$
     ytickv=[10,20,40,60,80,100,200,400,600,800,1000,2 ...

我想问一下ytickv=[10,20,40,60,80,100,200,400,600,800,1000,2000,4000,6000,8000,10000]是什么意思？运行出来报错

莎朗嘿哟YY

chenhongfengye 发表于 2014-5-25 15:12
axis,yaxis=0,/ylog,yticks=15,yticklen=1,$
     ytickv=[10,20,40,60,80,100,200,400,600,800,1000,2 ...

我知道了，你那个值是你纵坐标的值，我改成我的后，还是报错，是什么原因呢？

莎朗嘿哟YY

chenhongfengye 发表于 2014-5-25 15:12
axis,yaxis=0,/ylog,yticks=15,yticklen=1,$
     ytickv=[10,20,40,60,80,100,200,400,600,800,1000,2 ...

报错说： ytickv=[1000,900,800,700,600,500,400,300,200,100]   
                                                           ^
% CONTOUR: Incorrect number of arguments.

莎朗嘿哟YY

chenhongfengye 发表于 2014-5-25 15:12
axis,yaxis=0,/ylog,yticks=15,yticklen=1,$
     ytickv=[10,20,40,60,80,100,200,400,600,800,1000,2 ...

前辈，你看到我的回复了吗？求帮忙解决一下~~万分感谢！

spswarm

做了个例子，你看看有没有帮助。
;＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
;54511  116.47   39.80   55  156
;994    9999      23      22     135       2
;925      66      20      18     130       3
;850     138      16      15     240       1
;708    9999       9       5    9999    9999
;700     301       9       5     310      11
;572    9999      -3      -5    9999    9999
;526    9999      -6     -14    9999    9999
;516    9999      -6     -48    9999    9999
;500     571      -6     -48     280      14
;498    9999      -6     -48    9999    9999
;477    9999      -9     -36    9999    9999
;462    9999     -11     -15    9999    9999
;419    9999     -15     -17    9999    9999
;407    9999     -15     -19    9999    9999
;400     742     -16     -23     255      26
;382    9999     -19     -27    9999    9999
;328    9999     -27     -29    9999    9999
;319    9999     -23     -49    9999    9999
;300     954     -26     -61     245      40
;250    1083     -36     -68     245      39
;200    1234     -47     -76     255      41
;197    9999     -49     -78    9999    9999
;150    1420     -58     -84     255      30
;121    9999     -65    9999    9999    9999
;105    9999     -66    9999    9999    9999
;100    1669     -66    9999     275      15
;

pro draw_a_t_log_P

t = [23$
, 20$
, 16$
, 9$
, 9$
, -3$
, -6$
, -6$
, -6$
, -6$
, -9$
, -11$
, -15$
  ,-15$
  ,-16$
  ,-19$
  ,-27$
, -23$
, -26$
, -36$
, -47$
, -49$
, -58$
, -65$
, -66$
, -66$
]
p = [994$
, 925$
,850$
,708$
,700$
,572$
,526$
,516$
,500$
,498$
,477$
,462$
,419$
,407$
,400$
,382$
,328$
,319$
,300$
,250$
,200$
,197$
,150$
,121$
,105 $
,100]
h = press2alt(p)
help, p, t, h
plot, t, h, ytickv=press2alt([1000, 900, 850, 700, 600, 500, 400, 300, 200, 100]), yticks=10, ytickname=[1000, 900, 850, 700, 600, 500, 400, 300, 200, 100]
end

;＝＝＝＝＝＝＝＝＝＝＝
;那个　press2alt.pro网下载一下就行。

ytickvalues : An array of tick mark locations.
ytickname: A string array containing the tick labels.

我猜，那个气压的标注应该是在对应的海拔高度上，所以用press2alt转换了一下。