Minicode/cata_saturne.py

   1 # debut entete
   2
   3 import Accas
   4 from Accas import *
   5
   6 import ops
   7
   8 JdC = JDC_CATA(code='SATURNE',
   9                execmodul=None,
  10                regles = (AU_MOINS_UN('DEBUT','POURSUITE'),
  11                          AU_MOINS_UN('FIN'),
  12                          A_CLASSER(('DEBUT','POURSUITE'),'FIN')
  13                         )
  14               );
  15
  16 # P. RASCLE MMN
  17 # remarques diverses sur le catalogue Saturne
  18 # - dans les blocs, il faut au moins un mot clé de statut obligatoire
  19 # probleme de rafraichissement des blocs dépendants quand la valeur d'un mot cle global (ITURB) passe de 1 à 0
  20
  21 class sonde(ASSD):pass
  22 class varsca(ASSD):pass
  23 class flusca(ASSD):pass
  24 class varpre(ASSD):pass
  25 class varvitx(ASSD):pass
  26 class varvity(ASSD):pass
  27 class varvitz(ASSD):pass
  28 class eturb(ASSD):pass
  29 class dturb(ASSD):pass
  30 class tsr11(ASSD):pass
  31 class tsr22(ASSD):pass
  32 class tsr33(ASSD):pass
  33 class tsr12(ASSD):pass
  34 class tsr13(ASSD):pass
  35 class tsr23(ASSD):pass
  36 class resti(ASSD):pass
  37
  38 class maillage(ASSD):pass
  39 class listr8          (ASSD):pass
  40
  41 # fin entete
  42
  43 INCLUDE=MACRO(nom="INCLUDE",op=-1,docu="U4.13.01-e",
  44              fr="Débranchement vers un fichier de commandes secondaires",
  45              sd_prod=ops.INCLUDE,op_init=ops.INCLUDE_context,fichier_ini=1,
  46          UNITE = SIMP(statut='o',typ='I'),
  47          INFO  = SIMP(statut='f',typ='I',defaut=1,into=(1,2)),
  48 );
  49
  50 POURSUITE=MACRO(nom="POURSUITE",op=0,repetable='n',fr="Poursuite d une étude",
  51                 docu="U4.11.03-f1",sd_prod = ops.POURSUITE,
  52                 op_init = ops.POURSUITE_context,fichier_ini = 1,
  53          PAR_LOT         =SIMP(fr="mode de traitement des commandes",statut='f',typ='TXM',
  54                            into=("OUI","NON"),defaut="OUI"),
  55          BASE            =FACT(fr="définition des parmètres associés aux bases JEVEUX",
  56                                statut='f',min=1,max=3,
  57            FICHIER         =SIMP(fr="nom de la base",statut='o',typ='TXM'),
  58            TITRE           =SIMP(statut='f',typ='TXM'),
  59            CAS             =SIMP(statut='f',typ='TXM'),
  60            NMAX_ENRE       =SIMP(fr="nombre maximum d enregistrements",statut='f',typ='I'),
  61            LONG_ENRE       =SIMP(fr="longueur des enregistrements",statut='f',typ='I'),
  62            LONG_REPE       =SIMP(fr="longueur du répertoire",statut='f',typ='I'),
  63          ),
  64          IMPRESSION      =FACT(statut='f',min=1,max=3,
  65            FICHIER         =SIMP(statut='o',typ='TXM'),
  66            UNITE           =SIMP(statut='o',typ='I'),
  67          ),
  68          CATALOGUE       =FACT(statut='f',min=1,max=10,
  69            FICHIER         =SIMP(statut='o',typ='TXM'),
  70            TITRE           =SIMP(statut='f',typ='TXM'),
  71            UNITE           =SIMP(statut='f',typ='I'),
  72          ),
  73          DEBUG           =FACT(fr="option de déboggage reservée aux développeurs",
  74                                statut='f',min=1,max=1,
  75            JXVERI          =SIMP(fr="vérifie l intégrité de la segmentation mémoire",
  76                                  statut='f',typ='TXM',into=('OUI','NON'),defaut='NON'),
  77            JEVEUX          =SIMP(fr="force les déchargement sur disque",
  78                                  statut='f',typ='TXM',into=('OUI','NON'),defaut='NON'),
  79            ENVIMA          =SIMP(fr="imprime les valeurs définies dans ENVIMA",
  80                                  statut='f',typ='TXM',into=('TES',)),
  81          ),
  82          MEMOIRE         =FACT(fr="mode de gestion mémoire utilisé",statut='f',min=1,max=1,
  83            GESTION         =SIMP(statut='f',typ='TXM',into=('COMPACTE','RAPIDE'),defaut='RAPIDE'),
  84            TYPE_ALLOCATION =SIMP(statut='f',typ='I',into=(1,2,3,4),defaut=1),
  85            TAILLE          =SIMP(statut='f',typ='I'),
  86            TAILLE_BLOC     =SIMP(statut='f',typ='R',defaut=800.),
  87            PARTITION       =SIMP(statut='f',typ='R' ),
  88          ),
  89          CODE            =FACT("définition d un nom pour l'esemble d'une étude",
  90                                statut='f',min=1,max=1,
  91            NOM             =SIMP(statut='o',typ='TXM'),
  92            UNITE           =SIMP(statut='f',typ='I',defaut=15),
  93          ),
  94 )  ;
  95
  96 FORMULE = FORM( nom='FORMULE',op=-5,sd_prod=fonction,
  97                 fr="Définition d une fonction",reentrant = 'n',
  98                 regles=(UN_PARMI('REEL','ENTIER','COMPLEXE'),),
  99                 REEL = SIMP(typ = 'shell',max=1),
 100                 ENTIER = SIMP(typ = 'shell',max=1),
 101                 COMPLEXE = SIMP(typ = 'shell',max=1),
 102 ) ;
 103
 104 DEFI_SONDE = OPER(nom="DEFI_SONDE",op= 1,sd_prod=sonde,
 105      docu="U2D1",
 106      fr="définition d'une sonde historique avec ses coordonnées",
 107                   X = SIMP(statut ='o',typ='R',
 108                            fr="coordonnée X de la sonde"),
 109                   Y = SIMP(statut ='o',typ='R',
 110                            fr="coordonnée Y de la sonde"),
 111                   Z = SIMP(statut ='o',typ='R',
 112                            fr="coordonnée Z de la sonde")
 113                  );
 114
 115 def defi_scala_prod(**args):
 116    return varsca
 117
 118 DEFI_SCALA = OPER(nom="DEFI_SCALA",op=2,sd_prod=defi_scala_prod,
 119 fr="définition d'une inconnue scalaire avec ses paramètres physico numériques",
 120                   NUMERIC = FACT(statut='o',max=01,
 121                                 fr="propriétés numériques liées à l'inconnue",
 122                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 123                                                   fr="indicateur de convection"),
 124                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 125                                                   fr="présence terme instationnaire dans les matrices"),
 126                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 127                                                   fr="diffusion"),
 128                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 129                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 130                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 131                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 132                                                                           fr="pourcentage schéma convectif au second ordre"),
 133                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 134                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 135                                                             ),
 136                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 137                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 138                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 139                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 140                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 141                                                   fr="type de limitation des gradients"),
 142                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 143                                                   fr="reconstruction des gradients"),
 144                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 145                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 146                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 147                                                   fr="facteur de limitation des gradients"),
 148                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 149                                                   fr="précision relative pour la reconstruction itérative des gradients")
 150                                  ),
 151                   RESTITUE = FACT(statut='o',max=01,
 152                                   fr = "Description de la grandeur restituée",
 153                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 154                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 155                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 156                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 157                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 158                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 159                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 160                                   ),
 161                   MODELE = FACT(statut='o',max=01,
 162                                 fr="modélisation",
 163                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 164                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 165                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 166                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 167                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 168                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 169                                 )
 170                  );
 171
 172 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",
 173                   VARFL  = SIMP(statut='o',typ=varsca,fr="scalaire associé au calcul de la variance des fluctuations"),
 174                   NUMERIC = FACT(statut='o',max=01,
 175                                 fr="propriétés numériques liées à l'inconnue",
 176                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 177                                                   fr="indicateur de convection"),
 178                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 179                                                   fr="présence terme instationnaire dans les matrices"),
 180                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 181                                                   fr="diffusion"),
 182                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 183                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 184                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 185                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 186                                                                           fr="pourcentage schéma convectif au second ordre"),
 187                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 188                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 189                                                             ),
 190                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 191                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 192                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 193                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 194                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 195                                                   fr="type de limitation des gradients"),
 196                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 197                                                   fr="reconstruction des gradients"),
 198                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 199                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 200                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 201                                                   fr="facteur de limitation des gradients"),
 202                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 203                                                   fr="précision relative pour la reconstruction itérative des gradients")
 204                                  ),
 205                   RESTITUE = FACT(statut='o',max=01,
 206                                   fr = "Description de la grandeur restituée",
 207                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 208                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 209                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 210                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 211                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 212                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 213                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 214                                   ),
 215                   MODELE = FACT(statut='o',max=01,
 216                                 fr="modélisation",
 217                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 218                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 219                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 220                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 221                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 222                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable"),
 223                                 RVARFL = SIMP(statut='o',typ='R',val_min=0,defaut=0.8,
 224                                                        fr="nombre de Prandtl pour la variance des fluctuations du scalaire")
 225                                 )
 226                   );
 227
 228 DEFI_PRESSION = OPER(nom="DEFI_PRESSION",op=2,sd_prod=varpre,fr="définition del'inconnue pression avec ses paramètres physico numériques",
 229                    NUMERIC = FACT(statut='o',max=01,
 230                                   fr="propriétés de l'équation de pression",
 231                                   ICONV  = SIMP(statut='o',typ='I',defaut=0,into=(0,1),
 232                                                 fr="indicateur de convection"),
 233                                   ISTAT  = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 234                                                 fr="présence terme instationnaire dans les matrices"),
 235                                   IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 236                                                 fr="diffusion"),
 237                                   IDIRCL = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 238                                                 fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 239                                   SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 240                                                           BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 241                                                                         fr="pourcentage schéma convectif au second ordre"),
 242                                                           ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 243                                                                         fr ="type schéma convectif au second ordre 1 : centré")
 244                                                           ),
 245                                   NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 246                                                 fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 247                                   EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 248                                                 fr="précision relative pour la résolution des systèmes linéaires"),
 249                                   IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 250                                                 fr="type de limitation des gradients"),
 251                                   NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 252                                                 fr="reconstruction des gradients"),
 253                                   NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 254                                                 fr="nombre d'itérations pour la reconstruction des seconds membres"),
 255                                   CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 256                                                 fr="facteur de limitation des gradients"),
 257                                   EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 258                                                 fr="précision relative pour la reconstruction itérative des gradients")
 259                                   ),
 260                   RESTITUE = FACT(statut='o',max=01,
 261                                   fr = "Description de la grandeur restituée",
 262                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée",
 263                                                 defaut="Pression"),
 264                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 265                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 266                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 267                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 268                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 269                                                 fr = "inventaire des sondes utilisées pour la grandeur")
 270                                   ),
 271                   MODELE = FACT(statut='o',max=01,
 272                                 fr="modélisation",
 273                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 274                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour la grandeur"),
 275                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour la grandeur"),
 276                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour la grandeur"),
 277                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée à la grandeur"),
 278                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 279                                 )
 280                     );
 281
 282 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",
 283                    NUMERIC = FACT(statut='o',max=01,
 284                                   fr="propriétés de l'équation de vitesse X",
 285                                   ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 286                                                 fr="indicateur de convection"),
 287                                   ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 288                                                 fr="présence terme instationnaire dans les matrices"),
 289                                   IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 290                                                 fr="diffusion"),
 291                                   IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 292                                                 fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 293                                   SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 294                                                           BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 295                                                                         fr="pourcentage schéma convectif au second ordre"),
 296                                                           ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 297                                                                         fr ="type schéma convectif au second ordre 1 : centré")
 298                                                           ),
 299                                   NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 300                                                 fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 301                                   EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 302                                                 fr="précision relative pour la résolution des systèmes linéaires"),
 303                                   IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 304                                                 fr="type de limitation des gradients"),
 305                                   NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 306                                                 fr="reconstruction des gradients"),
 307                                   NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 308                                                 fr="nombre d'itérations pour la reconstruction des seconds membres"),
 309                                   CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 310                                                 fr="facteur de limitation des gradients"),
 311                                   EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 312                                                 fr="précision relative pour la reconstruction itérative des gradients")
 313                                   ),
 314                   RESTITUE = FACT(statut='o',max=01,
 315                                   fr = "Description de la grandeur restituée",
 316                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée",
 317                                                 defaut="vitesse_u1"),
 318                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 319                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 320                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 321                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 322                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 323                                                 fr = "inventaire des sondes utilisées pour la grandeur")
 324                                   ),
 325                   MODELE = FACT(statut='o',max=01,
 326                                 fr="modélisation",
 327                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 328                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour la grandeur"),
 329                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour la grandeur"),
 330                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour la grandeur"),
 331                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée à la grandeur"),
 332                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 333                                 )
 334                 );
 335
 336 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",
 337                    NUMERIC = FACT(statut='o',max=01,
 338                                   fr="propriétés de l'équation de vitesse Y",
 339                                   ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 340                                                 fr="indicateur de convection"),
 341                                   ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 342                                                 fr="présence terme instationnaire dans les matrices"),
 343                                   IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 344                                                 fr="diffusion"),
 345                                   IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 346                                                 fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 347                                   SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 348                                                           BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 349                                                                         fr="pourcentage schéma convectif au second ordre"),
 350                                                           ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 351                                                                         fr ="type schéma convectif au second ordre 1 : centré")
 352                                                           ),
 353                                   NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 354                                                 fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 355                                   EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 356                                                 fr="précision relative pour la résolution des systèmes linéaires"),
 357                                   IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 358                                                 fr="type de limitation des gradients"),
 359                                   NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 360                                                 fr="reconstruction des gradients"),
 361                                   NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 362                                                 fr="nombre d'itérations pour la reconstruction des seconds membres"),
 363                                   CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 364                                                 fr="facteur de limitation des gradients"),
 365                                   EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 366                                                 fr="précision relative pour la reconstruction itérative des gradients")
 367                                   ),
 368                   RESTITUE = FACT(statut='o',max=01,
 369                                   fr = "Description de la grandeur restituée",
 370                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée",
 371                                                 defaut="vitesse_v1"),
 372                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 373                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 374                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 375                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 376                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 377                                                 fr = "inventaire des sondes utilisées pour la grandeur")
 378                                   ),
 379                   MODELE = FACT(statut='o',max=01,
 380                                 fr="modélisation",
 381                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 382                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour la grandeur"),
 383                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour la grandeur"),
 384                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour la grandeur"),
 385                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée à la grandeur"),
 386                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 387                                 )
 388                 );
 389
 390 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",
 391                    NUMERIC = FACT(statut='o',max=01,
 392                                   fr="propriétés de l'équation de vitesse Z",
 393                                   ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 394                                                 fr="indicateur de convection"),
 395                                   ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 396                                                 fr="présence terme instationnaire dans les matrices"),
 397                                   IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 398                                                 fr="diffusion"),
 399                                   IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 400                                                 fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 401                                   SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 402                                                           BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 403                                                                         fr="pourcentage schéma convectif au second ordre"),
 404                                                           ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 405                                                                         fr ="type schéma convectif au second ordre 1 : centré")
 406                                                           ),
 407                                   NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 408                                                 fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 409                                   EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 410                                                 fr="précision relative pour la résolution des systèmes linéaires"),
 411                                   IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 412                                                 fr="type de limitation des gradients"),
 413                                   NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 414                                                 fr="reconstruction des gradients"),
 415                                   NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 416                                                 fr="nombre d'itérations pour la reconstruction des seconds membres"),
 417                                   CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 418                                                 fr="facteur de limitation des gradients"),
 419                                   EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 420                                                 fr="précision relative pour la reconstruction itérative des gradients")
 421                                   ),
 422                   RESTITUE = FACT(statut='o',max=01,
 423                                   fr = "Description de la grandeur restituée",
 424                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée",
 425                                                 defaut="vitesse_w1"),
 426                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 427                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 428                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 429                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 430                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 431                                                 fr = "inventaire des sondes utilisées pour la grandeur")
 432                                   ),
 433                   MODELE = FACT(statut='o',max=01,
 434                                 fr="modélisation",
 435                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 436                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour la grandeur"),
 437                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour la grandeur"),
 438                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour la grandeur"),
 439                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée à la grandeur"),
 440                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable"),
 441                                 )
 442                 );
 443
 444
 445 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",
 446                   NUMERIC = FACT(statut='o',max=01,
 447                                     fr="propriétés numériques liées à l'inconnue",
 448                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 449                                                   fr="indicateur de convection"),
 450                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 451                                                   fr="présence terme instationnaire dans les matrices"),
 452                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 453                                                   fr="diffusion"),
 454                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 455                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 456                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 457                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 458                                                                           fr="pourcentage schéma convectif au second ordre"),
 459                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 460                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 461                                                             ),
 462                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 463                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 464                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 465                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 466                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 467                                                   fr="type de limitation des gradients"),
 468                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 469                                                   fr="reconstruction des gradients"),
 470                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 471                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 472                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 473                                                   fr="facteur de limitation des gradients"),
 474                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 475                                                   fr="précision relative pour la reconstruction itérative des gradients")
 476                                  ),
 477                   RESTITUE = FACT(statut='o',max=01,
 478                                   fr = "Description de la grandeur restituée",
 479                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 480                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 481                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 482                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 483                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 484                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 485                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 486                                   ),
 487                   MODELE = FACT(statut='o',max=01,
 488                                 fr="modélisation",
 489                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 490                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 491                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 492                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 493                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 494                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 495                                 )
 496                  );
 497
 498 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",
 499                   NUMERIC = FACT(statut='o',max=01,
 500                                 fr="propriétés numériques liées à l'inconnue",
 501                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 502                                                   fr="indicateur de convection"),
 503                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 504                                                   fr="présence terme instationnaire dans les matrices"),
 505                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 506                                                   fr="diffusion"),
 507                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 508                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 509                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 510                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 511                                                                           fr="pourcentage schéma convectif au second ordre"),
 512                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 513                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 514                                                             ),
 515                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 516                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 517                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 518                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 519                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 520                                                   fr="type de limitation des gradients"),
 521                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 522                                                   fr="reconstruction des gradients"),
 523                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 524                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 525                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 526                                                   fr="facteur de limitation des gradients"),
 527                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 528                                                   fr="précision relative pour la reconstruction itérative des gradients")
 529                                  ),
 530                   RESTITUE = FACT(statut='o',max=01,
 531                                   fr = "Description de la grandeur restituée",
 532                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 533                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 534                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 535                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 536                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 537                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 538                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 539                                   ),
 540                   MODELE = FACT(statut='o',max=01,
 541                                 fr="modélisation",
 542                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 543                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 544                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 545                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 546                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 547                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 548                                 )
 549                  );
 550
 551 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",
 552                   NUMERIC = FACT(statut='o',max=01,
 553                                 fr="propriétés numériques liées à l'inconnue",
 554                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 555                                                   fr="indicateur de convection"),
 556                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 557                                                   fr="présence terme instationnaire dans les matrices"),
 558                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 559                                                   fr="diffusion"),
 560                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 561                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 562                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 563                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 564                                                                           fr="pourcentage schéma convectif au second ordre"),
 565                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 566                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 567                                                             ),
 568                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 569                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 570                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 571                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 572                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 573                                                   fr="type de limitation des gradients"),
 574                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 575                                                   fr="reconstruction des gradients"),
 576                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 577                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 578                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 579                                                   fr="facteur de limitation des gradients"),
 580                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 581                                                   fr="précision relative pour la reconstruction itérative des gradients")
 582                                  ),
 583                   RESTITUE = FACT(statut='o',max=01,
 584                                   fr = "Description de la grandeur restituée",
 585                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 586                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 587                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 588                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 589                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 590                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 591                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 592                                   ),
 593                   MODELE = FACT(statut='o',max=01,
 594                                 fr="modélisation",
 595                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 596                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 597                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 598                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 599                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 600                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 601                                 )
 602                  );
 603
 604 DEFI_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",
 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_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",
 658                   NUMERIC = FACT(statut='o',max=01,
 659                                 fr="propriétés numériques liées à l'inconnue",
 660                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 661                                                   fr="indicateur de convection"),
 662                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 663                                                   fr="présence terme instationnaire dans les matrices"),
 664                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 665                                                   fr="diffusion"),
 666                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 667                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 668                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 669                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 670                                                                           fr="pourcentage schéma convectif au second ordre"),
 671                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 672                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 673                                                             ),
 674                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 675                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 676                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 677                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 678                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 679                                                   fr="type de limitation des gradients"),
 680                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 681                                                   fr="reconstruction des gradients"),
 682                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 683                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 684                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 685                                                   fr="facteur de limitation des gradients"),
 686                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 687                                                   fr="précision relative pour la reconstruction itérative des gradients")
 688                                  ),
 689                   RESTITUE = FACT(statut='o',max=01,
 690                                   fr = "Description de la grandeur restituée",
 691                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 692                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 693                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 694                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 695                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 696                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 697                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 698                                   ),
 699                   MODELE = FACT(statut='o',max=01,
 700                                 fr="modélisation",
 701                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 702                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 703                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 704                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 705                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 706                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 707                                 )
 708                  );
 709
 710 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",
 711                   NUMERIC = FACT(statut='o',max=01,
 712                                 fr="propriétés numériques liées à l'inconnue",
 713                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 714                                                   fr="indicateur de convection"),
 715                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 716                                                   fr="présence terme instationnaire dans les matrices"),
 717                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 718                                                   fr="diffusion"),
 719                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 720                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 721                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 722                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 723                                                                           fr="pourcentage schéma convectif au second ordre"),
 724                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 725                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 726                                                             ),
 727                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 728                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 729                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 730                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 731                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 732                                                   fr="type de limitation des gradients"),
 733                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 734                                                   fr="reconstruction des gradients"),
 735                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 736                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 737                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 738                                                   fr="facteur de limitation des gradients"),
 739                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 740                                                   fr="précision relative pour la reconstruction itérative des gradients")
 741                                  ),
 742                   RESTITUE = FACT(statut='o',max=01,
 743                                   fr = "Description de la grandeur restituée",
 744                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 745                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 746                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 747                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 748                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 749                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 750                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 751                                   ),
 752                   MODELE = FACT(statut='o',max=01,
 753                                 fr="modélisation",
 754                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 755                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 756                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 757                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 758                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 759                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 760                                 )
 761                  );
 762
 763 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",
 764                   NUMERIC = FACT(statut='o',max=01,
 765                                 fr="propriétés numériques liées à l'inconnue",
 766                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 767                                                   fr="indicateur de convection"),
 768                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 769                                                   fr="présence terme instationnaire dans les matrices"),
 770                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 771                                                   fr="diffusion"),
 772                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 773                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 774                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 775                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 776                                                                           fr="pourcentage schéma convectif au second ordre"),
 777                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 778                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 779                                                             ),
 780                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 781                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 782                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 783                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 784                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 785                                                   fr="type de limitation des gradients"),
 786                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 787                                                   fr="reconstruction des gradients"),
 788                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 789                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 790                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 791                                                   fr="facteur de limitation des gradients"),
 792                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 793                                                   fr="précision relative pour la reconstruction itérative des gradients")
 794                                  ),
 795                   RESTITUE = FACT(statut='o',max=01,
 796                                   fr = "Description de la grandeur restituée",
 797                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 798                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 799                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 800                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 801                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 802                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 803                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 804                                   ),
 805                   MODELE = FACT(statut='o',max=01,
 806                                 fr="modélisation",
 807                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 808                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 809                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 810                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 811                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 812                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 813                                 )
 814                  );
 815
 816 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",
 817                   NUMERIC = FACT(statut='o',max=01,
 818                                 fr="propriétés numériques liées à l'inconnue",
 819                                     ICONV  = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 820                                                   fr="indicateur de convection"),
 821                                     ISTAT  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 822                                                   fr="présence terme instationnaire dans les matrices"),
 823                                     IDIFF  = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 824                                                   fr="diffusion"),
 825                                     IDIRCL = SIMP(statut='f',typ='I',defaut=1,into=(0,1),
 826                                                   fr="décalage du spectre des valeurs propres en l'absence de Dirichlet"),
 827                                     SCHEMA_CONVECTIF = BLOC(condition="ICONV == 1", fr ="ordre du schéma convectif",
 828                                                             BLENCV = SIMP(statut='o',typ='R',defaut=0,val_min=0,val_max=1,
 829                                                                           fr="pourcentage schéma convectif au second ordre"),
 830                                                             ISCHCV = SIMP(statut='f',typ='I',defaut=0,into=(0,1),
 831                                                                           fr ="type schéma convectif au second ordre 1 : centré")
 832                                                             ),
 833                                     NITMAX = SIMP(statut='f',typ='I',defaut=10000,val_min=0,
 834                                                   fr="nombre max d'itération pour la résolution des systèmes linéaires"),
 835                                     EPSILO = SIMP(statut='f',typ='R',defaut=1.e-4,val_min=0,
 836                                                   fr="précision relative pour la résolution des systèmes linéaires"),
 837                                     IMLIGR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,-1,0,1),
 838                                                   fr="type de limitation des gradients"),
 839                                     NSWRGR = SIMP(statut='f',typ='I',defaut=100,val_min=0,
 840                                                   fr="reconstruction des gradients"),
 841                                     NSWRSM = SIMP(statut='f',typ='I',defaut=2,val_min=0,
 842                                                   fr="nombre d'itérations pour la reconstruction des seconds membres"),
 843                                     CLIMGR = SIMP(statut='f',typ='R',defaut=1.5,val_min=0,
 844                                                   fr="facteur de limitation des gradients"),
 845                                     EPSRGR = SIMP(statut='f',typ='R',defaut=1.e-5,val_min=0,
 846                                                   fr="précision relative pour la reconstruction itérative des gradients")
 847                                  ),
 848                   RESTITUE = FACT(statut='o',max=01,
 849                                   fr = "Description de la grandeur restituée",
 850                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 851                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 852                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 853                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 854                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 855                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 856                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 857                                   ),
 858                   MODELE = FACT(statut='o',max=01,
 859                                 fr="modélisation",
 860                                 VALREF = SIMP(statut='o',typ='R',fr="valeur de référence (température, masse volumique...)"),
 861                                 SCAMIN = SIMP(statut='f',typ='R',defaut=1.e12,fr="valeur minimale pour le scalaire"),
 862                                 SCAMAX = SIMP(statut='f',typ='R',defaut=-1.e12,fr="valeur maximale pour le scalaire"),
 863                                 SIGMAS = SIMP(statut='f',typ='R',val_min=0,defaut=1,fr="nombre de Prandtl pour le scalaire"),
 864                                 VISLS0 = SIMP(statut='o',typ='R',fr="valeur de référence de la viscosité associée au scalaire"),
 865                                 IVISLS = SIMP(statut='f',typ='I',defaut=0,into=(0,1),fr="viscosité variable")
 866                                 )
 867                  );
 868
 869
 870 DEFI_RESTI = OPER(nom="DEFI_RESTI",op=2,sd_prod=resti,fr = "Description de la grandeur restituée",
 871                                   NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur restituée"),
 872                                   ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 873                                                 fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
 874                                   ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
 875                                                 fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
 876                                   IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
 877                                                 fr = "inventaire des sondes utilisées pour la grandeur scalaire")
 878                   );
 879
 880 CALCUL_SATURNE = PROC(nom = "CALCUL_SATURNE",op = 1,
 881                       fr = "définition des paramètres généraux pour un calcul Saturne",
 882                       docu = "néant",reentrant = 'n',
 883         NomsBibliotheque   = NUPL ( max      = '**',
 884                                     statut   = 'o',
 885                                     elements = (  SIMP (typ='TXM',fr="Identificateur Procedure Bibliotheque"),
 886                                                   SIMP (typ='TXM',fr="Identifiant de l'isotope dans la bibliotheque")
 887                                                )
 888                                   ),
 889         ENVELOPPE = FACT(statut='o',max=01,fr = "gestion de l'enveloppe",
 890                          IFOENV = SIMP(statut='o',typ='I',defaut=2,into=(0,1,2),position='global',
 891                                        fr = "mode de communication enveloppe solveur"),
 892                          ENVTOSOLV = BLOC(condition="IFOENV > 0",
 893                                           fr = "liaison enveloppe vers solveur",
 894                                           IMPEVI = SIMP(statut='o',typ='I',defaut=13,val_min=0,val_max=99,
 895                                                         fr = "numéro de fichier enveloppe vers solveur"),
 896                                           FICEVI = SIMP(statut='o',typ='TXM',defaut="enveloppe_vers_solveur          ",
 897                                                         fr = "nom de fichier enveloppe vers solveur")
 898                                          ),
 899                          SOLVTOENV = BLOC(condition="IFOENV > 0",
 900                                           fr = "liaison solveur vers enveloppe",
 901                                           IMPEVO = SIMP(statut='o',typ='I',defaut=14,val_min=0,val_max=99,
 902                                                         fr = "numéro de fichier solveur vers enveloppe"),
 903                                           FICEVO = SIMP(statut='o',typ='TXM',defaut="solveur_vers_enveloppe          ",
 904                                                         fr = "nom de fichier solveur vers enveloppe")
 905                                          )
 906                         ),
 907         FICHIERS_CALCUL = FACT(statut='f',max=01,
 908                                fr ="définition des fichiers géométrie, suite, stop",
 909                                GEOMETRIE = BLOC(condition="IFOENV == 0",
 910                                                 fr = "fichier géométrique (pas d'enveloppe)",
 911                                                 IMPGEO = SIMP(statut='o',typ='I',defaut=10,val_min=0,val_max=99,
 912                                                               fr = "numéro de fichier géométrique"),
 913                                                 FICGEO = SIMP(statut='o',typ='TXM',defaut="geomet",
 914                                                               fr = "nom de fichier géométrique")
 915                                                ),
 916                                SUITE_AMONT = BLOC(condition="ISUITE == 1",
 917                                                   fr = "fichier suite amont",
 918                                                   IMPAMO = SIMP(statut='o',typ='I',defaut=11,val_min=0,val_max=99,
 919                                                               fr = "numéro de fichier suite amont"),
 920                                                   FICGEO = SIMP(statut='o',typ='TXM',defaut="suiamo",
 921                                                               fr = "nom de fichier suite amont"),
 922                                                   IFOAMO = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 923                                                               fr = "format du fichier suite amont, 0 = binaire")
 924                                                  ),
 925                                FICHIER_STOP = FACT(fr = "fichier stop",
 926                                                    IMPSTP = SIMP(statut='o',typ='I',defaut=12,val_min=0,val_max=99,
 927                                                                  fr = "numéro de fichier stop"),
 928                                                    FICSTP = SIMP(statut='o',typ='TXM',defaut="ficstp",
 929                                                                  fr = "nom de fichier stop")
 930                                                   ),
 931                                SUITE_AVAL = FACT(fr = "fichier suite aval",
 932                                                  IMPAVA = SIMP(statut='o',typ='I',defaut=20,val_min=0,val_max=99,
 933                                                              fr = "numéro de fichier suite aval"),
 934                                                  FICAVA = SIMP(statut='o',typ='TXM',defaut="suiamo",
 935                                                              fr = "nom de fichier suite aval"),
 936                                                  IFOAVA = SIMP(statut='o',typ='I',defaut=1,into=(0,1),
 937                                                              fr = "format du fichier suite aval, 0 = binaire")
 938                                                 )
 939                               ),
 940         POST_PROC_ENSIGHT = FACT(statut='f',max=01,
 941                                  fr = "options de post processing par Ensight",
 942                                  IFOENS = SIMP(statut='f',typ='I',defaut=-2,into=(-2,-1,0,1),
 943                                                fr = "option type de fichier et processus générateur"),
 944                                  NTCHR  = SIMP(statut='f',typ='I',defaut=-1,
 945                                                fr = "périodicité de sortie des fichiers Ensight"),
 946                                  ITCHR  = SIMP(statut='f',typ='I',defaut=0,val_min=0,
 947                                                fr = "compteur des sorties des fichiers Ensight"),
 948                                  N3S_ASCII_NOY = BLOC(condition="(IFOENS == -1) and (IFOENV == 0)",
 949                                                       fr = "format Ensight N3S ASCII généré par le noyau",
 950                                                       IMPPST = SIMP(statut='o',typ='I',defaut=21,val_min=0,val_max=99,
 951                                                                     fr = "numéro de fichier Ensight"),
 952                                                       FICPST = SIMP(statut='o',typ='TXM',defaut="dessin",
 953                                                                     fr = "nom de fichier Ensight")
 954                                                       ),
 955                                  P0_NOY =        BLOC(condition="((IFOENS == 0) or (IFOENS == 1)) and (IFOENV == 0)",
 956                                                       fr = "format Ensight P0 généré par le noyau",
 957                                                       IMPEP0 = SIMP(statut='o',typ='I',defaut=22,val_min=0,val_max=99,
 958                                                                     fr = "numéro de fichier Ensight"),
 959                                                       EMPCHR = SIMP(statut='o',typ='TXM',defaut="./",
 960                                                                     fr = "répertoire de fichier Ensight"),
 961                                                       ENTCHR = SIMP(statut='o',typ='TXM',defaut="chr",
 962                                                                     fr = "préfixe nom de fichier Ensight")
 963                                                       )
 964                                 ),
 965         HISTORIQUE_PONCTUEL = FACT(statut='o',max=01,
 966                                    fr = "Sondes historiques",
 967                                    FICHIERS_HISTORIQUES = FACT(statut='f',max=01,
 968                                                                fr = "description des fichiers historiques",
 969                                                                EMPHIS = SIMP(statut='o',typ='TXM',defaut="./",
 970                                                                              fr="répertoire fichiers historiques"),
 971                                                                EXTHIS = SIMP(statut='o',typ='TXM',defaut="hst",
 972                                                                              fr="extension fichiers historiques")
 973                                                                ),
 974                                    NTHIST = SIMP(statut='f',typ='I',defaut=-999,
 975                                                  fr="fréquence de sortie des historiques en pas de temps"),
 976                                    NTHSAV = SIMP(statut='f',typ='I',defaut=-999,
 977                                                  fr="fréquence de sauvegarde des historiques en pas de temps")
 978                                   ),
 979         OPTIONS_TURBULENCE = FACT(statut='o',max=01,
 980                                   fr="modèle de turbulence",
 981                                   ITURB  = SIMP(statut='f',fr="laminaire : 0, k-epsilon :1,Rij-epsilon :2",
 982                                                 typ='I',into=(0,1,2),defaut=1,position='global'),
 983                                   MODTURB = BLOC(condition="ITURB == 1",
 984                                                  fr = "option k-epsilon",
 985                                                  IGRAKE = SIMP(statut='o',typ='I',
 986                                                                fr="prise en compte gravité dans k-epsilon",
 987                                                                into=(0,1),defaut=0),
 988                                                  IDEUCH = SIMP(statut='f',typ='I',
 989                                                                fr="prise en compte k-epsilon deux échelles",
 990                                                                into=(0,1),defaut=1),
 991                                                  IKEKOU = SIMP(statut='f',typ='I',
 992                                                                fr="prise en compte couplage en incréments sur k-epsilon",
 993                                                                into=(0,1),defaut=1)
 994                                                  ),
 995                                   TEMPTURB = BLOC(condition="ITURB == 1",
 996                                                   fr = "option k-epsilon",
 997                                                   ISCALT = SIMP(statut='o',typ=varsca,
 998                                                                 fr=" identificateur inconnue scalaire température ")
 999                                                   )
1000                                   ),
1001         MARCHE_TEMPS = FACT(statut='o',max=01,
1002                             fr = "définition de la marche en temps",
1003                             DTREF  = SIMP(statut='o',fr="pas de temps de référence",
1004                                                         typ='R',val_min=0),
1005                             IDTVAR = SIMP(statut='f',fr="pas de temps constant : 0, variable temps espace : 1, variable temps : 2",
1006                                           typ='I',into=(0,1,2),defaut=0,position='global'),
1007                                   # probleme trace eficas quand un mot cle position global change
1008                             PASVAR = BLOC(condition="IDTVAR != 0",fr="options pas de temps variable",
1009                                           XCFMAX = SIMP(statut='o',fr="nombre de Courant-Fourier cible",
1010                                                         typ='R',defaut=0.5,val_min=0),
1011                                           FMIN =   SIMP(statut='f',fr="rapport min pas calculé DTREF",
1012                                                         typ='R',defaut=0.1,val_min=0),
1013                                           FMAX =   SIMP(statut='f',fr="rapport max pas calculé DTREF",
1014                                                         typ='R',defaut=1000,val_min=0),
1015                                           VARRDT = SIMP(statut='f',fr="variation relative max pas calculé entre deux instants",
1016                                                         typ='R',defaut=0.1,val_min=0)
1017                                           )
1018                             ),
1019         OPTIONS_EQUATIONS = FACT(statut='o',max=01,
1020                                    fr = "propriétés des équations, inconnues principales",
1021                                    IMGR   = SIMP(statut='f',fr="utilisation du multigrille pour la résolution des systèmes linéaires",
1022                                                  typ='I',into=(0,1),defaut=0),
1023                                    IMRGRA = SIMP(statut='f',fr="type de reconstruction des gradients 1 : moindres carrés",
1024                                                  typ='I',into=(0,1),defaut=0),
1025                                    ),
1026         VARIABLES = FACT(statut='o',max=01,
1027                          fr = "Restitution des grandeurs principales",
1028                          NTLIST = SIMP(statut='f',typ='I',defaut=1,val_min=1,val_max='**',
1029                                        fr = "fréquence de sortie (en pas de temps) dans le compte rendu d'éxécution"),
1030                          IWARNI = SIMP(statut='f',typ='I',defaut=2,val_min=0,val_max='**',
1031                                        fr = "niveau de détail des impressions dans le compte rendu d'éxécution"),
1032                          MASVOL1 = FACT(statut='o',max=01,
1033                                          fr = "Description de la grandeur restituée : masse volumique",
1034                                          NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur : masse volumique",
1035                                                        defaut="Masse_vol1"),
1036                                          ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1037                                                        fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
1038                                          ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1039                                                        fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
1040                                          IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
1041                                                        fr = "inventaire des sondes utilisées pour la grandeur masse volumique")
1042                                         ),
1043                         TURB_KE = BLOC(condition="ITURB == 1",
1044                                        E_TURB = FACT(statut='o',max=01,
1045                                                      fr = "Description de la grandeur restituée : énergie turbulente",
1046                                                      NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur : energie turbulente",
1047                                                                    defaut="Energie_1"),
1048                                                      ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1049                                                                    fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
1050                                                      ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1051                                                                    fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
1052                                                      IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
1053                                                                    fr = "inventaire des sondes utilisées pour la grandeur énergie turbulente")
1054                                                     ),
1055                                        D_TURB = FACT(statut='o',max=01,
1056                                                      fr = "Description de la grandeur restituée : dissipation turbulente",
1057                                                      NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur : dissipation turbulente",
1058                                                                    defaut="Dissipation"),
1059                                                      ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1060                                                                    fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
1061                                                      ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1062                                                                    fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
1063                                                      IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
1064                                                                    fr = "inventaire des sondes utilisées pour la grandeur dissipation turbulente")
1065                                                     ),
1066                                        V_TURB = FACT(statut='o',max=01,
1067                                                      fr = "Description de la grandeur restituée : viscosité turbulente",
1068                                                      NOMVAR = SIMP(statut='o',typ='TXM',fr = "nom de la grandeur : viscosité turbulente",
1069                                                                    defaut="Visc_turb1"),
1070                                                      ICHRVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1071                                                                    fr = "post-traitement (-999 : le code décide ; 0 : non ; 1 : oui)"),
1072                                                      ILISVR = SIMP(statut='f',typ='I',defaut=-999,into=(-999,0,1),
1073                                                                    fr = "impression dans CR exécution (-999 : le code décide ; 0 : non ; 1 : oui)"),
1074                                                      IHISVR = SIMP(statut='f',typ=sonde,min=0,max='**',
1075                                                                    fr = "inventaire des sondes utilisées pour la grandeur viscosité turbulente")
1076                                                     ),
1077                                       )
1078                         ),
1079         GESTION_CALCUL = FACT(statut='o',max=01,
1080                               fr = "calcul suite et numéros de pas de temps début et fin",
1081                               ISUITE = SIMP(statut='o',typ='I',defaut=0,into=(0,1),position='global',
1082                                             fr = "indicateur calcul suite (1 = suite)"),
1083                               NTPABS = SIMP(statut='f',typ='I',defaut=0,
1084                                             fr = "numéro dernier pas de temps calcul précédent (initialisation automatique)"),
1085                               NTMABS = SIMP(statut='o',typ='I',
1086                                             fr = "numéro dernier pas de temps visé (absolu)"),
1087                               TTPABS = SIMP(statut='f',typ='R',defaut=0,
1088                                             fr = "temps simulation au dernier pas de temps précédent (initialisation automatique)")
1089                              ),
1090         CONSTANTES_PHYSIQUES=FACT(statut='o',max=01,
1091                                   fr = "Grandeurs physiques et modélisation",
1092                                   GRAVITE = FACT(statut='o',max=01,
1093                                                  fr = "composantes de la gravité",
1094                                                  GX = SIMP(statut ='o',typ='R',defaut=0, fr = "gravité selon X"),
1095                                                  GY = SIMP(statut ='o',typ='R',defaut=0, fr = "gravité selon Y"),
1096                                                  GZ = SIMP(statut ='o',typ='R',defaut=0, fr = "gravité selon Z")
1097                                                  ),
1098                                   FLUIDE = FACT(statut='o',max=01,
1099                                                 fr = "propriétés du fluide",
1100                                                 RO0 = SIMP(statut ='o',typ='R',val_min=0,fr="masse volumique de référence"),
1101                                                 VISCL0 = SIMP(statut ='o',typ='R',val_min=0,fr="viscosité dynamique de référence"),
1102                                                 P0 = SIMP(statut ='o',typ='R',val_min=0,fr="pression de référence")
1103                                                 ),
1104                                   TURBULENCE = BLOC(condition="ITURB >= 0",
1105                                                     XKAPPA = SIMP(statut ='o',typ='R',defaut=0.42,val_min=0,fr="constante de Karman"),
1106                                                     CSTLOG = SIMP(statut ='f',typ='R',defaut=5.2,val_min=0,fr="constante de la loi log"),
1107                                                     YPLULI = SIMP(statut ='f',typ='R',defaut=2/0.42,val_min=0,fr="valeur limite de y+ pour la sous couche visqueuse"),
1108                                                     CMU = SIMP(statut ='f',typ='R',defaut=0.009,val_min=0,fr="constante C mu"),
1109                                                     CE1 = SIMP(statut ='f',typ='R',defaut=1.44,val_min=0,fr="constante C epsilon 1"),
1110                                                     CE2 = SIMP(statut ='f',typ='R',defaut=1.92,val_min=0,fr="constante C epsilon 2"),
1111                                                     CE3 = SIMP(statut ='f',typ='R',defaut=1.0,val_min=0,fr="constante C epsilon 3"),
1112                                                     SIGMAK = SIMP(statut ='f',typ='R',defaut=1.0,val_min=0,fr="nombre de Prandtl pour k en k-epsilon"),
1113                                                     SIGMAE = SIMP(statut ='f',typ='R',defaut=1.3,val_min=0,fr="nombre de Prandtl pour epsilon en k-epsilon"),
1114                                                     ALMAX = SIMP(statut ='f',typ='R',val_min=0,fr="longueur macroscopique caractéristique du domaine"),
1115                                                     UREF = SIMP(statut ='f',typ='R',val_min=0,fr="vitesse caractéristique de l'écoulement pour l'initialisation du k-epsilon")
1116                                                     )
1117                                   )
1118                                  );
1119
1120 CALC_FONC_INTERP=OPER(nom="CALC_FONC_INTERP",op= 134,sd_prod=fonction,
1121                       docu="U4.32.01-d1",reentrant='f',
1122          regles=(UN_PARMI('VALE_R','LIST_PARA'),),
1123          FONCTION        =SIMP(statut='o',typ=fonction ),
1124          NOM_RESU        =SIMP(statut='f',typ='TXM',defaut="TOUTRESU"),
1125          VALE_R          =SIMP(statut='f',typ='R',max='**'),
1126          LIST_PARA       =SIMP(statut='f',typ=listr8 ),
1127          INTERPOL        =SIMP(statut='f',typ='TXM',max=2,defaut="LIN",into=("NON","LIN","LOG","INT") ),
1128          PROL_DROITE     =SIMP(statut='f',typ='TXM',defaut="EXCLU",into=("EXCLU","CONSTANT","LINEAIRE","INTERPRE") ),
1129          PROL_GAUCHE     =SIMP(statut='f',typ='TXM',defaut="EXCLU",into=("EXCLU","CONSTANT","LINEAIRE","INTERPRE") ),
1130          TITRE           =SIMP(statut='f',typ='TXM',max='**'),
1131          INFO            =SIMP(statut='f',typ='I',defaut= 1,into=( 1 , 2 ) ),
1132 )  ;
1133
1134 LIRE_MAILLAGE=OPER(nom="LIRE_MAILLAGE",op=   1,sd_prod=maillage,
1135                    fr="Lecture d'un fichier de maillage",
1136                    ang="Readings of a mesh file",
1137                    docu="U4.21.01-f",reentrant='n',
1138 #
1139          UNITE           =SIMP(statut='f',typ='I',defaut= 20 ),
1140 #
1141          FORMAT          =SIMP(statut='f',typ='TXM',defaut="ASTER",into=("ASTER","MED"),
1142                             fr="Format du fichier : ASTER ou MED.",
1143                             ang="Format of the file : ASTER or MED.",),
1144 #
1145          ABSC_CURV       =FACT(statut='f',min=00,max=01,
1146                TOUT          =SIMP(statut='f',typ='TXM',defaut="NON",into=("OUI","NON") ),
1147          ),
1148 #
1149          VERI_MAIL       =FACT(statut='d',min=01,max=01,
1150                APLAT         =SIMP(statut='f',typ='R',defaut= 1.0E-3 ),
1151          ),
1152 #
1153          b_format_med =BLOC( condition = " ( FORMAT == 'MED' ) " ,
1154                            fr="Informations complémentaires pour la lecture MED.",
1155                            ang="Further information for MED readings.",
1156 #
1157 # Pour une lecture dans un fichier MED, on peut préciser le nom sous lequel
1158 # le maillage y a été enregistré. Par défaut, on va le chercher sous le nom du concept à créer.
1159 #
1160               NOM_MED    = SIMP(statut='f',typ='TXM',
1161                             fr="Nom du maillage dans le fichier MED.",
1162                             ang="Name of the mesh into the MED file.",),
1163 #
1164               INFO_MED   = SIMP(statut='f',typ='I',defaut= 1,into=(1,2) ),
1165 #
1166                            ) ,
1167 #
1168          INFO            =SIMP(statut='f',typ='I',defaut= 1,into=(1,2) ),
1169 #
1170 )  ;
1171
1172 DEBUT=MACRO(nom="DEBUT",op=0 ,docu="U4.11.01-f1",repetable='n',
1173            fr="Ouverture d une étude. Allocation des ressources mémoire et disque",
1174           sd_prod=ops.DEBUT,
1175
1176          PAR_LOT         =SIMP(fr="mode de traitement des commandes",statut='f',typ='TXM',
1177                            into=("OUI","NON"),defaut="OUI"),
1178          BASE            =FACT(fr="définition des paramètres associés aux bases JEVEUX",
1179                                statut='f',min=01,max=03,
1180            FICHIER         =SIMP(fr="nom de la base",statut='o',typ='TXM',
1181                                  into=('GLOBALE','VOLATILE','LOCALE'),),
1182            TITRE           =SIMP(statut='f',typ='TXM'),
1183            CAS             =SIMP(statut='f',typ='TXM'),
1184            NMAX_ENRE       =SIMP(fr="nombre maximum d enregistrements",statut='f',typ='I'),
1185            LONG_ENRE       =SIMP(fr="longueur des enregistrements",statut='f',typ='I'),
1186            LONG_REPE       =SIMP(fr="longueur du répertoire",statut='f',typ='I'),
1187          ),
1188          IMPRESSION      =FACT(statut='f',min=01,max=03,
1189            FICHIER         =SIMP(statut='o',typ='TXM'),
1190            UNITE           =SIMP(statut='o',typ='I'),
1191          ),
1192          CATALOGUE       =FACT(statut='f',min=01,max=10,
1193            FICHIER         =SIMP(statut='o',typ='TXM'),
1194            TITRE           =SIMP(statut='f',typ='TXM'),
1195            UNITE           =SIMP(statut='f',typ='I'),
1196          ),
1197          CODE            =FACT("définition d un nom pour l'esemble d'une étude",
1198                                statut='f',min=01,max=01,
1199            NOM             =SIMP(statut='o',typ='TXM'),
1200            UNITE           =SIMP(statut='f',typ='I',defaut=15),
1201          ),
1202          DEBUG           =FACT(fr="option de déboggage reservée aux développeurs",
1203                                statut='f',min=01,max=01,
1204            JXVERI          =SIMP(fr="vérifie l intégrité de la segmentation mémoire",
1205                                  statut='f',typ='TXM',into=('OUI','NON'),defaut='NON'),
1206            JEVEUX          =SIMP(fr="force les déchargement sur disque",
1207                                  statut='f',typ='TXM',into=('OUI','NON'),defaut='NON'),
1208            ENVIMA          =SIMP(fr="imprime les valeurs définies dans ENVIMA",
1209                                  statut='f',typ='TXM',into=('TEST',)),
1210          ),
1211          MEMOIRE         =FACT(fr="mode de gestion mémoire utilisé",statut='f',min=01,max=01,
1212            GESTION         =SIMP(statut='f',typ='TXM',into=('COMPACTE','RAPIDE'),defaut='RAPIDE'),
1213            TYPE_ALLOCATION =SIMP(statut='f',typ='I',into=(1,2,3,4),defaut=1),
1214            TAILLE          =SIMP(statut='f',typ='I'),
1215            TAILLE_BLOC     =SIMP(statut='f',typ='R',defaut=800.),
1216            PARTITION       =SIMP(statut='f',typ='R'),
1217          ),
1218  );
1219
1220 FIN=PROC(nom="FIN",op=9999,repetable='n',fr="Fin d'une étude",
1221          docu="U4.11.02-f",
1222          RETASSAGE       =SIMP(fr="provoque le retassage de la base GLOBALE",
1223                                statut='f',typ='TXM',defaut="NON",into=("OUI","NON",) ),
1224          PERFORMANCE     =SIMP(fr="provoque l'impression d'un résumé des mesures de temps ",
1225                                statut='f',typ='TXM',defaut="OUI",into=("OUI","NON",) ),
1226          INFO_RESU       =SIMP(fr="provoque l'impression des informations sur les structures de données",
1227                                statut='f',typ='TXM',defaut="OUI",into=("OUI","NON",) ),
1228          FICHIER         =SIMP(statut='f',typ='TXM',defaut="MESSAGE"),
1229 )  ;
1230
1231