Saturne/cata_saturne.py

   1 # -*- coding: utf-8 -*-
   2 # debut entete
   3 import Accas
   4 from Accas import AU_MOINS_UN,UN_PARMI,EXCLUS,PRESENT_PRESENT
   5 from Accas import ENSEMBLE,PRESENT_ABSENT
   6 from Accas import JDC_CATA,NIVEAU
   7 from Accas import OPER,PROC,MACRO
   8 from Accas import SIMP,FACT,BLOC
   9 from Accas import ASSD
  10 #
  11 __version__="$Name:  $"
  12 __Id__="$Id: cata_saturne.py,v 1.1.1.1 2002/03/26 09:08:47 eficas Exp $"
  13 #
  14
  15 JdC = JDC_CATA(code='SATURNE',
  16                execmodul=None,
  17                regles = (AU_MOINS_UN('CALCUL_SATURNE'),),
  18                niveaux=(NIVEAU(nom='SONDES',label='Définition des sondes'),
  19                      NIVEAU(nom='EQUATIONS',label='Définition des équations'),
  20                      NIVEAU(nom='CALCUL',label ='Données du calcul')))
  21
  22 # P. RASCLE MMN
  23 # remarques diverses sur le catalogue Saturne
  24 # - dans les blocs, il faut au moins un mot clé de statut obligatoire
  25 # probleme de rafraichissement des blocs dépendants quand la valeur d'un mot cle global (ITURB) passe de 1 à 0
  26
  27 class sonde(ASSD):pass
  28 class varsca(ASSD):pass
  29 class flusca(ASSD):pass
  30 class varpre(ASSD):pass
  31 class varvitx(ASSD):pass
  32 class varvity(ASSD):pass
  33 class varvitz(ASSD):pass
  34 class eturb(ASSD):pass
  35 class dturb(ASSD):pass
  36 class tsr11(ASSD):pass
  37 class tsr22(ASSD):pass
  38 class tsr33(ASSD):pass
  39 class tsr12(ASSD):pass
  40 class tsr13(ASSD):pass
  41 class tsr23(ASSD):pass
  42 class resti(ASSD):pass
  43
  44 # fin entete
  45
  46 DEFI_SONDE = OPER(nom="DEFI_SONDE",op= 1,sd_prod=sonde,fr="définition d'une sonde historique avec ses coordonnées",
  47                   niveau = 'SONDES',
  48                   X = SIMP(statut ='o',typ='R',
  49                            fr="coordonnée X de la sonde"),
  50                   Y = SIMP(statut ='o',typ='R',
  51                            fr="coordonnée Y de la sonde"),
  52                   Z = SIMP(statut ='o',typ='R',
  53                            fr="coordonnée Z de la sonde")
  54                  );
  55
  56 DEFI_SCALA = OPER(nom="DEFI_SCALA",op=2,sd_prod=varsca,fr="définition d'une inconnue scalaire avec ses paramètres physico numériques",
  57                   niveau = 'EQUATIONS',
  58                   NUMERIC = FACT(statut='o',max=01,
  59                                 fr="propriétés numériques liées à l'inconnue",
  60                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
  61                                                   fr="indicateur de convection"),
  62                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
  63                                                   fr="présence terme instationnaire dans les matrices"),
  64                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
  65                                                   fr="diffusion"),
  66                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
  67                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
  68                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
  69                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
  70                                                                           fr="pourcentage schéma convectif au second ordre"),
  71                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
  72                                                                           fr ="type schéma convectif au second ordre 1 : centré")
  73                                                             ),
  74                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
  75                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
  76                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
  77                                                   fr="précision relative pour la résolution des systèmes linéaires"),
  78                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
  79                                                   fr="type de limitation des gradients"),
  80                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
  81                                                   fr="reconstruction des gradients"),
  82                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
  83                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
  84                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
  85                                                   fr="facteur de limitation des gradients"),
  86                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
  87                                                   fr="précision relative pour la reconstruction itérative des gradients")
  88                                  ),
  89                   RESTITUE = FACT(statut='o',max=01,
  90                                   fr = "Description de la grandeur restituée",
  91                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
  92                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
  93                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
  94                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
  95                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
  96                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
  97                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
  98                                   ),
  99                   MODELE = FACT(statut='o',max=01,
 100                                 fr="modélisation",
 101                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 102                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 103                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 104                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 105                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 106                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 107                                 )
 108                  );
 109
 110 DEFI_FLUSCA = OPER(nom="DEFI_FLUSCA",op=2,sd_prod=flusca,fr="définition des fluctuations d'une inconnue scalaire avec ses paramètres physico numériques",
 111                    niveau = 'EQUATIONS',
 112                   VARFL  = SIMP(statut='o',typ=varsca,fr="scalaire associé au calcul de la variance des fluctuations"),
 113                   NUMERIC = FACT(statut='o',max=01,
 114                                 fr="propriétés numériques liées à l'inconnue",
 115                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 116                                                   fr="indicateur de convection"),
 117                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 118                                                   fr="présence terme instationnaire dans les matrices"),
 119                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 120                                                   fr="diffusion"),
 121                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 122                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 123                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 124                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 125                                                                           fr="pourcentage schéma convectif au second ordre"),
 126                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 127                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 128                                                             ),
 129                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 130                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 131                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 132                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 133                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 134                                                   fr="type de limitation des gradients"),
 135                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 136                                                   fr="reconstruction des gradients"),
 137                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 138                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 139                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 140                                                   fr="facteur de limitation des gradients"),
 141                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 142                                                   fr="précision relative pour la reconstruction itérative des gradients")
 143                                  ),
 144                   RESTITUE = FACT(statut='o',max=01,
 145                                   fr = "Description de la grandeur restituée",
 146                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 147                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 148                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 149                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 150                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 151                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 152                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 153                                   ),
 154                   MODELE = FACT(statut='o',max=01,
 155                                 fr="modélisation",
 156                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 157                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 158                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 159                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 160                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 161                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable"),
 162                                 RVARFL = SIMP(statut='o',typ='R',val_min=0,defaut=0.8,
 163                                                        fr="nombre de Prandtl pour la variance des fluctuations du scalaire")
 164                                 )
 165                   );
 166
 167 DEFI_PRESSION = OPER(nom="DEFI_PRESSION",op=2,sd_prod=varpre,fr="définition del'inconnue pression avec ses paramètres physico numériques",
 168                      niveau = 'EQUATIONS',
 169                    NUMERIC = FACT(statut='o',max=01,
 170                                   fr="propriétés de l'équation de pression",
 171                                   ICONV  = SIMP(statut='o',typ='I',defaut=0,into=(0,1),
 172                                                 fr="indicateur de convection"),
 173                                   ISTAT  = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 174                                                 fr="présence terme instationnaire dans les matrices"),
 175                                   IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 176                                                 fr="diffusion"),
 177                                   IDIRCL = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 178                                                 fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 179                                   SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 180                                                           BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 181                                                                         fr="pourcentage schéma convectif au second ordre"),
 182                                                           ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 183                                                                         fr ="type schéma convectif au second ordre 1 : centré")
 184                                                           ),
 185                                   NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 186                                                 fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 187                                   EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 188                                                 fr="précision relative pour la résolution des systèmes linéaires"),
 189                                   IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 190                                                 fr="type de limitation des gradients"),
 191                                   NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 192                                                 fr="reconstruction des gradients"),
 193                                   NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 194                                                 fr="nombre d'itérations pour la reconstruction des seconds membres"),
 195                                   CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 196                                                 fr="facteur de limitation des gradients"),
 197                                   EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 198                                                 fr="précision relative pour la reconstruction itérative des gradients")
 199                                   ),
 200                   RESTITUE = FACT(statut='o',max=01,
 201                                   fr = "Description de la grandeur restituée",
 202                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée",
 203                                                 defaut="Pression"),
 204                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 205                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 206                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 207                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 208                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 209                                                 fr = "inventaire des sondes utilisées pour la grandeur")
 210                                   ),
 211                   MODELE = FACT(statut='o',max=01,
 212                                 fr="modélisation",
 213                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 214                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour la grandeur"),
 215                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour la grandeur"),
 216                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour la grandeur"),
 217                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée à la grandeur"),
 218                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 219                                 )
 220                     );
 221
 222 DEFI_VITX = OPER(nom="DEFI_VITX",op=2,sd_prod=varvitx,fr="définition de l'inconnue vitesse X avec ses paramètres physico numériques",
 223                  niveau = 'EQUATIONS',
 224                    NUMERIC = FACT(statut='o',max=01,
 225                                   fr="propriétés de l'équation de vitesse X",
 226                                   ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 227                                                 fr="indicateur de convection"),
 228                                   ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 229                                                 fr="présence terme instationnaire dans les matrices"),
 230                                   IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 231                                                 fr="diffusion"),
 232                                   IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 233                                                 fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 234                                   SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 235                                                           BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 236                                                                         fr="pourcentage schéma convectif au second ordre"),
 237                                                           ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 238                                                                         fr ="type schéma convectif au second ordre 1 : centré")
 239                                                           ),
 240                                   NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 241                                                 fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 242                                   EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 243                                                 fr="précision relative pour la résolution des systèmes linéaires"),
 244                                   IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 245                                                 fr="type de limitation des gradients"),
 246                                   NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 247                                                 fr="reconstruction des gradients"),
 248                                   NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 249                                                 fr="nombre d'itérations pour la reconstruction des seconds membres"),
 250                                   CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 251                                                 fr="facteur de limitation des gradients"),
 252                                   EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 253                                                 fr="précision relative pour la reconstruction itérative des gradients")
 254                                   ),
 255                   RESTITUE = FACT(statut='o',max=01,
 256                                   fr = "Description de la grandeur restituée",
 257                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée",
 258                                                 defaut="vitesse_u1"),
 259                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 260                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 261                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 262                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 263                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 264                                                 fr = "inventaire des sondes utilisées pour la grandeur")
 265                                   ),
 266                   MODELE = FACT(statut='o',max=01,
 267                                 fr="modélisation",
 268                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 269                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour la grandeur"),
 270                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour la grandeur"),
 271                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour la grandeur"),
 272                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée à la grandeur"),
 273                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 274                                 )
 275                 );
 276
 277 DEFI_VITY = OPER(nom="DEFI_VITY",op=2,sd_prod=varvity,fr="définition de l'inconnue vitesse Y avec ses paramètres physico numériques",
 278                  niveau = 'EQUATIONS',
 279                    NUMERIC = FACT(statut='o',max=01,
 280                                   fr="propriétés de l'équation de vitesse Y",
 281                                   ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 282                                                 fr="indicateur de convection"),
 283                                   ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 284                                                 fr="présence terme instationnaire dans les matrices"),
 285                                   IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 286                                                 fr="diffusion"),
 287                                   IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 288                                                 fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 289                                   SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 290                                                           BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 291                                                                         fr="pourcentage schéma convectif au second ordre"),
 292                                                           ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 293                                                                         fr ="type schéma convectif au second ordre 1 : centré")
 294                                                           ),
 295                                   NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 296                                                 fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 297                                   EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 298                                                 fr="précision relative pour la résolution des systèmes linéaires"),
 299                                   IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 300                                                 fr="type de limitation des gradients"),
 301                                   NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 302                                                 fr="reconstruction des gradients"),
 303                                   NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 304                                                 fr="nombre d'itérations pour la reconstruction des seconds membres"),
 305                                   CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 306                                                 fr="facteur de limitation des gradients"),
 307                                   EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 308                                                 fr="précision relative pour la reconstruction itérative des gradients")
 309                                   ),
 310                   RESTITUE = FACT(statut='o',max=01,
 311                                   fr = "Description de la grandeur restituée",
 312                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée",
 313                                                 defaut="vitesse_v1"),
 314                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 315                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 316                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 317                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 318                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 319                                                 fr = "inventaire des sondes utilisées pour la grandeur")
 320                                   ),
 321                   MODELE = FACT(statut='o',max=01,
 322                                 fr="modélisation",
 323                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 324                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour la grandeur"),
 325                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour la grandeur"),
 326                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour la grandeur"),
 327                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée à la grandeur"),
 328                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 329                                 )
 330                 );
 331
 332 DEFI_VITZ = OPER(nom="DEFI_VITZ",op=2,sd_prod=varvitz,fr="définition de l'inconnue vitesse Z avec ses paramètres physico numériques",
 333                  niveau = 'EQUATIONS',
 334                    NUMERIC = FACT(statut='o',max=01,
 335                                   fr="propriétés de l'équation de vitesse Z",
 336                                   ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 337                                                 fr="indicateur de convection"),
 338                                   ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 339                                                 fr="présence terme instationnaire dans les matrices"),
 340                                   IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 341                                                 fr="diffusion"),
 342                                   IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 343                                                 fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 344                                   SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 345                                                           BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 346                                                                         fr="pourcentage schéma convectif au second ordre"),
 347                                                           ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 348                                                                         fr ="type schéma convectif au second ordre 1 : centré")
 349                                                           ),
 350                                   NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 351                                                 fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 352                                   EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 353                                                 fr="précision relative pour la résolution des systèmes linéaires"),
 354                                   IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 355                                                 fr="type de limitation des gradients"),
 356                                   NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 357                                                 fr="reconstruction des gradients"),
 358                                   NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 359                                                 fr="nombre d'itérations pour la reconstruction des seconds membres"),
 360                                   CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 361                                                 fr="facteur de limitation des gradients"),
 362                                   EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 363                                                 fr="précision relative pour la reconstruction itérative des gradients")
 364                                   ),
 365                   RESTITUE = FACT(statut='o',max=01,
 366                                   fr = "Description de la grandeur restituée",
 367                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée",
 368                                                 defaut="vitesse_w1"),
 369                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 370                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 371                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 372                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 373                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 374                                                 fr = "inventaire des sondes utilisées pour la grandeur")
 375                                   ),
 376                   MODELE = FACT(statut='o',max=01,
 377                                 fr="modélisation",
 378                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 379                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour la grandeur"),
 380                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour la grandeur"),
 381                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour la grandeur"),
 382                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée à la grandeur"),
 383                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable"),
 384                                 )
 385                 );
 386
 387
 388 DEFI_ETURB = OPER(nom="DEFI_ETURB",op=2,sd_prod=eturb,fr="définition de l'inconnue energie turbulente k-eps avec ses paramètres physico numériques",
 389                   niveau = 'EQUATIONS',
 390                   NUMERIC = FACT(statut='o',max=01,
 391                                     fr="propriétés numériques liées à l'inconnue",
 392                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 393                                                   fr="indicateur de convection"),
 394                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 395                                                   fr="présence terme instationnaire dans les matrices"),
 396                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 397                                                   fr="diffusion"),
 398                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 399                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 400                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 401                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 402                                                                           fr="pourcentage schéma convectif au second ordre"),
 403                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 404                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 405                                                             ),
 406                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 407                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 408                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 409                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 410                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 411                                                   fr="type de limitation des gradients"),
 412                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 413                                                   fr="reconstruction des gradients"),
 414                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 415                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 416                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 417                                                   fr="facteur de limitation des gradients"),
 418                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 419                                                   fr="précision relative pour la reconstruction itérative des gradients")
 420                                  ),
 421                   RESTITUE = FACT(statut='o',max=01,
 422                                   fr = "Description de la grandeur restituée",
 423                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 424                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 425                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 426                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 427                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 428                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 429                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 430                                   ),
 431                   MODELE = FACT(statut='o',max=01,
 432                                 fr="modélisation",
 433                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 434                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 435                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 436                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 437                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 438                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 439                                 )
 440                  );
 441
 442 DEFI_DTURB = OPER(nom="DEFI_DTURB",op=2,sd_prod=dturb,fr="définition de l'inconnue dissipation turbulente k-eps avec ses paramètres physico numériques",
 443                   niveau = 'EQUATIONS',
 444                   NUMERIC = FACT(statut='o',max=01,
 445                                 fr="propriétés numériques liées à l'inconnue",
 446                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 447                                                   fr="indicateur de convection"),
 448                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 449                                                   fr="présence terme instationnaire dans les matrices"),
 450                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 451                                                   fr="diffusion"),
 452                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 453                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 454                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 455                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 456                                                                           fr="pourcentage schéma convectif au second ordre"),
 457                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 458                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 459                                                             ),
 460                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 461                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 462                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 463                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 464                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 465                                                   fr="type de limitation des gradients"),
 466                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 467                                                   fr="reconstruction des gradients"),
 468                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 469                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 470                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 471                                                   fr="facteur de limitation des gradients"),
 472                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 473                                                   fr="précision relative pour la reconstruction itérative des gradients")
 474                                  ),
 475                   RESTITUE = FACT(statut='o',max=01,
 476                                   fr = "Description de la grandeur restituée",
 477                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 478                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 479                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 480                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 481                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 482                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 483                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 484                                   ),
 485                   MODELE = FACT(statut='o',max=01,
 486                                 fr="modélisation",
 487                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 488                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 489                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 490                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 491                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 492                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 493                                 )
 494                  );
 495
 496 DEFI_TSR11 = OPER(nom="DEFI_TSR11",op=2,sd_prod=tsr11,fr="définition de l'inconnue tension Reynolds R11 Rij-eps avec ses paramètres physico numériques",
 497                   niveau = 'EQUATIONS',
 498                   NUMERIC = FACT(statut='o',max=01,
 499                                 fr="propriétés numériques liées à l'inconnue",
 500                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 501                                                   fr="indicateur de convection"),
 502                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 503                                                   fr="présence terme instationnaire dans les matrices"),
 504                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 505                                                   fr="diffusion"),
 506                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 507                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 508                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 509                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 510                                                                           fr="pourcentage schéma convectif au second ordre"),
 511                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 512                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 513                                                             ),
 514                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 515                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 516                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 517                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 518                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 519                                                   fr="type de limitation des gradients"),
 520                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 521                                                   fr="reconstruction des gradients"),
 522                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 523                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 524                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 525                                                   fr="facteur de limitation des gradients"),
 526                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 527                                                   fr="précision relative pour la reconstruction itérative des gradients")
 528                                  ),
 529                   RESTITUE = FACT(statut='o',max=01,
 530                                   fr = "Description de la grandeur restituée",
 531                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 532                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 533                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 534                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 535                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 536                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 537                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 538                                   ),
 539                   MODELE = FACT(statut='o',max=01,
 540                                 fr="modélisation",
 541                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 542                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 543                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 544                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 545                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 546                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 547                                 )
 548                  );
 549
 550 DEFI_TSR22 = OPER(nom="DEFI_TSR22",op=2,sd_prod=tsr11,fr="définition de l'inconnue tension Reynolds R22 Rij-eps avec ses paramètres physico numériques",
 551                   niveau = 'EQUATIONS',
 552                   NUMERIC = FACT(statut='o',max=01,
 553                                 fr="propriétés numériques liées à l'inconnue",
 554                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 555                                                   fr="indicateur de convection"),
 556                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 557                                                   fr="présence terme instationnaire dans les matrices"),
 558                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 559                                                   fr="diffusion"),
 560                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 561                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 562                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 563                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 564                                                                           fr="pourcentage schéma convectif au second ordre"),
 565                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 566                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 567                                                             ),
 568                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 569                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 570                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 571                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 572                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 573                                                   fr="type de limitation des gradients"),
 574                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 575                                                   fr="reconstruction des gradients"),
 576                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 577                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 578                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 579                                                   fr="facteur de limitation des gradients"),
 580                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 581                                                   fr="précision relative pour la reconstruction itérative des gradients")
 582                                  ),
 583                   RESTITUE = FACT(statut='o',max=01,
 584                                   fr = "Description de la grandeur restituée",
 585                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 586                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 587                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 588                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 589                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 590                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 591                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 592                                   ),
 593                   MODELE = FACT(statut='o',max=01,
 594                                 fr="modélisation",
 595                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 596                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 597                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 598                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 599                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 600                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 601                                 )
 602                  );
 603
 604 DEFI_TSR33 = OPER(nom="DEFI_TSR33",op=2,sd_prod=tsr11,fr="définition de l'inconnue tension Reynolds R33 Rij-eps avec ses paramètres physico numériques",
 605                   niveau = 'EQUATIONS',
 606                   NUMERIC = FACT(statut='o',max=01,
 607                                 fr="propriétés numériques liées à l'inconnue",
 608                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 609                                                   fr="indicateur de convection"),
 610                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 611                                                   fr="présence terme instationnaire dans les matrices"),
 612                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 613                                                   fr="diffusion"),
 614                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 615                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 616                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 617                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 618                                                                           fr="pourcentage schéma convectif au second ordre"),
 619                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 620                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 621                                                             ),
 622                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 623                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 624                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 625                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 626                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 627                                                   fr="type de limitation des gradients"),
 628                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 629                                                   fr="reconstruction des gradients"),
 630                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 631                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 632                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 633                                                   fr="facteur de limitation des gradients"),
 634                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 635                                                   fr="précision relative pour la reconstruction itérative des gradients")
 636                                  ),
 637                   RESTITUE = FACT(statut='o',max=01,
 638                                   fr = "Description de la grandeur restituée",
 639                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 640                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 641                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 642                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 643                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 644                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 645                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 646                                   ),
 647                   MODELE = FACT(statut='o',max=01,
 648                                 fr="modélisation",
 649                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 650                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 651                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 652                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 653                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 654                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 655                                 )
 656                  );
 657
 658 DEFI_TSR12 = OPER(nom="DEFI_TSR12",op=2,sd_prod=tsr11,fr="définition de l'inconnue tension Reynolds R12 Rij-eps avec ses paramètres physico numériques",
 659                   niveau = 'EQUATIONS',
 660                   NUMERIC = FACT(statut='o',max=01,
 661                                 fr="propriétés numériques liées à l'inconnue",
 662                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 663                                                   fr="indicateur de convection"),
 664                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 665                                                   fr="présence terme instationnaire dans les matrices"),
 666                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 667                                                   fr="diffusion"),
 668                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 669                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 670                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 671                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 672                                                                           fr="pourcentage schéma convectif au second ordre"),
 673                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 674                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 675                                                             ),
 676                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 677                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 678                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 679                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 680                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 681                                                   fr="type de limitation des gradients"),
 682                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 683                                                   fr="reconstruction des gradients"),
 684                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 685                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 686                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 687                                                   fr="facteur de limitation des gradients"),
 688                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 689                                                   fr="précision relative pour la reconstruction itérative des gradients")
 690                                  ),
 691                   RESTITUE = FACT(statut='o',max=01,
 692                                   fr = "Description de la grandeur restituée",
 693                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 694                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 695                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 696                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 697                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 698                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 699                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 700                                   ),
 701                   MODELE = FACT(statut='o',max=01,
 702                                 fr="modélisation",
 703                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 704                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 705                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 706                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 707                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 708                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 709                                 )
 710                  );
 711
 712 DEFI_TSR13 = OPER(nom="DEFI_TSR13",op=2,sd_prod=tsr11,fr="définition de l'inconnue tension Reynolds R13 Rij-eps avec ses paramètres physico numériques",
 713                   niveau = 'EQUATIONS',
 714                   NUMERIC = FACT(statut='o',max=01,
 715                                 fr="propriétés numériques liées à l'inconnue",
 716                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 717                                                   fr="indicateur de convection"),
 718                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 719                                                   fr="présence terme instationnaire dans les matrices"),
 720                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 721                                                   fr="diffusion"),
 722                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 723                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 724                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 725                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 726                                                                           fr="pourcentage schéma convectif au second ordre"),
 727                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 728                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 729                                                             ),
 730                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 731                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 732                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 733                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 734                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 735                                                   fr="type de limitation des gradients"),
 736                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 737                                                   fr="reconstruction des gradients"),
 738                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 739                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 740                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 741                                                   fr="facteur de limitation des gradients"),
 742                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 743                                                   fr="précision relative pour la reconstruction itérative des gradients")
 744                                  ),
 745                   RESTITUE = FACT(statut='o',max=01,
 746                                   fr = "Description de la grandeur restituée",
 747                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 748                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 749                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 750                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 751                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 752                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 753                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 754                                   ),
 755                   MODELE = FACT(statut='o',max=01,
 756                                 fr="modélisation",
 757                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 758                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 759                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 760                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 761                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 762                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 763                                 )
 764                  );
 765
 766 DEFI_TSR23 = OPER(nom="DEFI_TSR23",op=2,sd_prod=tsr11,fr="définition de l'inconnue tension Reynolds R23 Rij-eps avec ses paramètres physico numériques",
 767                   niveau = 'EQUATIONS',
 768                   NUMERIC = FACT(statut='o',max=01,
 769                                 fr="propriétés numériques liées à l'inconnue",
 770                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 771                                                   fr="indicateur de convection"),
 772                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 773                                                   fr="présence terme instationnaire dans les matrices"),
 774                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 775                                                   fr="diffusion"),
 776                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 777                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 778                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 779                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 780                                                                           fr="pourcentage schéma convectif au second ordre"),
 781                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 782                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 783                                                             ),
 784                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 785                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 786                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 787                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 788                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 789                                                   fr="type de limitation des gradients"),
 790                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 791                                                   fr="reconstruction des gradients"),
 792                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 793                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 794                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 795                                                   fr="facteur de limitation des gradients"),
 796                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 797                                                   fr="précision relative pour la reconstruction itérative des gradients")
 798                                  ),
 799                   RESTITUE = FACT(statut='o',max=01,
 800                                   fr = "Description de la grandeur restituée",
 801                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 802                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 803                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 804                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 805                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 806                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 807                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 808                                   ),
 809                   MODELE = FACT(statut='o',max=01,
 810                                 fr="modélisation",
 811                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 812                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 813                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 814                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 815                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 816                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 817                                 )
 818                  );
 819
 820
 821 DEFI_RESTI = OPER(nom="DEFI_RESTI",op=2,sd_prod=resti,fr = "Description de la grandeur restituée",
 822                   niveau = 'CALCUL',
 823                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 824                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 825                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 826                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 827                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 828                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 829                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 830                   );
 831
 832 CALCUL_SATURNE = PROC(nom = "CALCUL_SATURNE",op = 1,
 833                       fr = "définition des paramètres généraux pour un calcul Saturne",
 834                       docu = "néant",reentrant = 'n',
 835                       niveau = 'CALCUL',
 836         ENVELOPPE = FACT(statut='o',max=01,fr = "gestion de l'enveloppe",
 837                          IFOENV = SIMP(statut='o',typ='I',defaut=2,into=(0,1,2),position='global',
 838                                        fr = "mode de communication enveloppe solveur"),
 839                          ENVTOSOLV = BLOC(condition="IFOENV > 0",
 840                                           fr = "liaison enveloppe vers solveur",
 841                                           IMPEVI = SIMP(statut='o',typ='I',defaut=13,val_min=0,val_max=99,
 842                                                         fr = "numéro de fichier enveloppe vers solveur"),
 843                                           FICEVI = SIMP(statut='o',typ='TXM',defaut="enveloppe_vers_solveur          ",
 844                                                         fr = "nom de fichier enveloppe vers solveur")
 845                                          ),
 846                          SOLVTOENV = BLOC(condition="IFOENV > 0",
 847                                           fr = "liaison solveur vers enveloppe",
 848                                           IMPEVO = SIMP(statut='o',typ='I',defaut=14,val_min=0,val_max=99,
 849                                                         fr = "numéro de fichier solveur vers enveloppe"),
 850                                           FICEVO = SIMP(statut='o',typ='TXM',defaut="solveur_vers_enveloppe          ",
 851                                                         fr = "nom de fichier solveur vers enveloppe")
 852                                          )
 853                         ),
 854         FICHIERS_CALCUL = FACT(statut='f',max=01,
 855                                fr ="définition des fichiers géométrie, suite, stop",
 856                                GEOMETRIE = BLOC(condition="IFOENV == 0",
 857                                                 fr = "fichier géométrique (pas d'enveloppe)",
 858                                                 IMPGEO = SIMP(statut='o',typ='I',defaut=10,val_min=0,val_max=99,
 859                                                               fr = "numéro de fichier géométrique"),
 860                                                 FICGEO = SIMP(statut='o',typ='TXM',defaut="geomet",
 861                                                               fr = "nom de fichier géométrique")
 862                                                ),
 863                                SUITE_AMONT = BLOC(condition="ISUITE == 1",
 864                                                   fr = "fichier suite amont",
 865                                                   IMPAMO = SIMP(statut='o',typ='I',defaut=11,val_min=0,val_max=99,
 866                                                               fr = "numéro de fichier suite amont"),
 867                                                   FICGEO = SIMP(statut='o',typ='TXM',defaut="suiamo",
 868                                                               fr = "nom de fichier suite amont"),
 869                                                   IFOAMO = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 870                                                               fr = "format du fichier suite amont, 0 = binaire")
 871                                                  ),
 872                                FICHIER_STOP = FACT(fr = "fichier stop",
 873                                                    IMPSTP = SIMP(statut='o',typ='I',defaut=12,val_min=0,val_max=99,
 874                                                                  fr = "numéro de fichier stop"),
 875                                                    FICSTP = SIMP(statut='o',typ='TXM',defaut="ficstp",
 876                                                                  fr = "nom de fichier stop")
 877                                                   ),
 878                                SUITE_AVAL = FACT(fr = "fichier suite aval",
 879                                                  IMPAVA = SIMP(statut='o',typ='I',defaut=20,val_min=0,val_max=99,
 880                                                              fr = "numéro de fichier suite aval"),
 881                                                  FICAVA = SIMP(statut='o',typ='TXM',defaut="suiamo",
 882                                                              fr = "nom de fichier suite aval"),
 883                                                  IFOAVA = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 884                                                              fr = "format du fichier suite aval, 0 = binaire")
 885                                                 )
 886                               ),
 887         POST_PROC_ENSIGHT = FACT(statut='f',max=01,
 888                                  fr = "options de post processing par Ensight",
 889                                  IFOENS = SIMP(statut='f',typ='I',defaut=-2,into=(-2,-1,0,1),
 890                                                fr = "option type de fichier et processus générateur"),
 891                                  NTCHR  = SIMP(statut='f',typ='I',defaut=-1,
 892                                                fr = "périodicité de sortie des fichiers Ensight"),
 893                                  ITCHR  = SIMP(statut='f',typ='I',defaut=0,val_min=0,
 894                                                fr = "compteur des sorties des fichiers Ensight"),
 895                                  N3S_ASCII_NOY = BLOC(condition="(IFOENS == -1) and (IFOENV == 0)",
 896                                                       fr = "format Ensight N3S ASCII généré par le noyau",
 897                                                       IMPPST = SIMP(statut='o',typ='I',defaut=21,val_min=0,val_max=99,
 898                                                                     fr = "numéro de fichier Ensight"),
 899                                                       FICPST = SIMP(statut='o',typ='TXM',defaut="dessin",
 900                                                                     fr = "nom de fichier Ensight")
 901                                                       ),
 902                                  P0_NOY =        BLOC(condition="((IFOENS == 0) or (IFOENS == 1)) and (IFOENV == 0)",
 903                                                       fr = "format Ensight P0 généré par le noyau",
 904                                                       IMPEP0 = SIMP(statut='o',typ='I',defaut=22,val_min=0,val_max=99,
 905                                                                     fr = "numéro de fichier Ensight"),
 906                                                       EMPCHR = SIMP(statut='o',typ='TXM',defaut="./",
 907                                                                     fr = "répertoire de fichier Ensight"),
 908                                                       ENTCHR = SIMP(statut='o',typ='TXM',defaut="chr",
 909                                                                     fr = "préfixe nom de fichier Ensight")
 910                                                       )
 911                                 ),
 912         HISTORIQUE_PONCTUEL = FACT(statut='o',max=01,
 913                                    fr = "Sondes historiques",
 914                                    FICHIERS_HISTORIQUES = FACT(statut='f',max=01,
 915                                                                fr = "description des fichiers historiques",
 916                                                                EMPHIS = SIMP(statut='o',typ='TXM',defaut="./",
 917                                                                              fr="répertoire fichiers historiques"),
 918                                                                EXTHIS = SIMP(statut='o',typ='TXM',defaut="hst",
 919                                                                              fr="extension fichiers historiques")
 920                                                                ),
 921                                    NTHIST = SIMP(statut='f',typ='I',defaut=-999,
 922                                                  fr="fréquence de sortie des historiques en pas de temps"),
 923                                    NTHSAV = SIMP(statut='f',typ='I',defaut=-999,
 924                                                  fr="fréquence de sauvegarde des historiques en pas de temps")
 925                                   ),
 926         OPTIONS_TURBULENCE = FACT(statut='o',max=01,
 927                                   fr="modèle de turbulence",
 928                                   ITURB  = SIMP(statut='f',fr="laminaire : 0, k-epsilon :1,Rij-epsilon :2",
 929                                                 typ='I',into=(0,1,2),defaut=1,position='global'),
 930                                   MODTURB = BLOC(condition="ITURB == 1",
 931                                                  fr = "option k-epsilon",
 932                                                  IGRAKE = SIMP(statut='o',typ='I',
 933                                                                fr="prise en compte gravité dans k-epsilon",
 934                                                                into=(0,1),defaut=0),
 935                                                  IDEUCH = SIMP(statut='f',typ='I',
 936                                                                fr="prise en compte k-epsilon deux échelles",
 937                                                                into=(0,1),defaut=1),
 938                                                  IKEKOU = SIMP(statut='f',typ='I',
 939                                                                fr="prise en compte couplage en incréments sur k-epsilon",
 940                                                                into=(0,1),defaut=1)
 941                                                  ),
 942                                   TEMPTURB = BLOC(condition="ITURB == 1",
 943                                                   fr = "option k-epsilon",
 944                                                   ISCALT = SIMP(statut='o',typ=varsca,
 945                                                                 fr=" identificateur inconnue scalaire température ")
 946                                                   )
 947                                   ),
 948         MARCHE_TEMPS = FACT(statut='o',max=01,
 949                             fr = "définition de la marche en temps",
 950                             DTREF  = SIMP(statut='o',fr="pas de temps de référence",
 951                                                         typ='R',val_min=0),
 952                             IDTVAR = SIMP(statut='f',fr="pas de temps constant : 0, variable temps espace : 1, variable temps : 2",
 953                                           typ='I',into=(0,1,2),defaut=0,position='global'),
 954                                   # probleme trace eficas quand un mot cle position global change
 955                             PASVAR = BLOC(condition="IDTVAR != 0",fr="options pas de temps variable",
 956                                           XCFMAX = SIMP(statut='o',fr="nombre de Courant-Fourier cible",
 957                                                         typ='R',defaut=0.5,val_min=0),
 958                                           FMIN =   SIMP(statut='f',fr="rapport min pas calculé DTREF",
 959                                                         typ='R',defaut=0.1,val_min=0),
 960                                           FMAX =   SIMP(statut='f',fr="rapport max pas calculé DTREF",
 961                                                         typ='R',defaut=1000,val_min=0),
 962                                           VARRDT = SIMP(statut='f',fr="variation relative max pas calculé entre deux instants",
 963                                                         typ='R',defaut=0.1,val_min=0)
 964                                           )
 965                             ),
 966         OPTIONS_EQUATIONS = FACT(statut='o',max=01,
 967                                    fr = "propriétés des équations, inconnues principales",
 968                                    IMGR   = SIMP(statut='f',fr="utilisation du multigrille pour la résolution des systèmes linéaires",
 969                                                  typ='I',into=(0,1),defaut=0),
 970                                    IMRGRA = SIMP(statut='f',fr="type de reconstruction des gradients 1 : moindres carrés",
 971                                                  typ='I',into=(0,1),defaut=0),
 972                                    ),
 973         VARIABLES = FACT(statut='o',max=01,
 974                          fr = "Restitution des grandeurs principales",
 975                          NTLIST = SIMP(statut='f',typ='I',defaut=1,val_min=1,val_max='**',
 976                                        fr = "fréquence de sortie (en pas de temps) dans le compte rendu d'éxécution"),
 977                          IWARNI = SIMP(statut='f',typ='I',defaut=2,val_min=0,val_max='**',
 978                                        fr = "niveau de détail des impressions dans le compte rendu d'éxécution"),
 979                          MASVOL1 = FACT(statut='o',max=01,
 980                                          fr = "Description de la grandeur restituée : masse volumique",
 981                                          NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur : masse volumique",
 982                                                        defaut="Masse_vol1"),
 983                                          ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 984                                                        fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 985                                          ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 986                                                        fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 987                                          IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 988                                                        fr = "inventaire des sondes utilisées pour la grandeur masse volumique")
 989                                         ),
 990                         TURB_KE = BLOC(condition="ITURB == 1",
 991                                        E_TURB = FACT(statut='o',max=01,
 992                                                      fr = "Description de la grandeur restituée : énergie turbulente",
 993                                                      NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur : energie turbulente",
 994                                                                    defaut="Energie_1"),
 995                                                      ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 996                                                                    fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 997                                                      ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 998                                                                    fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 999                                                      IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
1000                                                                    fr = "inventaire des sondes utilisées pour la grandeur énergie turbulente")
1001                                                     ),
1002                                        D_TURB = FACT(statut='o',max=01,
1003                                                      fr = "Description de la grandeur restituée : dissipation turbulente",
1004                                                      NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur : dissipation turbulente",
1005                                                                    defaut="Dissipation"),
1006                                                      ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1007                                                                    fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
1008                                                      ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1009                                                                    fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
1010                                                      IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
1011                                                                    fr = "inventaire des sondes utilisées pour la grandeur dissipation turbulente")
1012                                                     ),
1013                                        V_TURB = FACT(statut='o',max=01,
1014                                                      fr = "Description de la grandeur restituée : viscosité turbulente",
1015                                                      NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur : viscosité turbulente",
1016                                                                    defaut="Visc_turb1"),
1017                                                      ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1018                                                                    fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
1019                                                      ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1020                                                                    fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
1021                                                      IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
1022                                                                    fr = "inventaire des sondes utilisées pour la grandeur viscosité turbulente")
1023                                                     ),
1024                                       )
1025                         ),
1026         GESTION_CALCUL = FACT(statut='o',max=01,
1027                               fr = "calcul suite et numéros de pas de temps début et fin",
1028                               ISUITE = SIMP(statut='o',typ='I',defaut=0,into=(0,1),position='global',
1029                                             fr = "indicateur calcul suite (1 = suite)"),
1030                               NTPABS = SIMP(statut='f',typ='I',defaut=0,
1031                                             fr = "numéro dernier pas de temps calcul précédent (initialisation automatique)"),
1032                               NTMABS = SIMP(statut='o',typ='I',
1033                                             fr = "numéro dernier pas de temps visé (absolu)"),
1034                               TTPABS = SIMP(statut='f',typ='R',defaut=0,
1035                                             fr = "temps simulation au dernier pas de temps précédent (initialisation automatique)")
1036                              ),
1037         CONSTANTES_PHYSIQUES=FACT(statut='o',max=01,
1038                                   fr = "Grandeurs physiques et modélisation",
1039                                   GRAVITE = FACT(statut='o',max=01,
1040                                                  fr = "composantes de la gravité",
1041                                                  GX = SIMP(statut ='o',typ='R',defaut=0, fr = "gravité selon X"),
1042                                                  GY = SIMP(statut ='o',typ='R',defaut=0, fr = "gravité selon Y"),
1043                                                  GZ = SIMP(statut ='o',typ='R',defaut=0, fr = "gravité selon Z")
1044                                                  ),
1045                                   FLUIDE = FACT(statut='o',max=01,
1046                                                 fr = "propriétés du fluide",
1047                                                 RO0 = SIMP(statut ='o',typ='R',val_min=0,fr="masse volumique de référence"),
1048                                                 VISCL0 = SIMP(statut ='o',typ='R',val_min=0,fr="viscosité dynamique de référence"),
1049                                                 P0 = SIMP(statut ='o',typ='R',val_min=0,fr="pression de référence")
1050                                                 ),
1051                                   TURBULENCE = BLOC(condition="ITURB >= 0",
1052                                                     XKAPPA = SIMP(statut ='o',typ='R',defaut=0.42,val_min=0,fr="constante de Karman"),
1053                                                     CSTLOG = SIMP(statut ='f',typ='R',defaut=5.2,val_min=0,fr="constante de la loi log"),
1054                                                     YPLULI = SIMP(statut ='f',typ='R',defaut=2/0.42,val_min=0,fr="valeur limite de y+ pour la sous couche visqueuse"),
1055                                                     CMU = SIMP(statut ='f',typ='R',defaut=0.009,val_min=0,fr="constante C mu"),
1056                                                     CE1 = SIMP(statut ='f',typ='R',defaut=1.44,val_min=0,fr="constante C epsilon 1"),
1057                                                     CE2 = SIMP(statut ='f',typ='R',defaut=1.92,val_min=0,fr="constante C epsilon 2"),
1058                                                     CE3 = SIMP(statut ='f',typ='R',defaut=1.0,val_min=0,fr="constante C epsilon 3"),
1059                                                     SIGMAK = SIMP(statut ='f',typ='R',defaut=1.0,val_min=0,fr="nombre de Prandtl pour k en k-epsilon"),
1060                                                     SIGMAE = SIMP(statut ='f',typ='R',defaut=1.3,val_min=0,fr="nombre de Prandtl pour epsilon en k-epsilon"),
1061                                                     ALMAX = SIMP(statut ='f',typ='R',val_min=0,fr="longueur macroscopique caractéristique du domaine"),
1062                                                     UREF = SIMP(statut ='f',typ='R',val_min=0,fr="vitesse caractéristique de l'écoulement pour l'initialisation du k-epsilon")
1063                                                     )
1064                                   )
1065                                  );
1066