参考GrADS软件基础教程115页例2,编写gs文件 A、计算1-12月多年平均(1951年-2010年60年) 500hPa气温场,并绘出1月和7月500hPa多年平均气温场,保存为.gmf图形; B、计算1951年1月至2010年12月500hPa气温距平,写入名为ta.grd的二进制数据文件中,注意确定写出数据的格点数;
时间: 2024-02-16 14:00:10 浏览: 26
A、计算1-12月多年平均(1951年-2010年60年)500hPa气温场并绘出1月和7月500hPa多年平均气温场,保存为.gmf图形
```
'reinit'
'sdfopen D:\data\air.2m.mon.mean.nc'
'set t 1 last'
'set z 16'
'set lon 70 140'
'set lat 0 60'
'define tave=ave(t, time=1, time=12)'
'define tave_jan=ave(t, time=1)'
'define tave_jul=ave(t, time=7)'
'gxout grfill'
'run gamitm'
'color white'
'draw string 1.5 55 500hPa Temperature (C) Multi-year Average'
'color navy'
'draw string 1.5 52 Jan'
'draw string 4.5 52 Jul'
'color'
'draw title 500hPa Multi-year Average Temperature'
'draw string 70.5 0.5 60N'
'draw string 139.5 0.5 60N'
'draw string 70.5 59 0N'
'draw string 139.5 59 0N'
'draw map'
'draw blockcont tave_jan'
'printim D:\output\jan_temp.gmf'
'draw blockcont tave_jul'
'printim D:\output\jul_temp.gmf'
```
B、计算1951年1月至2010年12月500hPa气温距平,写入名为ta.grd的二进制数据文件中
```
'reinit'
'sdfopen D:\data\air.2m.mon.mean.nc'
'set t 1 last'
'set z 16'
'set lon 70 140'
'set lat 0 60'
'define tave=ave(t, time=1, time=12)'
'define anom=t-tave'
'set fwrite -be ta.grd'
'define anom_short=anom[gt=ta]'
'define anom_short=anom_short[lt=ta]'
'fwrite anom_short'
```
以上代码中,`reinit` 表示重置所有设定,`sdfopen` 表示打开数据文件,`set` 表示设定维度范围,`define` 表示计算变量,`gxout` 表示输出设定,`run` 表示运行指定的脚本,`color` 表示设定颜色,`draw` 表示绘图,`printim` 表示输出图像。其中,`.gmf` 文件是 GrADS 自己的图形文件格式,`.grd` 文件是 GrADS 的二进制数据文件格式。