GrADS第4章 各类参数设置

站点数据变量或者表达式
x,y是窗口中的坐标位置
得到的结果：窗口中坐标为（x,y）最近站 点的站点号；经纬度；该站点物理量值 选择显示站点代码 例如：'set stid on' 'set gxout findstn' 'd rain;7.5;3.3'

'reinit' 'open e:\bs\rain\rain.qht7.ctl' 'set lon 73 136' 'set lat 17 54'



set lon -90 或set x 1 set lat 40 或set y 22 set lev 500 或set z 6 set time 01jan1985 或set t 5 set x 1 144 或 set lon 0 360 set y 1 73 或 set lat -90 90 set z 1 3 或 set lev 1000 850 set t 1 3 或 set time 01jan1985 03jan1985
vector: 矢量箭头形式绘二维风场（缺省设置） set gxout vector d u;v （u,v在同一个数据文件中可这样） 或 d u.1;v.2 （u,v在两个数据文件中） stream: 流线形式绘制二维风场 set gxout stream d u;v barb: 风向杆形式绘制二维风场 set gxout barb d u;v

linefill: 两单曲线之间填色 与命令set lfcols col1 col2合用。 用法：set gxout linefill

errbar: 单线图误差分布 set gxout errbar d line1;line2
'reinit' 'open D:\grads\data\u.ctl' 'enable print D:\grads\data\2.gmf' 'set grid off' 'set parea 1 10 1 8' 'set grads off' 'set xlopts 1 2 0.15' 'set ylopts 1 2 0.15' 'set t 1 48' 'set vrange -40 60'


fwrite: 将输出结果(数据)存入一个由‚set fwrite 文件名‛所指定的文件中
‘open c:\pcgrads\model\model.le.ctl’ ‘set gxout fwrite’ （设置输出类型） ‘set fwrite d:\result\u.grd’ it=1 while(it<=5) ‘set lev 500 ’ ‘set lon 120 270’ ‘set lat –10 10’ （设置维数环境） ‘set t ‘it’’ ‘d u ’ it=it+1 endwhile ‘disable fwrite’ （关闭文件u.grd）


symbol — 欲绘天气符号的代码 x — x 坐标（英寸） y — y 坐标 size — 天气符号的大小 color — 天气符号的颜色，-1为标准色(红色代表热带风 暴, 蓝色代表降雪等等) thickness — 天气符号的线宽
说明：


当图形输出类型为vector、stream和barb时, 在显示命令后需给出两种数据，两者间用 分号‚ ; ”隔开，前者理解为U分量，后者 理解为V分量 当图形输出类型为errbar、scatter、linefill 时，也须用分号将两种数据分开 例如： display u ；v 显示u、v 风场的合成矢量图

'set gxout line' 'set cstyle 1' 'set cmark 0' 'set ccolor 2' set lfcols col1 col2 'd v.4' 'set gxout line' d line1;line2 当line1<line2时，用颜色1（col1） 'set cstyle 1' 'set cmark 0' 当line1>line2时，用颜色2（col2） 'set ccolor 1' 'd v.5 ' 'set gxout linefill' 'set lfcols 4 10' 'd v.4; v.5'
'set stat off' 'reinit'
'set gxout stnmark 站点位置分布 'set stid on' 'set gxout stnmark' 'd rain'

对于wxsym类型的图形, 站点坐标的每一个数值都 假定为一种wx天气符号的代码，台风和热带风暴 的符号包含在wx天气符号集中。 draw wxsym symbol x y size <color <thickness>> （在一指定的坐标绘一指定的天气符号） 其中：
'set gxout line' 'set cmark 2' 'd u' 'set gxout errbar'
'set ccolor 2'
'd u-12.8;u+12.8'
'print' 'c' 'disable print' 'reinit' ;


scatter 散点图 'reinit'
set gxout scatter 'open d:\lon.ctl' 'open d:\lat.ctl‘ d v.1;v.2 'enable print D:\grads\data\1.gmf' 'set x 1 ';'set y 1';'set z 1 '; 'set t 1 24' 'set gxout scatter ' 'd v.1;v.2 ' 'print' 'c' 'disable print' 'reinit' ;
注意: grads只能显示两维或一维不固定的图！

当所有维数都固定时，得到的是一个单 值数据点

维数环境的设置一直保持到再次设定时 都有效。要查看当前维数情况，可以在 运行环境中输入命令： ga-> q dims
二维平面图

set lon 180 360 set lat 0 90 set lev 500 set t 2 d hgt set lon 115 set lat 15 set lev 850 set t 1 12 du
不固定 不固定

不固定
一维曲线图
4.2 图形类型的设置
当维数环境确定后，缺省情况下， 一维变量输出的图形为单线图，二 维变量为等值线图，改变缺省图形 输出类型的命令为： set gxout graphics_type
1. 格点数据
contour: 二维等值线图（缺省设置） set gxout contour shaded: 二维填色图(通常与contour并用) 先画填色图，再画等值线图；否则等值线会 被覆盖。 set gxout shaded d hgt set gxout contour d hgt

set lon 0 210 set lat 0 70 set lev 500 set t 2 set gxout contour d hgt/10
set t 1 set gxout shaded set clevs -1 1 set ccols 10 0 4 d (hgt(t=1)-ave(hgt,t=1,t=50) )/10 set gxout contour d (hgt(t=1)-ave(hgt,t=1,t=50))/10


注意： set gxout fwrite是一种特殊的输出格式，它 将向图形输出的数据转向一个二进制文件 输出。 由于地图坐标转换向格点坐标转换过程中 的误差，需要仔细的确定输出数据的维数。 确定输出数据的维数后，可利用Fortran程 序进行读取。
2. 站点数据
value:
在各站点标值（缺省设置） barb: 在各站点绘风向标（缺省设置） findstn: 搜索最近的站点 stnmark:站点位置分布 model: 以天气图形式将天气观测各 分量填放在站点四周 wxsym: 绘wx天气符号

grid: 二维场不绘图，以网格形式在各网格点中央标 出该点数值

fgrid: 用指定颜色填充二维格点场，对二维场不 绘制等值线图，只是将特定值的格点用指定的颜 色填充该网格。与命令set fgvals val col合用 用法：set gxout fgrid
set fgvals val1 col1
第4章来自百度文库各类参数设置
set是一项非常重要的命令，用于设置 各种运行环境的参数，包括维数环境、 图形类型、图形要素、屏幕显示等等。
4.1 维数环境的设置
虽然在数据描述文件中给出了各物理变 量数组的时空维数范围，但在GrADS 运 行环境中还需设定全数据集中参与操作 的部分或全部数据集的维数情况，以供 以后的表达式、显示命令等使用。这就 是维数环境的设置。
set fgvals val2 col2
set lon 70 130 set lat 10 30 set lev 1000 set t 1 set gxout grid du
set lon 70 130 set lat 10 30 set lev 1000 set t 1 set gxout fgrid set fgvals -1 4 Set fgvals -3 5 du

set gxout value 各站点标值
'reinit' 'open e:\bs\rain\rain.qht7.ctl' 'enable print D:\grads\data\rain.41.gmf' 'set grid off' 'set xlint 5' ‘set mpdset mres' 'set lat 15 55' 'set lon 100 135' 'set t 1'

line: 对一维场绘制单线图 set lon 120 set lat 20 set lev 200 set t 1 12 set gxout line du

bar: 对一维场不绘单线图，而绘制直方图 set lon 120 set lat 20 set lev 200 set t 1 12 set gxout bar set bargap 40 Set barbase 0 dv
3. 维数环境的定义
（1）一种是地球坐标（world coordinate）
以经纬度为度量单位。 形式如下： set lon|lat|lev|time val1 <val2> （2）一种是格点坐标（gridcoordinate） 以网格点数为度量单位。 形式如下： set x | y | z | t val1 <val2>
'set gxout value' 'd rain'
'run e:\bs\rain\test.gs' 'run e:\bs\rain\southsea.gs' 'print' 'c' 'disable' 'reinit‘ ;
set gxout findstn 搜索最近的站点 d rain;x;y
1. 维数环境的概念
GrADS系统视每一个物理变量（VAR） 场为一个四维数据集（4D data set）， 即包括三维空间（x，y，z）和一维时间 （t）。也可固定其中的一维或几维以获 得实际的低于四维的数据子集。
2.设置的作用
GrADS 中设置维数用以说明或指定 随后的分析或图形操作时参加操作的 原始数据集的维数范围。该工作数据 集可以是整个原始数据场，也可以是 原始数据场的一部分。
将findstn的结果赋 值到id，x,y,value
写出结果
'set stid on' 'set gxout findstn' 'd rain;7.5;3.3' id=subwrd(result,1) x=subwrd(result,2) y=subwrd(result,3) value=subwrd(result,4) say '========================' say 'stid: ' id say 'value: ' value say 'lon: ' x say 'lat: ' y