grads基础使用教程.ppt

台站型离散资料画图的一般原则 1、把台站资料r.dat写成二进制文件r.grd 2、编写r.grd数据的数据说明文件r.ctl 3、用r.ctl生成r.map文件 4、生成插值用的背景格点文件 5、编写grid.grd的.ctl文件：grid.ctl 6、把r.grd插值到grid.grd的网格点上，并画图 第一、二步：把台站资料*.dat写成二进制文件*.grd，并编写相应ctl文件 Fortran 程序片段 对应的数据描述文件 dset ^sample.grd title Sample Model Data dtype station stnmap ^sample.map undef 0.10000E+16 xdef 144 linear 0 2.5 ydef 73 linear -90 2.0 zdef 7 levels 1000 900 800 700 500 300 100 tdef 12 linear 00z01apr1985 6hr vars 3 slp 7 99 sea level pressure u 7 99 U component v 7 99 V component endvars parameter (mx=144, my=73, ml=7, mt=12) dimension slp(mx,my,mt), Dimension uwnd(mx,my,ml,mt), vwnd(mx,my,ml,mt) open (10, file=‘sample.grd’, form=‘unformatted’, access=‘direct’, recl=mx*my ) irec=1 do it=1,mt write (10, rec=irec ) (slp(i,j,it),i=1,mx),j=1,my) irec=irec+1 do ilev=1,ml write (10, rec=irec ) (( uwnd(I, j, ilev, it), i=1,mx), j=1,my) irec=irec+1 enddo do ilev=1,ml write (10, rec=irec ) (( vwnd(I,j, ilev, it), i=1,mx), j=1,my) irec=irec+1 enddo enddo close (10) end 第三步：用stnmap 外部命令生成.map文件 方法1、ga-!stnmap enter stn ctl filename : e:/data/sample.ctl 方法2、!stnmap – i e:/data/sample.ctl 特别注意：GrADS中的路径表达方法: 1,引用外部命令时，用unix表达法反斜杠‘/’， 例如用 !stnmap –i e:/data/sample.ctl， 否则要用 !stnmap –i e:\\ data\\sample.ctl 2,在grads的命令提示符下输入的内部命令，如open等，和ctl描述文件中的路径则不区分’/’,和‘\’, 例如：ga-open e:\data\sample.ctl 等同 ga-open e :/data/sample.ctl dset e:\data\sample.grd 等同 dset e:/data/sample.grd 第四步、生成插值用的背景格点文件 因为GrADS只识别格点数据，所以需要把站点上的资料通过插值函数插值到某个格点文件上才能被识别。 根据r.dat提供的信息来写格点文件grid.grd（二进制） 要求： （1） grid.grd文件的精度要高于或等于r.dat的精度 （2） grid.grd文件的范围要大于或等于r.dat的范围 （3）grid.grd文件的每个点上均赋值1（范围与r.dat相同） （4）当grid.grd文件的范围大于r.dat的范围时，除了相同范围的每个点赋值为1外，其余外围的点也都赋值为0，（也可以赋值为1，不过这样做最后在站点数据插值到网格点时，会插值到所期望的网格区域外） （5）插值网格文件的ctl文件的时间一定要与数据资料的时间一致！！！否则会出现Warning from OACRES: Less than two stations的错误信息提示 dset e:\data\china.grd title grid undef 999.9 xdef 620 linear 73.0 0.1 ydef 359 linear 18.0 0.1 zdef 1 levels 1000 tdef 1 linear jan1951 1yr vars 1 g 0 99 grid d