Sources Esope | Cast3M

movcen
C MOVCEN    SOURCE    CHAT      05/01/13    01:52:23     5004
      SUBROUTINE MOVCEN(PENTE,DEPS,SELAS,JA,IBOU,DALPHA,SIGMA,SPHER,
     1   SIGEL)
C
C         CALCULE LE MOUVEMENT DU CENTRE EN FLUAGE
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
      DIMENSION DALPHA(*),SIGMA(*),SPHER(*),SIGEL(*)
C
      JX=JA
      Z=PENTE*DEPS/SELAS
      IF(Z.LT.0.01D00) GO TO 920
      EZZ=1.-EXP(-Z)
      Z=1.-EZZ/Z
      DO 921 IB=1,IBOU
      DALPHA(IB)=EZZ*(SIGMA(JX)-SPHER(JX))+Z*(SIGEL(IB)-SIGMA(JX))
  921 JX=JX+1
      GO TO 913
  920 DO 922 IB=1,IBOU
      DALPHA(IB)=Z*(SIGMA(JX)-SPHER(JX)+0.5*(SIGEL(IB)-SIGMA(JX)) )
  922 JX=JX+1
  913 CONTINUE
      RETURN
      END