grads 对数横坐标[已解决]

xman · 发表于 2015-7-5 20:14:36

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 xman 于 2015-7-7 19:33 编辑

   各位大侠，云物理中冰晶粒子数浓度经常高达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 · 发表于 2015-7-6 09:49:42

没有大侠知道吗？

553458324@ · 发表于 2015-7-6 12:53:12

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

xman · 发表于 2015-7-7 19:36:55

已解决，庆祝一下

river · 发表于 2015-7-7 20:20:25

感谢分享，很厉害

desperado · 发表于 2015-12-21 21:33:39

先顶一个升级

shavian · 发表于 2017-10-11 11:32:43

先学习一下

webpie · 发表于 2017-10-14 10:52:37

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

-逆光年年- · 发表于 2017-12-3 15:33:01

先学习下，谢谢！

-逆光年年- · 发表于 2017-12-3 17:27:19

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

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

[图形美化] grads 对数横坐标[已解决]

登录后查看更多精彩内容~

评分

相关帖子