Télécharger gtsino.eso

Retour à la liste

Numérotation des lignes :

gtsino
  1. C GTSINO SOURCE GOUNAND 21/06/02 21:16:22 11022
  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. -INC TNLIN
  42. *-INC SPOGAU
  43. POINTEUR PGCOUR.POGAU
  44. *
  45. INTEGER NOPG,DIMSRF
  46. REAL*8 XCOR(DIMSRF+1)
  47. REAL*8 POIDS
  48. INTEGER IMPR,IRET
  49. *
  50. INTEGER ICOOR
  51. *
  52. * Executable statements
  53. *
  54. IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans gtsino'
  55. *
  56. * On ajoute le point
  57. *
  58. NOPG=NOPG+1
  59. DO 1 ICOOR=1,DIMSRF
  60. PGCOUR.XCOPG(ICOOR,NOPG)=XCOR(ICOOR+1)
  61. 1 CONTINUE
  62. PGCOUR.XPOPG(NOPG)=POIDS
  63. *
  64. * Normal termination
  65. *
  66. IRET=0
  67. RETURN
  68. *
  69. * Format handling
  70. *
  71. *
  72. * Error handling
  73. *
  74. 9999 CONTINUE
  75. IRET=1
  76. WRITE(IOIMP,*) 'An error was detected in subroutine gtsino'
  77. RETURN
  78. *
  79. * End of subroutine GTSINO
  80. *
  81. END
  82.  
  83.  
  84.  
  85.  
  86.  

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