Télécharger la plaquette Cast3M


fr en

Télécharger initfl.eso

Retour à la liste

Numérotation des lignes : 

initfl

  1. C INITFL    SOURCE    PV        22/04/15    13:20:09     11344          
  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.  
  25. -INC PPARAM
  26. -INC CCOPTIO
  27. -INC SMCHPOI
  28. -INC CCREEL
  29.       SEGMENT IDEMEN(0)
  30.       CHARACTER*(LOCOMP) MOC
  31. *
  32.       xspetl = xspeti
  33.       IFLU=0
  34. *
  35. *      write(6,*) ' initfl  nomdu ' , nomdu(15)
  36.       MCHPOI=IPMV1
  37.       SEGACT MCHPOI
  38.       NSOUPO = IPCHP(/1)
  39.       DO 10 ISOUPO=1,NSOUPO
  40.       MSOUPO = IPCHP(ISOUPO)
  41.       SEGACT MSOUPO
  42.       MPOVAL=IPOVAL
  43.       segact,MPOVAL
  44.       NC = NOCOMP(/2)
  45.       DO 20 IC = 1,NC
  46.       MOC=NOCOMP(IC)
  47. c       IF(MOC.EQ.NOMDU(15)) GOTO 110
  48.       IF(MOC.EQ.NOMDU(15)) IFLU=IFLU+VPOCHA(/1)
  49. 20    CONTINUE
  50.       segdes,MPOVAL
  51.       SEGDES MSOUPO
  52. 10    CONTINUE
  53.       SEGDES MCHPOI
  54. c       GO TO 100
  55.       if(IFLU.eq.0) GO TO 100
  56. 110   continue
  57. C
  58. C     IL Y A DES TERMES EN PI
  59. C
  60. *      write(6,*) ' on passe dans la mise à zero des termes en PI'
  61.       IPMV1L=IPMV1
  62.       CALL ANCHPO(IPMV1L,NOMDU(15),IPMV1)
  63.       CALL DTCHPO(IPMV1L)
  64.       CALL DTCHPO(IPV1)
  65.       SEGINI IDEMEN
  66.       IDEMEN(**)=IPMV1
  67.       NOID=0
  68.       NOEN=0
  69.       SEGDES IDEMEN
  70.       IDEM = IDEMEN
  71.       CALL RESOU1(IPRIGI,IDEM,NOID,NOEN,xspetl,0,0)
  72. *  recuperer nombre de modes d'ensemble et maillage pts contrainnts
  73.       CALL LIRENT(NBENT,1,IRETOU)
  74.       call lirobj('MAILLAGE',ipt8,1,iretou)
  75.       IDEMEN = IDEM
  76.       SEGACT IDEMEN
  77.       IPV1=IDEMEN(1)
  78.       SEGSUP IDEMEN
  79.       CALL MUCPRI(IPV1,IPMASS,IPMV1L)
  80.       CALL ANCHPO(IPMV1L,NOMDU(15),IPMV1)
  81.       CALL DTCHPO(IPMV1L)
  82. c       IFLU = 1
  83. *
  84. 100   CONTINUE
  85. *
  86. *     RETURN
  87.       END
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales