Sources Esope | Cast3M

scalt
C SCALT     SOURCE    CB215821  16/04/21    21:18:22     8920
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(ABS(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(ABS(V1(1)).LT.1.D-4) V1(1)=0.D0
      IF(ABS(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