投影程序运行的相关问题

麦田_smile

程序目的 ：将一时间段uv850资料投影到另一时间段特征向量上 得到相应的时间系数

麦田_smile

c~~~~~~~~~~~~~~~~~~~~
      parameter(m=252,n=97*37*2,kvt=2,n0=97*37)
        dimension x(n,m),ect(kvt,m),v(n,kvt),oy0(n0,m),oy(n0,m),
     $ ff0(n,m),ic(n0),x_ni(m,n)

        open(2,file='e:\uv+\result\ecx.grd',form='binary',
     $ access='direct',recl=n0*4)
        irec1=1
        do i=1,kvt
      read(2,rec=irec1) (v(j,i),j=1,n)
        irec1=irec1+1
        write(*,*) v(7178,i)
        enddo


      irec=1
        open(1,file='e:\uv-\u850_3.grd',
     $form='binary',access='direct',recl=n0*4)
        open(11,file='e:\uv-\v850_3.grd',
     $form='binary',access='direct',recl=n0*4)
        do  ky=1,m
      read(1,rec=irec) (oy0(i,ky),i=1,n0)
        read(11,rec=irec) (oy(i,ky),i=1,n0)
      irec=irec+1
        do i=1,n0
          if(abs(oy0(i,ky)).ge.999) then
            ff0(i,ky)=-999.0
            ic(i)=0
          elseif(abs(oy(i,ky)).ge.999) then
            ff0(i+n0,ky)=-999.0
            ic(i)=0
            else
            ic(i)=1
            ff0(i,ky)=oy0(i,ky)
            ff0(i+n0,ky)=oy(i,ky)
            endif
        enddo
             enddo
        write(*,*) ff0(15,12)
      
        call proj(n,m,kvt,v,ect,ff0)

        open(111,file='e:\3\ect.txt')
      do i=1,m
        write(111,100) (ect(j,i),j=1,kvt)
100   format(1x,f10.4)
        enddo
        end

c~~~~~~~~~~~~~~~~~~~~~~~~~~~projection obs by V & eigenvalue~~~~~~~~~~~~~~
        subroutine proj(nn,mm,kvt1,v1,ect1,x)
      dimension x(nn,mm),ect1(kvt1,mm),v1(nn,kvt1),
     $ v_ni(kvt1,nn)

      do i=1,nn
        do j=1,kvt1
        v_ni(j,i)=v1(i,j)
        enddo;enddo

      do 50 it=1,mm
      do 50 j=1,kvt1
        ect1(j,it)=0.0
      do 54 k=1,nn
  54  ect1(j,it)=ect1(j,it)+v1(j,k)*x(k,it)
  50  continue
      END

麦田_smile

运行总司出现 severe(157): Program Exception - access violation

麦田_smile

麦田_smile

找出错误了  是倒数第三行 把v1改为v_ni

river

麦田_smile 发表于 2013-5-4 10:25

                               
登录/注册后可看大图

找出错误了  是倒数第三行 把v1改为v_ni

自问自答的好啊

mofangbao

不知道说啥好了。。。

麦田_smile

mofangbao 发表于 2013-5-4 12:38

                               
登录/注册后可看大图

不知道说啥好了。。。

{:eb302:}

Miao

请问下楼主会算一个空间场到另一个空间场的投影系数么？

立夏shero

Miao 发表于 2013-11-23 09:40
请问下楼主会算一个空间场到另一个空间场的投影系数么？

我现在也在搞这个问题，请问你会了吗？