src/tool/Decision/deiarr.F

   1       subroutine deiarr ( nivmax,
   2      >                    decare,
   3      >                    hetare,
   4      >                    posifa, facare,
   5      >                    nivtri,
   6      >                    nivqua,
   7      >                    arsupp, arindi,
   8      >                    ulsort, langue, codret)
   9 c ______________________________________________________________________
  10 c
  11 c                             H O M A R D
  12 c
  13 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
  14 c
  15 c Version originale enregistree le 18 juin 1996 sous le numero 96036
  16 c aupres des huissiers de justice Simart et Lavoir a Clamart
  17 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
  18 c aupres des huissiers de justice
  19 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
  20 c
  21 c    HOMARD est une marque deposee d'Electricite de France
  22 c
  23 c Copyright EDF 1996
  24 c Copyright EDF 1998
  25 c Copyright EDF 2002
  26 c Copyright EDF 2020
  27 c ______________________________________________________________________
  28 c
  29 c traitement des DEcisions - Initialisation de l'indicateur entier
  30 c                --          -
  31 c                          - cas des ARetes - Raffinement
  32 c                                    --       -
  33 c ______________________________________________________________________
  34 c .        .     .        .                                            .
  35 c .  nom   . e/s . taille .           description                      .
  36 c .____________________________________________________________________.
  37 c . nivmax . e   .   1    . niveau max a ne pas depasser en raffinement.
  38 c . decare .  s  .0:nbarto. decisions des aretes                       .
  39 c . hetare . e   . nbarto . historique de l'etat des aretes            .
  40 c . posifa . e   . nbarto . pointeur sur tableau facare                .
  41 c . facare . e   . nbfaar . liste des faces contenant une arete        .
  42 c . nivtri . e   . nbtrto . niveau des triangles                       .
  43 c . nivqua . e   . nbquto . niveau des quadrangles                     .
  44 c . arsupp . e   . nbarto . support pour les aretes                    .
  45 c . arindi . e   . nbarto . valeurs entieres pour les aretes           .
  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 = 'DEIARR' )
  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 "impr02.h"
  75 c
  76 c 0.3. ==> arguments
  77 c
  78       integer nivmax
  79       integer decare(0:nbarto)
  80       integer hetare(nbarto)
  81       integer posifa(0:nbarto), facare(nbfaar)
  82       integer nivtri(nbtrto)
  83       integer nivqua(nbquto)
  84       integer arsupp(nbarto), arindi(nbarto)
  85 c
  86       integer ulsort, langue, codret
  87 c
  88 c 0.4. ==> variables locales
  89 c
  90       integer etat
  91       integer larete
  92       integer iaux, jaux, kaux, ideb, ifin
  93 c
  94       integer nbmess
  95       parameter (nbmess = 10 )
  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 c
 114       codret = 0
 115 c
 116 c====
 117 c 2. traitement des indicateurs portant sur les aretes
 118 c    pour le filtrage sur les niveaux, on tient compte du fait que
 119 c    le niveau d'une arete est identifie a celui de l'une quelconque de
 120 c    ses faces voisines quand elle en a. Sinon, on ne filtre pas.
 121 c====
 122 c
 123 #ifdef _DEBUG_HOMARD_
 124         write(ulsort,texte(langue,4)) mess14(langue,3,1)
 125 #endif
 126 c
 127 #ifdef _DEBUG_HOMARD_
 128         write(ulsort,texte(langue,5))
 129 #endif
 130 c
 131       iaux = 0
 132 c
 133       do 21 , larete = 1, nbarto
 134 c
 135         if ( arsupp(larete).ne.0 ) then
 136           etat = mod(hetare(larete),10)
 137           if ( etat.eq.0 ) then
 138             if ( arindi(larete).eq.1 ) then
 139               ideb = posifa(larete-1)+1
 140               ifin = posifa(larete)
 141               jaux = 0
 142               if ( ifin.ge.ideb .and. nivmax.ge.0 ) then
 143                 if ( facare(ideb).gt.0 ) then
 144                   kaux = nivtri(facare(ideb))
 145                 else
 146                   kaux = nivqua(-facare(ideb))
 147                 endif
 148                 if ( kaux.ge.nivmax ) then
 149                   jaux = 1
 150                 endif
 151               endif
 152               if ( jaux.eq.0 ) then
 153                 decare(larete) = 2
 154 cgn          write(ulsort,*) 'mise a 2 de decare pour arete', larete
 155               else
 156                 iaux = iaux + 1
 157               endif
 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,1)
 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