src/tool/Decision/deipei.F

   1       subroutine deipei ( decare, decfac,
   2      >                    aretri, pertri,
   3      >                    arequa, perqua,
   4      >                    facpen,
   5      >                    pesupp,
   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 PEntaedres - Initialisation
  30 c                                    --           -
  31 c ______________________________________________________________________
  32 c .        .     .        .                                            .
  33 c .  nom   . e/s . taille .           description                      .
  34 c .____________________________________________________________________.
  35 c . decare .  s  .0:nbarto. decisions des aretes                       .
  36 c . decfac .  s  . -nbquto. decision sur les faces (quad. + tri.)      .
  37 c .        .     . :nbtrto.                                            .
  38 c . aretri . e   .nbtrto*3. numeros des 3 aretes des triangles         .
  39 c . pertri . e   . nbtrto . pere des triangles                         .
  40 c . arequa . e   .nbquto*4. numeros des 4 aretes des quadrangles       .
  41 c . perqua . e   . nbquto . pere des quadrangles                       .
  42 c . facpen . e   .nbpecf*5. numeros des faces des pentaedres           .
  43 c . pesupp . e   . nbpeto . support pour les pentaedres                .
  44 c . ulsort . e   .   1    . numero d'unite logique de la liste standard.
  45 c . langue . e   .    1   . langue des messages                        .
  46 c .        .     .        . 1 : francais, 2 : anglais                  .
  47 c . codret . es  .    1   . code de retour des modules                 .
  48 c .        .     .        . 0 : pas de probleme                        .
  49 c .        .     .        . 2 : probleme dans le traitement            .
  50 c ______________________________________________________________________
  51 c
  52 c====
  53 c 0. declarations et dimensionnement
  54 c====
  55 c
  56 c 0.1. ==> generalites
  57 c
  58       implicit none
  59       save
  60 c
  61       character*6 nompro
  62       parameter ( nompro = 'DEIPEI' )
  63 c
  64 #include "nblang.h"
  65 c
  66 c 0.2. ==> communs
  67 c
  68 #include "envex1.h"
  69 #include "nombar.h"
  70 #include "nombtr.h"
  71 #include "nombqu.h"
  72 #include "nombpe.h"
  73 #include "impr02.h"
  74 c
  75 c 0.3. ==> arguments
  76 c
  77       integer decare(0:nbarto), decfac(-nbquto:nbtrto)
  78       integer aretri(nbtrto,3), pertri(nbtrto)
  79       integer arequa(nbquto,4), perqua(nbquto)
  80       integer facpen(nbpecf,5)
  81       integer pesupp(nbpeto)
  82 c
  83       integer ulsort, langue, codret
  84 c
  85 c 0.4. ==> variables locales
  86 c
  87       integer areloc, facloc
  88       integer laface, lepent, lamere
  89       integer iaux
  90 c
  91       integer nbmess
  92       parameter (nbmess = 10 )
  93       character*80 texte(nblang,nbmess)
  94 c ______________________________________________________________________
  95 c
  96 c====
  97 c 1. initialisation
  98 c====
  99 c
 100 c 1.1. ==> Les messages
 101 c
 102 #include "impr01.h"
 103 c
 104 #ifdef _DEBUG_HOMARD_
 105       write (ulsort,texte(langue,1)) 'Entree', nompro
 106       call dmflsh (iaux)
 107 #endif
 108 c
 109 #include "impr05.h"
 110 c
 111       codret = 0
 112 c
 113 c====
 114 c 2. traitement des indicateurs portant sur les pentaedres
 115 c====
 116 c
 117 #ifdef _DEBUG_HOMARD_
 118         write(ulsort,texte(langue,4)) mess14(langue,3,7)
 119 #endif
 120 c
 121       do 21 , lepent = 1, nbpeto
 122 c
 123         if ( pesupp(lepent).ne.0 ) then
 124 c
 125           do 210 , facloc = 1, 5
 126 c
 127             laface = facpen(lepent,facloc)
 128 c
 129 c 2.1. ==> Face triangulaire
 130 c
 131             if ( facloc.le.2 ) then
 132 c
 133 c 2.1.1. ==> Inhibition du raffinement par defaut : on garde la face
 134 c            designee
 135 c
 136               decfac(laface) = 0
 137               do 211 , areloc = 1, 3
 138                 decare(aretri(laface,areloc)) = 0
 139   211         continue
 140 c
 141 c 2.1.2. ==> Inhibition du deraffinement par defaut : on garde la mere
 142 c            de la face designee s'il existe
 143 c
 144               lamere = pertri(laface)
 145 c
 146               if ( lamere.gt.0 ) then
 147 c
 148                 decfac(lamere) = 0
 149                 do 212 , areloc = 1, 3
 150                   decare(aretri(lamere,areloc)) = 0
 151   212           continue
 152 c
 153               endif
 154 c
 155 c 2.2. ==> Face quadrangulaire
 156 c
 157             else
 158 c
 159 c 2.2.1. ==> Inhibition du raffinement par defaut : on garde la face
 160 c            designee
 161 c
 162               decfac(-laface) = 0
 163               do 221 , areloc = 1, 4
 164                 decare(arequa(laface,areloc)) = 0
 165   221         continue
 166 c
 167 c 2.2.2. ==> Inhibition du deraffinement par defaut : on garde la mere
 168 c            de la face designee s'il existe
 169 c
 170               lamere = perqua(laface)
 171 c
 172               if ( lamere.gt.0 ) then
 173 c
 174                 decfac(-lamere) = 0
 175                 do 222 , areloc = 1, 4
 176                   decare(arequa(lamere,areloc)) = 0
 177   222           continue
 178 c
 179               endif
 180 c
 181             endif
 182 c
 183   210     continue
 184 c
 185         endif
 186 c
 187    21 continue
 188 c
 189 c====
 190 c 3. la fin
 191 c====
 192 c
 193       if ( codret.ne.0 ) then
 194 c
 195 #include "envex2.h"
 196 c
 197       write (ulsort,texte(langue,1)) 'Sortie', nompro
 198       write (ulsort,texte(langue,2)) codret
 199 c
 200       endif
 201 c
 202 #ifdef _DEBUG_HOMARD_
 203       write (ulsort,texte(langue,1)) 'Sortie', nompro
 204       call dmflsh (iaux)
 205 #endif
 206 c
 207       end