Télécharger la plaquette Cast3M


fr en

Télécharger fimvf2.eso

Retour à la liste

Numérotation des lignes : 

fimvf2

  1. C FIMVF2    SOURCE    OF166741  24/12/13    21:15:48     12097          
  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**** Les variables
  46. C
  47.       IMPLICIT INTEGER(I-N)
  48.       INTEGER ILIINC,ICEN,IRO,IROVIT,IGRAV,IRES,N1,NLCE
  49.      &     ,IGEOM
  50.       CHARACTER*8 TYPE
  51.       REAL*8 RO, RVX, RVY, RVZ, GX, GY, GZ
  52. C
  53. C**** Les includes
  54. C
  55.  
  56. -INC PPARAM
  57. -INC CCOPTIO
  58. -INC SMCHPOI
  59. -INC SMELEME
  60. C
  61.       IPT1 = ICEN
  62.       SEGACT IPT1
  63.       N1 = IPT1.NUM(/2)
  64.       SEGDES IPT1
  65. C
  66. C**** Creation de CHPOINT IRES
  67. C
  68.       TYPE='        '
  69.       CALL KRCHP1(TYPE,ICEN,IRES,ILIINC)
  70. C
  71. C
  72. C**** Lecture de MPOVALs
  73. C
  74.       CALL LICHT(IRO,MPOVA1,TYPE,IGEOM)
  75.       CALL LICHT(IROVIT,MPOVA2,TYPE,IGEOM)
  76.       CALL LICHT(IGRAV,MPOVA3,TYPE,IGEOM)
  77.       CALL LICHT(IRES,MPOVA4,TYPE,IGEOM)
  78. C
  79. C**** LICHT active les MPOVALs en *MOD
  80. C
  81. C     i.e.
  82. C
  83. C     SEGACT MPOVA1*MOD
  84. C     SEGACT MPOVA2*MOD
  85. C     SEGACT MPOVA3*MOD
  86. C     SEGACT MPOVA4*MOD
  87. C
  88. C
  89.       GZ=0.0D0
  90.       RVZ=0.0D0
  91.       DO NLCE = 1, N1, 1
  92. C
  93. C******* Les differents variables a chaque centre
  94. C
  95.          RO = MPOVA1.VPOCHA(NLCE,1)
  96.          RVX = MPOVA2.VPOCHA(NLCE,1)
  97.          RVY = MPOVA2.VPOCHA(NLCE,2)
  98.          GX = MPOVA3.VPOCHA(NLCE,1)
  99.          GY = MPOVA3.VPOCHA(NLCE,2)
  100.          IF(IDIM .EQ. 3)THEN
  101.             RVZ = MPOVA2.VPOCHA(NLCE,3)
  102.             GZ = MPOVA3.VPOCHA(NLCE,3)
  103.          ENDIF
  104.          MPOVA4.VPOCHA(NLCE,1) = 0.0D0
  105.          MPOVA4.VPOCHA(NLCE,2) = RO*GX
  106.          MPOVA4.VPOCHA(NLCE,3) = RO*GY
  107.          IF(IDIM .EQ. 3)THEN
  108.             MPOVA4.VPOCHA(NLCE,4) = RO*GZ
  109.          ENDIF
  110.          MPOVA4.VPOCHA(NLCE,IDIM+2)=(RVX*GX)+(RVY*GY)+(RVZ*GZ)
  111.       ENDDO
  112. C
  113.       SEGDES MPOVA1
  114.       SEGDES MPOVA2
  115.       SEGDES MPOVA3
  116.       SEGDES MPOVA4
  117. C
  118.       RETURN
  119.       END
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  

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