佩服
李老师奉献精神着实令人钦佩，有大师风范。今后会向李老师多多请教。 下载看一下,楼主辛苦了! 支持 谢谢李老师 这真的要赞一下 ding 原帖由 dushuxiaogongzi 于 2013-6-1 22:28 发表 http://www.cfluid.com/bbs/images/common/back.gif李老师好：
我键入make后，显示mpif90 -fast -c opencfd2d.f90
f951:错误：无法识别的命令行选项“-fast”
make：***错误1
我用的Lin ...
有些编译器不支持 -fast 选项。 可改为 -O2 或 -O3 编译选项。 原帖由 ch06 于 2011-3-19 22:12 发表 http://www.cfluid.com/bbs/images/common/back.gif
在window xp下有编译成功的吗？执行build时，一下493个错误
都是类似这样的错误
error LNK2001: unresolved external symbol _MPI_BCAST@24
是mpi的库没找着吗？
我看external dependencies里有mpif.h，而且我将 ...
今天试着在win7下编译成功了。 # npx0 npy0 LAP
2 2 4
好像只能4核，9核，16核这样算。
8核的cpu，如果设定为24 8，好像不让算。 非常感谢李老师这样的奉献精神。 原帖由 莫口 于 2010-11-3 20:34 发表 http://www.cfluid.com/bbs/images/common/back.gif
不过DNS，没有机器啊……
自己的电脑也用不了……
4核的机器可以跑。 原帖由 shirazbj 于 2013-8-16 10:46 发表 http://www.cfluid.com/bbs/images/common/back.gif
# npx0 npy0 LAP
2 2 4
好像只能4核，9核，16核这样算。
8核的cpu，如果设定为24 8，好像不让算。
原来lap最大是5. 还以为是2*4呢，呵呵。 原帖由 shirazbj 于 2013-8-22 10:12 发表 http://www.cfluid.com/bbs/images/common/back.gif
原来lap最大是5. 还以为是2*4呢，呵呵。
没事写了个求lap最大值的代码，可以把lap设成找到的最大值maxLAP。
首先在opencfd2d.f90的subroutine read_parameter里声明新变量：
!------------------------------------
! bjStart
!------------------------------------
integer nLAP,aLAP(100),maxLAP
!------------------------------------
! bjEnd
!------------------------------------
然后在读完opencfd2d.in后，找最大值：
!------------------------------------
! bjStart
!------------------------------------
!------------------------------------
! Hoam-OpenCFD V1.7.2 使用手册
! 表4 差分方法描述符
!------------------------------------
nLAP=13
aLAP(1)=3
aLAP(2)=3
aLAP(3)=3
aLAP(4)=4
aLAP(5)=4
aLAP(6)=4
aLAP(7)=4
aLAP(8)=4
aLAP(9)=3
aLAP(10)=3
aLAP(11)=3
aLAP(12)=3
aLAP(13)=3
maxLAP=0
do j=1,3
do i=1,3
if aLAP(NUM_Methodx(i,j))>maxLAP then
maxLAP=aLAP(NUM_Methodx(i,j))
endif
enddo
enddo
do j=1,3
do i=1,3
if aLAP(NUM_Methody(i,j))>maxLAP then
maxLAP=aLAP(NUM_Methody(i,j))
endif
enddo
enddo
!------------------------------------
! bjEnd
!------------------------------------
最后把原来的 nparameter(7)=LAP语句，替换为：
!------------------------------------
! bjStart
!------------------------------------
nparameter(7)=maxLAP
!------------------------------------
! bjEnd
!------------------------------------
[ 本帖最后由 shirazbj 于 2013-8-23 07:42 编辑 ] 原帖由 lixl-imech 于 2011-4-12 16:12 发表 http://www.cfluid.com/bbs/images/common/back.gif
压缩文档中是Mach 6 钝锥绕流问题的完整算例及详细的步骤说明，供学习OpenCFD使用。 该算例使用Hoam-Opencfd2d 1.5.2 计算，该软件的源代码在本贴的主贴中下载。
下图为计算得到的压力分布图。
opencfd2d.in里Num_Method_x有0,20,21，但1.7.2使用手册里没介绍。不知道是什么？ if(T(i,j) .lt. 0) then
print*, " T < 0 !!! Computation Stop ! "
print*, i_offset(npx)+i-1, j_offset(npy)+j-1 , T(i,j)
出现(T(i,j)<0,怎么调呀？