src/tool/Decision/deiher.F

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