GrADS第2章数据处理技术总结.ppt

  1. 1、本文档共72页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
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起点坐标负数表示南纬; 在高斯映射中, Start起点坐标表示高斯纬度

文档评论(0)

yy558933 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档