src/tool/AV_Conversion/vcmmpo.F

   1       subroutine vcmmpo ( noempo, hetmpo, coexmp,
   2      >                    nnosca, nmpsho, nmpsca,
   3      >                    fameel, typele,
   4      >                    povoso, voisom,
   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    aVant adaptation - Conversion de Maillage - Mailles-POints
  27 c     -                 -             -          -       --
  28 c ______________________________________________________________________
  29 c .        .     .        .                                            .
  30 c .  nom   . e/s . taille .           description                      .
  31 c .____________________________________________________________________.
  32 c . noempo .  s  . nbmpto . numeros des noeuds associes aux mailles    .
  33 c . hetmpo .  s  . rbar00 . historique de l'etat des mailles-points    .
  34 c . coexmp .  s  . nbmpto*. codes externes sur les mailles-points      .
  35 c .        .     . nctfmp .   1 : famille MED                          .
  36 c .        .     .        .   2 : type de maille-point                 .
  37 c . nnosca . e   . rsnoto . numero des noeuds dans le calcul           .
  38 c . nmpsho .  s  . rsmpac . numero des mailles-points dans HOMARD      .
  39 c . nmpsca .  s  . rsmpto . numero des mailles-points du calcul        .
  40 c . fameel . e   . nbelem . famille med des elements                   .
  41 c . typele . e   . nbelem . type des elements pour le code de calcul   .
  42 c . povoso . e   .0:nbnoto. pointeur des voisins par sommet            .
  43 c . voisom . e   . nvosom . voisins des sommets en stockage morse      .
  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 .        .     .        . 3 : probleme                               .
  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 = 'VCMMPO' )
  63 c
  64 #include "nblang.h"
  65 c
  66 c 0.2. ==> communs
  67 c
  68 #include "envex1.h"
  69 c
  70 #include "dicfen.h"
  71 #include "nbutil.h"
  72 #include "nombno.h"
  73 #include "refert.h"
  74 #include "nombmp.h"
  75 #include "nombsr.h"
  76 #include "rftmed.h"
  77 c
  78 c 0.3. ==> arguments
  79 c
  80       integer noempo(nbmpto), hetmpo(nbmpto)
  81       integer nnosca(rsnoto)
  82       integer nmpsho(rsmpac), nmpsca(rsmpto)
  83       integer coexmp(nbmpto,nctfmp)
  84       integer fameel(nbelem), typele(nbelem)
  85       integer voisom(nvosom), povoso(0:nbnoto)
  86 c
  87       integer ulsort, langue, codret
  88 c
  89 c 0.4. ==> variables locales
  90 c
  91       integer iaux
  92       integer poinde, poinfi, point, nucode
  93       integer sommet
  94       integer lamapo, elem, typhom
  95 c
  96       integer nbmess
  97       parameter ( nbmess = 10 )
  98       character*80 texte(nblang,nbmess)
  99 c
 100 c 0.5. ==> initialisations
 101 c ______________________________________________________________________
 102 c
 103 c====
 104 c 1. preliminaires
 105 c====
 106 c
 107 c 1.1. ==> messages
 108 c
 109 #include "impr01.h"
 110 c
 111 #ifdef _DEBUG_HOMARD_
 112       write (ulsort,texte(langue,1)) 'Entree', nompro
 113       call dmflsh (iaux)
 114 #endif
 115 c
 116       texte(1,4) = '(''Nombre de mailles-points attendues : '',i10)'
 117       texte(1,5) = '(''Nombre de mailles-points trouvees  : '',i10)'
 118 c
 119       texte(2,4) = '(''Expected point-meshes : '',i10)'
 120       texte(2,5) = '(''Found point-meshes    : '',i10)'
 121 c
 122 c 1.2. ==> mise a zero
 123 c
 124       codret = 0
 125 c
 126       do 11 , lamapo = 1 , rsmpac
 127         nmpsho(lamapo) = 0
 128    11 continue
 129 c
 130       do 12 , nucode = 1 , nctfmp
 131         do 121 , lamapo = 1 , nbmpto
 132           coexmp(lamapo,nucode) = 0
 133   121   continue
 134    12 continue
 135 c
 136       lamapo = 0
 137 c
 138 c====
 139 c 2. on passe en revue chaque sommet
 140 c     remarque : l'exploration se fait dans la numerotation HOMARD
 141 c     ses elements voisins sont dans le tableau voisom, aux places
 142 c     comprises entre povoso(somm-1)+1 et povoso(somm), somm etant le
 143 c     numero dans le calcul correspondant au numero sommet dans
 144 c     homard
 145 c====
 146 c
 147       do 21 , sommet = 1 , nbnoto
 148 c
 149         poinde = povoso(nnosca(sommet)-1) + 1
 150         poinfi = povoso(nnosca(sommet))
 151 c
 152         do 22 , point = poinde , poinfi
 153 c
 154 c 2.1. ==> caracterisation de l'element
 155 c         elem : son numero global
 156 c         typhom : son type dans HOMARD
 157 c
 158           elem = voisom(point)
 159           typhom = medtrf(typele(elem))
 160 c
 161 c 2.2. ==> si l'element est une maille-point, il faut
 162 c              se souvenir de son type et de sa famille MED.
 163 c
 164           if ( typhom.eq.tyhmpo ) then
 165 c
 166              lamapo = lamapo + 1
 167              noempo(lamapo) = sommet
 168              coexmp(lamapo,1) = fameel(elem)
 169              coexmp(lamapo,2) = typele(elem)
 170              nmpsho(elem) = lamapo
 171              nmpsca(lamapo) = elem
 172            endif
 173 c
 174    22   continue
 175 c
 176    21 continue
 177 c
 178       if ( lamapo.ne.nbmpto ) then
 179         write(ulsort,texte(langue,4)) nbmpto
 180         write(ulsort,texte(langue,5)) lamapo
 181         codret = 3
 182       endif
 183 c
 184 c====
 185 c 3. consequences
 186 c====
 187 c
 188       do 32 , lamapo = 1 , nbmpto
 189         hetmpo(lamapo) = 0
 190    32 continue
 191 c
 192 c====
 193 c 4. la fin
 194 c====
 195 c
 196       if ( codret.ne.0 ) then
 197 c
 198 #include "envex2.h"
 199 c
 200       write (ulsort,texte(langue,1)) 'Sortie', nompro
 201       write (ulsort,texte(langue,2)) codret
 202 c
 203       endif
 204 c
 205 #ifdef _DEBUG_HOMARD_
 206       write (ulsort,texte(langue,1)) 'Sortie', nompro
 207       call dmflsh (iaux)
 208 #endif
 209 c
 210       end