|
发表于 2013-8-23 07:37:34
|
显示全部楼层
原帖由 shirazbj 于 2013-8-22 10:12 发表
原来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 编辑 ] |
|