ccgdiv
C CCGDIV SOURCE GOUNAND 21/06/02 21:15:11 11022 $ FC, $ IMPR,IRET) IMPLICIT REAL*8 (A-H,O-Z) IMPLICIT INTEGER (I-N) C*********************************************************************** C NOM : CCGDIV C DESCRIPTION : Calcul de la loi de comportement aux points de Gauss : C une divergence C 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 : C*********************************************************************** C ENTREES : C ENTREES/SORTIES : C SORTIES : - C TRAVAIL : C*********************************************************************** C VERSION : v1, 11/05/04, version initiale C HISTORIQUE : v1, 11/05/04, 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 SMCHAEL INTEGER NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM,N1 POINTEUR FC.MCHEVA POINTEUR LCOF.LCHEVA POINTEUR MYCOF.MCHEVA * INTEGER IMPR,IRET * * Executable statements * IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans ccgdiv' IF (.NOT.(IDIM.GE.1.AND.IDIM.LE.3)) THEN WRITE(IOIMP,*) 'IDIM=',IDIM,' ?' GOTO 9999 ENDIF NLFC=FC.WELCHE(/6) NPFC=FC.WELCHE(/5) DO IIDIM=1,IDIM MYCOF=LCOF.LISCHE(IIDIM) NLC=MYCOF.WELCHE(/6) NPC=MYCOF.WELCHE(/5) DO ILFC=1,NLFC IF (NLC.EQ.1) THEN ILC=1 ELSE ILC=ILFC ENDIF DO IPFC=1,NPFC IF (NPC.EQ.1) THEN IPC=1 ELSE IPC=IPFC ENDIF FC.WELCHE(1,1,1,1,IPFC,ILFC)= $ FC.WELCHE(1,1,1,1,IPFC,ILFC)+ $ MYCOF.WELCHE(1,1,1,IIDIM,IPC,ILC) ENDDO ENDDO ENDDO * * Normal termination * IRET=0 RETURN * * Format handling * * * Error handling * 9999 CONTINUE IRET=1 WRITE(IOIMP,*) 'An error was detected in subroutine ccgdiv' RETURN * * End of subroutine CCGDIV * END
© Cast3M 2003 - Tous droits réservés.
Mentions légales