Esope source files

   1 : C INITFL    SOURCE    BP208322  12/01/10    21:15:01     7244           
   2 :       SUBROUTINE INITFL(IPRIGI,IPMASS,IPMV1,IPV1,IFLU)                          
   3 : *                                                                               
   4 : *                                                                               
   5 : **********************************************************************          
   6 : *                                                                               
   7 : *     SBR REALISANT L INITAILISATION DES VECTEURS DE PROJECTION                 
   8 : *     POUR LA RECHERCHE DE MODES PAR LA METHODE SIMULTANEE                      
   9 : *                                                                               
  10 : *     UNE ATTENTION PARTICULIERE POUR LE CAS OU L ON A DU FLUIDE                
  11 : *     GENERATON D UN VECTEUR RESPECTANT LAPL.PI = 0                             
  12 : *                                                                               
  13 : *                                                                               
  14 : *     D. BROCHARD  (EMSI) 6994       13/07/88                                   
  15 : *     BP 10/01/2012 : IFLU = nombre d'inconnues PI                                                                          
  16 : *                                                                               
  17 : *                                                                               
  18 : **********************************************************************          
  19 : *                                                                               
  20 : *                                                                               
  21 :       IMPLICIT INTEGER(I-N)
  22 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  23 : -INC CCHAMP                                                                     
  24 : -INC CCOPTIO                                                                    
  25 : -INC SMCHPOI                                                                    
  26 :       SEGMENT IDEMEN(0)                                                         
  27 :       CHARACTER*4 MOC                                                           
  28 : *                                                                               
  29 :       IFLU=0                                                                    
  30 : *     
  31 : *      write(6,*) ' initfl  nomdu ' , nomdu(15)                                                                          
  32 :       MCHPOI=IPMV1                                                               
  33 :       SEGACT MCHPOI                                                             
  34 :       NSOUPO = IPCHP(/1)                                                        
  35 :       DO 10 ISOUPO=1,NSOUPO                                                     
  36 :       MSOUPO = IPCHP(ISOUPO)                                                    
  37 :       SEGACT MSOUPO                                                             
  38 :       MPOVAL=IPOVAL
  39 :       segact,MPOVAL
  40 :       NC = NOCOMP(/2)                                                           
  41 :       DO 20 IC = 1,NC                                                           
  42 :       MOC=NOCOMP(IC)                                                            
  43 : c       IF(MOC.EQ.NOMDU(15)) GOTO 110                                             
  44 :       IF(MOC.EQ.NOMDU(15)) IFLU=IFLU+VPOCHA(/1)
  45 : 20    CONTINUE                                                                  
  46 :       segdes,MPOVAL
  47 :       SEGDES MSOUPO                                                       
  48 : 10    CONTINUE                                                                  
  49 :       SEGDES MCHPOI
  50 : c       GO TO 100
  51 :       if(IFLU.eq.0) GO TO 100
  52 : 110   continue                                                       
  53 : C                                                                               
  54 : C     IL Y A DES TERMES EN PI                                                   
  55 : C
  56 : *      write(6,*) ' on passe dans la mise à zero des termes en PI'                                                                               
  57 :       IPMV1L=IPMV1                                                              
  57 :       CALL ANCHPO(IPMV1L,NOMDU(15),IPMV1)                                       
  57 :       CALL DTCHPO(IPMV1L)                                                       
  57 :       CALL DTCHPO(IPV1)                                                         
  58 :       SEGINI IDEMEN                                                             
  59 :       IDEMEN(**)=IPMV1                                                          
  60 :       NOID=0                                                                    
  61 :       NOEN=0                                                                    
  62 :       SEGDES IDEMEN                                                             
  63 :       IDEM = IDEMEN                                                            
  63 :       CALL RESOU1(IPRIGI,IDEM,NOID,NOEN,1D-18)
  64 : *  recuperer nombre de modes d'ensemble et maillage pts contrainnts
  64 :       CALL LIRENT(NBENT,1,IRETOU)
  64 :       call lirobj('MAILLAGE',ipt8,1,iretou)
  65 :       IDEMEN = IDEM                                                             
  66 :       SEGACT IDEMEN                                                             
  67 :       IPV1=IDEMEN(1)                                                            
  68 :       SEGSUP IDEMEN                                                          
  68 :       CALL MUCPRI(IPV1,IPMASS,IPMV1L)                                          
  68 :       CALL ANCHPO(IPMV1L,NOMDU(15),IPMV1)                                      
  68 :       CALL DTCHPO(IPMV1L)                                                       
  69 : c       IFLU = 1                                                                  
  70 : *                                                                               
  71 : 100   CONTINUE                                                                  
  72 : *                                                                               
  73 : *     RETURN                                                                    
  74 :       END                                                                       
  75 :  
  76 :  
  77 :  
  78 :  
  79 :  
  80 :  
  81 :  
  82 :  
  83 :