Télécharger initfl.eso

Retour à la liste

Numérotation des lignes :

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

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