第2章数据处理研究报告.pptVIP

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) 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插值函数将站点资料插值到该格点文件上然后使用了。 设置垂直网格点与气压面的对应关系 ZDEF number LINEAR Start increment 或ZDEF number LEVELS value-list zdef 7 levels 1000 850 700 500 300 200 100 列出具体坐标值 方向: 自下而上 单位: 层 hPa/其他 TDEF number LINEAR start-time increment 设置网格值与时间的映射关系 。 时间表示法 ：hh:mmZddmmmyyyy hour minute day mon year tdef 12 linear 00:00z21jan2000 1mo 时间步长 只有linear关键字 起始时间必须按照上述格式写 年份用4位数写 21日00:00时 按情况可省略 2000年1月 (不可省略) 3mn 分钟 6hr 小时 1dy 天 1mo 月 1yr 年 例如： 12Z1JAN1990 14:20Z22JAN1987 JUN1960 Z：世界时 12Z代表北京时间20点 VARS number 变量描述开始，同时给出变量个数number 。 abrev levs units description ……….. abrev，变量名称，为 1 到8 个字符组成的该变量的缩写名，画图或处理数据将用到，以字母开头 levs，数字，代表变量层数，0表示只有一层 units，单位，为GRIB预留，给99 description，对变量的文字描述，最多40个字符 用ENDVARS表示数据描述文件结束。 dtype 用于说明数据格式的类型 当数据为格点binary格式时省略 dtype的可选参数 – grib – station 站点数据 – hdfsds – netcdf v1.9版本新增 ? ctl书写注意事项 ??? 每行顶格写; 虽然大小写没关系, 推荐在一个关键字中都用 大写或都用小写; 如写成 dset, 而不写成 DsEt； 每行记录的各项以空格分开； 每行记录不超过80个字符； 注释行请用 * 开始, 并顶格写;不能出