Télécharger movcen.eso

Retour à la liste

Numérotation des lignes :

  1. C MOVCEN SOURCE CHAT 05/01/13 01:52:23 5004
  2. SUBROUTINE MOVCEN(PENTE,DEPS,SELAS,JA,IBOU,DALPHA,SIGMA,SPHER,
  3. 1 SIGEL)
  4. C
  5. C CALCULE LE MOUVEMENT DU CENTRE EN FLUAGE
  6. C
  7. IMPLICIT INTEGER(I-N)
  8. IMPLICIT REAL*8(A-H,O-Z)
  9. C
  10. DIMENSION DALPHA(*),SIGMA(*),SPHER(*),SIGEL(*)
  11. C
  12. JX=JA
  13. Z=PENTE*DEPS/SELAS
  14. IF(Z.LT.0.01D00) GO TO 920
  15. EZZ=1.-EXP(-Z)
  16. Z=1.-EZZ/Z
  17. DO 921 IB=1,IBOU
  18. DALPHA(IB)=EZZ*(SIGMA(JX)-SPHER(JX))+Z*(SIGEL(IB)-SIGMA(JX))
  19. 921 JX=JX+1
  20. GO TO 913
  21. 920 DO 922 IB=1,IBOU
  22. DALPHA(IB)=Z*(SIGMA(JX)-SPHER(JX)+0.5*(SIGEL(IB)-SIGMA(JX)) )
  23. 922 JX=JX+1
  24. 913 CONTINUE
  25. RETURN
  26. END
  27.  
  28.  

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