1 : C DFNFF1 SOURCE GOUNAND 06/08/04 21:15:24 5520 2 : SUBROUTINE DFNFF1(NDDL,IESREF,IESREL,NBPOGO,NBELEV,NBELEF,NBELFV, 3 : $ NPJMI, 4 : $ DFFPG,JMIJAC,SSFACT, 5 : $ JDFFPG, 6 : $ IMPR,IRET) 7 : IMPLICIT REAL*8 (A-H,O-Z) 8 : IMPLICIT INTEGER (I-N) 9 : C*********************************************************************** 10 : C NOM : DFNFF1 11 : C DESCRIPTION : 12 : C 13 : C 14 : C 15 : C LANGAGE : Fortran 77 (sauf E/S) 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 : DFNFR 21 : C*********************************************************************** 22 : C ENTREES : * NDDL (type entier) : nombre de ddl (par 23 : C élément) de la transformation géométrique. 24 : C * IESREF (type entier) : dimension de l'espace de 25 : C référence. 26 : C * IESREL (type entier) : dimension de l'espace 27 : C réel (i.e. géométrique). 28 : C * NBPOGO (type entier) : nombre de points 29 : C d'intégration. 30 : C * NBELEM (type entier) : nombre d'éléments du 31 : C maillage élémentaire courant. 32 : C * DFFPG (type réel) : valeurs des dérivées 33 : C premières des fonctions d'interpolation aux 34 : C points de gauss sur l'élément de référence. 35 : C * JMIJAC (type réel) : valeurs de l'inverse de 36 : C la matrice jacobienne aux points de Gauss sur 37 : C le maillage élémentaire. 38 : C ENTREES/SORTIES : * JDFFPG (type réel) : valeurs des dérivées 39 : C premières des fonctions d'interpolation 40 : C aux points de gauss sur l'élément réel. 41 : C SORTIES : - 42 : C*********************************************************************** 43 : C VERSION : v2, 03/10/03, refonte complète (modif SMPOUET) 44 : C VERSION : v1, 17/01/03, version initiale 45 : C HISTORIQUE : v1, 17/01/03, création 46 : C HISTORIQUE : 47 : C HISTORIQUE : 48 : C*********************************************************************** 49 : C Prière de PRENDRE LE TEMPS de compléter les commentaires 50 : C en cas de modification de ce sous-programme afin de faciliter 51 : C la maintenance ! 52 : C*********************************************************************** 53 : -INC CCOPTIO 54 : INTEGER NDDL,IESREF,IESREL,NBPOGO,NBELEV,NBELEF,NBELFV,NPJMI 55 : REAL*8 DFFPG (NDDL,IESREF,NBPOGO,NBELFV) 56 : REAL*8 JMIJAC(IESREF,IESREL,NPJMI,NBELEF) 57 : LOGICAL SSFACT(NBELFV,NBELEV) 58 : REAL*8 JDFFPG(NDDL,IESREL,NBPOGO,NBELEF) 59 : * 60 : INTEGER IMPR,IRET 61 : * 62 : INTEGER IDDL,ISREF,ISREL,IBPOGO,IBELEV,IBELEF,IBELFV,IPJMI 63 : * 64 : * Executable statements 65 : * 66 : IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans dfnff1' 67 : IBELEF=0 68 : DO 1 IBELEV=1,NBELEV 69 : DO IBELFV=1,NBELFV 70 : IF (SSFACT(IBELFV,IBELEV)) THEN 71 : IBELEF=IBELEF+1 72 : DO 12 IBPOGO=1,NBPOGO 73 : IF (NPJMI.EQ.1) THEN 74 : IPJMI=1 75 : ELSE 76 : IPJMI=IBPOGO 77 : ENDIF 78 : DO 122 ISREL=1,IESREL 79 : DO 1222 ISREF=1,IESREF 80 : DO 12222 IDDL=1,NDDL 81 : JDFFPG(IDDL,ISREL,IBPOGO,IBELEF)= 82 : $ JDFFPG(IDDL,ISREL,IBPOGO,IBELEF)+ 83 : $ (DFFPG(IDDL,ISREF,IBPOGO,IBELFV) 84 : $ *JMIJAC(ISREF,ISREL,IPJMI,IBELEF)) 85 : 12222 CONTINUE 86 : 1222 CONTINUE 87 : 122 CONTINUE 88 : 12 CONTINUE 89 : ENDIF 90 : ENDDO 91 : 1 CONTINUE 92 : 93 : * 94 : * Normal termination 95 : * 96 : IRET=0 97 : RETURN 98 : * 99 : * Format handling 100 : * 101 : * 102 : * Error handling 103 : * 104 : 9999 CONTINUE 105 : IRET=1 106 : WRITE(IOIMP,*) 'An error was detected in subroutine dfnff1' 107 : RETURN 108 : * 109 : * End of subroutine DFNFF1 110 : * 111 : END 112 : 113 :
© Cast3M 2003 - All rights reserved.
Disclaimer