GrADS第2章数据处理讲课.ppt

stnmap 用于将观测站点数据与图形显示时的站点图建立链接以加速显示。 语法为： ga->stnmap -i 站点数据描述文件路径/**.ctl 也可直接进入 stnmap 交互式操作：stnmap.exe位于pcgrads目录中win32e中， 双击进入后输入：站点数据描述文件路径/ **.ctl .map的生成 站点数据小结——格式是死的！ 1. 每个时次的站数随意; 2. 每组站点报是否包含地面报随意;是否包含高空报也随意; 3. 地面报是一个整体, 每次出现时要素及其顺序须完全相同; 4. 高空报也是一个整体, 每次出现时要素及其顺序、层次须完全相同; 5. time表示的是某组站点报相对于 .ctl文件中坐标时刻的时间偏差, 通常取值范围: time ? [-0.5, +0.5] 6. stnmap–i 出错, 排除.ctl文件问题外, 一定一定是数据格式不对! 格点文件的生成： 虽然GrADS系统开发了一些支持站点资料的功能，但其强大功能的运用主要体现在处理格点资料上，所以实际应用时，为了更广泛深入的使用站点资料，常常是将站点资料通过插值函数插值到某个格点文件上然后再使用。 这个格点文件（定为“grid.grd”）可以根据原有站点资料的信息经过编程生成。 在编辑grid.grd程序时，需要注意以下几点要求： grid.grd文件的精度要高于或等于站点数据的精度； grid.grd文件的范围要大于或等于站点数据的范围； grid.grd文件的每个点上均赋值为1； grid.grd文件的时间一定要与站点数据时间一致。 下面给出生成格点文件“grid.grd”的Fortran程序，编写如下： parameter(nx=71,ny=41) real lat(ny),lon(nx) real s(nx,ny) open(1,file='d:\data\grid.grd',form='binary') lat(1)=15.0 lon(1)=70.0 do j=1,ny-1 lat(j+1)=lat(j)+1.0 enddo do i=1,nx-1 lon(i+1)=lon(i)+1.0 enddo do i=1,nx do j=1,ny s(i,j)=1 enddo enddo write(1)((s(i,j),i=1,nx),j=1,ny) end 与该格点文件相对应的数据描述文件grid.ctl如下： dset d:\data\grid.grd undef -999.0 title Sample GRID Data xdef 71 linear 70 1 ydef 41 linear 15 1 zdef 1 linear 500 1 tdef 1 linear jul1951 1mo vars 1 g 0 99 grid data endvars 有了这个格点文件及与其相对应的数据描述文件后，就可以通过 GrADS插值函数将站点资料插值到该格点文件上然后使用了。 例如： 正确的替换为： %y2 两位数年 %y4 4 位数年 %m1 1 或2 位数月 %m2 2 位数月（用 0补齐 1 位数） %mc 3 字符月份缩写 %d1 1 或2 位数天 %d2 2 位数天 %h1 1 或2 位数小时 %h2 2 位小时 设置X方向格点与经度的对应关系 XDEF number <LINEAR Start increment> 或 XDEF number <LEVELS value-list> number（>=1）给定格点数，整形； LINEAR指明是等间隔分布格点， Start起点坐标，负数表示西经； increment网格间距。 LEVELS参数指明是不等间隔分布格点，列 出具体每个格点的坐标值（以空格分开）。 设置Y方向格点与纬度的对应关系 YDEF number <LINEAR Start increment> 或 YDEF number <LEVELS value-list> 或 YDEF number < mapping Start> mapping表示映射方式，有以下几种： GAUSR15 -- 高斯R15 纬度 (Gaussian) GAUSR20 -- 高斯R20 纬度 GAUSR30 -- 高斯R30 纬度 GAUSR40 -- 高斯R40 纬度 在linear中，Start起点坐标负数表示南纬； 在高斯映射中，