p5.for

关于simple算法

*=======================================================================
      SUBROUTINE USER(K)
*-----------------------------------------------------------------------
$INCLUDE:'SIMPLE.INC'
      DIMENSION TH(22),THU(22),THDIF(22),THCV(22),THCVS(22)
      EQUIVALENCE(X,TH),(XU,THU),(XDIF,THDIF),(XCV,THCV),
     + (XCVS,THCVS),(XL,THL)
*-----------------------------------------------------------------------
      COMMON/CNTL/LSTOP
      COMMON/SORC/SMAX,SSUM
      COMMON/COEF/FLOW,DIFF,ACOF
*-----------------------------------------------------------------------
*--------------------------   PROBLEM  5   -----------------------------
*-----------------------------------------------------------------------
	DIMENSION T(22,22)
	EQUIVALENCE (F(1,1,4),T(1,1))
*-----------------------------------------------------------------------
	GOTO (1,2,3,4,5,6),K
*-----------------------------------------------------------------------
* 	ENTRY GRID
1 	MODE=1
	XL=1.
	YL=1.
	L1=7
	M1=7
	CALL SUPPLY(1)
	RETURN
*-----------------------------------------------------------------------
*	ENTRY START
2	CONTINUE
 	LSOLVE(4)=.TRUE.
	DO 100 K=1,4
100	LPRINT(K)=.TRUE.
	TITLE(1)=' VEL U '
	TITLE(2)=' VEL V '
	TITLE(3)=' STR FN'
	TITLE(4)=' TEMP '
	LAST=5
	DO 110 J=1,M1
	DO 110 I=1,L1
	U(I,J)=10.*FLOAT(I-2)/FLOAT(L1-2)
	V(I,J)=-10.*FLOAT(J-2)/FLOAT(M1-2)
	T(I,J)=50.
	T(1,J)=100.
110	CONTINUE
	RETURN
*-----------------------------------------------------------------------
*	ENTRY DENSE
3	CONTINUE
	RETURN
*-----------------------------------------------------------------------
*	ENTRY BOUND
4	CONTINUE
	DO 300 J=2,M2
300	T(L1,J)=T(L2,J)
	DO 301 I=2,L2
301	T(I,1)=T(I,2)
	RETURN
*-----------------------------------------------------------------------
*	ENTRY OUTPUT
5	CONTINUE
	IF(ITER.NE.0)GOTO 400
	WRITE(*,401)
401	FORMAT('   ITER',11X,'T(4,5)',14X,'T(5,3)')
400	WRITE(*,402)ITER,T(4,5),T(5,3)
402	FORMAT(I6,2F20.5)
	IF(ITER.EQ.LAST)CALL SUPPLY(2)
	RETURN
*-----------------------------------------------------------------------
*	ENTRY GAMSOR
6	CONTINUE
	DO 500 J=1,M1
	DO 500 I=1,L1
	GAM(I,J)=1.
	GAM(I,1)=0.
	GAM(L1,J)=0.
500	CONTINUE
	RETURN
	END