Esope source files

   1 : C GTSINO    SOURCE    GOUNAND   05/12/21    21:29:31     5281           
   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 : -INC CCOPTIO
  39 : CBEGININCLUDE SPOGAU
  40 :       SEGMENT POGAU
  41 :       CHARACTER*(LNNPG)  NOMPG
  42 :       CHARACTER*(LNTPG)  TYPMPG
  43 :       CHARACTER*(LNFPG)  FORLPG
  44 :       INTEGER NORDPG
  45 :       REAL*8 XCOPG(NDLPG,NBPG)
  46 :       REAL*8 XPOPG(NBPG)
  47 :       ENDSEGMENT
  48 :       SEGMENT POGAUS
  49 :       POINTEUR LISPG(0).POGAU
  50 :       ENDSEGMENT
  51 : CENDINCLUDE SPOGAU
  52 :       POINTEUR PGCOUR.POGAU
  53 : *
  54 :       INTEGER NOPG,DIMSRF
  55 :       REAL*8 XCOR(DIMSRF+1)
  56 :       REAL*8 POIDS
  57 :       INTEGER IMPR,IRET
  58 : *
  59 :       INTEGER ICOOR
  60 : *     
  61 : * Executable statements
  62 : *
  63 :       IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans gtsino'
  64 : *
  65 : * On ajoute le point
  66 : *
  67 :       NOPG=NOPG+1
  68 :       DO 1 ICOOR=1,DIMSRF
  69 :          PGCOUR.XCOPG(ICOOR,NOPG)=XCOR(ICOOR+1)
  70 :  1    CONTINUE
  71 :       PGCOUR.XPOPG(NOPG)=POIDS
  72 : *
  73 : * Normal termination
  74 : *
  75 :       IRET=0
  76 :       RETURN
  77 : *
  78 : * Format handling
  79 : *
  80 : *
  81 : * Error handling
  82 : *
  83 :  9999 CONTINUE
  84 :       IRET=1
  85 :       WRITE(IOIMP,*) 'An error was detected in subroutine gtsino'
  86 :       RETURN
  87 : *
  88 : * End of subroutine GTSINO
  89 : *
  90 :       END
  91 : 
  92 :  
  93 :