找回密码
 注册
查看: 7432|回复: 17

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

  [复制链接]
发表于 2012-4-14 23:38:55 | 显示全部楼层 |阅读模式

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

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

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

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

程序计算结果

程序计算结果
文献计算结果.jpg

codes.rar

1.18 MB, 下载次数: 1088

发表于 2012-4-18 12:34:47 | 显示全部楼层

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

很好的学习例子
发表于 2012-4-18 17:31:25 | 显示全部楼层
好东西,谢谢分享
发表于 2012-4-18 19:32:54 | 显示全部楼层
支持
 楼主| 发表于 2012-5-27 13:39:34 | 显示全部楼层
:lol :lol :lol
发表于 2012-6-6 16:58:49 | 显示全部楼层
支持~~~~~
 楼主| 发表于 2012-8-19 23:02:57 | 显示全部楼层
自己顶自己一下
 楼主| 发表于 2012-12-1 16:13:28 | 显示全部楼层
:victory: :victory: :victory: :victory: :victory: :victory: :victory: :victory:
 楼主| 发表于 2013-2-26 12:24:54 | 显示全部楼层
再顶!
发表于 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,算的就不对了。
 楼主| 发表于 2013-10-7 07:34:12 | 显示全部楼层

回复 10# shirazbj 的帖子

我就是用F77编的程序,F90没用过,你可以仔细调试一下,看看能否找出问题,我现在已经不编程序了。
发表于 2013-10-7 23:55:50 | 显示全部楼层
楼主好人哪。
发表于 2013-10-8 11:30:08 | 显示全部楼层
好东西,顶起来
发表于 2017-7-15 08:10:54 | 显示全部楼层
谢谢楼主
发表于 2018-4-2 15:58:06 | 显示全部楼层
请问楼主能留下邮箱吗
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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