src/tool/Utilitaire/utflt0.F

   1       subroutine utflt0 ( somare, aretes,
   2      >                    champ, flux,
   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    UTilitaire - FLux a travers un Triangle - 0
  25 c    --           --                -          -
  26 c Calcule le flux d'une grandeur a travers un triangle par la somme
  27 c de la grandeur le long des aretes
  28 c Variante de utfltr quand le champ est deja une circulation
  29 c ______________________________________________________________________
  30 c .        .     .        .                                            .
  31 c .  nom   . e/s . taille .           description                      .
  32 c .____________________________________________________________________.
  33 c . somare . e   .2*nbaret. numeros des extremites d'arete             .
  34 c . aretes . e   .   3    . numeros des 3 aretes du triangle           .
  35 c . champ  . e   .   3    . champ sur les 3 aretes du triangle         .
  36 c . flux   .  s  . 1      . le flux                                    .
  37 c . ulsort . e   .   1    . numero d'unite logique de la liste standard.
  38 c . langue . e   .    1   . langue des messages                        .
  39 c .        .     .        . 1 : francais, 2 : anglais                  .
  40 c . codret . es  .    1   . code de retour des modules                 .
  41 c .        .     .        . 0 : pas de probleme                        .
  42 c .        .     .        . 1 : probleme                               .
  43 c ______________________________________________________________________
  44 c
  45 c====
  46 c 0. declarations et dimensionnement
  47 c====
  48 c
  49 c 0.1. ==> generalites
  50 c
  51       implicit none
  52       save
  53 c
  54       character*6 nompro
  55       parameter ( nompro = 'UTFLT0' )
  56 c
  57 #include "nblang.h"
  58 c
  59 c 0.2. ==> communs
  60 c
  61 #include "envex1.h"
  62 c
  63 c 0.3. ==> arguments
  64 c
  65       integer somare(2,*)
  66       integer aretes(3)
  67 c
  68       double precision champ(3)
  69       double precision flux
  70 c
  71       integer ulsort, langue, codret
  72 c
  73 c 0.4. ==> variables locales
  74 c
  75       integer iaux
  76       integer orient(3)
  77 c
  78       integer nbmess
  79       parameter ( nbmess = 10 )
  80       character*80 texte(nblang,nbmess)
  81 c
  82 c 0.5. ==> initialisations
  83 c ______________________________________________________________________
  84 c
  85 c====
  86 c 1. initialisations
  87 c====
  88 c
  89 #include "impr01.h"
  90 c
  91 #ifdef _DEBUG_HOMARD_
  92       write (ulsort,texte(langue,1)) 'Entree', nompro
  93       call dmflsh (iaux)
  94 #endif
  95 c
  96 #include "impr03.h"
  97 c
  98       codret = 0
  99 c
 100 c====
 101 c 2. Orientations des aretes
 102 c====
 103 c
 104       call utorat ( somare,
 105      >              aretes(1), aretes(2), aretes(3),
 106      >              orient(1), orient(2), orient(3) )
 107 #ifdef _DEBUG_HOMARD_
 108       write (ulsort,90002) '. Aretes      ', aretes
 109       write (ulsort,90002) '. Orientations', orient
 110 #endif
 111 c
 112 c====
 113 c 3. Parcours des aretes
 114 c====
 115 c
 116       flux = 0.d0
 117 c
 118       do 31 , iaux = 1 , 3
 119 c
 120 cgn      write (ulsort,90004) 'champ',champ(iaux),orient(iaux)*champ(iaux)
 121         flux = flux + dble(orient(iaux))*champ(iaux)
 122 c
 123    31 continue
 124 c
 125 #ifdef _DEBUG_HOMARD_
 126       write (ulsort,90004) '==> Flux', flux
 127 #endif
 128 c
 129 c====
 130 c 4. la fin
 131 c====
 132 c
 133       if ( codret.ne.0 ) then
 134 c
 135 #include "envex2.h"
 136 c
 137       write (ulsort,texte(langue,1)) 'Sortie', nompro
 138       write (ulsort,texte(langue,2)) codret
 139 c
 140       endif
 141 c
 142 #ifdef _DEBUG_HOMARD_
 143       write (ulsort,texte(langue,1)) 'Sortie', nompro
 144       call dmflsh (iaux)
 145 #endif
 146 c
 147       end