src/tool/Creation_Maillage/cmnosu.F

   1       subroutine cmnosu ( nohman,
   2      >                    indnoe,
   3      >                    lgetco, taetco,
   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    Creation du Maillage - les NOeuds SUpplementairs
  26 c    -           -              --     --
  27 c ______________________________________________________________________
  28 c
  29 c but : gestion de la creation des noeuds p2 et internes aux mailles
  30 c ______________________________________________________________________
  31 c .        .     .        .                                            .
  32 c .  nom   . e/s . taille .           description                      .
  33 c .____________________________________________________________________.
  34 c . nohman . e   .  ch8   . nom de l'objet contenant le maillage       .
  35 c . indnoe . es  . 1      . indice du dernier noeud cree               .
  36 c . lgetco . e   .   1    . longueur du tableau de l'etat courant      .
  37 c . taetco . e   . lgetco . tableau de l'etat courant                  .
  38 c . ulsort . e   .   1    . numero d'unite logique de la liste standard.
  39 c . langue . e   .    1   . langue des messages                        .
  40 c .        .     .        . 1 : francais, 2 : anglais                  .
  41 c . codret . es  .    1   . code de retour des modules                 .
  42 c .        .     .        . 0 : pas de 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 = 'CMNOSU' )
  56 c
  57 #include "nblang.h"
  58 c
  59 c 0.2. ==> communs
  60 c
  61 #include "envex1.h"
  62 #include "envca1.h"
  63 c
  64 #include "nombar.h"
  65 #include "nouvnb.h"
  66 #include "impr02.h"
  67 c
  68 c 0.3. ==> arguments
  69 c
  70       character*8 nohman
  71 c
  72       integer indnoe
  73       integer lgetco
  74       integer taetco(lgetco)
  75 c
  76       integer ulsort, langue, codret
  77 c
  78 c 0.4. ==> variables locales
  79 c
  80       integer codava
  81       integer nrosec
  82       integer nretap, nrsset
  83       integer iaux
  84       integer nuard1, nuarf1, nuard2, nuarf2
  85 c
  86       character*6 saux
  87 c
  88       integer nbmess
  89       parameter ( nbmess = 10 )
  90       character*80 texte(nblang,nbmess)
  91 c
  92 c 0.5. ==> initialisations
  93 c ______________________________________________________________________
  94 c
  95 c====
  96 c 1. messages
  97 c====
  98 c
  99       codava = codret
 100 c
 101 c=======================================================================
 102       if ( codava.eq.0 ) then
 103 c=======================================================================
 104 c
 105 c 1.1. ==> le debut des mesures de temps
 106 c
 107       nrosec = taetco(4)
 108       call gtdems (nrosec)
 109 c
 110 c 1.3. ==> les messages
 111 c
 112 #include "impr01.h"
 113 c
 114       texte(1,4) = '(/,a6,'' CREATION DES NOEUDS SUPPLEMENTAIRES'')'
 115       texte(1,5) = '(42(''=''),/)'
 116       texte(1,6) =
 117      > '(5x,''On attendait un nombre total de noeuds de : '',i10)'
 118       texte(1,7) =
 119      > '(5x,''Mais on en a cree un nombre total de      : '',i10)'
 120 c
 121       texte(2,4) = '(/,a6,'' CREATION OF ADDITIONAL NODES'')'
 122       texte(2,5) = '(35(''=''),/)'
 123       texte(2,6) = '(5x,i10,'' nodes were expected.'')'
 124       texte(2,7) = '(5x,i10,'' nodes were created.'')'
 125 c
 126 #include "impr03.h"
 127 c
 128 c 1.4. ==> le numero de sous-etape
 129 c
 130       nretap = taetco(1)
 131       nrsset = taetco(2) + 1
 132       taetco(2) = nrsset
 133 c
 134       call utcvne ( nretap, nrsset, saux, iaux, codret )
 135 c
 136 c 1.5. ==> le titre
 137 c
 138 #ifdef _DEBUG_HOMARD_
 139       write (ulsort,texte(langue,1)) 'Entree', nompro
 140       call dmflsh (iaux)
 141 #endif
 142 c
 143       write ( ulsort,texte(langue,4)) saux
 144       write ( ulsort,texte(langue,5))
 145 c
 146 #ifdef _DEBUG_HOMARD_
 147       write (ulsort,90002) 'indnoe', indnoe
 148       write (ulsort,90002) 'permar', permar
 149       write (ulsort,90002) 'nouvtr', nouvtr
 150 #endif
 151 c
 152 c====
 153 c 2. Prealables
 154 c    Il faut trier pour creer d'abord les noeuds supplementaires
 155 c    attaches aux entites permanentes, puis ceux attaches aux entites
 156 c    provisoires. Cela est indispensable pour traiter correctement
 157 c    la future suppression de la conformite.
 158 c====
 159 c
 160       if ( codret.eq.0 ) then
 161 c
 162       nuard1 = nbarma + 1
 163 c
 164       if ( mod(mailet,2).eq.0 .or.
 165      >     mod(mailet,3).eq.0 .or.
 166      >     mod(mailet,5).eq.0 ) then
 167         nuarf1 = permar
 168         nuard2 = nuarf1 + 1
 169         nuarf2 = nouvar
 170       else
 171         nuarf1 = nouvar
 172         nuard2 = nuarf1 + 1
 173         nuarf2 = -1
 174       endif
 175 c
 176       endif
 177 c
 178 c====
 179 c 2. Creation des noeuds P2 sur les aretes :
 180 c    . permanentes tout le temps
 181 c    . provisoires si pas de noeuds internes
 182 c====
 183 #ifdef _DEBUG_HOMARD_
 184       write (ulsort,*) '2. noeuds P2 permanents ; codret = ', codret
 185       write (ulsort,90002) 'nuard1, nuarf1', nuard1, nuarf1
 186 #endif
 187 c
 188       if ( codret.eq.0 ) then
 189 c
 190 #ifdef _DEBUG_HOMARD_
 191       write (ulsort,texte(langue,3)) 'CMNO22 permanent', nompro
 192 #endif
 193       call cmno22 ( nohman,
 194      >              indnoe, nuard1, nuarf1,
 195      >              ulsort, langue, codret )
 196 c
 197       endif
 198 c
 199 c====
 200 c 3. Creation des noeuds internes aux mailles permanentes
 201 c====
 202 #ifdef _DEBUG_HOMARD_
 203       write (ulsort,*) '3. internes permanents ; codret = ', codret
 204 #endif
 205 c
 206       if ( mod(mailet,2).eq.0 .or.
 207      >     mod(mailet,3).eq.0 .or.
 208      >     mod(mailet,5).eq.0 ) then
 209 c
 210         if ( codret.eq.0 ) then
 211 c
 212         iaux = 0
 213 #ifdef _DEBUG_HOMARD_
 214       write (ulsort,texte(langue,3)) 'CMNOIM permanent', nompro
 215 #endif
 216         call cmnoim ( nohman,
 217      >                indnoe, iaux,
 218      >                ulsort, langue, codret )
 219 c
 220         endif
 221 c
 222       endif
 223 c
 224 c====
 225 c 4. Creation des noeuds P2 sur les aretes :
 226 c    . provisoires s'il existe des noeuds internes
 227 c====
 228 #ifdef _DEBUG_HOMARD_
 229       write (ulsort,*) '4. noeuds P2 provisoires ; codret = ', codret
 230       write (ulsort,90002) 'nuard2, nuarf2', nuard2, nuarf2
 231 #endif
 232 c
 233       if ( nuarf2.ge.nuard2 ) then
 234 c
 235         if ( codret.eq.0 ) then
 236 c
 237 #ifdef _DEBUG_HOMARD_
 238       write (ulsort,texte(langue,3)) 'CMNO22 provisoire', nompro
 239 #endif
 240         call cmno22 ( nohman,
 241      >                indnoe, nuard2, nuarf2,
 242      >                ulsort, langue, codret )
 243 c
 244         endif
 245 c
 246       endif
 247 c
 248 c====
 249 c 5. Creation des noeuds internes aux mailles provisoires
 250 c====
 251 #ifdef _DEBUG_HOMARD_
 252       write (ulsort,*) '5. internes provisoires ; codret = ', codret
 253 #endif
 254 c
 255       if ( mod(mailet,2).eq.0 .or.
 256      >     mod(mailet,3).eq.0 .or.
 257      >     mod(mailet,5).eq.0 ) then
 258 c
 259         if ( codret.eq.0 ) then
 260 c
 261         iaux = 1
 262 #ifdef _DEBUG_HOMARD_
 263       write (ulsort,texte(langue,3)) 'CMNOIM provisoire', nompro
 264 #endif
 265         call cmnoim ( nohman,
 266      >                indnoe, iaux,
 267      >                ulsort, langue, codret )
 268 c
 269         endif
 270 c
 271       endif
 272
 273 c====
 274 c 6. verifications et impressions
 275 c====
 276 c
 277 #ifdef _DEBUG_HOMARD_
 278       write (ulsort,90002) 'indnoe', indnoe
 279 #endif
 280 c
 281       if ( codret.eq.0 ) then
 282 c
 283       if ( indnoe.ne.nouvno ) then
 284         write(ulsort,texte(langue,6)) nouvno
 285         write(ulsort,texte(langue,7)) indnoe
 286         codret = 4
 287       endif
 288 c
 289       endif
 290 c
 291 c====
 292 c 7. la fin
 293 c====
 294 c
 295 c 7.1. ==> message si erreur
 296 c
 297       if ( codret.ne.0 ) then
 298 c
 299 #include "envex2.h"
 300 c
 301       write (ulsort,texte(langue,1)) 'Sortie', nompro
 302       write (ulsort,texte(langue,2)) codret
 303 c
 304       endif
 305 c
 306 c 7.2. ==> fin des mesures de temps de la section
 307 c
 308       call gtfims (nrosec)
 309 c
 310 #ifdef _DEBUG_HOMARD_
 311       write (ulsort,texte(langue,1)) 'Sortie', nompro
 312       call dmflsh (iaux)
 313 #endif
 314 c
 315 c=======================================================================
 316       endif
 317 c=======================================================================
 318 c
 319       end