src/tool/Decision/deiqud.F

   1       subroutine deiqud ( nivmin,
   2      >                    decare, decfac,
   3      >                    arequa, hetqua, filqua, nivqua,
   4      >                    qusupp, quindi,
   5      >                    ulsort, langue, codret)
   6 c ______________________________________________________________________
   7 c
   8 c                             H O M A R D
   9 c
  10 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
  11 c
  12 c Version originale enregistree le 18 juin 1996 sous le numero 96036
  13 c aupres des huissiers de justice Simart et Lavoir a Clamart
  14 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
  15 c aupres des huissiers de justice
  16 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
  17 c
  18 c    HOMARD est une marque deposee d'Electricite de France
  19 c
  20 c Copyright EDF 1996
  21 c Copyright EDF 1998
  22 c Copyright EDF 2002
  23 c Copyright EDF 2020
  24 c ______________________________________________________________________
  25 c
  26 c traitement des DEcisions - Initialisation de l'indicateur entier
  27 c                --          -
  28 c                          - cas des QUadrangles - Deraffinement
  29 c                                    --            -
  30 c ______________________________________________________________________
  31 c .        .     .        .                                            .
  32 c .  nom   . e/s . taille .           description                      .
  33 c .____________________________________________________________________.
  34 c . nivmin . e   .   1    . niveau min a ne pas depasser en deraffinemt.
  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 . arequa . e   .nbquto*4. numeros des 4 aretes des quadrangles       .
  39 c . hetqua . e   . nbquto . historique de l'etat des quadrangles       .
  40 c . filqua . e   . nbquto . premier fils des quadrangles               .
  41 c . nivqua . e   . nbquto . niveau des quadrangles                     .
  42 c . qusupp . e   . nbquto . support pour les quadrangles               .
  43 c . quindi . e   . nbquto . valeurs entieres pour les quadrangles      .
  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 = 'DEIQUD' )
  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 "impr02.h"
  73 c
  74 c 0.3. ==> arguments
  75 c
  76       integer nivmin
  77       integer decare(0:nbarto), decfac(-nbquto:nbtrto)
  78       integer arequa(nbquto,4), hetqua(nbquto), filqua(nbquto)
  79       integer nivqua(nbquto)
  80       integer qusupp(nbquto), quindi(nbquto)
  81 c
  82       integer ulsort, langue, codret
  83 c
  84 c 0.4. ==> variables locales
  85 c
  86       integer areloc, etat
  87       integer lequad
  88       integer fille1
  89       integer iaux
  90 c
  91       integer nbmess
  92       parameter (nbmess = 30 )
  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 #include "derco1.h"
 111 c
 112       codret = 0
 113 c
 114 c====
 115 c 2. traitement des indicateurs portant sur les quadrangles
 116 c====
 117 c
 118 #ifdef _DEBUG_HOMARD_
 119         write(ulsort,texte(langue,4)) mess14(langue,3,4)
 120 cgn      write (ulsort,*) 'quadrangle 7'
 121 cgn      write (ulsort,*) 'decfac(q7) =',decfac(-7)
 122 cgn      write (ulsort,*) arequa(7,1),arequa(7,2),
 123 cgn     >arequa(7,3),arequa(7,4)
 124 cgn      write (ulsort,*) decare(arequa(7,1)),decare(arequa(7,2)),
 125 cgn     >decare(arequa(7,3)),decare(arequa(7,4))
 126 cgn      write (ulsort,*) hetare(arequa(7,1)),hetare(arequa(7,2)),
 127 cgn     >hetare(arequa(7,3)),hetare(arequa(7,4))
 128 cgn      write (ulsort,*) ' '
 129 #endif
 130 cgn      write(ulsort,*)'quindi :'
 131 cgn      write(ulsort,1789)(lequad, quindi(lequad),lequad = 1, nbquto)
 132 cgn 1789 format(5(i3,' : ',i4,', '))
 133 c
 134 #ifdef _DEBUG_HOMARD_
 135         write(ulsort,texte(langue,6))
 136 #endif
 137 c
 138       iaux = 0
 139 c
 140       do 21 , lequad = 1, nbquto
 141         etat = mod(hetqua(lequad),100)
 142 #ifdef _DEBUG_HOMARD_
 143       write (ulsort,*) 'quad ',lequad,', etat = ',etat
 144 #endif
 145         if ( etat.eq.4  ) then
 146           fille1 = filqua(lequad)
 147           if ( qusupp(fille1)   .ne. 0 .and.
 148      >         qusupp(fille1+1) .ne. 0 .and.
 149      >         qusupp(fille1+2) .ne. 0 .and.
 150      >         qusupp(fille1+3) .ne. 0 ) then
 151             if ( quindi(fille1)   .eq. -1  .and.
 152      >           quindi(fille1+1) .eq. -1 .and.
 153      >           quindi(fille1+2) .eq. -1 .and.
 154      >           quindi(fille1+3) .eq. -1 ) then
 155               if ( nivqua(lequad).lt.nivmin ) then
 156                 iaux = iaux + 4
 157               else
 158                 decfac(-lequad) = -1
 159 #ifdef _DEBUG_HOMARD_
 160       write (ulsort,texte(langue,30))'decfac',
 161      >-lequad,decfac(-lequad),' '
 162 #endif
 163                 do 22 , areloc = 1, 4
 164                   decare(arequa(lequad,areloc)) = -1
 165 #ifdef _DEBUG_HOMARD_
 166       write (ulsort,texte(langue,30))'decare',arequa(lequad,areloc),
 167      >  decare(arequa(lequad,areloc)),' (une de ses aretes)'
 168 #endif
 169    22           continue
 170               endif
 171             endif
 172           endif
 173         endif
 174    21 continue
 175 cgn      write (ulsort,*) 'apres 312'
 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,8)) nivmin
 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