src/tool/Decision/deited.F

   1       subroutine deited ( nivmin,
   2      >                    decare, decfac,
   3      >                    aretri, nivtri,
   4      >                    tritet, hettet, filtet,
   5      >                    tesupp, teindi,
   6      >                    ulsort, langue, codret)
   7 c ______________________________________________________________________
   8 c
   9 c                             H O M A R D
  10 c
  11 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
  12 c
  13 c Version originale enregistree le 18 juin 1996 sous le numero 96036
  14 c aupres des huissiers de justice Simart et Lavoir a Clamart
  15 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
  16 c aupres des huissiers de justice
  17 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
  18 c
  19 c    HOMARD est une marque deposee d'Electricite de France
  20 c
  21 c Copyright EDF 1996
  22 c Copyright EDF 1998
  23 c Copyright EDF 2002
  24 c Copyright EDF 2020
  25 c ______________________________________________________________________
  26 c
  27 c traitement des DEcisions - Initialisation de l'indicateur entier
  28 c                --          -
  29 c                          - cas des TEtraedres - Deraffinement
  30 c                                    --           -
  31 c ______________________________________________________________________
  32 c .        .     .        .                                            .
  33 c .  nom   . e/s . taille .           description                      .
  34 c .____________________________________________________________________.
  35 c . nivmin . e   .   1    . niveau min a ne pas depasser en deraffinemt.
  36 c . decare .  s  .0:nbarto. decisions des aretes                       .
  37 c . decfac .  s  . -nbquto. decision sur les faces (quad. + tri.)      .
  38 c .        .     . :nbtrto.                                            .
  39 c . aretri . e   .nbtrto*3. numeros des 3 aretes des triangles         .
  40 c . nivtri . e   . nbtrto . niveau des triangles                       .
  41 c . tritet . e   .nbtecf*4. numeros des 4 triangles des tetraedres     .
  42 c . hettet . e   . nbteto . historique de l'etat des tetraedres        .
  43 c . filtet . e   . nbteto . premier fils des tetraedres                .
  44 c . tesupp . e   . nbteto . support pour les tetraedres                .
  45 c . teindi . e   . nbteto . valeurs entieres pour les tetraedres       .
  46 c . ulsort . e   .   1    . numero d'unite logique de la liste standard.
  47 c . langue . e   .    1   . langue des messages                        .
  48 c .        .     .        . 1 : francais, 2 : anglais                  .
  49 c . codret . es  .    1   . code de retour des modules                 .
  50 c .        .     .        . 0 : pas de probleme                        .
  51 c .        .     .        . 2 : probleme dans le traitement            .
  52 c ______________________________________________________________________
  53 c
  54 c====
  55 c 0. declarations et dimensionnement
  56 c====
  57 c
  58 c 0.1. ==> generalites
  59 c
  60       implicit none
  61       save
  62 c
  63       character*6 nompro
  64       parameter ( nompro = 'DEITED' )
  65 c
  66 #include "nblang.h"
  67 c
  68 c 0.2. ==> communs
  69 c
  70 #include "envex1.h"
  71 #include "nombar.h"
  72 #include "nombtr.h"
  73 #include "nombqu.h"
  74 #include "nombte.h"
  75 #include "impr02.h"
  76 c
  77 c 0.3. ==> arguments
  78 c
  79       integer nivmin
  80       integer decare(0:nbarto), decfac(-nbquto:nbtrto)
  81       integer aretri(nbtrto,3), nivtri(nbtrto)
  82       integer tritet(nbtecf,4), hettet(nbteto), filtet(nbteto)
  83       integer tesupp(nbteto), teindi(nbteto)
  84 c
  85       integer ulsort, langue, codret
  86 c
  87 c 0.4. ==> variables locales
  88 c
  89       integer areloc, facloc, etat
  90       integer letria, letetr
  91       integer fils1
  92       integer iaux, jaux
  93 c
  94       integer nbmess
  95       parameter (nbmess = 30 )
  96       character*80 texte(nblang,nbmess)
  97 c ______________________________________________________________________
  98 c
  99 c====
 100 c 1. initialisation
 101 c====
 102 c
 103 c 1.1. ==> Les messages
 104 c
 105 #include "impr01.h"
 106 c
 107 #ifdef _DEBUG_HOMARD_
 108       write (ulsort,texte(langue,1)) 'Entree', nompro
 109       call dmflsh (iaux)
 110 #endif
 111 c
 112 #include "impr05.h"
 113 #include "derco1.h"
 114 c
 115       codret = 0
 116 c
 117 c====
 118 c 2. traitement des indicateurs portant sur les tetraedres
 119 c    pour le filtrage sur les niveaux, on tient compte du fait que
 120 c    le niveau d'un tetraedre est identifie a celui de l'une quelconque
 121 c    de ses triangles.
 122 c====
 123 c
 124 #ifdef _DEBUG_HOMARD_
 125         write(ulsort,texte(langue,4)) mess14(langue,3,3)
 126 #endif
 127 c
 128 #ifdef _DEBUG_HOMARD_
 129         write(ulsort,texte(langue,6))
 130 #endif
 131 c
 132       iaux = 0
 133 c
 134       do 21 , letetr = 1, nbteto
 135 c
 136         etat = mod( hettet(letetr) , 100 )
 137         if ( etat.ge.85 .and. etat.le.87 ) then
 138           fils1 = filtet(letetr)
 139           if ( tesupp(fils1)   .ne. 0 .and.
 140      >         tesupp(fils1+1) .ne. 0 .and.
 141      >         tesupp(fils1+2) .ne. 0 .and.
 142      >         tesupp(fils1+3) .ne. 0 .and.
 143      >         tesupp(fils1+4) .ne. 0 .and.
 144      >         tesupp(fils1+5) .ne. 0 .and.
 145      >         tesupp(fils1+6) .ne. 0 .and.
 146      >         tesupp(fils1+7) .ne. 0 ) then
 147             if ( teindi(fils1)   .eq. -1 .and.
 148      >           teindi(fils1+1) .eq. -1 .and.
 149      >           teindi(fils1+2) .eq. -1 .and.
 150      >           teindi(fils1+3) .eq. -1 .and.
 151      >           teindi(fils1+4) .eq. -1 .and.
 152      >           teindi(fils1+5) .eq. -1 .and.
 153      >           teindi(fils1+6) .eq. -1 .and.
 154      >           teindi(fils1+7) .eq. -1 ) then
 155               jaux = tritet(letetr,1)
 156               if ( nivtri(jaux).lt.nivmin ) then
 157                 iaux = iaux + 8
 158               else
 159                 do 22 , facloc = 1, 4
 160                   letria = tritet(letetr,facloc)
 161                   decfac(letria) = -1
 162 #ifdef _DEBUG_HOMARD_
 163       write (ulsort,texte(langue,30))'decfac',
 164      >letria,decfac(letria),' '
 165 #endif
 166                   do 23 , areloc = 1, 3
 167                     decare(aretri(letria,areloc)) = -1
 168 #ifdef _DEBUG_HOMARD_
 169       write (ulsort,texte(langue,30))'decare',aretri(letria,areloc),
 170      >  decare(aretri(letria,areloc)),' (une de ses aretes)'
 171 #endif
 172    23             continue
 173    22           continue
 174               endif
 175             endif
 176           endif
 177         endif
 178 c
 179    21 continue
 180 c
 181       if ( iaux.ne.0 ) then
 182         write(ulsort,texte(langue,10))
 183         write(ulsort,texte(langue,4)) mess14(langue,3,3)
 184         write(ulsort,texte(langue,8)) nivmin
 185         write(ulsort,texte(langue,9)) iaux
 186       endif
 187 c
 188 c====
 189 c 3. la fin
 190 c====
 191 c
 192       if ( codret.ne.0 ) then
 193 c
 194 #include "envex2.h"
 195 c
 196       write (ulsort,texte(langue,1)) 'Sortie', nompro
 197       write (ulsort,texte(langue,2)) codret
 198 c
 199       endif
 200 c
 201 #ifdef _DEBUG_HOMARD_
 202       write (ulsort,texte(langue,1)) 'Sortie', nompro
 203       call dmflsh (iaux)
 204 #endif
 205 c
 206       end