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