老师讲课速记：
打开ctl文件后，可以用q file 查看文件信息
ctl文件可以用记事本打开编辑，把dat文件的路径设置好。
若要画一个这样的：1987年1月1号，高度为地面的，经度纬度为全球的，变量为纬向的风场
先设定时间 set time 01jan1987
再设置高度 set lev 1000
设置纬度 set lat -90 90
设置经度 set lon 0 360
显示图形 d u
Ok了
 
若要画1987年01月01日-到01月05日全球纬向的平均风场：
set lev 1000
set lat -90 90
set lon 0 360
d ave(u,t=1,t=5) 二维的等值线图
 
c 把图形擦掉，再继续
 
画全球的地面纬向风场填色图（u即表示是纬向的，v是经向的）
set lev 1000
set lat -90 90
set lon 0 360
set gxout shaded
d ave(u,t=1,t=5)
 
cbarc 在图的右上角给出了标尺
cbarn 在图的下方给出标尺
 
填色图的基础上叠加等值线图，不要清掉，再输入一个命令。Set gxout contour
全部步骤：
set lev 1000
set lat -90 90
set lon 0 360
set gxout shaded
d ave(u,t=1,t=5)
cbarc
cbarn
Set gxout contour
--------------------------
复习题：
一、填空
1、  等值线图中强迫所有等值线标记：set clab forced
在grads中定义图形类型：二维等值线set gxout contour二维填图set gxout shaded
2、  Grads地图投影方式：Set mproj scaled（不固定的投影角）
3、  设置y坐标轴标记：set ylevs 1000 700 500 300 200……
4、  设置纵坐标使用对数压力坐标：set zlog on（ 对z方向取对数）
等值线绘图中，人为的设置不等间距绘图，并设定特定的等值线：Set clevs clev1 clev2 …
5、  一般变量的表达方式：temp 3 （lev=500）解释：表示第三个文件中高度为500hpa的temp变量
6、  时间的格式：hh mm z dd mmm yyyy 如 12 Z 01 jan 1996   (实际应用中间无空格12z01jan1996)
        Z前面表示是小时和分钟。小时分钟缺省值为00分，日期缺省为01日。年份和月份不能缺省。
7、  通用求平均函数的用法：ave （expr，time=a，time=b，flag ）a表示时间起点，b表示时间终点，flag表示取值间隔
例如：ave (z, t=1, t=20)     表示，时间从t=1时次到t=20时次对变量Z求时间平均。
9、*.nc文件的打开方式：sdfopen *.Nc
10、grads中数据的存放顺序：X→Y→Z→Var→T
11、常用函数 abs→绝对值 mag→平方根 hcurl→相对涡度 hdivg→水平辐散
12、用grads制作gif和gmf图形：
输出gmf图形的命令（三个命令）：1、enable print *.gmf （创建一个gmf文件） 2、print（当前窗口中的图形保存到创建的文件中 ）3、disable print（关闭创建的文件）
输出gif图形的命令（仅一个命令）：printim *.gif
（Gmf文件用grads安装目录下的gv32.exe打开。

二、
1、看程序说明：（括号内为说明文字）这是一个ctl文件，可用记事本编辑
dset timedata.dat （表示在ctl文件中调用的dat文件名）
xdef 145 linear2.0 2.5 (表示全球水平方向0—360°，格距为2.5 )
ydef 73 linear -90.0 2.5 (表示纬度为 -90—90)
zdef 7 levels 1000 925 850 700 600 500 400（ 表示高度有7层，分别为1000，925…）
tdef 366 linear 01jan2000 1dy    （表示366个时次）
vars 2 （表示有2个变量）
temp 7 99 （表示温度）
vhum 0 99（表示相对湿度）
endvars    (结束)

2、作一全球二月平均的地面温度等值线
Open *.ctl
Set lat -90 90
Set lon 0 360
Set lev 1000
Set cint 4 (设置等值线间隔为4°)
D ave(temp,time=01feb2000,time=29feb2000)
D ave(temp,time=01feb2000,time=29feb2000)-273 (这句表示转化为华氏温度)
Draw title 需要添加的标题
 
2、作一个夏季6-8月北半球500hpa位势高度场（zg）的极地投影。保存到d:\date\daiy.gmf
Open *.ctl
Set mproj nps (设置为北半球)
set lev 500
Set lat 0 90
D ave(zg,time=01jun2000,time=31aug2000) 6月份的正规缩写为june
Enable print d:\date\daiy.gmf
Print
Disable print
 
3、7月北半球，35°N纬向平均的地面温度逐日变化图，保存到d:\date\temp.gmf
Open *.ctl
Set lat 35
set lev 1000
Set lon 0
Set time 01jul2000 31jul2000
D ave(temp lon=0,lon=360)
存图步骤同上题加粗部分，修改保存的文件名即可。
 
三 1、d:\dat目录下有mmmyyyy.dat，缺省值为2.56E33（表示有缺测数据时，统一用该数值代替，可自定义）
时间为1979年1月到2003年12月共300个文件，144*73格距，1000，850...17个层次，6个要素u,v,t,ps,ts,q
Ctl文件如下
Dset d:\dat\%m3%y4.dat（注释：%m3表示三位数的月份，%y4表示四位数年份）
Undef 2.56 E33
xdef 144 （x方向网格点数）
Ydef 73 （y方向网格点数）
Zdef 17 （表示Z方向网格点数，共17层等压面）
Vars 6 （表示变量描述开始，且6个变量）
Temp 7 99 temperature
u 17  99  u wind  
v 17 99  v wind  
ps 0  99  surface pressure  英语的注释可以免写
ts 0 99  surface temperature 英语的注释可以免写
q 17 99  spocific humidity     英语的注释可以免写
Endvars
 
2、多年平均夏季7月份北半球500hpa的风场
Open *.ctl
Set lat 0 90
Set lev 500
D ave (u,t=7,t=300,12)
D ave (v,t=7,t=300,12)

3、全球多年平均的7月份地面气压场（填色图）和温度场（等值线图）
Open *.ctl
Set lat -90 90
Set lon 0 360
Set lev 1000
Set gxout shaded （填色图）
D ave（ps，t=7,t=300,12）（ps表示变量为气压）
Set gxout contour（等值线图）
D ave (ts, t=7,t=300,12) （ts表示变量为温度）

有不明白的地方，可以查grads中文手册

第二页是中长期天气预报复习。在成都呆了这么久，就学了这么些东西，也许很快就忘记了，保存一下，就当纪念吧。