- 积分
- 939
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-7-17
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 云一林 于 2019-11-14 19:34 编辑
数据是1979-2017年39年的5月组成的大气热源q1值grads打开如下
ga-> open F:\lunwen\q1nwns05.ctl
Scanning description file: F:\lunwen\q1nwns05.ctl
Data file f:\lunwen\q1nwns05.dat is open as file 1
LON set to 0 360
LAT set to -85 85
LEV set to 1000 1000
Time values set: 1979:5:1:0 1979:5:1:0
E set to 1 1
ctl如下:
ga-> q ctlinfo
dset f:\lunwen\q1nwns05.dat
title
undef -9.99e+33
xdef 144 linear 0 2.5
ydef 69 linear -85 2.5
zdef 1 linear 1000 1000
tdef 39 linear 00Z01MAY1979 1mo
vars 1
q1 0 99
endvars
reof代码:
! program reof written by prof.wuhongbao.
! nn---时间;mm---空间;np---特征向量个数
parameter(nn=39,mm=9936)
parameter(np=10)
dimension lon(nn,mm),lat(nn,mm)
dimension rh(nn,mm),b(nn,np),reco(nn,mm)
dimension a(mm,mm),v(mm,mm),ramda(mm),rc(np,np),x(nn),y(nn),h(mm),mx(mm),xw(mm),d(mm),rr(np,np),st(np),mmx(mm)
dimension u(mm),vr(mm),wt(mm),wk(mm),wbt(nn),wbk(nn)
dimension vrlv(50)
integer t,k
real mx,mmx
!***************************************************************
!打开数据
open(5,file='F:\lunwen\q1nwns05.dat',form='binary') !若为夏季则读入夏季数据
do t=1,nn
do k=1,mm
read(5)rh(t,k)
enddo;enddo
close(5)
!************************************************************************************8
!结果输出
! 此文件用来写计算误差
open(11,file='F:\lunwen\reof\eof_wuchajisuan05.dat',form='binary')
! 此文件用来写eof分析得到的空间模态
open(21,file='F:\lunwen\reof\eof_kongjian05.dat',form='binary')
! 此文件用来写eof分析得到的时间系数
open(22,file='F:\lunwen\reof\eof_time05.dat',form='binary')
! 此文件用来写reof分析得到主成分荷载(相关系数)
open(10,file='F:\lunwen\reof\correlation_coefficients_05.dat',status='replace')
! 此文件用来写reof分析得到主成分(标准化值)
open(12,file='F:\lunwen\reof\biaozhunhua_05.dat',status='replace')
! open(222,file='G:\reof\1_ioee.txt',status='replace')
! do i=1,nn
! do j=1,mm
! write(222,*) lon(i,j),lat(i,j),rh(i,j)
! enddo;enddo
!
write(*,*)'read data ok'
383 format(1x,10f8.3)
15 format(583f6.0)
reof运行问题:
--------------------Configuration: REOF - Win32 Debug--------------------
Linking...
Debug/REOF.exe : warning LNK4084: total image size 827506688 exceeds max (268435456); image may not run
REOF.exe - 0 error(s), 1 warning(s)
这个到底应该怎么改啊,运行不出来,附件是所有数据从原始到处理成39年的5月的q1数据到画reof全部都有,有谁能看看哪里出问题了吗?
|
|