1 : C SCALT SOURCE CHAT 05/01/13 03:12:33 5004 2 : C SCALT SOURCE INSL 24/10/96 3 : SUBROUTINE SCALT(SIGM,SIGR,S1,V1,SCT,NSTRS,DTAU,TU, 4 : & AA,BB,DK1,DK2,ALPHA,RB,DTR1,DTR2,TETA) 5 : C 6 : C ************************************************************** 7 : C ** CALCUL DU SCALAIRE (t) ** 8 : C ** CE SOUS-PROGRAMME CALCUL LE SCALAIRE t PROPORTION DE ** 9 : C ** DEFORMATION POUR LAQUELLE L'ECOULEMENT EST RENCONTRE ** 10 : C ** POUR LA PREMIERE FOIS ** 11 : C ************************************************************** 12 : C 13 : IMPLICIT INTEGER(I-N) 14 : IMPLICIT REAL*8(A-H,O-Z) 15 : DIMENSION SIGM(NSTRS),SIGR(NSTRS),S1(NSTRS),V1(4) 16 : C------------------------------------------------------------------ 17 : DO 9 I=1,NSTRS 18 : S1(I)=SIGR(I)+SIGM(I) 19 : 9 CONTINUE 20 : C------------------------------------------------------------------ 20 : CALL CRIOTO(SIGR,SEQ0,FCRI0,NSTRS,TU,AA,BB,DK1,DK2) 20 : CALL CRIOTO(S1,SEQ1,FCRI1,NSTRS,TU,AA,BB,DK1,DK2) 21 : C 22 : DTAU=SEQ1-SEQ0 23 : IF(DABS(DTAU).LT.1.D-4) THEN 24 : SCT=1.D0 25 : ELSE 26 : SCT=(TU-SEQ0)/DTAU 27 : ENDIF 28 : * 29 : DO 2 I=1,NSTRS 30 : S1(I)=SIGR(I)+SCT*SIGM(I) 31 : 2 CONTINUE 32 : C 33 : IF(SCT.LT.0.D0) SCT=0.D0 34 : IF(SCT.GT.1.D0) SCT=1.D0 35 : SCT1=SCT 36 : C 36 : CALL CAPRIN(S1,V1,NSTRS,ALPHA,RB,DTR1,DTR2,TETA) 37 : IF(DABS(V1(1)).LT.1.D-4) V1(1)=0.D0 38 : IF(DABS(V1(2)).LT.1.D-4) V1(2)=0.D0 39 : C------------------------------------------------------------------ 40 : * WRITE (*,*) 'DANS SCALT DTR1=',DTR1,'DTR2=',DTR2,' TETA=',TETA, 41 : * &' SCT=',SCT 42 : C------------------------------------------------------------------ 43 : 1991 FORMAT(18(1X,E12.5)) 44 : RETURN 45 : END 46 : 47 :
© Cast3M 2003 - All rights reserved.
Disclaimer