From 0d4f43d355de79178b1142e9735902cf641670b6 Mon Sep 17 00:00:00 2001 From: Dimitri Sokolyuk Date: Mon, 11 May 2009 00:27:49 +0000 Subject: Xfoil 6.97 --- orrs/src/fscorr.f | 87 +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 87 insertions(+) create mode 100755 orrs/src/fscorr.f (limited to 'orrs/src/fscorr.f') diff --git a/orrs/src/fscorr.f b/orrs/src/fscorr.f new file mode 100755 index 0000000..f6392a1 --- /dev/null +++ b/orrs/src/fscorr.f @@ -0,0 +1,87 @@ + PROGRAM FSCORR + PARAMETER (NHX=200) + REAL H(NHX), M(NHX), LSQ(NHX), FUN(NHX), FCORR(NHX) + REAL X(NHX), Y(NHX), Z(NHX) +C + IDEV = 6 + SIZE = 8.0 + CH = 0.02 +C + HMAX = 10.0 + DH = 1.0 +C + FMAX = 0.5 + DF = 0.1 +C + PAR = 0.75 +C + HWT = 1.0/HMAX + FWT = PAR/FMAX +C + OPEN(7,FILE='hfun.fs',STATUS='OLD') + DO 10 I=1, NHX + READ(7,*,END=11) H(I),M(I),LSQ(I),FUN(I) + 10 CONTINUE + 11 CONTINUE + N = I-1 + CLOSE(7) +C + DO 20 I=1, N + HB = 1.0/(H(I)-1.0) +ccc F = 0.22*(1.0 - (5.0*HB-1.0)**2) + 0.5*HB - 0.05 + 3.0*HB**3 +C + F = -0.05 + 2.7*HB - 5.5*HB**2 + 3.0*HB**3 +C + HK = H(I) +ccc TFS = (6.54*HK - 14.07 )/HK**2 +c + TFS = 4.70*HB - 8.45*HB**2 + 3.41*HB**3 + AM = 2.0*F/TFS - 1.0 +C + Z(I) = M(I) +C + Y(I) = AM + X(I) = 10.0*HB + + FUN(I) = M(I) + FCORR(I) = AM +C +CC F = 0.5*(BUH + 1.0)*TFS +C +C X(I) = 10.0*HB +CC +C F = 0.395*(1.0 - 5.8*(HB-0.485)**2) +C Y(I) = F +C Z(I) = FUN(I)*4.0/H(I) +CC +C FCORR(I) = 0.25*F*H(I) +C + 20 CONTINUE +C +C + CALL PLOTS(0,-999,IDEV) + CALL FACTOR(SIZE) +C + CALL PLOT(8.0*CH,8.0*CH,-3) +C + CALL PLOTON +C + CALL XAXIS(0.0,0.0,1.0,DH*HWT,0.0,DH,CH,1) + CALL YAXIS(0.0,0.0,PAR,DF*FWT,0.0,DF,CH,1) +C + CALL XYPLOT(N,H,FUN ,0.0,HWT,0.0,FWT,1,0.3*CH,+1) + CALL XYPLOT(N,H,FCORR,0.0,HWT,0.0,FWT,1,0.3*CH, 0) +C + CALL XYPLOT(N,X,Z,0.0,HWT,0.0,FWT,1,0.3*CH,+1) + CALL XYPLOT(N,X,Y,0.0,HWT,0.0,FWT,1,0.3*CH, 0) +C + CALL PLOTOF +C + WRITE(*,*) 'Hit ' + READ(*,1000) ANS + 1000 FORMAT(A4) +C + CALL PLOT(0.0,0.0,+999) + STOP + END + -- cgit v1.2.3