Télécharger initfl.eso

Retour à la liste

Numérotation des lignes :

  1. C INITFL SOURCE PV 15/11/25 21:15:09 8707
  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
  58. CALL ANCHPO(IPMV1L,NOMDU(15),IPMV1)
  59. CALL DTCHPO(IPMV1L)
  60. CALL DTCHPO(IPV1)
  61. SEGINI IDEMEN
  62. IDEMEN(**)=IPMV1
  63. NOID=0
  64. NOEN=0
  65. SEGDES IDEMEN
  66. IDEM = IDEMEN
  67. CALL RESOU1(IPRIGI,IDEM,NOID,NOEN,1D-18,0)
  68. * recuperer nombre de modes d'ensemble et maillage pts contrainnts
  69. CALL LIRENT(NBENT,1,IRETOU)
  70. call lirobj('MAILLAGE',ipt8,1,iretou)
  71. IDEMEN = IDEM
  72. SEGACT IDEMEN
  73. IPV1=IDEMEN(1)
  74. SEGSUP IDEMEN
  75. CALL MUCPRI(IPV1,IPMASS,IPMV1L)
  76. CALL ANCHPO(IPMV1L,NOMDU(15),IPMV1)
  77. CALL DTCHPO(IPMV1L)
  78. c IFLU = 1
  79. *
  80. 100 CONTINUE
  81. *
  82. * RETURN
  83. END
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  

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