Minicode/cata_saturne2.py

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