作业要求:
作图内容是(选作两个,并简单文字分析):
NCEP再分析资料下载地址:ftp://ftp.cdc.noaa.gov/Datasets
FTP下:海平面气压资料在目录ncep.reanalysis.derived/surface下,多层资料在ncep.reanalysis.derived/pressure下
总结的画图流程:
)的一月份和七月份的平均值的等值线图。所以最后应该画出四张图。
下载的数据如果没有注明是从1971年开始,那就可能是从更早年份开始,数据需要自己做处理,不能直接来画图。
准备工作:
(1)安装grads
(2)下载资料:部分资料(文件大小较小)的名称:olr.mon.mean.nc,sst.ltm.1971-2000.nc
说明:mon.mean.nc和mon.ltm.nc的区别:mon.mean里的数据是从某一年(未必是1971年)开始,1月至12月的月平均值逐年依次排列下来的数据;mon.ltm是从某一年开始到某一年结束这几十年的每年一月至12月分别平均后的12个平均值。也就是mon.ltm是由mon.mean处理得来的。
绘制部分:
(1)最好将下载的数据放在C:或f:这样比较直接的目录下,路径里最好不要有中文。
(2)打开grads,终端显示:“是否显示绘图框”的询问,输入y
(3)终端绘制
分别以sst.ltm.1971-2000.nc 和olr.mon.mean.nc为例:
(I)sst.ltm.1971-2000.nc
Sdfopen
q file
可以看到,x360格,y180格,z一层,时间t有12个时次(即12个月),变量是sst
q dims
set
绘制填色图
set gxout shaded
set cint 5
d sst
cbarn
绘制窗口的结果:
c
绘制等值线图(等值线没色标)
set gxout contour
set cint 5
d sst
绘制窗口的结果:
printim
c 清除图片
(II)olr.mon.mean.nc
reinit
打开方式和查询资料信息类似(I)
注意到,olr资料始于1974.6.01.00,因此,无法画出1971-2000三十年平均的,调整一下,就画1974-2003三十年平均。
以等值线图为例(填色图画法类似)
set
d ave(olr,t+0,t+359,12)
define olrave= ave(olr,t+0,t+359,12)
d olrave
绘制窗口的结果:
q dims
(手工计算一下当前的t+359,此情况时是361)
set
q dims
可以看出,刚好加到2004年6月,恰好1974-2003,30个七月的平均。
printim
quite
(4).gs文件批量绘制
可以发现,用终端得不断地输入,作为练习或尝试比较好,但是已经能完成绘制后,开始大量处理绘图时却很不方便。于是,可以把刚才在终端的输入整理好,写在文本文件里,然后把.txt,改呈.gs,直接用grads运行.gs文件。
例如sst七月平均的.gs文件可写为sst.gs,内容如下(汉字部分为注释,写.gs时不要写;这个.gs的示例是图片格式,还是请手动输入编写一下.gs文件)
终端输入部分:
绘制窗口的输出:
更多的入门和进阶练习:
参考【GrADS绘图软件使用手册】(网上可以下载到),将model.ctl中的链接model.dat的路径设订对,就可以在终端下运行open model.ctl即可,里面有较多的模版和绘图练习。
其他:
此流程主要适合.nc格式的数据,第一次作业的站点数据以及模式计算生成的格点的画法与此处的.nc格式数据不同,需要写.ctl文件,具体请查看grads的指导手册。
老师的建议:
1)我一般是双击grads后,直接就sdfopen *.nc
2)假设t=277是1971年1月, 那么set t 277等价于set time 00z01jan1971;
3)图像输出我是用"enable print **.gmf" 回车 “print" 回车 "disable print"
4)重新打开一个文件 一种是"reinit"表示重新初始化就是不保存前面第一个nc文件 还有种“c" ,那么后面一个变量名要变为'prec.2'
5)如果画图更美观,如“set grads off" "set xlopts