Télécharger la plaquette Cast3M


fr en

Télécharger gcpasy.eso

Retour à la liste

Numérotation des lignes : 

gcpasy

  1. C GCPASY    SOURCE    GOUNAND   21/06/02    21:16:03     11022          
  2.       SUBROUTINE GCPASY(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         : GCPASY
  8. C PROJET      : Noyau linéaire NLIN
  9. C DESCRIPTION : Rajoute un point et son image par une symétrie par
  10. C               rapport à l'axe des x en 2D dans une méthode
  11. C               d'intégration type Gauss (PGCOUR).
  12. C               Domaine de type [ Fully symmetric region ]
  13. C                       i.e. segment, carré et cube pour nous.
  14. C               Générateur de type [ Partial symmetry ].
  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       : INGAQU
  22. C***********************************************************************
  23. C ENTREES            : DIMSRF, XCOR, POIDS
  24. C ENTREES/SORTIES    : PGCOUR (actif en *MOD), NOPG
  25. C SORTIES            : MYPGS
  26. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  27. C***********************************************************************
  28. C VERSION    : v1, 20/10/99, version initiale
  29. C HISTORIQUE : v1, 20/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. * Executable statements
  50. *
  51.       IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans gcpasy'
  52.       IF (DIMSRF.NE.2) THEN
  53.          WRITE(IOIMP,*) 'DIMSRF doit etre égal à 2...'
  54.          GOTO 9999
  55.       ENDIF
  56. *
  57. * On ajoute le point...
  58. *
  59.       NOPG=NOPG+1
  60.       PGCOUR.XCOPG(1,NOPG)= XCOR(1)
  61.       PGCOUR.XCOPG(2,NOPG)= XCOR(2)
  62.       PGCOUR.XPOPG(NOPG)=POIDS
  63. *
  64. * ...et son image par la symétrie / à l'axe des x
  65. *
  66.       NOPG=NOPG+1
  67.       PGCOUR.XCOPG(1,NOPG)= XCOR(1)
  68.       PGCOUR.XCOPG(2,NOPG)=-XCOR(2)
  69.       PGCOUR.XPOPG(NOPG)=POIDS
  70. *
  71. * Normal termination
  72. *
  73.       IRET=0
  74.       RETURN
  75. *
  76. * Format handling
  77. *
  78. *
  79. * Error handling
  80. *
  81.  9999 CONTINUE
  82.       IRET=1
  83.       WRITE(IOIMP,*) 'An error was detected in subroutine gcpasy'
  84.       RETURN
  85. *
  86. * End of subroutine GCPASY
  87. *
  88.       END
  89.  
  90.  
  91.  
  92.  
  93.  

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