1 : C KPRIOR SOURCE CHAT 05/01/13 01:04:48 5004 2 : SUBROUTINE KPRIOR(NES,NSA,NSB,A,B,U,V,IVU) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8 (A-H,O-Z) 5 : C 6 : DIMENSION A(NES,NSA),B(NES,NSB),U(1),V(1) 7 : C 8 : EPS = 1E-5 9 : C 10 : LS1 = 0 11 : DO 1 IS=1,NSB 12 : F=0. 13 : DO 11 IES = 1,NES 14 : F = F+B(IES,IS)*U(IES) 15 : 11 CONTINUE 16 : F = F + U(NES+1) 17 : LS2=KSIG(F,EPS) 18 : IF (IS.GE.2.AND.LS1.NE.LS2.AND.LS2.NE.0.AND.LS1.NE.0) GOTO 4 19 : IF(LS2.NE.0) LS1=LS2 20 : 1 CONTINUE 21 : C 22 : IF(LS1.EQ.1) THEN 23 : LS1 = 0 24 : DO 2 IS=1,NSA 25 : F=0. 26 : DO 21 IES = 1,NES 27 : F = F+A(IES,IS)*V(IES) 28 : 21 CONTINUE 29 : F = F + V(NES+1) 30 : LS2=KSIG(F,EPS) 31 : IF (IS.GE.2.AND.LS1.NE.LS2.AND.LS2.NE.0.AND.LS1.NE.0) THEN 32 : IVU = 2 33 : RETURN 34 : ENDIF 35 : IF(LS2.NE.0) LS1=LS2 36 : 2 CONTINUE 37 : C 38 : IF(LS1.EQ.1) THEN 39 : IVU=1 40 : RETURN 41 : ELSE 42 : IVU = 0 43 : RETURN 44 : ENDIF 45 : C 46 : ELSE 47 : IVU = 0 48 : RETURN 49 : ENDIF 50 : C 51 : 4 CONTINUE 52 : F=0. 53 : DO 41 IES = 1,NES 54 : F = F+A(IES,1)*V(IES) 55 : 41 CONTINUE 56 : F = F + V(NES+1) 57 : LS2 = KSIG(F,EPS) 58 : IF (LS2.EQ.-1) THEN 59 : IVU = 0 60 : RETURN 61 : ELSE 62 : IVU = 3 63 : RETURN 64 : ENDIF 65 : END 00000070 66 :
© Cast3M 2003 - All rights reserved.
Disclaimer