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