GrADS站点文件作图详细解决方案.pdf

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