菜鸟再次求助各位大虾！关于eof用grads作图！

realzly

用魏凤英老师的eof程序做出了结果，如何用grads做出特征向量的站点分布图？

realzly

是不是如grads做站点等值线一样，要写成经度、纬度、数值？？？
还想到一个问题：特征向量是一列取值的话，从上到下数值的排列方式是对应的原始数据中列的排列方式吗？

言深深

1.参看论坛很多的帖子，写一个gs文件；http://bbs.06climate.com/forum.p ... &extra=page%3D1
2.用surfer直接绘图，surfer版面或者fortran版面有http://bbs.06climate.com/forum.php?mod=viewthread&tid=4753

realzly

我还是想问这个：特征向量是以列取值的话，从上到下数值的排列方式是对应的原始数据中列的排列方式吗？因为我这里有一个画图的小程序，基于surfer的但需要手工输入每个站点的数值，我需要知道站点对应的数值。谢谢！

mofangbao

realzly 发表于 2012-2-14 23:14

                               
登录/注册后可看大图

我还是想问这个：特征向量是以列取值的话，从上到下数值的排列方式是对应的原始数据中列的排列方式吗？因为 ...

要看数据怎么放的，看源程序的输出部分，先把源程序理解了，内部运算暂时不理解的先看输入输出部分

mofangbao

realzly 发表于 2012-2-14 23:14

                               
登录/注册后可看大图

我还是想问这个：特征向量是以列取值的话，从上到下数值的排列方式是对应的原始数据中列的排列方式吗？因为 ...

要看数据怎么放的，看源程序的输出部分，先把源程序理解了，内部运算暂时不理解的先看输入输出部分

realzly

嗯，谢谢各位热心的大虾！半路出家搞气象，数理计算机知识还不过关，哎！难啊！多亏了有大家的帮助，不甚感激啊！  又遇到一个问题想要求助！
我用魏凤英老师的小波分析程序：
C        THIS IS A PROGRAM OF WAVELET ANALYSIS BY USING
C       GAUSS REAL OF MCXICAN HAT
        PROGRAM WA
        PARAMETER (N=40,M=2,NMAX=N*3,NS=0,KS=10)  注：n,m值我自己改过        
        DIMENSION F(M,N),X(N),XM(NMAX),IT(M)
        REAL D,Y,Z
C        *****************************************************
C       * N:  SAMPLE SIZE                                   *
C       * M:  NUMBER OF TIME SCALE                          *
C        * NS: STARTING SCALE                                *
C        * KS: INTERVAL OF SCALE                             *
C        * X:  PRIMARY TIME SERIES                           *
C        * F:  WAVELET TRANSFORMATION COEFFICIENTS           *
C        * IT: EXPANDED AND CONTRACTED  SCALES               *
C        *****************************************************
        WAVE(Y)=(1.0-16.0*Y*Y)*EXP(-8.0*Y*Y)
        OPEN(6,FILE='E:\WA\scdata.txt')
        READ(6,*)(X(I),I=1,N)
        DO 10 I=1,N
        XM(I)=X(N+1-I)
        XM(N+I)=X(I)
        XM(2*N+I)=X(N+1-I)
  10        CONTINUE
        DO 20 I=1,M
        IT(I)=I*KS+NS
  20        CONTINUE
        DO 30 J=1,N
        F(I,J)=0.0   到这里老提示：array bounds exeeded  
30        CONTINUE
我自己读程序也觉得不太对劲，这个F(I,J)定义为二维数组，这样赋值是怎么回事啊？

梦想425

同问啊，F(I,J)=0.0,出现同样的错误啊

wode_q_x

很理解楼主的心情，我也是半路出来搞气象，很多东西都是自己摸索，不过，幸好有论坛，可以学到很多。支持一下你，只要勇于尝试，不怕挫折，一定会达到和别人一样的水平！至于小波分析，我没搞过，一直在用功率谱分析。。。。。。加油！

山峰崖渊路

realzly 发表于 2012-2-15 15:34

                               
登录/注册后可看大图

嗯，谢谢各位热心的大虾！半路出家搞气象，数理计算机知识还不过关，哎！难啊！多亏了有大家的帮助，不甚感 ...

        do I=1,M
           dO J=1,N
            F(I,J)=0.0   !此处为替换下面的语句
           enddo
        enddo
C        DO 30 J=1,N
C        F(I,J)=0.0
C  30        CONTINUE
这样就好了，我也算过，是在这里出的问题，可能魏老师不亲自算，都是我们给算得，有时候得自己改一下。