src/tool/Decision/deinst.F

   1       subroutine deinst ( typenh,
   2      >                    seuihe, seuibe,
   3      >                    nbenti, suppor, indire, indien,
   4      >                    ulsort, langue, codret)
   5 c ______________________________________________________________________
   6 c
   7 c                             H O M A R D
   8 c
   9 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
  10 c
  11 c Version originale enregistree le 18 juin 1996 sous le numero 96036
  12 c aupres des huissiers de justice Simart et Lavoir a Clamart
  13 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
  14 c aupres des huissiers de justice
  15 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
  16 c
  17 c    HOMARD est une marque deposee d'Electricite de France
  18 c
  19 c Copyright EDF 1996
  20 c Copyright EDF 1998
  21 c Copyright EDF 2002
  22 c Copyright EDF 2020
  23 c ______________________________________________________________________
  24 c
  25 c traitement des DEcisions - INitialisation des Seuils - Tri
  26 c                --          --                 -        -
  27 c ______________________________________________________________________
  28 c .        .     .        .                                            .
  29 c .  nom   . e/s . taille .           description                      .
  30 c .____________________________________________________________________.
  31 c . typenh . e   .   1    . type d'entites concernees                  .
  32 c .        .     .        . 0 : noeuds                                 .
  33 c .        .     .        . 1 : aretes                                 .
  34 c .        .     .        . 2 : triangles                              .
  35 c .        .     .        . 3 : tetraedres                             .
  36 c .        .     .        . 4 : quadrangles                            .
  37 c .        .     .        . 5 : pyramides                              .
  38 c .        .     .        . 6 : hexaedres                              .
  39 c .        .     .        . 7 : pentaedres                             .
  40 c . seuihe . e   .   1    . borne superieure absolue de l'erreur entite.
  41 c . seuibe . e   .   1    . borne inferieure absolue de l'erreur entite.
  42 c . nbenti . e   .   1    . nombre d'entites pour les entites          .
  43 c . suppor . e   . nbenti . support pour les entites                   .
  44 c . indire . e   . nbenti . valeurs reelles pour les entites           .
  45 c . indien .  s  . nbenti . valeurs entieres filtrees pour les entites .
  46 c . ulsort . e   .   1    . numero d'unite logique de la liste standard.
  47 c . langue . e   .    1   . langue des messages                        .
  48 c .        .     .        . 1 : francais, 2 : anglais                  .
  49 c . codret . es  .    1   . code de retour des modules                 .
  50 c .        .     .        . 0 : pas de probleme                        .
  51 c .        .     .        . 2 : probleme dans le traitement            .
  52 c ______________________________________________________________________
  53 c
  54 c====
  55 c 0. declarations et dimensionnement
  56 c====
  57 c
  58 c 0.1. ==> generalites
  59 c
  60       implicit none
  61       save
  62 c
  63       character*6 nompro
  64       parameter ( nompro = 'DEINST' )
  65 c
  66 #include "nblang.h"
  67 c
  68 c 0.2. ==> communs
  69 c
  70 #include "envex1.h"
  71 #include "impr02.h"
  72 c
  73 c 0.3. ==> arguments
  74 c
  75       integer typenh
  76       integer nbenti
  77       integer suppor(nbenti), indien(nbenti)
  78 c
  79       double precision seuibe, seuihe
  80       double precision indire(nbenti)
  81 c
  82       integer ulsort, langue, codret
  83 c
  84 c 0.4. ==> variables locales
  85 c
  86       integer iaux
  87 c
  88       integer nbmess
  89       parameter (nbmess = 10 )
  90       character*80 texte(nblang,nbmess)
  91 c ______________________________________________________________________
  92 c
  93 c====
  94 c 1. initialisation
  95 c====
  96 c
  97 #include "impr01.h"
  98 c
  99 #ifdef _DEBUG_HOMARD_
 100       write (ulsort,texte(langue,1)) 'Entree', nompro
 101       call dmflsh (iaux)
 102 #endif
 103 c
 104       texte(1,4) = '(''Transfert de reel a entier pour les '',a))'
 105       texte(1,5) = '(''Seuil haut = '',g13.5)'
 106       texte(1,6) = '(''Seuil bas  = '',g13.5)'
 107 c
 108       texte(2,4) = '(''Transfert de reel a entier pour les '',a))'
 109       texte(2,5) = '(''High threshold = '',g13.5)'
 110       texte(2,6) = '(''Low threshold  = '',g13.5)'
 111 c
 112 #ifdef _DEBUG_HOMARD_
 113       write (ulsort,texte(langue,4)) mess14(langue,3,typenh)
 114       write (ulsort,texte(langue,5)) seuihe
 115       write (ulsort,texte(langue,6)) seuibe
 116 #endif
 117 c
 118 #include "impr03.h"
 119 c
 120 c====
 121 c 2. transfert de reel a entier
 122 c====
 123 c
 124       do 21 , iaux = 1, nbenti
 125 cgn      write (ulsort,90012) 'support pour', iaux, suppor(iaux)
 126 c
 127         if ( suppor(iaux).ne.0 ) then
 128 cgn      write (ulsort,90024) 'indire pour', iaux, indire(iaux)
 129 c
 130           if ( indire(iaux).le.seuibe ) then
 131             indien(iaux) = -1
 132           elseif ( indire(iaux).ge.seuihe ) then
 133             indien(iaux) = 1
 134           else
 135             indien(iaux) = 0
 136           endif
 137 cgn      write (ulsort,90012) '==> indien pour', iaux, indien(iaux)
 138 c
 139         endif
 140 c
 141    21  continue
 142 c
 143 c====
 144 c 3. la fin
 145 c====
 146 c
 147       if ( codret.ne.0 ) then
 148 c
 149 #include "envex2.h"
 150 c
 151       write (ulsort,texte(langue,1)) 'Sortie', nompro
 152       write (ulsort,texte(langue,2)) codret
 153 c
 154       endif
 155 c
 156 #ifdef _DEBUG_HOMARD_
 157       write (ulsort,texte(langue,1)) 'Sortie', nompro
 158       call dmflsh (iaux)
 159 #endif
 160 c
 161       end