grads 对数横坐标[已解决]

xman

       各位大侠，云物理中冰晶粒子数浓度经常高达10^5次方量级以上，因此当横坐标为数浓度，纵坐标为海拔高度时，作图很难看。
用别的工具可以将横坐标做成对数的，但是我看GRADS里面，只有纵坐标有对数表示  'set zlog on'，不知道是否也有对横坐标如此设置的命令。
       还望大家赐教

————————————————————————————————————————————————————
升级到“小雨”后，下载了清风大神的axis.gs，根据其小刻度绘图方法，设计出了对数横坐标的绘图程序。

function DrawSAxislog10(lastx,lasty,nx,ny,size,sinterv,posi,type,color,cthick,hoffset,voffset,nowlbl,suff)
sidx = 1
while(sidx < 9)
  if(lastx=nx)
    nowx=lastx
    dy=(ny-lasty)*math_log10(sidx+1.0)
    nowy=lasty+dy
  else
    nowy=lasty
    dx=(nx-lastx)*math_log10(sidx+1.0)
    nowx=lastx+dx
  endif
  DrawAxis(posi,type,nowx,nowy,size/2.0,color,cthick/2.0,hoffset,voffset,nowlbl,suff,'off')
  sidx=sidx+1
endwhile
return

基本思路依然是，找到小刻度的起始[x，y]坐标，然后画竖线。当然，文字标示也要相应修改：  nowlbl='10`a'now'`n'
下面给出粒子数浓度随高度变化的对数坐标图片。

GRADS中原有的命令 ‘set zlog on', 仅针对ctl文件指定的气压或者高度坐标有效，对于其他物理量无效。但是上面的办法完美的解决了这个问题，

xman

没有大侠知道吗？

553458324@

有很多有用的东西，来迟了

xman

已解决，庆祝一下

river

感谢分享，很厉害

desperado

先顶一个升级

shavian

先学习一下

webpie

感谢分享。还是毛毛雨，好多无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助）呢。不过没关系，继续学习中

-逆光年年-

先学习下，谢谢！

-逆光年年-

你好！我是grads菜鸟一只，对数坐标的问题一直困扰我，想请问下，您这个脚本是的加入清风的脚本里一起调用吗？能看下您画图的完整脚本吗？万分感谢，真的不知道怎么用，尝试了好久都不对