Esope source files

   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 :