资源描述:
《grads常用命令.doc》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、第一个命令:这样将得到一个动画显示的序列,现在输入: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现在已经知道怎样使用选定的部分数据绘图了,下面我们学习怎样操作数据。设置维数
2、环境为Z,Y变化:clear set lon -180 0 set lat 0 90 set lev 500 set t 1现在假设我们想显示华氏温度代替开尔文温度。可以这样转换:display (t-273.16)*9/532任何由标准的运算符,-,*和/组成的表达式都是允许的,并且操作数可以包含常量,变量,或者函数。一个包含函数的例子:clear d sqrt(u*uv*v)有一个函数用来计算风的级数。d mag(u,v)另一个内建函数计算平均值:cleardave(a,t=1,t=5)这种情况我们可以计算5天的平均。我们也可以从数据中移除平均值(距平值):dz-ave(z,t=
3、1,t=5)也可以在x方向作平均并求距平:clear d z-ave(z,x=1,x=72)也可以做时间差分:clear d z(t=2)-z(t=1)完整规范的变量名是:name.file(dim
4、-
5、=value,…)如果我们打开了两个文件,也许一个是模式输出,另一个是分析,我们应该区分用如下方法二者:display z.2-z.1另一个内置的函数通过有线差分计算水平涡度相关cleard hcurl(u,v)还有另外一个计算数值方向的质量积分:clear d vint(ps,q,275)这儿我们计算了可降水量(单位mm)现在来讨论控制图形输出的话题。目前为止,我们一直让GrADS选择默认
6、的画线间隔。我们可以这样自己设定: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使用注解:drawtitle500mbHeightsandVorticity风矢量:clear set gxout vec
7、tor 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
8、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文档并且尝试此处描述的其它选项
9、。英文原文并没有完全按原文翻译。reinit 初始化、关闭文件q/queryfile 查询文件,xyz格点数,时间,变量,(变量层<小于z>,)qctlinfo路径名字缺色值xyz格点数,初始值,stepsopen 打开文件d/display(temp) 绘图输出c/clear清除drawtitle(name)给图添加名字,如果加上反斜杠,名字分两行保存文件enableprint