src/tool/Decision/deiqur.F

   1       subroutine deiqur ( nivmax,
   2      >                    decare, decfac,
   3      >                    hetare,
   4      >                    arequa, hetqua, nivqua,
   5      >                    qusupp, quindi,
   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 QUadrangles - Raffinement
  30 c                                    --            -
  31 c ______________________________________________________________________
  32 c .        .     .        .                                            .
  33 c .  nom   . e/s . taille .           description                      .
  34 c .____________________________________________________________________.
  35 c . nivmax . e   .   1    . niveau max a ne pas depasser en raffinement.
  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 . hetare . e   . nbarto . historique de l'etat des aretes            .
  40 c . arequa . e   .nbquto*4. numeros des 4 aretes des quadrangles       .
  41 c . hetqua . e   . nbquto . historique de l'etat des quadrangles       .
  42 c . nivqua . e   . nbquto . niveau des quadrangles                     .
  43 c . qusupp . e   . nbquto . support pour les quadrangles               .
  44 c . quindi . e   . nbquto . valeurs entieres pour les quadrangles      .
  45 c . ulsort . e   .   1    . numero d'unite logique de la liste standard.
  46 c . langue . e   .    1   . langue des messages                        .
  47 c .        .     .        . 1 : francais, 2 : anglais                  .
  48 c . codret . es  .    1   . code de retour des modules                 .
  49 c .        .     .        . 0 : pas de probleme                        .
  50 c .        .     .        . 2 : probleme dans le traitement            .
  51 c ______________________________________________________________________
  52 c
  53 c====
  54 c 0. declarations et dimensionnement
  55 c====
  56 c
  57 c 0.1. ==> generalites
  58 c
  59       implicit none
  60       save
  61 c
  62       character*6 nompro
  63       parameter ( nompro = 'DEIQUR' )
  64 c
  65 #include "nblang.h"
  66 c
  67 c 0.2. ==> communs
  68 c
  69 #include "envex1.h"
  70 #include "nombar.h"
  71 #include "nombtr.h"
  72 #include "nombqu.h"
  73 #include "impr02.h"
  74 c
  75 c 0.3. ==> arguments
  76 c
  77       integer nivmax
  78       integer decare(0:nbarto), decfac(-nbquto:nbtrto)
  79       integer hetare(nbarto)
  80       integer arequa(nbquto,4), hetqua(nbquto)
  81       integer nivqua(nbquto)
  82       integer qusupp(nbquto), quindi(nbquto)
  83 c
  84       integer ulsort, langue, codret
  85 c
  86 c 0.4. ==> variables locales
  87 c
  88       integer areloc
  89       integer lequad
  90       integer iaux, kaux
  91 c
  92       integer nbmess
  93       parameter (nbmess = 30 )
  94       character*80 texte(nblang,nbmess)
  95 c ______________________________________________________________________
  96 c
  97 c====
  98 c 1. initialisation
  99 c====
 100 c
 101 #include "impr01.h"
 102 c
 103 #ifdef _DEBUG_HOMARD_
 104       write (ulsort,texte(langue,1)) 'Entree', nompro
 105       call dmflsh (iaux)
 106 #endif
 107 c
 108 #include "impr03.h"
 109 c
 110 #include "impr05.h"
 111 #include "derco1.h"
 112 c
 113       codret = 0
 114 c
 115 c====
 116 c 2. traitement des indicateurs portant sur les quadrangles
 117 c====
 118 c
 119 #ifdef _DEBUG_HOMARD_
 120         write(ulsort,texte(langue,4)) mess14(langue,3,4)
 121 cgn      write (ulsort,*) 'quadrangle 7'
 122 cgn      write (ulsort,*) 'decfac(q7) =',decfac(-7)
 123 cgn      write (ulsort,*) arequa(7,1),arequa(7,2),
 124 cgn     >arequa(7,3),arequa(7,4)
 125 cgn      write (ulsort,*) decare(arequa(7,1)),decare(arequa(7,2)),
 126 cgn     >decare(arequa(7,3)),decare(arequa(7,4))
 127 cgn      write (ulsort,*) hetare(arequa(7,1)),hetare(arequa(7,2)),
 128 cgn     >hetare(arequa(7,3)),hetare(arequa(7,4))
 129 cgn      write (ulsort,*) ' '
 130 #endif
 131 c
 132 #ifdef _DEBUG_HOMARD_
 133       write(ulsort,texte(langue,5))
 134 #endif
 135 c
 136       iaux = 0
 137 c
 138       do 21 , lequad = 1, nbquto
 139 c
 140 cgn      write (ulsort,90112) 'qusupp/quindi',
 141 cgn     >                     lequad,qusupp(lequad),quindi(lequad)
 142         if ( qusupp(lequad).ne.0 ) then
 143           if ( quindi(lequad).eq.1 ) then
 144             if ( mod(hetqua(lequad),100).eq.0 ) then
 145               if ( nivmax.gt.0 .and. nivqua(lequad).ge.nivmax ) then
 146                 iaux = iaux + 1
 147               else
 148                 decfac(-lequad) = 4
 149 #ifdef _DEBUG_HOMARD_
 150       write (ulsort,texte(langue,30))'decfac',
 151      >-lequad,decfac(-lequad),' '
 152 #endif
 153                 do 22 , areloc = 1, 4
 154                   kaux = arequa(lequad,areloc)
 155                   if ( mod(hetare(kaux),10).eq.0 ) then
 156                     decare(kaux) = 2
 157 #ifdef _DEBUG_HOMARD_
 158       write (ulsort,texte(langue,30))'decare',kaux,
 159      >  decare(kaux),' (une de ses aretes)'
 160 #endif
 161                   elseif ( mod(hetare(kaux),10).eq.2 ) then
 162                     decare(kaux) = 0
 163 #ifdef _DEBUG_HOMARD_
 164       write (ulsort,texte(langue,30))'decare',kaux,
 165      >  decare(kaux),' (une de ses aretes)'
 166 #endif
 167                   endif
 168    22           continue
 169               endif
 170             endif
 171           endif
 172         endif
 173 c
 174    21 continue
 175 cgn      write (ulsort,*) 'apres 322'
 176 cgn      write (ulsort,*) 'quadrangle 7'
 177 cgn      write (ulsort,*) 'decfac(q7) =',decfac(-7)
 178 cgn      write (ulsort,*) arequa(7,1),arequa(7,2),
 179 cgn     >arequa(7,3),arequa(7,4)
 180 cgn      write (ulsort,*) decare(arequa(7,1)),decare(arequa(7,2)),
 181 cgn     >decare(arequa(7,3)),decare(arequa(7,4))
 182 cgn      write (ulsort,*) hetare(arequa(7,1)),hetare(arequa(7,2)),
 183 cgn     >hetare(arequa(7,3)),hetare(arequa(7,4))
 184 cgn      write (ulsort,*) ' '
 185 cgn        print 1789,(lequad, decfac(lequad),lequad = -nbquto,0)
 186 c
 187       if ( iaux.ne.0 ) then
 188         write(ulsort,texte(langue,10))
 189         write(ulsort,texte(langue,4)) mess14(langue,3,2)
 190         write(ulsort,texte(langue,7)) nivmax
 191         write(ulsort,texte(langue,9)) iaux
 192       endif
 193 c
 194 c====
 195 c 3. la fin
 196 c====
 197 c
 198       if ( codret.ne.0 ) then
 199 c
 200 #include "envex2.h"
 201 c
 202       write (ulsort,texte(langue,1)) 'Sortie', nompro
 203       write (ulsort,texte(langue,2)) codret
 204 c
 205       endif
 206 c
 207 #ifdef _DEBUG_HOMARD_
 208       write (ulsort,texte(langue,1)) 'Sortie', nompro
 209       call dmflsh (iaux)
 210 #endif
 211 c
 212       end