grads站点文件作图详细解决方案.docVIP

Fortran+GrADS站点作图详解 mofangbao@气象家园 帖子导读： 以前论坛也有过相关的帖子，但是到目前为止，仍然不断的有朋友看了那些帖子之后还 是不知道如何下手，于是我打算为纯粹的新手写一份说明，这份文档会非常详细，因此适合 刚接触 GrADS 的朋友来参考。 帖子介绍了 fortran 文件中读写数据的基本知识，grads 的站点文件的存放格式，如何通 过 fortran 来写出这个格式的文件，如何建立站点文件的 CTL 文件，如何为站点数据生成映 射文件，站点映射文件的作用，如何通过站点 CTL 文件来预览站点数据，如何用 fortran 生 成 grads 支持的格点文件，如何给格点配置 CTL 文件，如何编写 GS 文件来进行插值，如何 绘制等值线以及阴影图，如何屏蔽区域外的图形等内容。如果你对以上问题存在疑惑，希望 这个帖子能够对你有所帮助，因此，看完这份文档之后，您应该可以做到这些： （会的就跳过吧，我也是犹豫了很久要不要写的，算是新手礼包的一部分吧） 开始之前，先了解站点数据的基本内容 一般情况下，我们看到的站点数据是这样的： 50353 126.39 51.43 100 89 50632 121.55 48.46 100 54 50527 119.45 49.13 100 29 50434 121.41 50.29 100 46 50557 125.14 49.10 100 83 50745 123.55 47.23 100 80 50756 126.58 47.26 100 74 50788 131.59 47.14 100 36 50873 130.17 46.49 100 31 50978 130.57 45.17 100 59 第一列为站号，第二、三列为经纬度，第四、五列为其他数据。当然，也可能是把 经纬度单独存放，然后数据单独存放，比如气候中心免费开放的 160 站的月平均温 度降水值。如果是 micaps 资料可以查看 micaps 用户手册。这个例子中使用的就是 micaps 的第三类数据格式，示例文件可以从这里下载： （资料下载） 除了了解这些，还必须了解该数据的缺测值用什么来表示，这个在后面的数据描述 文件中要用到。 1、在 fortran中读入这些资料 在使用 fortran 之前，请确保你已经安装了 fortran 编译器，如果你是 xp 系统， 建议安装 CFV6，如果是 win7 系统，那么可以安装 microsoft fortran powerstation4.0 （以下简称 4.0），这两个软件在气象资料站均有下载，安装完成后，请在你磁盘的 某个文件夹下面新建一个文件夹用于本次作图，如果使用的是 CVF6，该文件夹的路 径中不要包含中文名。新建文件夹完成后，从开始菜单打开 fortran 编程窗口，然后 点击 file-new，新建一个 free format 的自由格式 fortran 文件,4.0 的直接建立一个 TextFile 即可： CVF 中新建 fortran 文件 4.0 中新建 TextFile 即可 注意：CVF 中先选路径和文件名，然后会建立该文件，4.0 中直接建立一个临 时文件，当你按保存按钮时才会提示你选择保存的路径，如果你之前用 4.0 的 fortran 编译过其他文件，那么请先关闭工作区再新建，从 file-close workspace 可以关闭。 当然也可以直接建立一个文本文件，然后重命名为 sta2grd.f90，记住，要修改后缀 为.f90，表示该文件为自由格式的 fortran 源文件。现在就可以开始着手读取数据了。 （1）先把 fortran 的开头和结尾写了，并且用 implicit none 声明不能使用未定义的变量， 同时添加一个注释，自由格式中可以用一个英文的!来表示该行该符号往后均为注释内容， 不参与编译： Program sta2grd Implicit none !这里是程序的变量声明 !变量声明结束 !程序开始 !程序结束 End （2）定义几个变量 Fortran 中常用的变量类型有整型：integer，单精度实型 real，双精度实型 double，字 符型 character，关于 fortran 变量的更多内容请参考彭国伦的 fortran95 程序设计。 通过观察我们要读取的文件： diamond 3 2011 年 05 月降水量值 0 0 0 0 0 0 1 0 0 1 1000 50353 126.39 51.43 100 89 50632 121.55 48.46 100 54 … 第一行和第二行对我们转换没有什么作用，因此我们只需要定义必须要用到的变量，站 号 stid，类型为字符型，长