爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 27240|回复: 30

[作图] 求问:NCL和Fortran做EOF分析的结果存在明显差别的原因

[复制链接]
发表于 2015-7-8 12:05:24 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 hjwen 于 2015-7-8 15:03 编辑

用NCL做EOF分析得到结果,感觉明显有问题。之后用Fotran程序做了下,结果跟预期的差不多,是正确的。两者结果存在较大差别,这可能是什么原因?困惑很久了额,求家园内大神们解答。。
附件说明:分别是NCL和Fortran脚本,输入的是一个指数资料Index(36,365),表示36年365天的指数,把365天当做测站点,36年为时间序列进行EOF分析。

EOF_Index_APO.for

8.68 KB, 下载次数: 49, 下载积分: 金钱 -5

Fortran程序

EOF_Index_APO_1979-2014(DOEd).ncl

3.31 KB, 下载次数: 67, 下载积分: 金钱 -5

NCL脚本

密码修改失败请联系微信:mofangbao
发表于 2015-7-8 12:38:49 | 显示全部楼层

回帖奖励 +1 金钱

做EOF的话注意输入数据的处理,是否去除季节循环之类的

请教把全年逐日逐年的APO指数进行EOF分解有什么物理意义吗?APO不同的季节模态差异还是挺大的啊
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2015-7-8 13:15:30 | 显示全部楼层
niuda 发表于 2015-7-8 12:38
做EOF的话注意输入数据的处理,是否去除季节循环之类的

请教把全年逐日逐年的APO指数进行EOF分解有什么 ...

输入数据两个都没去除,两种方式标准化后得到的EOF结果还是差很多。这个有影响吗?
对多年逐日的APO指数做EOF,是想把多年逐日变化的主模态提取出来,方便看转变时间
密码修改失败请联系微信:mofangbao
发表于 2015-7-8 14:38:24 | 显示全部楼层
上个图看看嘛!
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2015-7-8 15:04:11 | 显示全部楼层

重新编辑貌似加不了图,附件也不行==
密码修改失败请联系微信:mofangbao
发表于 2015-7-9 23:06:35 | 显示全部楼层
hjwen 发表于 2015-7-8 13:15
输入数据两个都没去除,两种方式标准化后得到的EOF结果还是差很多。这个有影响吗?
对多年逐日的APO指数 ...

请教你的APO指数怎么定义的? 是PC1还是选择区域平均?
因为不同的季节对应的区域范围是不同的,而且pc1也是不同的
那么你的每年的365天的apo指数是怎么定义出来的呢
密码修改失败请联系微信:mofangbao
发表于 2015-7-10 07:02:07 | 显示全部楼层
hjwen 发表于 2015-7-8 13:15
输入数据两个都没去除,两种方式标准化后得到的EOF结果还是差很多。这个有影响吗?
对多年逐日的APO指数 ...

好像一般EOF前都要是距平场输入的,不过用ncl验证好像跟原始场结果是一样的,不过跟标准场不一样
密码修改失败请联系微信:mofangbao
发表于 2015-7-10 07:13:17 | 显示全部楼层
hjwen 发表于 2015-7-8 13:15
输入数据两个都没去除,两种方式标准化后得到的EOF结果还是差很多。这个有影响吗?
对多年逐日的APO指数 ...

楼主ncl脚本里用eofcov和eofunc这两个出来的结果有区别吗?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2015-7-10 10:57:31 | 显示全部楼层
niuda 发表于 2015-7-9 23:06
请教你的APO指数怎么定义的? 是PC1还是选择区域平均?
因为不同的季节对应的区域范围是不同的,而且pc1 ...

区域平均的。不同季节对应范围差不太多,就按照赵平老师定义的选的
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2015-7-10 10:59:26 | 显示全部楼层
niuda 发表于 2015-7-10 07:02
好像一般EOF前都要是距平场输入的,不过用ncl验证好像跟原始场结果是一样的,不过跟标准场不一样

我这个问题后来解决了。是发现NCL做EOF会自动处理为距平,而我Fortran是做的原场。NCL可以对原场做EOF么??
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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