src/tool/Suivi_Frontiere/sfmop2.F

   1       subroutine sfmop2 ( coonoe, hetnoe, arenoe,
   2      >                    somare,
   3      >                    ulsort, langue, codret)
   4 c ______________________________________________________________________
   5 c                             H O M A R D
   6 c
   7 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
   8 c
   9 c Version originale enregistree le 18 juin 1996 sous le numero 96036
  10 c aupres des huissiers de justice Simart et Lavoir a Clamart
  11 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
  12 c aupres des huissiers de justice
  13 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
  14 c
  15 c    HOMARD est une marque deposee d'Electricite de France
  16 c
  17 c Copyright EDF 1996
  18 c Copyright EDF 1998
  19 c Copyright EDF 2002
  20 c Copyright EDF 2020
  21 c ______________________________________________________________________
  22 c
  23 c   Suivi de Frontiere - MOdification des noeuds P2
  24 c   -        -           --                      --
  25 c ______________________________________________________________________
  26 c .        .     .        .                                            .
  27 c .  nom   . e/s . taille .           description                      .
  28 c .____________________________________________________________________.
  29 c . coonoe . es  . nbnoto . coordonnees des noeuds                     .
  30 c .        .     . *sdim  .                                            .
  31 c . hetnoe . e   . nbnoto . historique de l'etat des noeuds            .
  32 c . arenoe . e   . nbnoto . 0 pour un sommet, le numero de l'arete pour.
  33 c .        .     .        . un noeud milieu                            .
  34 c . somare . es  .2*nbarto. numeros des extremites d'arete             .
  35 c . ulsort . e   .   1    . numero d'unite logique de la liste standard.
  36 c . langue . e   .    1   . langue des messages                        .
  37 c .        .     .        . 1 : francais, 2 : anglais                  .
  38 c . codret . es  .    1   . code de retour des modules                 .
  39 c .        .     .        . 0 : pas de probleme                        .
  40 c .        .     .        . x : probleme                               .
  41 c ______________________________________________________________________
  42 c
  43 c====
  44 c 0. declarations et dimensionnement
  45 c====
  46 c
  47 c 0.1. ==> generalites
  48 c
  49       implicit none
  50       save
  51 c
  52       character*6 nompro
  53       parameter ( nompro = 'SFMOP2' )
  54 c
  55 #include "nblang.h"
  56 c
  57 c 0.2. ==> communs
  58 c
  59 #include "envex1.h"
  60 #include "envca1.h"
  61 c
  62 #include "nombno.h"
  63 #include "nombar.h"
  64 c
  65 c 0.3. ==> arguments
  66 c
  67       double precision coonoe(nbnoto,sdim)
  68 c
  69       integer hetnoe(nbnoto), arenoe(nbnoto)
  70       integer somare(2,nbarto)
  71 c
  72       integer ulsort, langue, codret
  73 c
  74 c 0.4. ==> variables locales
  75 c
  76       integer iaux
  77 #ifdef _DEBUG_HOMARD_
  78       integer jaux
  79 #endif
  80 c
  81       integer lenoeu, lenoe1, lenoe2
  82       integer larete
  83 c
  84       integer nbmess
  85       parameter ( nbmess = 10 )
  86       character*80 texte(nblang,nbmess)
  87 c
  88 c 0.5. ==> initialisations
  89 c ______________________________________________________________________
  90 c
  91 c====
  92 c 1. messages
  93 c====
  94 c
  95       codret = 0
  96 c
  97 c 1.3. ==> les messages
  98 c
  99 #include "impr01.h"
 100 c
 101 #ifdef _DEBUG_HOMARD_
 102       write (ulsort,texte(langue,1)) 'Entree', nompro
 103       call dmflsh (iaux)
 104 #endif
 105 c
 106       texte(1,4) = '(/5x,''Retouche des coordonnees des noeuds P2'')'
 107       texte(1,5) = '(7x,''Nombre de noeuds modifies : '',i10)'
 108       texte(1,6) = '(7x,''==> Arete '',i10)'
 109       texte(1,7) = '(7x,''==> Repositionnement du noeud '',i10)'
 110 c
 111       texte(2,4) = '(/5x,''Updating of P2 node coordinates'')'
 112       texte(2,5) = '(7x,''Number of modified nodes : '',i10)'
 113       texte(2,6) = '(7x,''==> Edge # '',i10)'
 114       texte(2,7) = '(7x,''==> Relocalization of node # '',i10,)'
 115 c
 116       codret = 0
 117 #ifdef _DEBUG_HOMARD_
 118       write (ulsort,texte(langue,1)) 'Entree', nompro
 119 #endif
 120 c
 121 c====
 122 c 2. retouche des coordonnees des noeuds P2
 123 c    remarque : on repositionne tous les nouveaux noeuds P2
 124 c              on pourrait faire plus malin en ne bougeant que ceux
 125 c              qui doivent l'etre mais c'est assez complique. Il faut
 126 c              tous ceux sur une arete dont l'extremite a bouge par
 127 c              suite de raffinement, et tous ceux qui etaient des P1
 128 c              et qui redeviennent P2 par deraffinement.
 129 c              Au final, cela suppose de reperer des voisinages et
 130 c             donc assez complique et pas forcement economique en CPU
 131 c    remarque : on en fait pas assez. il faudrait faire les noeuds
 132 c               au milieu des quadrangles.
 133 c====
 134 c
 135 #ifdef _DEBUG_HOMARD_
 136       write (ulsort,*) '2. noeuds P2 ; codret = ', codret
 137 #endif
 138 c
 139 #ifdef _DEBUG_HOMARD_
 140       jaux = 0
 141 #endif
 142 c
 143       do 21 , lenoeu = 1 , nbnoto
 144 c
 145         if ( mod(hetnoe(lenoeu),10).eq.2 ) then
 146 c
 147 c 2.1. ==> Reperage
 148 c
 149           larete = arenoe(lenoeu)
 150           lenoe1 = somare(1,larete)
 151           lenoe2 = somare(2,larete)
 152 c
 153 #ifdef _DEBUG_HOMARD_
 154       write(ulsort,texte(langue,6)) larete
 155       write(ulsort,texte(langue,7)) lenoeu
 156 #endif
 157 c
 158 c 2.2. ==> Calcul des coordonnees milieu
 159 c
 160          do 221 , iaux = 1 , sdim
 161            coonoe(lenoeu,iaux) =
 162      >           (coonoe(lenoe1,iaux)+coonoe(lenoe2,iaux))*0.5d0
 163   221 continue
 164 #ifdef _DEBUG_HOMARD_
 165          jaux = jaux + 1
 166 #endif
 167 c
 168         endif
 169 c
 170    21 continue
 171 c
 172 #ifdef _DEBUG_HOMARD_
 173       write (ulsort,texte(langue,5)) jaux
 174       write (ulsort,*)
 175 #endif
 176 c
 177 c====
 178 c 3. la fin
 179 c====
 180 c
 181       if ( codret.ne.0 ) then
 182 c
 183 #include "envex2.h"
 184 c
 185       write (ulsort,texte(langue,1)) 'Sortie', nompro
 186       write (ulsort,texte(langue,2)) codret
 187 c
 188       endif
 189 c
 190 #ifdef _DEBUG_HOMARD_
 191       write (ulsort,texte(langue,1)) 'Sortie', nompro
 192       call dmflsh (iaux)
 193 #endif
 194 c
 195       end