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. -INC CCOPTIO
  39. CBEGININCLUDE SPOGAU
  40. SEGMENT POGAU
  41. CHARACTER*(LNNPG) NOMPG
  42. CHARACTER*(LNTPG) TYPMPG
  43. CHARACTER*(LNFPG) FORLPG
  44. INTEGER NORDPG
  45. REAL*8 XCOPG(NDLPG,NBPG)
  46. REAL*8 XPOPG(NBPG)
  47. ENDSEGMENT
  48. SEGMENT POGAUS
  49. POINTEUR LISPG(0).POGAU
  50. ENDSEGMENT
  51. CENDINCLUDE SPOGAU
  52. POINTEUR PGCOUR.POGAU
  53. *
  54. INTEGER NOPG,DIMSRF
  55. REAL*8 XCOR(DIMSRF+1)
  56. REAL*8 POIDS
  57. INTEGER IMPR,IRET
  58. *
  59. INTEGER ICOOR
  60. *
  61. * Executable statements
  62. *
  63. IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans gtsino'
  64. *
  65. * On ajoute le point
  66. *
  67. NOPG=NOPG+1
  68. DO 1 ICOOR=1,DIMSRF
  69. PGCOUR.XCOPG(ICOOR,NOPG)=XCOR(ICOOR+1)
  70. 1 CONTINUE
  71. PGCOUR.XPOPG(NOPG)=POIDS
  72. *
  73. * Normal termination
  74. *
  75. IRET=0
  76. RETURN
  77. *
  78. * Format handling
  79. *
  80. *
  81. * Error handling
  82. *
  83. 9999 CONTINUE
  84. IRET=1
  85. WRITE(IOIMP,*) 'An error was detected in subroutine gtsino'
  86. RETURN
  87. *
  88. * End of subroutine GTSINO
  89. *
  90. END
  91.  
  92.  
  93.  
  94.  

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