请教关于谱分析和grb2d等资料该如何用

空空空

如题，并没有破坏版规呀，想请教关于谱分析和grb2d等资料该如何用？
这是谱分析的源代码
      SUBROUTINE POW(XS,N,M,PR,P11,PN1)
DIMENSION C11(500),X2(30)
data X2/3.841,6.0,7.8,9.5,11.070,
$      12.592,14.067,15.507,17,18.307,19.675,21.026,22.362,23.685,
$      25,26.292,27.6,28.87,30.144,31.410,32.671,33.924,35.172,
$      31.415,37.652,38.885,40.113,41.337,42.551,43.773
X=0.
DO 20 I=1,N
X=X+XS(I)/FLOAT(N)
20    CONTINUE
      DO 21 I=1,N
XS(I)=XS(I)-X
21    CONTINUE
      CO1=0.
DO 22 I=1,N
CO1=CO1+XS(I)**2/FLOAT(N)
22    CONTINUE
40    FORMAT(/4F8.4)
      M1=M+1
DO 12 J=1,M1
XX=0.
N1=N-J+1
DO 11 I=1,N1
XX=XX+XS(I)*XS(I+J-1)
11    CONTINUE
      C11(J)=XX/FLOAT(N-J+1)/CO1
12    CONTINUE
      PI=4.*ATAN(1.)
DO 13 J=1,M1
CP1=0.
M2=M-1
DO 14 I=1,M2
CS=COS(PI*(J-1)*I/M)
SST=SIN(PI*(J-1)*I/M)
FC=1+COS(PI*I/M)
CP1=CP1+C11(I+1)*CS*FC
14    CONTINUE
      P11(J)=(C11(1)+CP1)/FLOAT(M)
13    CONTINUE
      P11(1)=P11(1)/2.
P11(M+1)=P11(M+1)/2.
DO 15 I=1,M1
IF(I-1) 16,16,17
17    PR(I)=2.*FLOAT(M)/FLOAT(I-1)
      GOTO 15
16    PR(I)=0.
15    CONTINUE
      CP1=0.
DO 37 I=1,M1
CP1=CP1+P11(I)/(M+1)
37    CONTINUE
99    FORMAT(10F7.3)
      V=(2.*N-M*1.0/2.)/M
      IV=NINT(V)
IF(C11(2).LT.0.1)GOTO 66
39    FORMAT(2(F16.4))
      AK1=1-C11(2)**2
BY1=1+C11(2)**2
DO 34 I=1,M1
CS=COS(PI*(I-1)/M)
PNO=CP1*AK1/(BY1-2.*C11(2)*CS)
PN1(1)=PNO*X2(IV)/IV
34    CONTINUE
      GOTO 67
66    DO 68 I=1,M1
      PN1(I)=CP1*X2(IV)/IV
68    CONTINUE
67    PM=0.0
      DO 18 I=1,M1
IF(P11(I).LE.PM)GOTO 18
PM=P11(I)
PRM=PR(I)
18    CONTINUE
      WRITE(7,*)'N,M,FREEDOM,SPECTRUM,PERIOD FOR MAXINMUM'
      WRITE(*,*)'N,M,FREEDOM,SPECTRUM,PERIOD FOR MAXINMUM'
      WRITE(7.*)N,M,V,PM,PRM
      WRITE(*.*)N,M,V,PM,PRM
WRITE(7.*)'AUTOCORRELATION,SPECTRUM,NOISE spe.,PERIOD'
WRITE(*.*)'AUTOCORRELATION,SPECTRUM,NOISE spe.,PERIOD'
      WRITE(*,19)(I,C11(I).P11(I),PN1(I),PR(I),I=1,M+1)
      WRITE(7,19)(I,C11(I).P11(I),PN1(I),PR(I),I=1,M+1)
19    FORMAT(I6,3F7.3,F6.1)
      RETURN
END

mofangbao

可能是你按错了，写成存草稿了，我看没有内容就删了。
楼主这种提问方式别人应该如何回答呢？你自己可以想象一下别人这么给你一个大段的程序和一句话，你觉得有没有必要去看完这一大段程序，然后自己在找个grib数据帮你算功率谱，而且还是在不知道你要干嘛的情况下，太不现实了，没有谁会有这么多的时间的。
在这里讨论需要有自己的观点，需要自己对问题已经有了一定的认识，而不是贴一大段代码加上一句这个该怎么做呢？
你可以先看看这里，提问的智慧：http://bbs.06climate.com/forum.php?mod=viewthread&tid=4571

如有冒犯，请谅解！

空空空

版主mofangbao 说的很对，是我太着急了，我最开始发的帖子讲明了问题，但是真是被我存成草稿了，我现在重新发一下试试，非常感谢

njzqxt

楼主，成了吗？

daerhaosidaer

grb2文件怎么读取阿