1 : C VERMAI SOURCE MAUGIS 05/08/05 21:15:14 5145 2 : SUBROUTINE VERMAI 3 : C---------------------------------------------------- 4 : C 5 : C Vérification de la conformité d'un maillage 6 : C 7 : C---------------------------------------------------- 8 : C 9 : C Création : ??? 10 : C 11 : C Modification : Pascal Maugis 04/08/2005 12 : C ajout commentaires 13 : C Identification des doublons/triplons,... 14 : C si mot-clef 'CORR'IGE, vire les doublons 15 : C 16 : C---------------------------------------------------- 17 : C 18 : C Remarques 19 : C 20 : C---------------------------------------------------- 21 : C 22 : C Variables 23 : C 24 : C---------------------------------------------------- 25 : IMPLICIT INTEGER(I-N) 26 : IMPLICIT REAL*8(A-H,O-Z) 27 : -INC CCOPTIO 28 : -INC SMELEME 29 : CHARACTER*4 MCLE(1) 30 : DATA MCLE/'CORR'/ 31 : 32 : * Lecture maillage 32 : CALL LIROBJ ('MAILLAGE',MELEME,1,IRETOU) 33 : IF(IERR.NE.0) RETURN 34 : 35 : C Y a-t-il des mailles en double ? 35 : CALL LIRMOT(MCLE,1,ICLE,0) 36 : IF (IERR.NE.0) RETURN 37 : IF (ICLE.EQ.0) THEN 38 : C diagnostique seul 38 : CALL DOUBL1(MELEME,NMESH) 39 : ELSE 40 : C on vire les doublons en mode bavard 40 : CALL DOUBL2(MELEME,.TRUE.) 41 : ENDIF 42 : 43 : C vérifie qu'il n'y a pas d'éléments de degré un accolé à un 44 : C élément de degré 2. 45 : IVO=MELEME 45 : CALL VTOP2D(IVO) 46 : 47 : C y a-t-il des éléments connectés à des éléments d'une autre forme. 48 : IF((ICLE.EQ.0).AND.(NMESH.NE.0)) THEN 49 : write(IOIMP,*) 'Presence de doublons, on saute le 3eme test' 50 : ELSE 51 : C on travaille sur le maillage corrigé sans doublon 51 : CALL VMDISC(MELEME) 52 : ENDIF 53 : 54 : RETURN 55 : END 56 : 57 :
© Cast3M 2003 - All rights reserved.
Disclaimer