- 积分
- 1581
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-3-20
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我做的是高原地表潜热通量的eof,用的是月平均资料处理成的距平年平均资料,这是我的第一模态图,感觉不对,看了相关的潜热的图,做出来并不是这个样子的,希望各位大神能帮我看看到底是哪里出了问题,下面是我的ctl,一共30个模态,把它写成30个时次了。
dset D:\gyw\slhf\eof\result\mo1\egvt_1.bin
undef -9.99E33
xdef 37 linear 77.25 0.75
ydef 20 linear 27.00 0.75
ZDEF 1 Levels 1000
tdef 30 linear 00z01jun1986 1yr
vars 1
eg 0 99 1000hpa
endvars
我用的是言深深大神的贴子里面李建平老师的eof程序,这是我提取第一模态用的Fortran,拜托大家帮忙看看。
parameter(m=37,n=20,it=30,mt=30)
real egvt(m,n,mt),ecof(it,mt)
open(1,file="D:\gyw\slhf\eof\result\egvt.txt")
do j=1,n
do i=1,m
read(1,*)(egvt(i,j,k),k=1,mt)
enddo
enddo
open(10,file="D:\gyw\slhf\eof\result\mo1\egvt_1.bin",form="binary")
do j=1,n
do i=1,m
write(10)(egvt(i,j,k),k=1,mt)
enddo
enddo
|
|