Question

fortran程序求解释

CCC SUB_COIFV : ITEM -- DENL*g ---- IS DELETED

C TURBULENT FLOW
C ADD AP=-SP ITEM

C THE SUBROUTINE FOR CALCULATING V-VELOCITY
C --------------------------------------------------
SUBROUTINE VELOSV
INCLUDE 'DEFIN.DIM'

DIMENSION RRR(IMAX,JJMAX)

c DO NTIME=1,30
CALL COIFV

c ----- at 4 lines of meniscus, V=0 ------
c DO NNTIME=1,3
DO NTIME=1,30
CALL SOLVER(V,2,IMAX-1,3,JJMAX-1)
c END DO
RRRM=0.0
RRMM=0.0
DO J=3,JJMAX-1
DO I=2,IMAX-1
RRR(I,J)=AP(I,J)*V(I,J)-AE(I,J)*V(I+1,J)-AW(I,J)
+ *V(I-1,J)-AN(I,J)*V(I,J+1)-AS(I,J)*V(I,J-1)
+ -BB(I,J)
RRR(I,J)=ABS(RRR(I,J))
if (RRR(I,J).gt.1.e+20) then
aaa=1.
end if

END DO
END DO

DO 112 J=3,JJMAX-1
DO 112 I=2,IMAX-1
IF(RRRM.LT.RRR(I,J)) THEN
RRRM=RRR(I,J) ! max
MAXI=I
MAXJ=J
END IF
RRMM=RRMM+RRR(I,J)
112 CONTINUE

Answer 1

但从这段程序里来看，这是个计算流速的子程序，开头引用了一个dim文件，里面可能有些参数，中间还引用了两个另外的子程序，不知道这两个子程序的内容的情况下，确实不太好理解到底是算啥的。不过里面其他语句除了循环都是单纯的计算，应该是直接从公式转化来的。