Grads实习报告.docVIP

实习五 1、实习资料： “data”文件夹下二进制月平均风场文件：uv.grd 数据描述文件：uv.ctl 时间范围：2002.1-2005.12 共48个月 水平范围：37×17个网格点，分辨率为2.5°×2.5° 经纬度范围：60-150°E， 0-40°N 垂直层次：850hPa、200hPa 物理量：u，v 2、gs文件 reinit open D:\data\uv.ctl set grid off set xlopts 3 5 0.16 set ylopts 3 5 0.16 set parea 1 10.4 1.2 7.9 set lon 120 set lat 0 40 set lev 200 set t 1 12 set annot 5 8 enable print D:\data\uv200.gmf set xyrev on set gxout shaded set cmin 0 d u cbarn 1 0 5.5 0.3 set grads off set gxout contour set cmax 0 set grads off d u set clevs 0 set clskip 1 set ccolor 2 set cthick 10 set grads off d u draw xlab time draw ylab lat draw title u 120E print disable print ; 3、 实习六 实习资料： “data”文件夹下有全球海表温度资料sst.mnmean.nc。 注意：海温时间从1854年1月开始。 2、shixi6.gs文件（用于提取数据资料） reinit sdfopen D:/data/sst.mnmean.nc set gxout fwrite set fwrite D:/data/sst.grd it=1405 while(it=1764) set lon 90 200 set lat –20 60 set t it d sst it=it+1 endwhile disable fwrite ; 3、Fortran程序 parameter(nx=56,ny=41,nz=1,nt=360) dimension sst(nx,ny,nz,nt) real sst open(15,file=D:\data\sst.grd,form=binary) open(16,file=D:\data\sst.txt) do it=1,nt do iz=1,nz read(15) ((sst(i,j,iz,it),i=1,nx),j=1,ny) enddo enddo do it=1,nt do iz=1,nz write(16,*) ((sst(i,j,iz,it),i=1,nx),j=1,ny) enddo enddo close(15) close(16) end 4、sst.ctl dset C:\data\sst.grd undef -9.99E+33 title sea surface tempture xdef 56 linear 90 2 ydef 41 linear -20 2 tdef 360 linear jan1971 1mo zdef 1 levels 0 vars 1 sst 0 99 surface sea tempture endvars 图片 总结 通过本次实习我掌握了许多知识，比如怎样处理特殊图形，怎样画阴影图形，怎样处理特别等值线与特殊字符。还学会了怎样设置色标，怎样从nc文件中提取数据资料及二进制与十进制数字之间的相互转化。我相信这对我未来的工作是非常有帮助的。