p1.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  1   ----------------------------
*-----------------------------------------------------------------------
	DIMENSION T(22,22)
	EQUIVALENCE (F(1,1,4),T(1,1))
*-----------------------------------------------------------------------
	GOTO (1,2,3,4,5,6),K
*-----------------------------------------------------------------------
*	ENTRY GRID
1	XL=1.
	YL=2.
	L1=7
	M1=7
	MODE=1
	CALL SUPPLY(1)
	RETURN
*-----------------------------------------------------------------------
*	ENTRY START
2	CONTINUE
 	LSOLVE(4)=.TRUE.
	LPRINT(4)=.TRUE.
	TITLE(4)='   TEMP   '
	LAST=10
	DO 100 J=1,M1
	DO 100 I=1,L1
	T(I,J)=0.
	IF(I.EQ.1.OR.I.EQ.L1)GOTO 101
	IF(J.EQ.1.OR.J.EQ.M1)GOTO 101
	GOTO 100
101	T(I,J)=(X(I)+Y(J)+X(I)*Y(J))
100	CONTINUE
	RETURN
*-----------------------------------------------------------------------
*	ENTRY DENSE
3	CONTINUE
	RETURN
*-----------------------------------------------------------------------
*	ENTRY BOUND
4	CONTINUE
	RETURN
*-----------------------------------------------------------------------
*	ENTRY OUTPUT
5	CONTINUE
	IF(ITER.NE.0)GOTO 400
	WRITE(*,401)
401	FORMAT('  ITER',13X,'T(4,4)',14X,'T(5,3)')
400	WRITE(*,403)ITER,T(4,4),T(5,3)
403	FORMAT(I5,2F20.6)
	IF(ITER.EQ.LAST)CALL SUPPLY(2)
	RETURN
*-----------------------------------------------------------------------
*	ENTRY GAMSOR
6	CONTINUE
	IF(ITER.GT.0)GOTO 501
	DO 500 J=1,M1
	DO 500 I=1,L1
500	GAM(I,J)=1.
501	CONTINUE
	RETURN
	END