src/tool/Decision/deiped.F

   1       subroutine deiped ( nivmin,
   2      >                    decare, decfac,
   3      >                    aretri, nivtri,
   4      >                    arequa,
   5      >                    facpen, hetpen, filpen,
   6      >                    pesupp, peindi,
   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 PEntaedres - Deraffinement
  31 c                                    --           -
  32 c ______________________________________________________________________
  33 c .        .     .        .                                            .
  34 c .  nom   . e/s . taille .           description                      .
  35 c .____________________________________________________________________.
  36 c . nivmin . e   .   1    . niveau min a ne pas depasser en deraffinemt.
  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 . aretri . e   .nbtrto*3. numeros des 3 aretes des triangles         .
  41 c . nivtri . e   . nbtrto . niveau des triangles                       .
  42 c . arequa . e   .nbquto*4. numeros des 4 aretes des quadrangles       .
  43 c . facpen . e   .nbpecf*5. numeros des faces des pentaedres           .
  44 c . hetpen . e   . nbpeto . historique de l'etat des pentaedres        .
  45 c . filpen . e   . nbpeto . premier fils des pentaedres                .
  46 c . pesupp . e   . nbpeto . support pour les pentaedres                .
  47 c . peindi . e   . nbpeto . valeurs entieres pour les pentaedres       .
  48 c . ulsort . e   .   1    . numero d'unite logique de la liste standard.
  49 c . langue . e   .    1   . langue des messages                        .
  50 c .        .     .        . 1 : francais, 2 : anglais                  .
  51 c . codret . es  .    1   . code de retour des modules                 .
  52 c .        .     .        . 0 : pas de probleme                        .
  53 c .        .     .        . 2 : probleme dans le traitement            .
  54 c ______________________________________________________________________
  55 c
  56 c====
  57 c 0. declarations et dimensionnement
  58 c====
  59 c
  60 c 0.1. ==> generalites
  61 c
  62       implicit none
  63       save
  64 c
  65       character*6 nompro
  66       parameter ( nompro = 'DEIPED' )
  67 c
  68 #include "nblang.h"
  69 c
  70 c 0.2. ==> communs
  71 c
  72 #include "envex1.h"
  73 #include "nombar.h"
  74 #include "nombtr.h"
  75 #include "nombqu.h"
  76 #include "nombpe.h"
  77 #include "impr02.h"
  78 c
  79 c 0.3. ==> arguments
  80 c
  81       integer nivmin
  82       integer decare(0:nbarto), decfac(-nbquto:nbtrto)
  83       integer aretri(nbtrto,3), nivtri(nbtrto)
  84       integer arequa(nbquto,4)
  85       integer facpen(nbpecf,5), hetpen(nbpeto), filpen(nbpeto)
  86       integer pesupp(nbpeto), peindi(nbpeto)
  87 c
  88       integer ulsort, langue, codret
  89 c
  90 c 0.4. ==> variables locales
  91 c
  92       integer areloc, facloc, etat
  93       integer lepent
  94       integer iaux, jaux
  95       integer letria, lequad, fils1
  96 c
  97       integer nbmess
  98       parameter (nbmess = 10 )
  99       character*80 texte(nblang,nbmess)
 100 c ______________________________________________________________________
 101 c
 102 c====
 103 c 1. initialisation
 104 c====
 105 c
 106 c 1.1. ==> Les messages
 107 c
 108 #include "impr01.h"
 109 c
 110 #ifdef _DEBUG_HOMARD_
 111       write (ulsort,texte(langue,1)) 'Entree', nompro
 112       call dmflsh (iaux)
 113 #endif
 114 c
 115 #include "impr05.h"
 116 c
 117       codret = 0
 118 c
 119 c====
 120 c 2. traitement des indicateurs portant sur les pentaedres
 121 c    pour le filtrage sur les niveaux, on tient compte du fait que
 122 c    le niveau d'un pentaedre est identifie a celui de n'importe lequel
 123 c    de ses faces.
 124 c====
 125 c
 126 #ifdef _DEBUG_HOMARD_
 127         write(ulsort,texte(langue,4)) mess14(langue,3,7)
 128 #endif
 129 c
 130 #ifdef _DEBUG_HOMARD_
 131         write(ulsort,texte(langue,6))
 132 #endif
 133 c
 134       iaux = 0
 135 c
 136       do 21 , lepent = 1, nbpeto
 137 c
 138         etat = mod( hetpen(lepent) , 100 )
 139         if ( etat.eq.80 ) then
 140           fils1 = filpen(lepent)
 141           if ( pesupp(fils1)   .ne. 0 .and.
 142      >         pesupp(fils1+1) .ne. 0 .and.
 143      >         pesupp(fils1+2) .ne. 0 .and.
 144      >         pesupp(fils1+3) .ne. 0 .and.
 145      >         pesupp(fils1+4) .ne. 0 .and.
 146      >         pesupp(fils1+5) .ne. 0 .and.
 147      >         pesupp(fils1+6) .ne. 0 .and.
 148      >         pesupp(fils1+7) .ne. 0 ) then
 149             if ( peindi(fils1)   .eq. -1 .and.
 150      >           peindi(fils1+1) .eq. -1 .and.
 151      >           peindi(fils1+2) .eq. -1 .and.
 152      >           peindi(fils1+3) .eq. -1 .and.
 153      >           peindi(fils1+4) .eq. -1 .and.
 154      >           peindi(fils1+5) .eq. -1 .and.
 155      >           peindi(fils1+6) .eq. -1 .and.
 156      >           peindi(fils1+7) .eq. -1 ) then
 157               jaux = facpen(lepent,1)
 158               if ( nivtri(jaux).lt.nivmin ) then
 159                 iaux = iaux + 8
 160               else
 161                 do 22 , facloc = 1, 2
 162                   letria = facpen(lepent,facloc)
 163                   decfac(letria) = -1
 164                   do 23 , areloc = 1, 3
 165                     decare(aretri(letria,areloc)) = -1
 166 cgn       write(ulsort,*) 'reactivation de arete', aretri(letria,areloc)
 167    23             continue
 168    22           continue
 169                 do 24 , facloc = 3, 5
 170                   lequad = facpen(lepent,facloc)
 171                   decfac(-lequad) = -1
 172                   do 25 , areloc = 1, 4
 173                     decare(arequa(lequad,areloc)) = -1
 174 cgn       write(ulsort,*) 'reactivation de arete', arequa(lequad,areloc)
 175    25             continue
 176    24           continue
 177               endif
 178             endif
 179           endif
 180         endif
 181    21 continue
 182 c
 183       if ( iaux.ne.0 ) then
 184         write(ulsort,texte(langue,10))
 185         write(ulsort,texte(langue,4)) mess14(langue,3,7)
 186         write(ulsort,texte(langue,8)) nivmin
 187         write(ulsort,texte(langue,9)) iaux
 188       endif
 189 c
 190 c====
 191 c 3. la fin
 192 c====
 193 c
 194       if ( codret.ne.0 ) then
 195 c
 196 #include "envex2.h"
 197 c
 198       write (ulsort,texte(langue,1)) 'Sortie', nompro
 199       write (ulsort,texte(langue,2)) codret
 200 c
 201       endif
 202 c
 203 #ifdef _DEBUG_HOMARD_
 204       write (ulsort,texte(langue,1)) 'Sortie', nompro
 205       call dmflsh (iaux)
 206 #endif
 207 c
 208       end