流线曲率法解S1流面问题计算程序.docVIP

流线曲率法解S1流面问题计算程序 ******MAIN****** EXTERNAL WG,GW DIMENSION XA(40),YP(40),YS(40) COMMON/HH1/JM,KN,K0,JF,JT,W1,AI,T1,P1,IA,IB,C,D,F,WJ(40),AL3 /HH2/P1S,T1S.B0,B1,B2,GM2,PAI,PT0,AD0,R1,CP,G0,RR,R0,T0,U0 /HH3/OM0,OM,NB,V01,W01,X(40),Y,AL,X0,Y0,Z0,X00,Y00,Z00 ,IM(40),FX /HH4/X1,X2,X3,Y1,Y2,Y3/HH12/A15 /HH5/AMJ(40),RJ(40),OPJ(40),OSJ(40),YPJ(40),YSJ(40)MDPWJ(40), NN Y00=Y00/R0 Z00=Z00/R0 XL=X1/R0 YL=Y1/R0 XE=X2/R0 YE=Y2/R0 W1=W1/U0 R1=R1/R0 388 RH1=RH1*OM0*R0**3/G0 T1S=T1S/T0 P1S=P1S/(G0*OM0*R0) OM=OM/OM0 RR=RR/(U0*U0/T0) CP=RR*GM2 FX=0. CALL QUGEN(WJ(1),1,GM) A15=2*PAI/NB JF1=JF-1 DO 75 J=1,JF1 CALL XY(XL,YL,B2,AMJ(J),Y) YPJ(J)=Y YSJ(J)=Y CONTINUE CALL Y2Y3(XL,YL,Z0,AMJ(JF-1),YPJ(JF-1),X2,Y2,X3,Y3) CALL ABC(AMJ(JF-1),X3,AMJ(JF+1),Y3,YPJ(JF+1), AMJ(JF),Y) YPJ(JF)=Y CALL ABC(AMJ(JF-1),X2,AMJ(JF+1), YPJ(JF-1),Y2,YPJ(JF+1), AMJ(JF),Y) YSJ(JF)=Y CALL Y2Y3(XE,YE,Z00,AMJ(JT+1),YPJ(JT+1),X2,Y2,X3,Y3) CALL ABC(AMJ(JT-1),X3,AMJ(JT+1),YPJ(JT-1),Y3,YPJ(JT+1 ),AMJ(JT),Y) YPJ(JT)=Y CALL ABC(AMJ(JT-1),X2,AMJ(JT+1),YSJ(JT-1),Y2,YSJ(JT+1 ),AMJ(JT),Y) YSJ(JT)=Y DO 85 J=1,JM OPJ(J)=YPJ(J)/RJ(J) OSJ(J)=YSJ(J)/RJ(J)+A15 AD0=(OSJ(J)-OPJ(J))/(KN-1) DO 85 K=1,KN O(J,K)=OPJ(J)+(K-1)*AD0 CONTINUE CALL SCB(1,KN,1) DO 87 J=1,JM DO 87 K=1,KN DWODM(J,K)=0. WRITE(*,450) FORMAT(1X,2HJ=,7X,3HYPJ,11X,3HYSJ,11X,2HAMJ) DO 88 J=1,JM AM1=AMJ(J)*R0 YP1=YPJ(J)*R0 YS1=YSJ(J)*R0 389 WRITE(*,350) J, YP1, YS1,AM1 CONTINUE WRITE(*,411) FORMAT(/5X,’----Go To The Last Program -------‘/) IF(N007.EQ.1) GOTO 95 IF(N008.EQ.1) GOTO 90 CALL TXY2(B1,B2,0,0.2*PT0,N006) GOTO 105 CALL TXY1(B1,B2,1,0.1*PT0,N006) GOTO 105 IF(N008.EQ.1) GOTO 100 CALL TXY1(B1,B2,1,0.1*PT0,N006) GOTO 105 CALL CON(B1,B2,F0) CONTINUE DO 554 J=1,JM AMJ(J)=AM(J)*R0 DO 554 K=1, KN P(J,K)=P(J,K)*(Go*OM0/R0) O(J,K)=O(J,K)*RJ(J)*R0 T(J,K)-=T(J,K)*T0 W(J,K)=W(J,K)*U0 RH(J,K)=RH(J,K)/(OM0*R0* *3/Go) CONTINUE DS1= -(O(JF,KN)-A15*RJ(JF)R0)+O(JF,1) DS2= O(JT,1) - (O(JT,KN)-A15*RJ(JT)*R0) OPEN(8,FILE = ‘YR,DAT’) WRITE(8,*)DS1,DS2 WRITE(8,*)