grads常用命令.doc

第一个命令： 这样将得到一个动画显示的序列，现在输入： claer??set?lon?-90?set?lat?-90?90?set?lev?1000?100?set?t?1??d?t??d?u 该情况下我们设置y（纬度）和z（高度）变化，所以我们得到一个竖直切面图。我们显示了两个变量，他们重叠在一起。你可以让任意 多的变量重叠在一起，只要你不使用clear命令就行了。另一个例子，这种情况X和T是变化的（哈莫图）： c??set?lon?-180?0??set?lat?40??set?lev?500??set?t?1?5??d?z 现在已经知道怎样使用选定的部分数据绘图了，下面我们学习怎样操作数据。设置维数环境为Z,Y变化： clear??set?lon?-180?0??set?lat?0?90??set?lev?500??set?t?1 现在假设我们想显示华氏温度代替开尔文温度。可以这样转换： display?(t-273.16)*9/5 32 任何由标准的运算符 ,-,* 和 / 组成的表达式都是允许的，并且操作数可以包含常量，变量，或者函数。一个包含函数的例子： clear??d?sqrt(u*u v*v) 有一个函数用来计算风的级数。 d?mag(u,v) 另一个内建函数计算平均值： clear d ave(a,t=1,t=5)这种情况我们可以计算5天的平均。我们也可以从数据中移除平均值（距平值）：d z-ave(z,t=1,t=5)也可以在x方向作平均并求距平： clear?d?z-ave(z,x=1,x=72) 也可以做时间差分： clear?d?z(t=2)-z(t=1) 完整规范的变量名是： name.file(dim |-|=value,…) 如果我们打开了两个文件，也许一个是模式输出，另一个是分析，我们应该区分用如下方法二者： display?z.2-z.1 另一个内置的函数通过有线差分计算水平涡度相关 cleard?hcurl(u,v) 还有另外一个计算数值方向的质量积分： clear??d?vint(ps,q,275) 这儿我们计算了可降水量（单位mm） 现在来讨论控制图形输出的话题。目前为止，我们一直让GrADS选择默认的画线间隔。我们可以这样自己设定： clear??set?cint?30d?z 也可以控制线的颜色：clear??set?ccolor?3? d?z 可以使用交替的方式显示数据： clear??set?gxout?shaded??d?hcurl(u,v) 这看起来不很平滑，可以使用三次平滑： clear??set?csmooth?on??d?hcurl(u,v) 我们也可以使用不同的图形样式 set?gxout?contour??set?ccolor?0??set?cint?30??d?z 使用注解：draw title 500mb Heights and Vorticity 风矢量： clear??set?gxout?vector??d?u,v 这里我们显示两个表达式，第一个是U分量矢，第二个是V分量矢。我们也可以使用第三个字段使图像彩色化 d?u;v;q 或者： d?u;v;hcurl(u,v) 你也可以显示任何区域以显示伪矢量： clear??d?mag(u,v);q*10000 此处U分量是风速和V分量是湿度。 我们也可以查看流线（并且着色） clear??set?gxout?stream??d?u;v;hcurl(u,v) 或者我们可以查看实际的格点值 clear?set?gxout?grid??d?u 我们可能希望改变地图背景： clear??set?lon?-110?-70?set?lat?30?45??set?mpdset?nam??set?digsize?0.2??set?dignum?2??d?u 改变投影： set?lon?-140?-40??set?lat?15?80??set?mpvals?-120?-75?25?65??set?mpproj?nps??set?gxout?contour??set?cint?30?d?z 我们让grads可以操作纬度从140W到40W之间，经度从15N到80N之间的数据。但是我们告诉它使用从120W到75W和25N到65N极坐标投影绘图。多余的部分被投影程序修剪掉了。该简单教程到此结束。此时，你可能希望深入查看数据集，或者你可能希望浏览GrADS文档并且尝试此处描述的其它选项。 HYPERLINK /grads/gadoc/tutorial.html \t _