DO 3 J=1,JV
DI(1)=DL
DO 2 I=1,ID
RR=R/DI(I)
RNI(I,J,K)=DI(I)*VI(J)*RHO/DV
RN=RNI(I,J,K)
F=.094*(RR**.225)+(.53*RR)+(88.*(RR**.44))/(RN**(1.62*(RR**
*.134)))
FI(I,J,K)=F
A1=A1+LOG(F)
A2=A2+LOG(RR)
A3=A3+LOG(RN)
A4=A4+LOG(F)*LOG(RR)
A5=A5+((LOG(RR))*(LOG(RR)))
A6=A6+LOG(RN)*LOG(RR)
A7=A7+LOG(RN)*LOG(F)
A8=A8+((LOG(RN))**2.)
* PRINT 100,RNI(J),DI(I),FI(I,J)
DI(I+1)=DI(I)+(DELTAD*I)
2 CONTINUE
VI(J+1)=VI(J)+VDELTA
3 CONTINUE
T(K+1)=T(K)+TDELTA
4 CONTINUE
A9=A6-(A2*A3/N)
A10=A5-(A2*A2/N)
A11=A6-(A2*A3/N)
A12=A4-(A1*A2/N)
A13=A8-(A3*A3/N)
A14=A7-(A1*A3/N)
A15=A13-(A11*A9/A10)
A16=A14-(A12*A9/A10)
C=A16/A15
B=A12/A10-C*A11/A10
ALN=A1/N-C*A3/N-B*A2/N
A=EXP(ALN)
* PRINT 200,A1,A2,A3,A4,A5,A6,A7,A8
PRINT 300,A,B,C
RN=RNL
DO 7 K=1,KT
DO 6 J=1,JV
D=DL
DO 5 I=1,ID
RN=RNI(I,J,K)
FFCALC(I,J,K)=A*((R/D)**B)*(RN* *C)
ERROR(I,J,K)=(FFCALC(I,J,K)-FI(I,J,K))/FI(I,J,K)
* PRINT 400,T(K),VI(J),D,RN,FI(I,J,K),FFCALC(I,J,K),ERROR(I,J,K)
D=D+(DELTAD*I)
EABS=ABS(ERROR(I,J,K))
84
Previous Page
