Télécharger gccesy.eso

Retour à la liste

Numérotation des lignes :

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

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