p4.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  4   -----------------------------
*-----------------------------------------------------------------------
	GOTO (1,2,3,4,5,6),K
*-----------------------------------------------------------------------
* 	ENTRY GRID
1	MODE=3
	PI=3.14159
	THL=15.*PI/180.
	YL=1.
	R(1)=1.
	L1=7
	M1=7
	CALL SUPPLY(1)
	RETURN
*-----------------------------------------------------------------------
*	ENTRY START
2	CONTINUE
 	LSOLVE(4)=.TRUE.
	LPRINT(4)=.TRUE.
	LPRINT(5)=.TRUE.
	TITLE(4)=' W/WBAR '
	TITLE(5)=' THETA  '
	LAST=6
	NTIMES(4)=4
	NTIMES(5)=4
	TW=1.
	DO 100 J=1,M1
 	DO 100 I=1,L1
	F(I,J,4)=0.
	F(I,J,5)=TW
100	CONTINUE
	AMU=1.
	DPDZ=-2000.
	COND=1.
	RHOCP=1.
	DTDZ=100.
	RETURN
*-----------------------------------------------------------------------
*	ENTRY DENSE
3	CONTINUE
	RETURN
*-----------------------------------------------------------------------
*	ENTRY BOUND
4	CONTINUE
	WSUM=0.
	ASUM=0.
	TSUM=0.
	DO 300 J=2,M2
	DO 300 I=2,L2
	IF(I.EQ.2.AND.J.LE.4)GOTO 300
	AR=YCVR(J)*THCV(I)
	WSUM=WSUM+F(I,J,4)*AR
	TSUM=TSUM+AR*F(I,J,5)*F(I,J,4)
	ASUM=ASUM+AR
300	CONTINUE
	WBAR=WSUM/ASUM
	WP=(R(1)+R(M1))*THL+(1.+THCV(2))*(RMN(5)-R(1))
	DH=4.*ASUM/WP
	RE=RHOCON*WBAR*DH/AMU
	FRE=-2.*DPDZ*DH/(RHOCON*WBAR**2+1.E-30)*RE
	TBULK=TSUM/(WSUM+1.E-30)
	HTP=WP-R(M1)*THL
	HTC=RHOCP*WSUM*DTDZ/((TW-TBULK+1.E-30)*HTP)
	ANU=HTC*DH/COND
	IF(ITER.LE.3)RETURN
	LSOLVE(5)=.TRUE.
	LSOLVE(4)=.FALSE.
	RETURN
*-----------------------------------------------------------------------
*	ENTRY OUTPUT
5	CONTINUE
	IF(ITER.NE.0)GOTO 400
	WRITE(*,401)
401	FORMAT('   ITER',12X,'F*RE',18X,'NU')
400	WRITE(*,402)ITER,FRE,ANU
402	FORMAT(I6,1P2E20.4)
	IF(ITER.NE.LAST)RETURN
	DO 410 J=1,M1
	DO 410 I=1,L1
	F(I,J,4)=F(I,J,4)/WBAR
410	F(I,J,5)=(F(I,J,5)-TW)/(TBULK-TW+1.E-30)
	CALL SUPPLY(2)
	RETURN
*-----------------------------------------------------------------------
*	ENTRY GAMSOR
6	CONTINUE
	DO 500 J=1,M1
	DO 500 I=1,L1
	GAM(I,J)=AMU
	IF(NF.NE.5)GAM(I,J)=COND
	GAM(1,J)=0.
	GAM(L1,J)=0.
	IF(NF.EQ.5)GAM(I,M1)=0.
	IF(J.LE.4)GAM(2,J)=1.E+10
	IF(NF.EQ.5.AND.J.LE.4)GAM(2,J)=10.*COND
500	CONTINUE
	DO 510 J=2,M2
	DO 510 I=2,L2
	CON(I,J)=-DPDZ
	IF(NF.EQ.5)CON(I,J)=-DTDZ*F(I,J,4)*RHOCP
510	CONTINUE
	RETURN
	END