爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5676|回复: 10

[求助] 投影程序运行的相关问题

[复制链接]

新浪微博达人勋

发表于 2013-5-4 10:15:38 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
程序目的 :将一时间段uv850资料投影到另一时间段特征向量上 得到相应的时间系数
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-5-4 10:16:27 | 显示全部楼层
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
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-5-4 10:16:54 | 显示全部楼层
运行总司出现 severe(157): Program Exception - access violation
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-5-4 10:18:21 | 显示全部楼层
360桌面截图20130504101612.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-5-4 10:25:38 | 显示全部楼层
找出错误了  是倒数第三行 把v1改为v_ni
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-5-4 12:07:00 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2013-5-4 12:38:32 | 显示全部楼层
不知道说啥好了。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-5-4 18:56:26 | 显示全部楼层
mofangbao 发表于 2013-5-4 12:38
不知道说啥好了。。。

{:eb302:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-23 09:40:13 | 显示全部楼层
请问下楼主会算一个空间场到另一个空间场的投影系数么?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-10-21 15:01:33 | 显示全部楼层
Miao 发表于 2013-11-23 09:40
请问下楼主会算一个空间场到另一个空间场的投影系数么?

我现在也在搞这个问题,请问你会了吗?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表