新GrADS使用笔记

新GrADS使用笔记
由wrfchina 于星期六, 2012-03-17 22:27 提交
求导———
define dv = cdiff(v,x)
define dx = cdiff(lon,x)*3.1416/180
define du = cdiff(u*cos(lat*3.1416/180),y)
define dy = cdiff(lat,y)*3.1416/180
display (dv/dx-du/dy)/(6.37e6*cos(lat*3.1416/180))
矢量图
之————控制矢量图中矢量标志的位置
'set gxout vector'
'set arrlab off' ＊关掉缺省的图标
'set arrscl 0.5 0.5' ＊表示在图中0.5英寸的长度代表0.5m/s 'd u;v'
'draw string 0.95 3.6 `3"'＊画箭头，使用set font 3的特殊图形'draw string 5.75 0.8 0.5'＊写箭头下文字
或
'set arrlab off' ＊关掉缺省的图标
'set arrscl 0.5 0.5' ＊表示在图中0.5英寸的长度代表0.5m/s 'd u;v'
'draw line 5.5 1 6 1' ＊画线
'draw line 5.95 0.95 6 1'＊画箭头
'draw line 5.95 1.05 6 1'
'draw string 5.75 0.8 0.5'＊写箭头下文字
之————规定矢量图输出的阈值，使用set cmin/cmax
但要注意的是，display时，必须是
'd u;v;mag(u,v)'
或
'd u;v;sqrt(u*u,v*v)'
而不能仅仅'd u;v'
之————控制矢量密度与颜色
'set ccolor 1'*控制颜色
'd skip(u,2);v'*控制密度
function
之————gs重复使用到的语句，打包成function
以function main(track)开头
以return结尾
函数
smth9对于资料分辨率过低引起的不平滑起作用，但如果对于模式输出等高分辨率的不平滑不会有任何效果
多幅图输出
之————多幅图输出时，边缘的经纬度为了避免重叠，必须使得第n张图的最后一个坐标不输出，
如set lon 132 144
set xlint 3
set xlopts color <thinkness<size>>
*因为144是（144-132）/3的公倍数，144必定会在最后一个经度坐标上出现，这样就会与第n 1个图的第一个横坐标：132重叠，解决方法有:
1.在最初的计算区域，就是分页前，就规定所选区域为[132,144)
2.分页前，依然选取[132,144]，在分页后，每幅图的display前，将区域选为[132,144) 经纬度对应图上位置——
之1.
lon = '133.5 '
lat = '33.8'
'q w2xy '%subwrd(lon,1)%' '%subwrd(lat,1)
x = subwrd(result,3)
y = subwrd(result,6)
之2.
lon='xxx.x xxx.x xxx.x'
lat='xx.x xx.x xx.x'
j=1
cn=999
while(j<cn)
'q w2xy '%subwrd(lon,i)%' '%subwrd(lat,i)
x = subwrd(result,3)
y = subwrd(result,6)
'draw mark 5 'x' 'y' 0.1'
j=j 1
endwhile
之3.
i=1
cnt=999
while(i<=cnt)
aa=read('.txt')
aa1=sublin(aa,2)
lat=subwrd(aa1,3)
lon=subwrd(aa1,4)
'q w2xy 'lon' 'lat
x = subwrd(result,3)
y = subwrd(result,6)
········
i = i 1
endwhile
ret=close('.txt')
后处理
之——色标
cbarc 在图的右上角给出了标尺
cbarn 在图的下方给出标尺
之——设置特殊字体与符号
用set font n设置字体
font 3时，可以输出很多symbols，
具体写法是:
'set font 3' *actually，this setting can be deleted，and it's strongly recommended not to write this command so that you can avoid unwanted character change
之字符：'draw string 9.09 3.03 `52900141' *then you can see bold roman numbers 2900141
之符号：'draw string xxx xxx `3#' then you can see a arrow,remember this the symbol must come together with the number "3"
关于画水汽通量，水汽通量为一个向量，方法如下：
set gxout shaded
d mag(q*u,q*v)
set gxout vector
d q*u;q*v （mag(q*u,q*v)只是画通量的绝对值，所以要加上d q*u;q*v才能看出通量的方向）
关于画水汽通量散度，应该为：
hdivg（q*u,q*v）
grads将数据写入不同类型文件：
之——存成ASCII码：
file='output.txt'
'set gxout print'
'd var'
rc=write(file,result)
rc=close(file)
之——存成binary：
'set gxout fwrite'
'set fwrite filename'
'd var'
'disable fwrite'。