GrADS学习资料：实习七变量与函数的设置.doc

实习七 变量与函数的设置1实习资料：现在“data”文件夹下有全球海表温度资料sst.mnmean.nc。2实习要求：1）根据所提供的数据资料，计算1-12月多年平均（1971年-2000年）海温场，并绘出1月海温多年平均图。2）进一步计算1971年至2000年1月海温距平场，绘制1998年1月和1999年1月海温距平场，海温正距平区填色，画色标，并分别写标题Jan1998和Jan1999。要求设置虚页，将两幅图并排放在一起。最终将图形保存。所有命令编写于.gs文件中。3. 实习目的：掌握函数的使用，进一步练习阴影图形与等值线图形的叠加画图方法。4. 实习步骤：reinitopen d:datasst2.ctlset t 1 12define tem=ave(sst,t+0,t=360,12)set t 1d temdraw title averageprintim d:data1yuepingjun.gif gif white;dset d:datasst2.grdundef 32767.00xdef 180 linear 0 2ydef 89 linear -88 2 zdef 1 linear 1 1 tdef 360 linear jan1971 1movars 1sst 0 -999 Monthly Means of Sea Surface Temperatureendvars4.1 在写字板或记事本程序中按要求编写*.gs文件；reinitsdfopen d:datasst.mnmean.ncset gxout fwriteset fwrite d:datasst2.grdit=1405while(it=1764)set x 1 180set lat -88 88set t itd sstit=it+1endwhiledisable fwrite;reinitopen d:datasst2.ctlset t 1 12define climate=ave(sst,t+0,t=360,12)modify climate seasonalset t 1 360define anom=sst-climateset vpage 0 8.5 5 10.5set gxout shadedset t 325set cmin 0d anomset gxout contourd anomdraw title jan1998cbarn 1 0 4.5 0.35set vpage 0 8.5 0 5.5set gxout shadedset t 337set cmin 0set rbcols autod anomset gxout contourd anomdraw title jan1999cbarn 1 0 4.5 0.35printim d:datashiyan7.gif gif white;4.2启动GrADS绘图软件，运行*.gs文件4.3完成实习报告1） 说明所用资料使用了sst.mnmean.nc文件，用第一个gs文件将其中1971年-2000年的海温场的数据取出来卸载一个grd文件里面，编写对应的CTL文件，写出对应要求的gs文件。2）给出所编写的*.gs文件3）给出所绘图形