奉献自己用lvevl set方法计算水下爆炸问题的Fortran源代码。

20040522 · 发表于 2012-4-14 23:38:55

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

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

x

本程序计算自由液面附近的水下爆炸问题，具体可参考相应文献，希望能对各位仁兄有所帮助。

[ 本帖最后由 20040522 于 2013-2-26 04:24 编辑 ]

imchengpeng · 发表于 2012-4-18 12:34:47

很好的学习例子

hilllike · 发表于 2012-4-18 17:31:25

好东西，谢谢分享

030301029 · 发表于 2012-4-18 19:32:54

支持

20040522 · 发表于 2012-5-27 13:39:34

:lol :lol :lol

hua123191492 · 发表于 2012-6-6 16:58:49

支持~~~~~

20040522 · 发表于 2012-8-19 23:02:57

自己顶自己一下

20040522 · 发表于 2012-12-1 16:13:28

:victory: :victory: :victory: :victory: :victory: :victory: :victory: :victory:

20040522 · 发表于 2013-2-26 12:24:54

再顶！

shirazbj · 发表于 2013-7-26 14:29:57

程序有点小问题，求教一下：

算出来的FI都是NaN.

看了一下在SOLVEFI.F的WENOF5子程序里
从124句开始：
   IF(VV(I,J).GT.0) THEN
   SKYY(1)=(FI(I,J-2)-FI(I,J-3))/DY
   SKYY(2)=(FI(I,J-1)-FI(I,J-2))/DY
   SKYY(3)=(FI(I,J)-FI(I,J-1))/DY
   SKYY(4)=(FI(I,J+1)-FI(I,J))/DY
   SKYY(5)=(FI(I,J+2)-FI(I,J+1))/DY
   ENDIF
   IF(VV(I,J).LT.0) THEN
   SKYY(1)=(FI(I,J+3)-FI(I,J+2))/DY
   SKYY(2)=(FI(I,J+2)-FI(I,J+1))/DY
   SKYY(3)=(FI(I,J+1)-FI(I,J))/DY
   SKYY(4)=(FI(I,J)-FI(I,J-1))/DY
   SKYY(5)=(FI(I,J-1)-FI(I,J-2))/DY
   ENDIF

定义了VV(I,J)大于和小于0的情况。可初值给的恰巧是0.
程序随便给了SKYY的值，造成溢出。按F77编译可能没问题，改成F90，算的就不对了。

20040522 · 发表于 2013-10-7 07:34:12

我就是用F77编的程序，F90没用过，你可以仔细调试一下，看看能否找出问题，我现在已经不编程序了。

xiehang199228 · 发表于 2013-10-7 23:55:50

楼主好人哪。

sunco · 发表于 2013-10-8 11:30:08

好东西，顶起来

太子 · 发表于 2017-7-15 08:10:54

谢谢楼主

傻太子 · 发表于 2018-4-2 15:58:06

请问楼主能留下邮箱吗

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

奉献自己用lvevl set方法计算水下爆炸问题的Fortran源代码。

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

奉献一个水下爆炸的Fortran源代码

回复 10# shirazbj 的帖子

站点信息