Esope source files

   1 : C ZXFIMP    SOURCE    CHAT      06/06/01    21:23:30     5450           
   2 :       SUBROUTINE ZXFIMP(NEL,K0,NP,IPADL,LE,DRR,G,FLUX,IKA)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
   5 : C************************************************************************       
   6 : C                                                                               
   7 : C     SYNTAXE :                                                                 
   8 : C                                                                               
   9 : C       FIMP (FLUX)
  10 : C                                                                               
  11 : C       CAS FORMULATION EFM1
  12 : C                                                                               
  13 : C************************************************************************       
  14 :       DIMENSION G(*),FLUX(*),DRR(NP,NEL)
  15 :       DIMENSION LE(NP,NEL),IPADL(*)                                                        
  16 : C                                                                               
  17 : C                                                                               
  18 : C                     **************
  19 : C                     *   2D & 3D  *
  20 : C                     **************
  21 :                                                                                 
  22 :                                                                                 
  23 : C     write(6,*)' NEL,NP=',NEL,NP
  24 :       DO 50 K=1,NEL                                                             
  25 :       NK=K0+K                                                                   
  26 :       KA=1+(1-IKA)*(NK-1)
  27 : 
  28 :       DO 50 I=1,NP
  29 :       NF=IPADL(LE(I,K))
  30 :       G(NF)=G(NF)-FLUX(KA)*DRR(I,K)
  31 :                                                                                 
  32 :  50   CONTINUE                                                                                                                                      
  33 :       RETURN                                                                    
  34 :                                                                                 
  35 :  1002 FORMAT(10(1X,1PE11.4))                                                    
  36 :  1001 FORMAT(20(1X,I5))                                                         
  37 :       END                                                                       
  38 :  
  39 :  
  40 :  
  41 :  
  42 :