Grads实习报告4.docVIP

实习七 实习资料： “data”文件夹下有全球海表温度资料sst.mnmean.nc。 gs文件: （1）转换为grd文件 reinit sdfopen d:\data\sst.mnmean.nc set fwrite d:\data\sst.grd set gxout fwrite it=1045 while(it=1764) set lat -88 88 set lon 0 360 set z 0 set t it d sst it=it+1 endwhile disable fwrite reinit ; （2）转换为ctl描述文件 dset d:\data\sst.grd title sst undef 32767 xdef 181 linear 0 2.0 ydef 89 linear -88 2.0 zdef 1 levels 0 tdef 360 linear jan1971 1mo vars 1 sst 0 0 sst Endvars （3）1月海温多年平均图 reinit open d:\mean.ctl enable print d:\sst.mnmean.gmf set gxout contour set lat -88 88 set lon 0 360 set t 1 12 define uclim=ave(sst,t+0,t=360,12) set t 1 d uclim print disable print reinit ; （4）1998年1月和1999年1月海温距平场 reinit open D:\data\sst.ctl set grid off enable print d:\data\shixi7.gmf set lat -90 90 set lon 0 357.5 set t 1 12 asst=ave(sst,t+0,t=360,12) modify asst seasonal set time jan1971 set t 1 360 nasst=sst-asst d asst set vpage 0 5.5 0 8.5 set gxout shaded set time JAN1998 set cmin 0 d nasst cbarn 0.5 0 2.5 1.5 draw title JAN 1998 set vpage 5.5 11 0 8.5 set gxout shaded set time JAN1999 set cmin 0 d nasst cbarn 0.5 0 2.5 1.5 draw title JAN 1999 print disable print； 图片 实习八 实习资料： “data”文件夹下有全球海表温度资料sst.mnmean.nc。 (1)gs文件 reinit open d:\data\mean8.ctl set fwrite d:\data\mean3.grd set gxout fwrite it=6 while(it=354) set lat -88 88 set lon 0 360 define uclim=ave(sst,t=it,t=it+2) d uclim it=it+12 endwhile disable fwrite reinit ; (2)转换为ctl文件 dset d:\mean3.grd title sst undef 32767 xdef 181 linear 0 2.0 ydef 89 linear -88 2.0 zdef 1 levels 0 tdef 30 linear jun1971 12mo vars 1 sst 0 0 sst Endvars （3）gs文件 1971年至2000年夏季海温距平图 reinit open d:\data\tiqu.ctl enable print d:\shixi8..gmf i=1 temp=1971 while(i=30) define a=ave(sst,t+0,t=30) define anom=sst-a set gxout shaded set lat -88 88 set lon 0 360 set cmin 0 d anom cbar 1 0 5.5 4 set gxout contour set lat -88 88 set lon 0 360 set cmax 0 d anom draw title temp temp=temp+1 print c i=i+1 endwhile disable print reinit ;