Télécharger fimvf2.eso

Retour à la liste

Numérotation des lignes :

  1. C FIMVF2 SOURCE CHAT 05/01/13 00:01:06 5004
  2. SUBROUTINE FIMVF2(ILIINC,ICEN,IRO,IROVIT,IGRAV,
  3. & IRES)
  4. C************************************************************************
  5. C
  6. C PROJET : CASTEM 2000
  7. C
  8. C NOM : FIMVF2
  9. C
  10. C DESCRIPTION : VOIR FIMVF1
  11. C
  12. C Gaz ideal mono-espece:
  13. C residu lié à la gravité.
  14. C
  15. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  16. C
  17. C AUTEUR : A. BECCANTINI, DEN/DM2S/SFME/LTMF
  18. C
  19. C************************************************************************
  20. C
  21. C ENTREES :
  22. C
  23. C ILIINC : LISTMOTS, noms des inconnues
  24. C
  25. C ICEN : SPG geometrique
  26. C
  27. C IRO : CHPOINT, densité
  28. C
  29. C IROVIT : CHPOINT qdm
  30. C
  31. C IGRAV : CHPOINT gravité
  32. C
  33. C
  34. C SORTIES : IRES : CHPOINT residu
  35. C
  36. C
  37. C************************************************************************
  38. C
  39. C HISTORIQUE (Anomalies et modifications éventuelles)
  40. C
  41. C HISTORIQUE : Créée le 24.1.03
  42. C
  43. C************************************************************************
  44. C
  45. C
  46. C**** Variables de COOPTIO
  47. C
  48. C INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  49. C & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ
  50. C & ,IOPER, IOSGB, IOGRA, IOSAU, IORES
  51. C & ,IECHO, IIMPI, IOSPI
  52. C & ,IDIM
  53. C & ,MCOORD
  54. C & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  55. C & ,NGMAXY, IZROSF, ISOTYP, IOSCR, LTEXLU
  56. C & ,NORINC, NORVAL, NORIND, NORVAD
  57. C & ,NUCROU, IPSAUV, IFICLE, IPREFI
  58. C
  59. C**** Les variables
  60. C
  61. IMPLICIT INTEGER(I-N)
  62. INTEGER ILIINC,ICEN,IRO,IROVIT,IGRAV,IRES,N1,NLCE
  63. & ,IGEOM
  64. CHARACTER*8 TYPE
  65. REAL*8 RO, RVX, RVY, RVZ, GX, GY, GZ
  66. C
  67. C**** Les includes
  68. C
  69. -INC CCOPTIO
  70. -INC SMCHPOI
  71. -INC SMELEME
  72. C
  73. IPT1 = ICEN
  74. SEGACT IPT1
  75. N1 = IPT1.NUM(/2)
  76. SEGDES IPT1
  77. C
  78. C**** Creation de CHPOINT IRES
  79. C
  80. TYPE=' '
  81. CALL KRCHP1(TYPE,ICEN,IRES,ILIINC)
  82. C
  83. C
  84. C**** Lecture de MPOVALs
  85. C
  86. CALL LICHT(IRO,MPOVA1,TYPE,IGEOM)
  87. CALL LICHT(IROVIT,MPOVA2,TYPE,IGEOM)
  88. CALL LICHT(IGRAV,MPOVA3,TYPE,IGEOM)
  89. CALL LICHT(IRES,MPOVA4,TYPE,IGEOM)
  90. C
  91. C**** LICHT active les MPOVALs en *MOD
  92. C
  93. C i.e.
  94. C
  95. C SEGACT MPOVA1*MOD
  96. C SEGACT MPOVA2*MOD
  97. C SEGACT MPOVA3*MOD
  98. C SEGACT MPOVA4*MOD
  99. C
  100. C
  101. GZ=0.0D0
  102. RVZ=0.0D0
  103. DO NLCE = 1, N1, 1
  104. C
  105. C******* Les differents variables a chaque centre
  106. C
  107. RO = MPOVA1.VPOCHA(NLCE,1)
  108. RVX = MPOVA2.VPOCHA(NLCE,1)
  109. RVY = MPOVA2.VPOCHA(NLCE,2)
  110. GX = MPOVA3.VPOCHA(NLCE,1)
  111. GY = MPOVA3.VPOCHA(NLCE,2)
  112. IF(IDIM .EQ. 3)THEN
  113. RVZ = MPOVA2.VPOCHA(NLCE,3)
  114. GZ = MPOVA3.VPOCHA(NLCE,3)
  115. ENDIF
  116. MPOVA4.VPOCHA(NLCE,1) = 0.0D0
  117. MPOVA4.VPOCHA(NLCE,2) = RO*GX
  118. MPOVA4.VPOCHA(NLCE,3) = RO*GY
  119. IF(IDIM .EQ. 3)THEN
  120. MPOVA4.VPOCHA(NLCE,4) = RO*GZ
  121. ENDIF
  122. MPOVA4.VPOCHA(NLCE,IDIM+2)=(RVX*GX)+(RVY*GY)+(RVZ*GZ)
  123. ENDDO
  124. C
  125. SEGDES MPOVA1
  126. SEGDES MPOVA2
  127. SEGDES MPOVA3
  128. SEGDES MPOVA4
  129. C
  130. RETURN
  131. END
  132.  
  133.  
  134.  
  135.  

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