求助有关SVD的问题请帮忙

匆匆 · 发表于 2012-8-20 14:54:50

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

本人求算极涡（500hPa高度距平场）和海冰的SVD，出现如下错误：

run-time error M6201:MATH-sqrt:DOMAIN error

对程序中的sqrt函数进行了检查，但未能发现怎样改正，请大家帮帮忙！！

其中：用北纬30度以北的500hPa高度距平场jwjp.bin代表极涡（2.5*2.5）

jwjp.bin (843.75 KB, 下载次数: 1)

568679112 · 发表于 2013-3-28 12:42:01

可能有缺测值要先剔除

568679112 · 发表于 2013-3-28 12:43:32

program tcundef
integer,parameter::NY=121*61
integer,parameter::NT=732
real,parameter::ZMV=-9.99e+33
real ::Y(NY,NT),f(ny,nt)
open(unit=1,file='sst.grd',form='unformatted',access='direct',recl=ny)
   do j=1,nt
      read(1,rec=j)(Y(i,j),i=1,ny)
   enddo
close(1)

l1=0
do i=1,ny
if(y(i,1).ne.zmv)then
      l1=l1+1
do k=1,nt

      f(l1,k)=y(i,k)
enddo
   endif
enddo
write(*,*) l1

open(unit=11,file='tcsst.grd',form='unformatted',access='direct',recl=l1)
   do j=1,nt
      write(11,rec=j)(f(i,j),i=1,l1)
   enddo
close(11)
end

回忆笔记 · 发表于 2013-3-29 12:13:29

数据具体是咋处理的？

云中月影 · 发表于 2013-3-30 10:36:45

我也出现同样的问题了

风儿飘飘 · 发表于 2013-7-23 15:53:27

同问，解决的话，把方法贴上来

风儿飘飘 · 发表于 2013-7-23 16:17:19

检查数据没有缺测值，也出现同样的问题，如何处理呀

越后之龙 · 发表于 2013-11-17 19:19:39

应该是缺测的问题

zjchris · 发表于 2015-12-25 18:53:54

我也出现同样的问题了怎么解决的

		自动登录	找回密码
密码			立即注册

[求助] 求助 有关SVD的问题 请帮忙

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

[求助] 求助有关SVD的问题请帮忙