src/tool/Modification/mmag91.F

   1       subroutine mmag91 ( larete, ordre, nujois,
   2      >                    nbduno, tbau30,
   3      >                    somare,
   4      >                    aredup,
   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    Modification de Maillage - AGRegat - phase 9.1
  27 c    -               -          --              - -
  28 c    Reperage des aretes liees aux joints simples formant
  29 c    le joint multiple
  30 c ______________________________________________________________________
  31 c .        .     .        .                                            .
  32 c .  nom   . e/s . taille .           description                      .
  33 c .____________________________________________________________________.
  34 c . larete . e   .   1    . arete support du joint multiple            .
  35 c . ordre  . e   .   1    . ordre du joint a explorer                  .
  36 c . nujois . e   . ordre  . numeros des joints simples associes        .
  37 c . nbduno . e   .   1    . nombre de duplication de noeuds            .
  38 c . tbau30 . e   .8*nbduno. Pour la i-eme duplication de noeud :       .
  39 c .        .     .        . (1,i) : noeud a dupliquer                  .
  40 c .        .     .        . (2,i) : arete construite sur le noeud      .
  41 c .        .     .        . (3,i) : noeud cree cote min(fammed)        .
  42 c .        .     .        . (4,i) : noeud cree cote max(fammed)        .
  43 c .        .     .        . (5,i) : numero du joint simple cree        .
  44 c .        .     .        . (6,i) : arete entrant dans le cote 1       .
  45 c .        .     .        . (7,i) : arete entrant dans le cote 2       .
  46 c .        .     .        . (8,i) : ordre de multiplicite              .
  47 c . somare . e   .2*nbarto. numeros des extremites d'arete             .
  48 c . aredup .  s  .2*nbduno. aretes issues de la duplication            .
  49 c . ulsort . e   .   1    . numero d'unite logique de la liste standard.
  50 c . langue . e   .    1   . langue des messages                        .
  51 c .        .     .        . 1 : francais, 2 : anglais                  .
  52 c . codret . es  .    1   . code de retour des modules                 .
  53 c .        .     .        . 0 : pas de probleme                        .
  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 = 'MMAG91' )
  67 c
  68 #include "nblang.h"
  69 c
  70 c 0.2. ==> communs
  71 c
  72 #include "envex1.h"
  73 c
  74 #include "nombar.h"
  75 #include "impr02.h"
  76 c
  77 c 0.3. ==> arguments
  78 c
  79       integer larete
  80       integer ordre, nujois(ordre)
  81       integer nbduno, tbau30(8,nbduno)
  82       integer somare(2,nbarto)
  83       integer aredup(2*ordre)
  84 c
  85       integer ulsort, langue, codret
  86 c
  87 c 0.4. ==> variables locales
  88 c
  89       integer iaux, jaux, kaux
  90       integer nbarlo
  91 c
  92       integer nbmess
  93       parameter ( nbmess = 40 )
  94       character*80 texte(nblang,nbmess)
  95 c
  96 c 0.5. ==> initialisations
  97 c ______________________________________________________________________
  98 c
  99 c====
 100 c 1. initialisations
 101 c====
 102 c 1.1. ==> messages
 103 c
 104 #include "impr01.h"
 105 c
 106 #ifdef _DEBUG_HOMARD_
 107       write (ulsort,texte(langue,1)) 'Entree', nompro
 108       call dmflsh (iaux)
 109 #endif
 110 c
 111 #include "impr03.h"
 112 #include "mmag01.h"
 113 #include "mmag02.h"
 114 c
 115       codret = 0
 116 c
 117 c====
 118 c 2. On explore tous les noeuds qui ont ete dupliques
 119 c    On est bon quand :
 120 c    . le noeud duplique est un des sommets de l'arete triple
 121 c    . la duplication a lieu pour un des joints simples associes
 122 c
 123 c    On stocke les aretes qui partent de chacun des noeuds.
 124 c====
 125 c
 126       nbarlo = 0
 127 c
 128       do 21 , iaux = 1 , nbduno
 129 c
 130         do 211 , jaux = 1 , 2
 131 c
 132 #ifdef _DEBUG_HOMARD_
 133       write (ulsort,texte(langue,4)) '. ', mess14(langue,1,-1),
 134      >     somare(jaux,larete)
 135 #endif
 136 c
 137           if ( tbau30(1,iaux).eq.somare(jaux,larete) ) then
 138 c
 139             do 2111 , kaux = 1 , ordre
 140 c
 141 c      write (ulsort,*) 'nujois(kaux) =', nujois(kaux)
 142 c
 143               if ( tbau30(5,iaux).eq.nujois(kaux) ) then
 144                 aredup(ordre*(jaux-1)+kaux) = tbau30(2,iaux)
 145 #ifdef _DEBUG_HOMARD_
 146       write (ulsort,texte(langue,40)) mess14(langue,1,-1),
 147      >          tbau30(1,iaux),mess14(langue,1,1),tbau30(2,iaux)
 148 #endif
 149                 nbarlo = nbarlo + 1
 150                 if ( nbarlo.eq.2*ordre ) then
 151                   goto 210
 152                 endif
 153                 goto 2110
 154               endif
 155  2111       continue
 156 c
 157  2110       continue
 158 c
 159           endif
 160 c
 161   211   continue
 162 c
 163    21 continue
 164 c
 165 c     Si on arrive ici, c'est que les aretes n'ont pas ete trouvees
 166 c
 167       write (ulsort,texte(langue,4)) ' ', mess14(langue,1,1), larete
 168       write (ulsort,texte(langue,38)) mess14(langue,3,-1)
 169       codret = 21
 170 c
 171   210 continue
 172 c
 173 c====
 174 c 3. la fin
 175 c====
 176 c
 177       if ( codret.ne.0 ) then
 178 c
 179 #include "envex2.h"
 180 c
 181       write (ulsort,texte(langue,1)) 'Sortie', nompro
 182       write (ulsort,texte(langue,2)) codret
 183 c
 184       endif
 185 c
 186 #ifdef _DEBUG_HOMARD_
 187       write (ulsort,texte(langue,1)) 'Sortie', nompro
 188       call dmflsh (iaux)
 189 #endif
 190 c
 191       end