爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 61361|回复: 93

[分享资料] 中国160站降水的画图

  [复制链接]

新浪微博达人勋

发表于 2016-5-6 08:28:58 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 洪hong 于 2016-5-15 20:33 编辑

1.在文件夹“fortran_160_1”中是计算中国夏季降水距平百分率和降水距平场(ps:降水距平百分率是短期气候实习老师编写),以及指数与中国降水距平百分率的相关系数子程序(我编写)。
2.在文件夹“fortran_zhandian”中是站点格式输出的程序(ps:短期实习老师编写的程序,我在其基础上做了修改)。
3.直接在文件“zother”是可以绘制中国降水距平的图,或是与中国相关系数场的图。有6个文件,一个是'grid.grd',是格点的二进制文件,这样才能插值绘图;'grid.ctl'是格点二进制文件'grid.grd'的描述文件;'rrSIM.grb'是自己通过(文件夹“fortran_160_1”和)文件夹“fortran_zhandian”处理以后生成的文件,站点格式输出肯定得有文件夹“fortran_zhandian”中处理,文件夹“fortran_160_1”的处理要根据自己的需要,可以进行修改;在'rrSIM.ctl'是自己文件'rrSIM.grb'的描述文件;'china.map'是映射文件,可以通过文件'rrSIM.ctl'在GRADS生成,具体命令看GRADS书第8章站点资料的使用;'rrSIM.gs'是绘制自己图形的gs文件了。有了这6个文件就可以绘制中国站点资料的图了。
4.为了简便,可以直接用我的程序,不用修改名称,只需要在读取文件的时候替换成你自己的数据即可。比如在文件夹“fortran_160_1”中,把你的数据直接放在文件夹“fortran_160_1”中,就和160.f90程序在同一目录下,不用修改路径,修改160.f90程序“open(4,file='indx1.txt')”读取文件名的语句,就替换成你自己的数据了,会生成'rrSMI.txt'的文件,然后直接把生成的文件'rrSMI.txt'直接放在文件夹“fortran_zhandian”中,把原来的'rrSMI.txt'替换掉,也就和zhandian_rrSIM.f90程序在同一目录,因此也不用修改路径,由于zhandian_rrSIM.f90程序中语句“open(2,file='RRSMI.txt')”和生成的数据是相同名字,所以不用修改读取文件的名称,直接运行,会生成'rrSIM.grb'的文件,把生成'rrSIM.grb'的文件放到文件夹“zother”下,把原来旧的'rrSIM.grb'替换掉,直接把'rrSIM.gs'在GRADS中运行即可出图。不会画中国160站要素的图形的同学,我只能帮你到这儿了。
ps.图形的美观自己去学习其他的语句或是软件,我就不教了。

相关系数图

相关系数图

zother.zip

401.4 KB, 下载次数: 191, 下载积分: 金钱 -5

cnriver等底图.zip

233.97 KB, 下载次数: 198, 下载积分: 金钱 -5

评分

参与人数 4金钱 +26 贡献 +5 收起 理由
高小一 + 5 很给力!
pubuzhaxi + 5 洪敏666
mofangbao + 15 + 5
星期六小妹妹 + 1 很给力!

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-6 09:21:19 | 显示全部楼层
弱弱的问一下,indx1.txt这里面是站点经纬度数据?在哪里可以下这个呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-6 09:26:20 | 显示全部楼层
楼主真细致
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-6 09:59:13 | 显示全部楼层
请问index.txt是什么?在哪里下的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-6 12:52:55 | 显示全部楼层
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2016-5-6 15:49:11 | 显示全部楼层
星期六小妹妹 发表于 2016-5-6 09:59
请问index.txt是什么?在哪里下的

那个是要修改成你自己的文件名字,你也可以不用修改,把你文件名修改成index.txt也可以。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-6 15:52:24 | 显示全部楼层
星期六小妹妹 发表于 2016-5-6 09:21
弱弱的问一下,indx1.txt这里面是站点经纬度数据?在哪里可以下这个呢?

这个文件是你自己通过计算以后得出的文件,是你自己的文件哦。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-6 15:54:42 | 显示全部楼层
洪hong 发表于 2016-5-6 15:49
那个是要修改成你自己的文件名字,你也可以不用修改,把你文件名修改成index.txt也可以。

按照您的我已经画出了相关系数。
请问如果单独画降水百分率还需要底下的子程序么?我想单独画百分率或者距平场,,嘻嘻
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-6 16:02:13 | 显示全部楼层
星期六小妹妹 发表于 2016-5-6 15:54
按照您的我已经画出了相关系数。
请问如果单独画降水百分率还需要底下的子程序么?我想单独画百分率或者 ...

那就不需要了,因为子程序是计算相关系数的。如果需要降水距平场与指数(或是一列数)求相关,调用程序的时候,修改一下调用数组的名称就好,如:“call relation(ra,SMI,rrSMI,m,n)”,把ra改成rain。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-6 16:08:40 | 显示全部楼层
洪hong 发表于 2016-5-6 16:02
那就不需要了,因为子程序是计算相关系数的。如果需要降水距平场与指数(或是一列数)求相关,调用程序的 ...

恩恩明白了。那index里面是算什么?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表