Télécharger gcsino.eso

Retour à la liste

Numérotation des lignes :

gcsino
  1. C GCSINO SOURCE GOUNAND 21/06/02 21:16:06 11022
  2. SUBROUTINE GCSINO(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 : GCSINO
  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 [ Fully symmetric region ]
  12. C i.e. segment, carré et cube pour nous
  13. C ou autres (on ne fait qu'ajouter un point)
  14. C Générateur de type [ Single node ].
  15. C
  16. C LANGAGE : ESOPE
  17. C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  18. C mél : gounand@semt2.smts.cea.fr
  19. C***********************************************************************
  20. C APPELES : -
  21. C APPELE PAR : INGASE, INGJ10, INGJ20, INGAQU, INGACU
  22. C***********************************************************************
  23. C ENTREES : DIMSRF, XCOR, POIDS
  24. C ENTREES/SORTIES : PGCOUR (actif en *MOD), NOPG
  25. C SORTIES : -
  26. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  27. C***********************************************************************
  28. C VERSION : v1, 19/10/99, version initiale
  29. C HISTORIQUE : v1, 19/10/99, création
  30. C HISTORIQUE :
  31. C HISTORIQUE :
  32. C***********************************************************************
  33. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  34. C en cas de modification de ce sous-programme afin de faciliter
  35. C la maintenance !
  36. C***********************************************************************
  37.  
  38. -INC PPARAM
  39. -INC CCOPTIO
  40. -INC TNLIN
  41. *-INC SPOGAU
  42. POINTEUR PGCOUR.POGAU
  43. *
  44. INTEGER NOPG,DIMSRF
  45. REAL*8 XCOR(DIMSRF)
  46. REAL*8 POIDS
  47. INTEGER IMPR,IRET
  48. *
  49. INTEGER ICOOR
  50. *
  51. * Executable statements
  52. *
  53. IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans gcsino'
  54. *
  55. * On ajoute le point
  56. NOPG=NOPG+1
  57. DO 1 ICOOR=1,DIMSRF
  58. PGCOUR.XCOPG(ICOOR,NOPG)=XCOR(ICOOR)
  59. 1 CONTINUE
  60. PGCOUR.XPOPG(NOPG)=POIDS
  61. *
  62. * Normal termination
  63. *
  64. IRET=0
  65. RETURN
  66. *
  67. * Format handling
  68. *
  69. *
  70. * Error handling
  71. *
  72. 9999 CONTINUE
  73. IRET=1
  74. WRITE(IOIMP,*) 'An error was detected in subroutine gcsino'
  75. RETURN
  76. *
  77. * End of subroutine GCSINO
  78. *
  79. END
  80.  
  81.  
  82.  
  83.  
  84.  

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