《GRADS内部呕心沥血版》由会员分享，可在线阅读，更多相关《GRADS内部呕心沥血版（7页珍藏版）》请在金锄头文库上搜索。
1、一.主要命令：启动和推出 GRADS：1. Grads启动 GRADS，gaGRADS 系统提示符；2. help给出 GRADS 的主要命令；3. quit退出 GRADS。最基本的 GRADS 操作：1. Open打开一个网格点或者站点数据文件供 GRADS 使用（sdfopen打开netCDF 数据） ；2. query（q）查看当前 GRADS 环境的各种参数；3. display（d）显示绘制一幅 GRADS“表达式”的图形；4. set设置绘图的操作方式（包括绘什么图，在哪里绘图，怎样绘图）5. clear（c）清除当前显示的图形。基本命令与参数：1. Lon经度；2. lat纬度
2、；3. lev海拔高度（单位 hpa） ；4. t时次；Z位势高度（返回值为一个数） ；5. mag（u，v）=aqrt（u*u+v*v）开平方和；6. are（Z，t=1,t=5）求位势高度 Z 的时间平均；7. Z- are（Z，t=1,t=5）距平，起始点 0，间距 5；8. Z- are（Z，x=1,x=5）=(lon=5,lon=360)x 表示格点数，距平；9. Z（t-2）-Z(t-1)=Z(t+1)-Z后一时次的 Z 减去当前时次的 Z（槽东移了） ；10.Hcurl（u，v）涡度；11.Vint积分；12.Ps地面气压；13.Cint30等值线间隔为 30；14.Ccolor
3、3等值线颜色设为第三个颜色；15.Gxout图形类型；16.Shaded填充图；17.Contour等值线图；18.Hires高精度地图；19.Csmooth on把平滑功能打开；20.Ccolor0（白色）将前面图完全覆盖（因为前面设成了填充图） ，如果想将两图叠加则先画填充图，后画等值线图；21.Draw title写题目（编写文件名称） ；22.Rector矢量场，右下角有个标尺，表示矢量大小，可以调节；23.Grid方格场；24.q比湿，颜色深浅表示 q 大小，箭头表示风场大小；25.stream流线（只给出运动方向，无大小） ，颜色表示涡度大小q*10000越平表示风速越大，q 越小
4、。mag26.mpdest调用地图；27.nam地图名称；28.digsiz 0.2设置字符大小为 0.2；29.dignum 2设置小数点后 2 位；30.mpral极射赤面范围；31.mproj设置地图投影；32.latlon兰伯特正形圆锥投影；33.scaled墨卡托投影；34.nps北半球极射赤面投影，sps南半球极射赤面投影。35. 一月，Jan;二月，Feb; 三月，Mar.;四月，Apr.;五月，May;六月，Jun;七月，Jul.;八月, Aug;九月，Sep; 十月, Oct;十一月 Nov; 十二月，Dec36. 一月，1-31;二月，32-60; 三月，61-91.;四月
5、，92-121.;五月，122-152;六月，153-182;七月，183-213.;八月,214-244;九月，245-274; 十月,275-305;十一月 306-335; 十二月，336-366二简答题（解释 ctl 文件，修改 ctl 文件）数据表数文件（ctl）：ASCII（二进制）文本文件，描述原始数据的基本信息，包括数据集的文件名、数据类型、数据结构、变量描述、时间空间映射关系等，可以是格点数据，也可以是站点数据，是从其他气象数据转换生成的。格点数据：格式为二进制无格式直接或顺序记录格式，是一个大的五维数据场，包括三维物理空间，一维物理变量，一维时间变量，存放以三维数组片的形式
6、按水平、垂直、物理变量、时间序列顺序排放。例：数据文件 DAILY2000.DAT 的描述文件如下，如果描述文件的文件名为DAILY2000.ctl，并且与数据文件存放在同一文件夹中。DEST E:dataDAILY2000.DAT */路径与文件名/*TITLE Daily data in 2000 */本数据集说明标题/*UNDEF -2.56E33 */数据集的数据类型、格式、选项和缺省值定义（dtype，format，options，undef）/*X 上格点 144 个 XDEF 144 LINEAR 0.0 2.5 */等压面上格点的水平分辨率Y 上格点 73 个 YDEF 73
7、LINEAR -90.0 2.5 是 2.5*2.5 经纬度网格/*Z 上 7 个层面 ZDEF 7 LEVELS 1000 850 700 500 300 200 100时次 366 天每年 TDEF 366 LINEAR 01JAN2000 1dy /*2000 年 1 年逐日的/*变量 9 个 VARS 9 */水平网格，严格对于 fortran 数组存放顺序，第一维从西到东，第二维从南到北/*等压面 UG 7 99 zonal wind （m/s） */纬向风速/*等压面 VG 7 99 meridional wind （m/s） */径向风速/*等压面 WG 7 99 vertica
8、l velocity (Pa/s) */垂直风度/*等压面 ZG 7 99 Gertical height (gpm) */位势高度/*等压面 TG 7 99 Air temperature (K) */温度/*只有下面 5 层 QG 5 99 specific humidity (kg/kg) */比湿/*地面上 PS 0 99 surface pressure （Pa） */气压/*地面上 TS 0 99 surface temperature （K） */气温/*地面上 PP 0 99 rainfall （mm） */降水量/*问：说明数据文件包含资料内容，如包含哪些要素，资料区域，资料
9、水平分布率，垂直层次。答：1.2000 年 1 年逐日的 7 个不同等压面上或地面上得共 9 个要素的全年范围格点资料。2.7 个等压面分别是 1000、850、100hpa。3.等压面上格点的水平分辨率是2.5*2.5 经纬度网格。4.等压面上有 6 个要素，分别是：纬向风速（东西风分量） 、径向风速（南北风分量） 、垂直速度、位势高度、温度和比湿，其中比湿只有下面 5 层。5.地面上有 3 个要素，分别是：气压、气温、降水量。三按要求绘图（写出完整的操作步骤，包括完整的维数环境设置等）1.启动 GRADS 后，打开上述文件，作一张一月份平均的全球地面温度（TS）图，等值线的间距为 4.Op
10、en E:dataDATA2000.CTLSet lon 0 360Set lat -90 90Set lev 1000 */地面/*Set t 1Set cint 4d ave(TS,t=1,t=31)2.使用极射赤面投影，做一张二月平均的北半球 500hpa 位势高度的图，等值线间距 40位势米。Set lon 0 360Set lat 0 90Set lev 500Ser t 32Set mproj npsSet cint 40d (ZG,t=32,t=60)3.做一张区域为 30E-180E，30S-40N 的 5 月 1 日的 850hpa 温度和流场图。Set lon 30 180
11、Set lat -30 40Set lev 850Set t 122d TGd UG,VG4.做一张 6 月 1 日，在 120E 处，90S-90N、1000-100hpa 范围的径向垂直剖面图，包括两个要素：温度，纬向风。Set lon 120Set lat -90 90Ser lev 1000 100Ser t 153Set cint 4 换行 d TGSet cint 5 换行 d UG5.作一张 1 月平均的再 90S-90N、1000-100hpa 范围的纬向平均径向垂直剖面图，包括两个要素，温度 T 和纬向风 USet lon 120Set lat -90 90Ser lev 1
12、000 100Set t 152Ser cint 4d ave(ave(TG,lon=0,lon=360),t=1,t=31)Ser cint 5d ave(ave(UG,lon=0,lon=360),t=1,t=31)6.同一张图上作位于 105E、30N 处一个站的 1 月份地面温度、500hpa 温度和 300hpa与 850hpa 之间的温度差的时间变化曲线。set lon 105set lat 30set lev 1000set t 1 31d TSset lev 500d TGd TG(lev=300)-TG(lev=850)7.作一张区域（70E-130E、15N-55N）夏季
13、6-8 月的总降水量。set lon 70 130set lat 15 55set lev 1000set t 153 d sum(PP,t=153,244)8.作一张 1 月份平均的北半球 500hpa 中纬度（35N-55N）地区纬向环流指数（平均纬向风）和径向环流指数（平均径向风）的变化曲线。set lon 0 360set lat 35 55set lev 500set t 1d ave(ave(UG,lat=35,lat=55),t=1,t=31)d ave(ave(VG,lat=35,lat=55),t=1,t=31)9.作一张 7 月份平均的北半球 500hpa 垂直速度和位势高
14、度图。要求垂直速度用二维填色等值线，位势高度用一般等值线。set lon 0 360set lat 0 90set lev 500set t 1ser gxout shadedd ave(WG,t=183,t=213)set gxout contourset cint 40d ave(ZG,t=183,t=213)示意图：100500100040S EQ 90N四。读程序（理解程序目的，按要求修改程序，完善程序功能）reinitopen E:dataDAILY2000.CTLset lon 0 180set lat 0 90set t 10set lev 500while(1)say say input variable name(zZ):view heightsay input variable name(tT):view temperaturesay input (tT)to quitpull vnamclearif(vnam=qvnam=Q)；break;endifif(vnam=zvnam=Z) set cint 40 d ZG endifif(vnam=tvnam=T) set cint 4 d TG endifendwhile问：1.这段程序的功能是什么？2.修改程序段，使程
