C SCALT SOURCE CHAT 05/01/13 03:12:33 5004 C SCALT SOURCE INSL 24/10/96 SUBROUTINE SCALT(SIGM,SIGR,S1,V1,SCT,NSTRS,DTAU,TU, & AA,BB,DK1,DK2,ALPHA,RB,DTR1,DTR2,TETA) C C ************************************************************** C ** CALCUL DU SCALAIRE (t) ** C ** CE SOUS-PROGRAMME CALCUL LE SCALAIRE t PROPORTION DE ** C ** DEFORMATION POUR LAQUELLE L'ECOULEMENT EST RENCONTRE ** C ** POUR LA PREMIERE FOIS ** C ************************************************************** C IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) DIMENSION SIGM(NSTRS),SIGR(NSTRS),S1(NSTRS),V1(4) C------------------------------------------------------------------ DO 9 I=1,NSTRS S1(I)=SIGR(I)+SIGM(I) 9 CONTINUE C------------------------------------------------------------------ CALL CRIOTO(SIGR,SEQ0,FCRI0,NSTRS,TU,AA,BB,DK1,DK2) CALL CRIOTO(S1,SEQ1,FCRI1,NSTRS,TU,AA,BB,DK1,DK2) C DTAU=SEQ1-SEQ0 IF(DABS(DTAU).LT.1.D-4) THEN SCT=1.D0 ELSE SCT=(TU-SEQ0)/DTAU ENDIF * DO 2 I=1,NSTRS S1(I)=SIGR(I)+SCT*SIGM(I) 2 CONTINUE C IF(SCT.LT.0.D0) SCT=0.D0 IF(SCT.GT.1.D0) SCT=1.D0 SCT1=SCT C CALL CAPRIN(S1,V1,NSTRS,ALPHA,RB,DTR1,DTR2,TETA) IF(DABS(V1(1)).LT.1.D-4) V1(1)=0.D0 IF(DABS(V1(2)).LT.1.D-4) V1(2)=0.D0 C------------------------------------------------------------------ * WRITE (*,*) 'DANS SCALT DTR1=',DTR1,'DTR2=',DTR2,' TETA=',TETA, * &' SCT=',SCT C------------------------------------------------------------------ 1991 FORMAT(18(1X,E12.5)) RETURN END