Télécharger la plaquette Cast3M


fr en

Télécharger gccesy.eso

Retour à la liste

Numérotation des lignes : 

gccesy

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

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