EOF分解程序,附测试数据,图

米农

刚好正在用，谢谢啦~

siciggl

学习啦，感谢

karst

好好。好好好好好好好好好好好好好好好好好好好好好好好好好好好好好好好好

duoduoduo

大家好，我按帖子的内容计算了EOF，用GRADs作图。计算和经纬度设置应该都没问题，想问一下EOF的特征向量这个空间分布合理么？谢谢！

言深深

duoduoduo 发表于 2015-4-28 10:41
大家好，我按帖子的内容计算了EOF，用GRADs作图。计算和经纬度设置应该都没问题，想问一下EOF的特征向量这 ...

看不出来合理与否。
如果你这个是地面要素，估计不太合理，地面要素和地形的关系一般比较大。
如果是高层资料，那么是有可能的，呈现一个类似于波列的传播···

hare

真心不错，下载啦，谢谢啊

hare

吧程序弄了弄，能算，结果没出图呢

学而不厌

感谢大师

星星阿雅

请问楼主，我现在在做毕业论文，搞得焦头烂额。程序运行不了，不知道怎么改，望楼主不吝赐教！！！

我把楼主open文件语句改了一下，变成了open(21,file='F:\890\890\test.txt',status='old'),
F:\890\890\是我存储test的路径。
然后下面对应的读取数据和输出数据为：
open(21,file='F:\890\890\test.txt',STATUS='old')
        do i=1,m
        read(21,*) (x(i,j),j=1,n)
        enddo
        close(21)
        call eof(m,n,mnl,x,ks,er,egvt,ecof)
        open(22,file='F:\890\890\er.txt',STATUS='old')
        do i=1,mnl
        write(22,"(4f30.8)") (er(i,j),j=1,4)
        enddo
        close(22)
        open(23,file='F:\890\890\ecof.txt',STATUS='old')
        do j=1,n
        write(23,"(<mnl>f20.6)") (ecof(i,j),i=1,mnl)
        enddo
        close(23)
        open(24,file='F:\890\890\egvt.txt',STATUS='old')
        do i=1,m
        write(24,"(<mnl>f20.6)") (egvt(i,j),j=1,mnl)
        enddo
        close(24)

执行时出现错误，显示forrt1:severe<59>:list-directed I/O syntax error,unit 21,file F:\890\890\test.txt
Image            PC                Routine            Line                  Source
890.exe         004111C9      Unknown         Unknown           Unkown
Press  any  key to continue

感觉打开文件读取数据就出问题了，但是不知道怎么改，因为以前没怎么接触过，现在做毕业论文又比较着急，所以求指教！！！谢谢！！！1

言深深

看看你的文件里面，是不是有字符，比如null什么的···