计算EU遥相关指数

TWiligHT

计算资料是1948年至2012年1月的标准化的EU遥相关指数并绘图，公式为
[img]file:///C:/Users/Administrator.R3LYAH2QFQYL51B/AppData/Roaming/Tencent/Users/673475008/QQ/WinTemp/RichOle/0$R%25NUR%7B(9%25(O_OD[4IGMKF.jpg[/img]



[imgfile:///C:/Users/Administrator.R3LYAH2QFQYL51B/AppData/Roaming/Tencent/Users/673475008/QQ/WinTemp/RichOle/%25J1QK%25T@ZJUGHCXCL@WJ77E.jpg[/img]
程序也出现了正确的grd程序如下program ex1118
real h(144,73,12,65),eu1(65),eu2(65),eu3(65),eu(65),euu(65),&
     d(65)
real mean,s,total
open(10,file='e:\dqsx\33\h500.grd',form='binary')
    do it=1,65
do k=1,12
do j=1,73
do i=1,144
read(10)h(i,j,k,it)
enddo;enddo;enddo;enddo
      write(*,*)'read data ok'
  do it=1,65
       eu1(it)=-0.25*h(9,23,1,it)
       eu2(it)=0.5*h(31,23,1,it)
    eu3(it)=-0.25*h(59,17,1,it)
    eu(it)=eu1(it)+eu2(it)+eu3(it)
    total=total+eu(it)
    write(*,*) eu(it)
    enddo
mean=total/65.0  !65年遥相关指数平均值
do it=1,65
d(it)=(eu(it)-mean)**2
s=s+d(it)
s=sqrt(s/65.0)
euu(it)=(eu(it)-mean)/s  !标准化后的1月份遥相关指数时间序列
      write(*,*) euu(it)
end do

open(12,file='e:\dqsx\33\wx.grd',form='binary')
do it=1,65
write(12) euu(it)
end do
end
ctl文件如下dset e:\dqsx\33\wx.grd
undef -9.99E+33
title  Jan
xdef   144    linear    0.0      2.5
ydef   73     linear   -90.0     2.5
zdef     1  levels   500
tdef    65   linear   jan1948  1yr
vars 1
euu    0   99  ******************
endvars
;  
gs文件如下
'reinit'
'open e:\dqsx\33\wx1.ctl'
'enable print e:\dqsx\33\wx1.gmf'
'd euu'
'draw title Jan standardized teleconnection index'
'print'
'disable print'
;
[img]file:///C:/Users/Administrator.R3LYAH2QFQYL51B/AppData/Roaming/Tencent/Users/673475008/QQ/WinTemp/RichOle/0$R%25NUR%7B(9%25(O_OD[4IGMKF.jpg[/img]

但是在grads里就出现了这样的内容
这是啥意思啊？我到底哪里错了啊？

╯①路霺笑

我也出现这问题了，想问问楼主最后如何解决的？

TWiligHT

╯①路霺笑 发表于 2014-12-24 18:28
我也出现这问题了，想问问楼主最后如何解决的？

我咨询了一个博士哥哥