Télécharger la plaquette Cast3M


fr en

Télécharger fimax2.eso

Retour à la liste

Numérotation des lignes : 

fimax2

  1. C FIMAX2    SOURCE    OF166741  24/12/13    21:15:44     12097          
  2.       SUBROUTINE FIMAX2(ICEN,IX2D,IY2D,IS2D,IVOL,IPN,IGPN,ILGPN,IRES)
  3. C************************************************************************
  4. C
  5. C PROJET            :  CASTEM 2000
  6. C
  7. C NOM               :  FIMAX2
  8. C
  9. C DESCRIPTION       :  Voir FIMAXI
  10. C
  11. C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  12. C
  13. C AUTEUR            :  A. BECCANTINI, LTMF
  14. C
  15. C************************************************************************
  16. C
  17. C
  18. C APPELES (Outils)  :  AUCUN
  19. C
  20. C APPELES (Calcul)  :  AUCUN
  21. C
  22. C
  23. C************************************************************************
  24. C
  25. C ENTREES
  26. C
  27. C
  28. C     1) Pointeurs de MELEMEs et de CHPOINTs "CENTRE" de la table DOMAINE
  29. C
  30. C     ICEN      : MELEME de 'POI1' SPG des CENTRES
  31. C
  32. C     IX2D,IY2D : coordonnés des centres 2D
  33. C
  34. C     IS2D      : surface 2D
  35. C     IVOL      : volume
  36. C
  37. C     2) CHPOINTs "CENTRE"
  38. C
  39. C     IPN,IGPN,ILGPN : pression, gradient, limiteur
  40. C
  41. C
  42. C SORTIES
  43. C
  44. C     IRES    : contribution au residu
  45. C
  46. C
  47. C************************************************************************
  48. C
  49. C HISTORIQUE (Anomalies et modifications éventuelles)
  50. C
  51. C HISTORIQUE :  Créée le 25.2.2004
  52. C
  53. C
  54. C
  55. C************************************************************************
  56. C
  57.       IMPLICIT INTEGER(I-N)
  58. C
  59.       INTEGER ICEN,IX2D,IY2D,IS2D,IVOL,IPN,IGPN,ILGPN,IRES
  60.      &     , NBEL, IELEM, NP, INDP, IGEOM
  61.       REAL*8 DX, DY, GX, GY, GEOM, AL, P
  62.       CHARACTER*8 TYPE
  63. C
  64. C**** Les Includes
  65. C
  66. -INC SMCOORD
  67.  
  68. -INC PPARAM
  69. -INC CCOPTIO
  70. -INC SMCHPOI
  71.       POINTEUR MP.MPOVAL,     MGRP.MPOVAL,
  72.      &         MALP.MPOVAL,   MPXC2D.MPOVAL,
  73.      &         MPYC2D.MPOVAL, MPSU2D.MPOVAL,
  74.      &         MPVOLU.MPOVAL, MPRES.MPOVAL
  75. -INC SMELEME
  76. C
  77.       CALL LICHT(IPN,MP,TYPE,IGEOM)
  78.       CALL LICHT(IGPN,MGRP,TYPE,IGEOM)
  79.       CALL LICHT(ILGPN,MALP,TYPE,IGEOM)
  80.       CALL LICHT(IX2D,MPXC2D,TYPE,IGEOM)
  81.       CALL LICHT(IY2D,MPYC2D,TYPE,IGEOM)
  82.       CALL LICHT(IS2D,MPSU2D,TYPE,IGEOM)
  83.       CALL LICHT(IVOL,MPVOLU,TYPE,IGEOM)
  84.       CALL LICHT(IRES,MPRES,TYPE,IGEOM)
  85. C
  86. C     SEGACT   MP
  87. C     SEGACT   MGRP
  88. C     SEGACT   MALP
  89. C     SEGACT   MPXC2D
  90. C     SEGACT   MPYC2D
  91. C     SEGACT   MPSU2D
  92. C     SEGACT   MPVOLU
  93. C     SEGACT   MPRES
  94. C
  95.       MELEME=ICEN
  96.       SEGACT MELEME
  97.       NBEL=MELEME.NUM(/2)
  98.       DO IELEM=1,NBEL,1
  99.          NP=MELEME.NUM(1,IELEM)
  100.          INDP = (NP-1)*3
  101.          DX=MPXC2D.VPOCHA(IELEM,1) - XCOOR(INDP+1)
  102.          DY=MPYC2D.VPOCHA(IELEM,1) - XCOOR(INDP+2)
  103.          GX=MGRP.VPOCHA(IELEM,1)
  104.          GY=MGRP.VPOCHA(IELEM,2)
  105.          GEOM=MPSU2D.VPOCHA(IELEM,1)/MPVOLU.VPOCHA(IELEM,1)
  106.          AL=MALP.VPOCHA(IELEM,1)
  107.          P=MP.VPOCHA(IELEM,1)
  108.          MPRES.VPOCHA(IELEM,1)=(P+(AL*((GX*DX)+(GY*DY))))*GEOM
  109.       ENDDO
  110. C
  111.       SEGDES  MELEME
  112.       SEGDES  MP
  113.       SEGDES  MGRP
  114.       SEGDES  MALP
  115.       SEGDES  MPXC2D
  116.       SEGDES  MPYC2D
  117.       SEGDES  MPSU2D
  118.       SEGDES  MPVOLU
  119.       SEGDES  MPRES
  120. C
  121.  9999 CONTINUE
  122. C
  123.       RETURN
  124.       END
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  

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