Télécharger scalt.eso

Retour à la liste

Numérotation des lignes :

  1. C SCALT SOURCE CB215821 16/04/21 21:18:22 8920
  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------------------------------------------------------------------
  21. CALL CRIOTO(SIGR,SEQ0,FCRI0,NSTRS,TU,AA,BB,DK1,DK2)
  22. CALL CRIOTO(S1,SEQ1,FCRI1,NSTRS,TU,AA,BB,DK1,DK2)
  23. C
  24. DTAU=SEQ1-SEQ0
  25. IF(ABS(DTAU).LT.1.D-4) THEN
  26. SCT=1.D0
  27. ELSE
  28. SCT=(TU-SEQ0)/DTAU
  29. ENDIF
  30. *
  31. DO 2 I=1,NSTRS
  32. S1(I)=SIGR(I)+SCT*SIGM(I)
  33. 2 CONTINUE
  34. C
  35. IF(SCT.LT.0.D0) SCT=0.D0
  36. IF(SCT.GT.1.D0) SCT=1.D0
  37. SCT1=SCT
  38. C
  39. CALL CAPRIN(S1,V1,NSTRS,ALPHA,RB,DTR1,DTR2,TETA)
  40. IF(ABS(V1(1)).LT.1.D-4) V1(1)=0.D0
  41. IF(ABS(V1(2)).LT.1.D-4) V1(2)=0.D0
  42. C------------------------------------------------------------------
  43. * WRITE (*,*) 'DANS SCALT DTR1=',DTR1,'DTR2=',DTR2,' TETA=',TETA,
  44. * &' SCT=',SCT
  45. C------------------------------------------------------------------
  46. 1991 FORMAT(18(1X,E12.5))
  47. RETURN
  48. END
  49.  
  50.  
  51.  
  52.  
  53.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales