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

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

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

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

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

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

很好的学习例子

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

好东西，谢谢分享

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

支持

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

:lol :lol :lol :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: :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

回复 10# shirazbj 的帖子

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

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

楼主好人哪。

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

好东西，顶起来:lol

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

谢谢楼主

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

请问楼主能留下邮箱吗
页: [1] 2
查看完整版本: 奉献自己用lvevl set方法计算水下爆炸问题的Fortran源代码。