[程序见GrADS Script Library]
1. basemap.gs:
basemap L | O | U <fill_color> <out_color> <hi/lo>
在低分辨率海岸廓线范围内用颜色覆盖陆地/海洋。适用于各种投影方式,需lpoly.asc, lpoly_hires.asc, lpoly_US.asc, opoly.asc, opoly_hires.asc文件。
其中:L(l):覆盖陆地,O(o):覆盖海洋,
U(u):覆盖20N-50N的墨西哥和加拿大领土(低分辨率,适用美国),
fill_color:填充色号,缺省为15,out_color:廓线的颜色号,缺省为15,
hi/lo:高分辨率('set mpdset hires',仅对15N-53N, 130W-60W区域)/低分辨率。
[
basemap.gs补充:
对高分辨率,新的
lpoly_hires.asc和
opoly_hires.asc文件可覆盖北美洲(美、加、格陵兰)极其沿岸。
ftp://grads.iges.org/grads/scripts/lpoly_hires.asc
ftp://grads.iges.org/grads/scripts/opoly_hires.asc
重要补充:
另外有一个类似basemap.gs的程序
makebg.gs可下载:
ftp://grads.iges.org/grads/scripts/makebg.gs
这个程序用到DODS资料,需高版本GrADS支持,另外还需联网使用,且因数据量原因,速度极慢,各位可根据自己的条件对其中的代码进行改进。
]
2. cbar.gs、cbarn.gs、cbarc.gs、cbar_l.gs和cbar_line.gs:
cbarn sf vert xmid ymid
cbarc center_x center_y back_color
绘制'set gxout shaded'图形的填色标尺。sf:色标尺寸,1为全尺寸,0.5为半尺寸;
vert:0为水平,1为垂直;xmid,ymid:色标中心点的位置。
cbar_l -x X -y Y -n number -t text -p
cbar_line -x X -y Y -c color -m mark -l linestyle -t text -p
加'set gxout line'的图例说明。其中:-x,-y:图中x和y的位置,-n:线条的数目(最多可为10条),-t:文字说明的内容(最多10条,需双引号括起),-c:线和标记的颜色,-m:定义标记; -l:定义线型,-p:用户可在图中点击给定图例的放置位置。
3. colors.gs:
为雪盖(颜色序号40~45)、降水(颜色序号50~59)及温度(颜色序号64~85)资料的shaded图设置填充色。
4. cmap.gs 和 colortab.gs:
run cmap: 生成新的grads.gct文件,
cmap filename.gct: 修改filename.gct文件.
rc = colortab ( filename ): 调用filename.gct文件的颜色设置。
Note: cmap以交互方式生成一个颜色列表。输入的颜色数量n需在1~84间。
5. showcol.gs:
showcol <palette> [<color_numbers>]
显示各调色板d的色卡,并生成palette.gmf和palette.gif文件。其中:
palette:调色板文件名,color_numbers:显示的颜色数目。
附:调用调色板时,直接在‘ga>’提示符下直接键入以下调色板的名称:
调色板名称 颜色号 说明
C16a、C16b、C16c: 16~29号 14阶'DKRZ historical'色,
C32a、C32b: 16~45号 32阶'DKRZ historical'色,
C256-1,-2,-3,-4: 16~79号 256色,每64种颜色一个文件,
Crainbow-1: 16~95号 80种五彩色(蓝~红),
Crainbow-2: 16~95号 80种五彩色(红~绿),
Crainbow-3: 16~95号 80种五彩色(绿~蓝),
BYR-02: 16~80号 蓝~黄~红(无绿色)65色,
BYR-03: 16~80号 蓝~黄~红(无绿色)65色,
BYG-01: 16~80号 蓝~黄~绿,65色,
GYR-01: 16~80号 绿~黄~红,65色,
RED-65: 16~80号 深红~浅橙,65色,
GREEN-65: 16~80号 深绿~浅绿,65色,
BLUE-65: 16~80号 深蓝~浅蓝,65色,
PalGrey: 16~99号 84阶灰度色。
6. 定义颜色的有关Script程序:
define_colors.gs、rgbset.gs、rbtab.gs、rgbmap.gs、showcol.gs、colors.gs
7. set font n 设置字体
n = 0:helvetica font(缺省)n = 1:roman font
n = 2:italics roman fontn = 3:symbols font
n = 4:bold helvetica fontn = 5:bold roman
特别推荐:font 3,可以输出很多symbols,让图例说明更完美,具体字符对照可用font.gs显示,用法: run font.gs n
8. GrADS 支持的投影方式 (I):
set mproj latlon|scaled|nps|sps|robinson|orthogr|mollweide|lambert|off
爱心小帖士:设置set mpvals lon1 lon2 lat1 lat2 可以帮你实现完美。
map.gs:帮你设定多种投影方式,不满足的自己加!用法:
run map.gs显示提示
run map.gs map_type直接设置,用熟了就这样。
GrADS 精致绘图说(四)
9. GrADS 支持的投影方式 (II)—<.ctl>文件中用pdef关键字 :
pdef isize jsize proj ...
如:
* pdef isize jsize nps|sps ipole jpole lonref gridinc
对Polar Stereo投影数据的说明。其中:
isize,jsize:X、Y方向的点数;ipole,jpole:极点的(i,j)位置;
lonref:参考经度; gridinc:以KM为单位的网格距。
* pdef isize jsize lcc latref lonref iref jref Slat Nlat standard_lon dx dy
对Lambert投影数据的说明。其中:
isize, jsize:X、Y方向的总点数; latref, lonref:参考点的经、纬度;
iref, jref:参考点的(i,j)位置; Slat, Nlat:Lambert投影m=1时南、北两个纬度;
standard_lon:标准经度; dx, dy:以M为单位的X和Y方向的网格距。
爱心小帖士:对lcc投影方式的数据,运用 set mproj lambert会有美妙的效果。
* pdef isize jsize eta.u latref lonref dlon dlat
对NMC Eta模式(unstaggered grids)数据的说明;
* pdef isize jsize pse slat slon polei polej dx dy sgn
对NMC high accuracy polar stereo for SSM/I data进行说明;
* pdef isize jsize ops latref lonref xref_offset yref_offset iref jref dx dy
对CSU RAMS Oblique Polar Stereo Grids数据进行说明;
* pdef isize jsize BILIN <STREAM> <BINARY> fname
<SEQUENTIAL> <BINARY-BIG>
<BINARY-LITTLE>
对除以上特殊投影方式外的数据的说明。
10. 为等值线图加注 H/L 中心标记 clhilo.gs、philo.gs和philo2.gs:
用法:
run clhilo.gs var <ptype <fmt <cint <rad>>>>
run philo.gs var <fmt <crad>>
run philo2.gs var <rad <cint <ptype <fmt>>>>
参数说明:
var:要素场;
ptype:maxmin—标记H/L,max—标记H,min—标记L;
fmt:fortran数据格式说明;
cint:等值线间隔,一般不需设;
rad:半径参量,可控制标记的多少;clhilo.gs和philo2.gs中rad以KM为单位,缺省值分别为500和1000;philo.gs中缺省为0.25。
这几个程序需要用到GrADS的"user defined functions",fortran执行程序名分别为:
clhilo.exe:clhilo.gs 和 philo2.gs 用
grhilo.exe:philo.gs用
在GrADS调用这两个用户自定义函数时,须先定义环境变量 GAUDFT 。具体做法如下:
UNIX/LINUX:setenv GAUDFT /usr/local/grads/udft
PC WINDOWS:GAUDFT=c:/pcgrads/udf/udft
实际路径因人而异。
注意:标记了彩色的"udft"是一个用户自定义函数列表文件。
设置是否成功,可启动GrADS后,用 'q udft'命令查看,系统将给出一个列表。
附:用户自定义函数说明:http://grads.iges.org/grads/gadoc/udf.html
11. Some math functions
rc = math_trigfunc(angle) (三角函数)
math_trigfunc可为:sin, cos, tan, asin, acos, atan, sinh, cosh, tanh, asinh, acosh, 或atanh。angle:需为弧度,rc:结果。
rc = math_format( format , num )
format:C语言形式的输出格式, num:需转换格式的数据, rc:转换后的数据。
【附】C语言输出格式说明符:
%d以带符号的十进制形式输出整数(正整数无符号),
%o以8进制无符号形式输出整数(不输出前导符0),
%x以16进制无符号形式输出整数(不输出前导符0x),
%u以无符号10进制形式输出整数,
%c以字符形式输出,只输出一个字符,
%s输出字符串,
%f以小数形式输出单、双精度数,隐含输出6位小数,
%e以标准指数形式输出单、双精度数,数字部分小数位数为6位,
%g选用%f或%e格式中输出宽度较短的一种格式,不输出无意义的0。
在%后可插入以下附加符号:
l(L)用于长整型整数,可加在d、o、x、u前面,
m(正整数m)数据最小的宽度,
.n(正整数n)对实数输出n位小数,对字符串截取n个字符,
—输出的数字或字符串在域内向左靠。
rc = math_nint(num)
num:小数形式的实型数,rc:与num最接近的整型数。
rc = math_int(num)
num:小数形式的实型数,rc:num的整数部分(小数部分截去)。
rc = math_pow(num,exponent)
num,exponent:任意数,rc:num的exponent次方。
rc = math_exp(num)
num:任意数,rc:e的num次方。
rc = math_fmod(num1,num2)
num1,num2:任意数,rc:num1除以num2的余数。
rc = math_strlen(string)
string:任意字符变量,rc:字符串长度。
rc = valnum(string)
string:任意字符串变量,rc:0为非数据,1为整数,2为非整数。
rc = wordpos(string,int)
string:任意字符串,常包括一个以上单词,
int:整数,rc:自第int个字符开始的单词数。
注意:
a) 仅对GrADS v1.8以上版本有效。用法示例参见 "script_math_demo.gs"。
12. Some script functions
panels.gsf 根据给定的数值,将实页分成形成若干行和列。
用法:panels(rows cols),示例见panels_demo.gs。
注意:
a) 仅对GrADS v1.8以上版本有效。
b) 要加载Script语言编写的函数,请在脚本文件开始时写上:
rc = gsfallow("on")
c) 通过gsfpath()函数可指定Script语言函数(*.gsf)的路径。用法:
rc = gsfpath("dir_list")
grafana gauge 颜色 grass是什么颜色
转载本文章为转载内容,我们尊重原作者对文章享有的著作权。如有内容错误或侵权问题,欢迎原作者联系我们进行内容更正或删除文章。
提问和评论都可以,用心的回复会被更多人看到
评论
发布评论
相关文章
-
MapReduce分布式软件框架属于国产化吗
&nbs
mapreduce hadoop 数据 应用程序 HDFS