IF(NAQ(I) .EQ. NT1
$
.AND.
NAQ(J) .EQ. NT2
$
.AND.
NT3.EQ.0) THEN
BACKSPACE(4)
READ(4,*) NT1,NT2,NT3,(TETM(I,J,K),K=1,6)
REWIND(4)
ELSE
GO TO 103
END IF
7001
CONTINUE
7000
CONTINUE
DO 8000 I=NC1,NN-1
DO 8001 J=I+1,NN
DO 8002 IA=1,NCAT
106
READ(4,*) NT1,NT2,NT3
IF(NAQ(I) .EQ. NT1
$
.AND.
NAQ(J) .EQ. NT2
$
.AND.
NAQ(IA) .EQ. NT3) THEN
BACKSPACE(4)
READ(4,*)NT1,NT2,NT3,(PSIM(I,J,IA,K),K=1,6)
REWIND(4)
ELSE
GO TO 106
END IF
8002
CONTINUE
8001
CONTINUE
8000
CONTINUE
104
CLOSE(4)
C..........................
END OF READING
......................
RETURN
1
CONTINUE
C.. CALCULATION OF VALUES OF THE PARAMETERES FOR CURRENT TEMPERATURE
A0= PF(T,A0M(1),A0M(2),A0M(3),A0M(4),A0M(5),A0M(6))
DO 9000 I=1,NCAT
DO 9001 J=NC1,NN
BET0(I,J)= PF(T,BET0M(I,J,1),BET0M(I,J,2),
$BET0M(I,J,3),BET0M(I,J,4),BET0M(I,J,5),BET0M(I,J,6) )
BET1(I,J)= PF(T,BET1M(I,J,1),BET1M(I,J,2),
$BET1M(I,J,3),BET1M(I,J,4),BET1M(I,J,5),BET1M(I,J,6) )
BET2(I,J)= PF(T,BET2M(I,J,1),BET2M(I,J,2),
$BET2M(I,J,3),BET2M(I,J,4),BET2M(I,J,5),BET2M(I,J,6) )
C0(I,J)= PF(T,C0M(I,J,1),C0M(I,J,2),C0M(I,J,3),
21
Previous Page
