1 : C INILAG SOURCE GOUNAND 05/12/21 21:32:55 5281 2 : SUBROUTINE INILAG(ELCOUR,IMPR,IRET) 3 : IMPLICIT REAL*8 (A-H,O-Z) 4 : IMPLICIT INTEGER (I-N) 5 : C*********************************************************************** 6 : C NOM : INILAG 7 : C DESCRIPTION : Initialise les tableaux QUENOD et ORDER reliant 8 : C ddl <-> Noeud + ordre de dérivation de la variable 9 : C pour un élément de Lagrange. 10 : C 11 : C LANGAGE : ESOPE 12 : C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF) 13 : C mél : gounand@semt2.smts.cea.fr 14 : C*********************************************************************** 15 : C APPELES : - 16 : C APPELE PAR : INELSE, INELTR, INELQU, INELTE, INELPR, INELCU 17 : C*********************************************************************** 18 : C ENTREES : - 19 : C ENTREES/SORTIES : ELCOUR (supposé actif en *MOD) 20 : C SORTIES : - 21 : C CODE RETOUR (IRET) : = 0 si tout s'est bien passé 22 : C*********************************************************************** 23 : C VERSION : v1, 09/05/2000, version initiale 24 : C HISTORIQUE : v1, 09/05/2000, création 25 : C HISTORIQUE : 26 : C HISTORIQUE : 27 : C*********************************************************************** 28 : C Prière de PRENDRE LE TEMPS de compléter les commentaires 29 : C en cas de modification de ce sous-programme afin de faciliter 30 : C la maintenance ! 31 : C*********************************************************************** 32 : -INC CCOPTIO 33 : CBEGININCLUDE SELREF 34 : SEGMENT ELREF 35 : CHARACTER*(LNNOM) NOMLRF 36 : CHARACTER*(LNFORM) FORME 37 : CHARACTER*(LNTYPL) TYPEL 38 : CHARACTER*(LNESP) ESPACE 39 : INTEGER DEGRE 40 : REAL*8 XCONOD(NDIMEL,NBNOD) 41 : INTEGER NPQUAF(NBDDL) 42 : INTEGER NUMCMP(NBDDL) 43 : INTEGER QUENOD(NBDDL) 44 : INTEGER ORDDER(NDIMEL,NBDDL) 45 : POINTEUR MBPOLY.POLYNS 46 : ENDSEGMENT 47 : SEGMENT ELREFS 48 : POINTEUR LISEL(0).ELREF 49 : ENDSEGMENT 50 : CENDINCLUDE SELREF 51 : POINTEUR ELCOUR.ELREF 52 : * 53 : INTEGER IMPR,IRET 54 : * 55 : INTEGER NDDL,NDEL 56 : INTEGER IDDL,IDEL 57 : * 58 : * Executable statements 59 : * 60 : IF (IMPR.GT.2) WRITE(IOIMP,*) 'Entrée dans inilag.eso' 61 : NDDL=ELCOUR.QUENOD(/1) 62 : NDEL=ELCOUR.ORDDER(/1) 63 : * Chaque ddl est associé à un noeud+la valeur de la variable non dérivée 64 : DO 1 IDDL=1,NDDL 65 : ELCOUR.QUENOD(IDDL)=IDDL 66 : DO 12 IDEL=1,NDEL 67 : ELCOUR.ORDDER(IDEL,IDDL)=0 68 : 12 CONTINUE 69 : 1 CONTINUE 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 inilag' 84 : RETURN 85 : * 86 : * End of subroutine INILAG 87 : * 88 : END 89 : 90 : 91 :
© Cast3M 2003 - All rights reserved.
Disclaimer