Télécharger gtsino.eso

Retour à la liste

Numérotation des lignes :

  1. C GTSINO SOURCE GOUNAND 05/12/21 21:29:31 5281
  2. SUBROUTINE GTSINO(PGCOUR,NOPG,DIMSRF,XCOR,POIDS,
  3. $ IMPR,IRET)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. IMPLICIT INTEGER (I-N)
  6. C***********************************************************************
  7. C NOM : GTSINO
  8. C PROJET : Noyau linéaire NLIN
  9. C DESCRIPTION : Rajoute un point dans une méthode d'intégration
  10. C type Gauss (PGCOUR).
  11. C Domaine de type [ Simplex ]
  12. C i.e. triangle et tétrahèdre pour nous
  13. C Générateur de type [ Single node ].
  14. C
  15. C XCOR sont les coordonnées barycentriques...
  16. C
  17. C LANGAGE : ESOPE
  18. C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  19. C mél : gounand@semt2.smts.cea.fr
  20. C***********************************************************************
  21. C APPELES : -
  22. C APPELE PAR : INGATR, INGATE
  23. C***********************************************************************
  24. C ENTREES : DIMSRF, XCOR, POIDS
  25. C ENTREES/SORTIES : PGCOUR (actif en *MOD), NOPG
  26. C SORTIES : -
  27. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  28. C***********************************************************************
  29. C VERSION : v1, 21/10/99, version initiale
  30. C HISTORIQUE : v1, 21/10/99, création
  31. C HISTORIQUE :
  32. C HISTORIQUE :
  33. C***********************************************************************
  34. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  35. C en cas de modification de ce sous-programme afin de faciliter
  36. C la maintenance !
  37. C***********************************************************************
  38.  
  39. -INC PPARAM
  40. -INC CCOPTIO
  41. CBEGININCLUDE SPOGAU
  42. SEGMENT POGAU
  43. CHARACTER*(LNNPG) NOMPG
  44. CHARACTER*(LNTPG) TYPMPG
  45. CHARACTER*(LNFPG) FORLPG
  46. INTEGER NORDPG
  47. REAL*8 XCOPG(NDLPG,NBPG)
  48. REAL*8 XPOPG(NBPG)
  49. ENDSEGMENT
  50. SEGMENT POGAUS
  51. POINTEUR LISPG(0).POGAU
  52. ENDSEGMENT
  53. CENDINCLUDE SPOGAU
  54. POINTEUR PGCOUR.POGAU
  55. *
  56. INTEGER NOPG,DIMSRF
  57. REAL*8 XCOR(DIMSRF+1)
  58. REAL*8 POIDS
  59. INTEGER IMPR,IRET
  60. *
  61. INTEGER ICOOR
  62. *
  63. * Executable statements
  64. *
  65. IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans gtsino'
  66. *
  67. * On ajoute le point
  68. *
  69. NOPG=NOPG+1
  70. DO 1 ICOOR=1,DIMSRF
  71. PGCOUR.XCOPG(ICOOR,NOPG)=XCOR(ICOOR+1)
  72. 1 CONTINUE
  73. PGCOUR.XPOPG(NOPG)=POIDS
  74. *
  75. * Normal termination
  76. *
  77. IRET=0
  78. RETURN
  79. *
  80. * Format handling
  81. *
  82. *
  83. * Error handling
  84. *
  85. 9999 CONTINUE
  86. IRET=1
  87. WRITE(IOIMP,*) 'An error was detected in subroutine gtsino'
  88. RETURN
  89. *
  90. * End of subroutine GTSINO
  91. *
  92. END
  93.  
  94.  
  95.  
  96.  

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