Download vermai.eso

Back to the list

   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