一些常用的grads命令与ncl命令的对应关系（续）

游子

4.grads中的图形类型设置与ncl的图形接口csm函数对应

说明	grads图形类型设置'set gxout grahics-style	ncl图形类型设置gsn_csm_choice[_choise[_choice[……]]]
一般等值线图	set gxout contour	gsn_csm_contour(无地图底图）gsn_csm_map_ce,之前都要设置cnLinesOn=True
等值线 + robinson投影 （cylindrical equidistant map）	set gxout contour set mproj robinson	gsn_csm_contour_map_ce
等值线+极射赤面投影 （Polar Stereographic Projections）	set gxout contour set mproj nps/sps	gsn_csm_contour_map_polar(之前设置： gsnPolar = "NH"/gsnPolar = "SH")
填色等值线图	set gxout shaded	gsn_csm_contour,之前设置 gsn_define_colormap(wks,colors),cnFillOn=True
填色等值线+等值线	set gxout shaded d .. set gxout contour d .. 注意：顺序不能反	同上，之前设置cnFillOn = True, cnLinesOn=True
仅仅标记网格点	set gxout grid	比较灵活，可以使用gsn_add_test做到，没有对应的csm接口，参照http://www.ncl.ucar.edu/Applications/Scripts/station_3.ncl
矢量图	set gxout vector	gsn_vector
流线图	set gxout streamine	gsn_streamline
一维数据图	set gxout bar/line/linefill	gsn_xy与各种属性联合使用达到目的，其中line是最简单的。
单线line图	set gxout line	gsn_xy
直方图	set gxout bar	gsn_xy,之前设置gsnScale = True（效果见图1），还可以使用gsn_histogram来统计一组数据的分布情况，分别落在什么区间，也就是说x轴是数值，y轴是发生频率或次数（有点概率密度的意思）(效果见图2）
两单线着色图	set gxout linefill（效果见图3）	gsn_xy，之前设置属性（效果见图4）
二维场方格子填充图	set gxout fgrid	同等值线图画法，设置相应的属性
误差盒子图	set gxout errbar	gsn_xy与相应属性联合使用达到目的

（欢迎大家补充，精力有限）





图1 ncl使用gsn_xy绘制直方图效果
-------------------------------------------------------------------------------------------------------------------------------------------------------

图2 ncl使用gsn_histogram绘制“概率密度”（其实是总个数*概率密度）的效果
------------------------------------------------------------------------------------------------------------------------------------------------

图3

图3 grads绘制填色单线图效果
----------------------------------------------------------------------------------------------------------------------------------------------------------------

图4 ncl使用gsn_xy绘制填色line图效果
---------------------------------------------------------------------------------------------------------------------------------------------
5.grads与ncl的函数对应关系

说明	grads函数	ncl函数
平均函数	ave（var,t/lon/lat=?,t/lon/lat=?,inc,flag）,求平均	avg, dim_avg, dim_avg_n, dim_avg_n_Wrap, dim_avg_Wrap
纬向平均	ave(var,x=x1,x=x2）	同上
时间平均	ave(var,t=t1,t=t2,dt)	同上，外加clmMonLLLT, clmMonLLT, clmMonTLL clmMonTLLL
高度权重平均	vint(ps,exp,top)	dim_avg_wgt, dim_avg_wgt_n, dim_avg_wgt_n_Wrap dim_avg_wgt_Wrap
区域面积平均	aave(u,x1=x,x=x2,y=y1,y=y2)	同上，注意一下参数的设定技巧,还可以直接使用： wgt_areaave, wgt_areaave2, wgt_areaave_Wrap
垂直涡度	hcurl(u,v)	uv2vr_cfd, uv2vrF, uv2vrf, uv2vrF_Wrap uv2vrG, uv2vrg, uv2vrG_Wrap
水平散度	hdivg(u,v)	uv2dv_cfd, uv2dvF, uv2dvf, uv2dvF_Wrap uv2dvG, uv2dvg, uv2dvG_Wrap
取样稀疏化	skip(u,dx,dy)	没有函数，可以直接引用u(...,::dx,::dy)
差分函数	cdiff(u),为中央差分函数	center_finite_diff, center_finite_diff_n
九点滑动平均	smth9(u)	smath9, smath9_Wrap
三角函数	sin,cos,tan	sin,cos,tan,此外还有反三角函数atan,acos,asin不记得 Pi=3.141592653的人一般这样表示之，pi=atan(1.0)*4， 经常可以在ncl脚本中看得到
对数指数函数	exp, log, log10	exp, log
求和	sumg	sum , dim_sum, dim_sum_n, dim_sum_n_Wrap dim_sum__Wrap
权重求和	sum	与权重平均一致，avg对应sum就好了

官网上各种函数，大家可以使用关键词搜索，http://www.ncl.ucar.edu/Document/Functions/list_alpha.shtml
分享总是快乐的

太多了，我看我的题目写的太大了。
有时间再补充，写一写grads的精细绘图与ncl属性的对应关系。

上一篇：http://bbs.06climate.com/forum.php?mod=viewthread&tid=10541
下一篇写了四种主流气象绘图软件的绘图命令对应一览表，欢迎大家讨论、交流。
   http://bbs.06climate.com/forum.php?mod=viewthread&tid=11453

善人/jw

顶起顶起

Aires

楼主好强悍！我的水平只限于看懂grads……

游子

Aires 发表于 2012-10-20 20:08

                               
登录/注册后可看大图

楼主好强悍！我的水平只限于看懂grads……

实在是太大的题目了，够我写一段时间的

Aires

游子 发表于 2012-10-20 20:12

                               
登录/注册后可看大图

实在是太大的题目了，够我写一段时间的

楼主慢慢写，谢谢分享好东西~

随缘

楼主辛苦了！致敬！

游子

随缘 发表于 2012-10-20 20:35

                               
登录/注册后可看大图

楼主辛苦了！致敬！

多谢捧场啊，写的多错的也可能越多

zysdg

Grads用了多年，也想学ncl来着，有了这个，应该会方便很多，谢谢楼主啊！

7788

给力！！！！！！！
补充下几个空白：
涡度：        vr = uv2vr_cfd（无缺测）
散度：        div = uv2dv_cfd（无缺测）
区域平均： wgt_areaave
纬向，经向，高度平均等，均可用dim_avg_n（wrap），n设定相应的维即可。

ayzqs

厉害呀，学习了