Download chaini.eso

Back to the list

   1 : C CHAINI    SOURCE    CHAT      05/01/12    21:54:23     5004           
   2 :       SUBROUTINE CHAINI(A,NMAX,IBOU,NPLEIN,E,DEPSUR,ICAS)                       
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 :       DIMENSION E(*),A(NMAX,*)                                                  
   6 :       DO 1 I=1,IBOU                                                             
   7 :       DO 1 J=1,IBOU                                                             
   8 :   1   A(I,J)=0.                                                                 
   9 :       IBOU2=IBOU                                                                
  10 :       IF(ICAS.EQ.2) IBOU2=IBOU/2                                                
  11 :       DO 6 IC=1,ICAS                                                            
  12 :       IP=(IC-1)*IBOU2                                                           
  13 :       DUM=E(1)*DEPSUR                                                           
  14 :       DO 2 I=1,NPLEIN                                                           
  15 :       IK=I+IP                                                                   
  16 :    2  A(IK,IK)=DUM                                                              
  17 :       IF(NPLEIN.EQ.1) GO TO 6                                                   
  18 :       DUM=E(2)*DEPSUR                                                           
  19 :       DO 3 I=2,NPLEIN                                                           
  20 :       I1=I-1                                                                    
  21 :       IK=I+IP                                                                   
  22 :       DO 3 J=1,I1                                                               
  23 :       JK=J+IP                                                                   
  24 :       A(IK,JK)=DUM                                                              
  25 :    3  A(JK,IK)=DUM                                                              
  26 :       NP1=NPLEIN+1                                                              
  27 :       DUM=E(3)*DEPSUR                                                           
  28 :       DO 4 I=NP1,IBOU2                                                          
  29 :       IK=I+IP                                                                   
  30 :    4  A(IK,IK)=DUM                                                              
  31 :    6  CONTINUE                                                                  
  32 :       RETURN                                                                    
  33 :       END                                                                       
  34 :  

© Cast3M 2003 - All rights reserved.
Disclaimer