23 lines
641 B
Fortran
23 lines
641 B
Fortran
SUBROUTINE TINT
|
|
C
|
|
C LOGARITHMIC INTERPOLATION COEFFICIENTS FOR INTERPOLATION OF
|
|
C TEMP(ID) TO THE VALUES 5000,10000,20000,40000
|
|
C
|
|
INCLUDE 'PARAMS.FOR'
|
|
INCLUDE 'MODELP.FOR'
|
|
DIMENSION TT(4)
|
|
DATA TT /3.699, 4.000, 4.301, 4.602/
|
|
C
|
|
DO ID=1,ND
|
|
T=LOG10(TEMP(ID))
|
|
J=3
|
|
IF(T.GT.TT(3)) J=4
|
|
JT(ID)=J
|
|
X=(TT(J)-TT(J-1))*(TT(J)-TT(J-2))*(TT(J-1)-TT(J-2))
|
|
TI0(ID)=(T-TT(J-2))*(T-TT(J-1))*(TT(J-1)-TT(J-2))/X
|
|
TI1(ID)=(T-TT(J-2))*(TT(J)-T)*(TT(J)-TT(J-2))/X
|
|
TI2(ID)=(T-TT(J-1))*(T-TT(J))*(TT(J)-TT(J-1))/X
|
|
ENd dO
|
|
RETURN
|
|
END
|