爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 24871|回复: 56

[分享资料] [已解决]求助关于利用fortran和grads组合进行小波分析

[复制链接]
发表于 2012-3-18 19:59:50 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 陌小默的 于 2012-3-24 11:38 编辑

我用墨西哥帽小波的程序对一个指数进行分析,指数是1950到2007年的逐年值;利用网上找到的fortran处理原始数据后,生成的dat文件无法绘制出图形,请问该怎么处理啊?下面是我用到的fortran程序,gs脚本还有对应ctl描述文件,刚刚接触grads,不太熟悉,不知道ctl文件是不是描述的不对导致画图出错,我把程序打包上传了,请教一下高手了!
墨西哥帽小波.zip (3.26 KB, 下载次数: 101)
密码修改失败请联系微信:mofangbao
发表于 2012-3-19 08:09:21 | 显示全部楼层
楼主程序里面有些小问题,由于时间仓促我大概改了一下,加上了注释,楼主可以仔细看看,能够出图。
for程序中,wlet的三维定义是可以的,但是我改成了2维的,方便一些;最重要的是输出小波系数那段,将recl改成了nb,先写1~58,外层循环是1~41;
描述文件中,xdef写成了58 (也就是横轴)ydef写成41;要加上zdef 1 linear 1 1
其余的就是grads操作了,打开描述文件wave.ctl  然后  d s1.。。细致的作图技巧我就不多赘述了

index.txt (580 Bytes, 下载次数: 47)

点评

小菜鸟正好要用,拿来学习学习,谢了  发表于 2012-12-20 12:53
厉害  发表于 2012-3-22 15:52

评分

参与人数 2金钱 +30 贡献 +7 收起 理由
mofangbao + 20 + 5
言深深 + 10 + 2

查看全部评分

密码修改失败请联系微信:mofangbao
 成长值: 0
发表于 2012-3-18 20:45:03 | 显示全部楼层
我记得我搞过一个dat转txt的程序,论坛有http://bbs.06climate.com/forum.php?mod=viewthread&tid=2001··通过txt文件可以用surfer绘图http://bbs.06climate.com/forum.php?mod=viewthread&tid=946··上手很快的
祝你成功
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-3-18 21:12:57 | 显示全部楼层

深深版主,我用了,画图时后有问题,才使用墨西哥帽的,
密码修改失败请联系微信:mofangbao
发表于 2012-3-18 22:49:21 | 显示全部楼层
你把你错误的提示发上来看看
密码修改失败请联系微信:mofangbao
 成长值: 0
发表于 2012-3-18 23:19:54 | 显示全部楼层
陌小默的 发表于 2012-3-18 21:12
深深版主,我用了,画图时后有问题,才使用墨西哥帽的,

看了真的需要传图看看··
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-3-19 08:51:55 | 显示全部楼层
言深深 发表于 2012-3-18 23:19
看了真的需要传图看看··

深深,我真的很想用你的那个方法的,但是我用了之后不知道是不是哪里出来问题,实部,模方的文件都生成,但是surfer做图之哟mod的那个能,而且图好像还是有问题的,剩下的两个文件我看了下数据都太小了,转化后基本都是0,没法做图,不知道是什么原因,我把程序和生成文件给您看看 mywave.zip (90.08 KB, 下载次数: 16)
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-3-19 09:08:42 | 显示全部楼层
edwardli 发表于 2012-3-19 08:09
楼主程序里面有些小问题,由于时间仓促我大概改了一下,加上了注释,楼主可以仔细看看,能够出图。
for程序 ...

非常感谢,我回来就研究下
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-3-19 09:15:31 | 显示全部楼层
善人/jw 发表于 2012-3-18 22:49
你把你错误的提示发上来看看

程序运行没有错误,就是绘图出错,生成的图片打不开
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-3-19 09:16:05 | 显示全部楼层
微博评论 发表于 2012-3-18 21:01
啧啧~这里哪有人会

来自 流年某年 的新浪微博

怎么会咩人会嘛,你看啊。高手这么多的
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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