gtsino
C GTSINO SOURCE GOUNAND 21/06/02 21:16:22 11022 $ IMPR,IRET) IMPLICIT REAL*8 (A-H,O-Z) IMPLICIT INTEGER (I-N) C*********************************************************************** C NOM : GTSINO C PROJET : Noyau linéaire NLIN C DESCRIPTION : Rajoute un point dans une méthode d'intégration C type Gauss (PGCOUR). C Domaine de type [ Simplex ] C i.e. triangle et tétrahèdre pour nous C Générateur de type [ Single node ]. C C XCOR sont les coordonnées barycentriques... C C LANGAGE : ESOPE C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF) C mél : gounand@semt2.smts.cea.fr C*********************************************************************** C APPELES : - C APPELE PAR : INGATR, INGATE C*********************************************************************** C ENTREES : DIMSRF, XCOR, POIDS C ENTREES/SORTIES : PGCOUR (actif en *MOD), NOPG C SORTIES : - C CODE RETOUR (IRET) : = 0 si tout s'est bien passé C*********************************************************************** C VERSION : v1, 21/10/99, version initiale C HISTORIQUE : v1, 21/10/99, création C HISTORIQUE : C HISTORIQUE : C*********************************************************************** C Prière de PRENDRE LE TEMPS de compléter les commentaires C en cas de modification de ce sous-programme afin de faciliter C la maintenance ! C*********************************************************************** -INC PPARAM -INC CCOPTIO -INC TNLIN *-INC SPOGAU POINTEUR PGCOUR.POGAU * INTEGER NOPG,DIMSRF REAL*8 XCOR(DIMSRF+1) REAL*8 POIDS INTEGER IMPR,IRET * INTEGER ICOOR * * Executable statements * IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans gtsino' * * On ajoute le point * NOPG=NOPG+1 DO 1 ICOOR=1,DIMSRF PGCOUR.XCOPG(ICOOR,NOPG)=XCOR(ICOOR+1) 1 CONTINUE PGCOUR.XPOPG(NOPG)=POIDS * * Normal termination * IRET=0 RETURN * * Format handling * * * Error handling * 9999 CONTINUE IRET=1 WRITE(IOIMP,*) 'An error was detected in subroutine gtsino' RETURN * * End of subroutine GTSINO * END
© Cast3M 2003 - Tous droits réservés.
Mentions légales