src/tool/Suivi_Frontiere/sfnnfl.F

   1       subroutine sfnnfl ( hetare,
   2      >                    cfaare, famare,
   3      >                    ulsort, langue, codret)
   4 c ______________________________________________________________________
   5 c
   6 c                             H O M A R D
   7 c
   8 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
   9 c
  10 c Version originale enregistree le 18 juin 1996 sous le numero 96036
  11 c aupres des huissiers de justice Simart et Lavoir a Clamart
  12 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
  13 c aupres des huissiers de justice
  14 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
  15 c
  16 c    HOMARD est une marque deposee d'Electricite de France
  17 c
  18 c Copyright EDF 1996
  19 c Copyright EDF 1998
  20 c Copyright EDF 2002
  21 c Copyright EDF 2020
  22 c ______________________________________________________________________
  23 c
  24 c but : compter le nombre de noeuds du maillage de calcul qui sont sur
  25 c       la frontiere
  26 c ______________________________________________________________________
  27 c
  28 c   Suivi de Frontiere - Nombre de Noeuds sur la Frontiere - Lignes
  29 c   -        -           -  -                    -           -
  30 c ______________________________________________________________________
  31 c .        .     .        .                                            .
  32 c .  nom   . e/s . taille .           description                      .
  33 c .____________________________________________________________________.
  34 c . hetare . e   . nbarto . historique de l'etat des aretes            .
  35 c . cfaare . e   . nctfar*. codes des familles des aretes              .
  36 c .        .     . nbfare .   1 : famille MED                          .
  37 c .        .     .        .   2 : type de segment                      .
  38 c .        .     .        .   3 : orientation                          .
  39 c .        .     .        .   4 : famille d'orientation inverse        .
  40 c .        .     .        .   5 : numero de ligne de frontiere         .
  41 c .        .     .        .  > 0 si concernee par le suivi de frontiere.
  42 c .        .     .        . <= 0 si non concernee                      .
  43 c .        .     .        .   6 : famille frontiere active/inactive    .
  44 c .        .     .        .   7 : numero de surface de frontiere       .
  45 c .        .     .        . + l : appartenance a l'equivalence l       .
  46 c . famare . e   . nbarto . famille des aretes                         .
  47 c . ulsort . e   .   1    . numero d'unite logique de la liste standard.
  48 c . langue . e   .    1   . langue des messages                        .
  49 c .        .     .        . 1 : francais, 2 : anglais                  .
  50 c . codret . es  .    1   . code de retour des modules                 .
  51 c .        .     .        . 0 : pas de probleme                        .
  52 c .        .     .        . x : probleme                               .
  53 c ______________________________________________________________________
  54 c
  55 c====
  56 c 0. declarations et dimensionnement
  57 c====
  58 c
  59 c 0.1. ==> generalites
  60 c
  61       implicit none
  62       save
  63 c
  64       character*6 nompro
  65       parameter ( nompro = 'SFNNFL' )
  66 c
  67 #include "nblang.h"
  68 #include "cofaar.h"
  69 c
  70 c 0.2. ==> communs
  71 c
  72 #include "envex1.h"
  73 #include "envca1.h"
  74 #include "front2.h"
  75 #include "dicfen.h"
  76 #include "nbfami.h"
  77 #include "nombar.h"
  78 c
  79 c 0.3. ==> arguments
  80 c
  81       integer hetare(nbarto)
  82       integer famare(nbarto), cfaare(nctfar,nbfare)
  83 c
  84       integer ulsort, langue, codret
  85 c
  86 c 0.4. ==> variables locales
  87 c
  88       integer iaux
  89       integer lig
  90       integer nbarfr
  91 c
  92       integer nbmess
  93       parameter ( nbmess = 10 )
  94       character*80 texte(nblang,nbmess)
  95 c
  96 c 0.5. ==> initialisations
  97 c ______________________________________________________________________
  98 c
  99 c====
 100 c 1. messages
 101 c====
 102 c
 103 c 1.3. ==> les messages
 104 c
 105 #include "impr01.h"
 106 c
 107 #ifdef _DEBUG_HOMARD_
 108       write (ulsort,texte(langue,1)) 'Entree', nompro
 109       call dmflsh (iaux)
 110 #endif
 111 c
 112       texte(1,4) =
 113      > '(''. Max du nombre de noeuds sur la frontiere :'',i10)'
 114       texte(1,5) = '(''Arete '',i10,'' de ligne'',i10)'
 115 c
 116       texte(2,4) = '(''. Max of number of nodes over boundary :'',i10)'
 117       texte(2,5) = '(''Edge '',i10,'' in line'',i10)'
 118 c
 119 c====
 120 c 2. boucle sur les aretes : situer les sommets sur la ligne
 121 c====
 122 #ifdef _DEBUG_HOMARD_
 123       write (ulsort,*) '2. boucle sur les aretes ; codret = ', codret
 124 #endif
 125 c
 126       nbarfr = 0
 127 c
 128       if ( codret.eq.0 ) then
 129 c
 130       do 21 , iaux = 1 , nbarto
 131 c
 132         lig = cfaare(cosfli,famare(iaux))
 133 #ifdef _DEBUG_HOMARD_
 134         if ( lig.gt.0 ) then
 135           write (ulsort,texte(langue,5)) iaux, lig
 136         endif
 137 #endif
 138 c
 139         if ( lig.gt.0 ) then
 140 c
 141           if ( hetare(iaux).ne.50 ) then
 142 c
 143             nbarfr = nbarfr + 1
 144 c
 145           endif
 146 c
 147         endif
 148 c
 149    21 continue
 150 c
 151       endif
 152 c
 153       if ( degre.eq.1 ) then
 154         mcnxnf = 2*nbarfr
 155       else
 156         mcnxnf = 3*nbarfr
 157       endif
 158 c
 159 #ifdef _DEBUG_HOMARD_
 160       write (ulsort,texte(langue,4)) mcnxnf
 161 #endif
 162 c
 163 c====
 164 c 3. la fin
 165 c====
 166 c
 167       if ( codret.ne.0 ) then
 168 c
 169 #include "envex2.h"
 170 c
 171       write (ulsort,texte(langue,1)) 'Sortie', nompro
 172       write (ulsort,texte(langue,2)) codret
 173 c
 174       endif
 175 c
 176 #ifdef _DEBUG_HOMARD_
 177       write (ulsort,texte(langue,1)) 'Sortie', nompro
 178       call dmflsh (iaux)
 179 #endif
 180 c
 181       end
 182
 183