Télécharger la plaquette Cast3M


fr en

Télécharger clim3dj.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  clim3dj.dgibi
  2. *********************************************************************
  3. * VF, CLIM 3D                                                       *
  4. * Jacobians in a rotated mesh                                       *
  5. *                                                                   *
  6. * BECCANTINI A., DM2S/SFME/LTMF,  JANVIER 2002                      *
  7. *********************************************************************
  8.  
  9.  TYEL = 'CU20' ;
  10.  
  11.  'OPTION'  'DIME' 3  'ELEM' TYEL  'ISOV' 'SULI' 'ECHO' 1 'TRAC' 'X' ;
  12.  
  13.  GRAPH = VRAI ;
  14.  GRAPH = FAUX ;
  15.  
  16. ******************
  17. **** MAILLAGE ****
  18. ******************
  19. *
  20. *  P5    P4
  21. *
  22. *  P2    P3
  23. *
  24.  
  25.  DX = 0.5 ;
  26.  P2 = 0.0 0.0 0.0 ;
  27.  P3 = DX  0.0 0.0 ;
  28.  P4 = DX  DX  0.0 ;
  29.  P5 = 0.0 DX  0.0 ;
  30.  P2P3 = P2 'DROIT' 1 P3 ;
  31.  P5P4 = P5 'DROIT' 1 P4 ;
  32.  ENTREE = P2P3 'REGLER' 1 P5P4 ;
  33.  DOMINT = ENTREE 'VOLU' 'TRAN' 4 (0.0 0.0 (4 '*' DX)) ;
  34.  
  35.  ENTREE = ENTREE 'TOURNER' 15 (-1 0.0 0.0) (2.0 0.0 0.0) ;
  36.  DOMINT = DOMINT 'TOURNER' 15 (-1 0.0 0.0) (2.0 0.0 0.0) ;
  37.  
  38.  'ELIMINATION' ENTREE DOMINT (DX '/' 100) ;
  39.  
  40.  DOMINT = 'CHANGER' DOMINT 'QUAF' ;
  41.  ENTREE = 'CHANGER' ENTREE 'QUAF' ;
  42.  
  43.  'ELIMINATION' DOMINT 0.1D-3 ENTREE ;
  44.  
  45. *
  46. **** Creation of DOMAINE tables via the MODEL object
  47. *
  48.  
  49.  MDOMINT  = 'MODELISER' DOMINT 'NAVIER_STOKES' 'LINE' ;
  50.  MENTREE    = 'MODELISER' ENTREE 'NAVIER_STOKES' 'LINE' ;
  51.  
  52.  $DOMINT = 'DOMA' MDOMINT 'TABLE' ;
  53.  $ENTREE   = 'DOMA' MENTREE   'TABLE' ;
  54.  
  55.  DOMINT = 'DOMA' MDOMINT 'MAILLAGE' ;
  56.  ENTREE   = 'DOMA' MENTREE 'MAILLAGE' ;
  57.  
  58.  PCEN = ($DOMINT . 'CENTRE') 'POIN' 'PROC'
  59.    ((DX '/' 2) (DX '/' 2) (DX '/' 2)) ;
  60.  
  61.  'SI' GRAPH ; 
  62.     'TRACER' (DOMINT 'ET' PCEN) ;
  63.  'FINSI' ;
  64.  
  65.  
  66. **********************************************************************************
  67. **********************************************************************************
  68. ***************************  Riemann invariants, inlet ***************************
  69. **********************************************************************************
  70. **********************************************************************************
  71.  
  72.  RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.15 ;
  73.  PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.011 ;
  74.  GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  75.  VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 1.19 'UX' 1.5
  76.     'UZ' 1.11 ;
  77.  
  78.  GN = RN '*' VN ;
  79.  RETN = (('INVERSE' (GAMN '-' 1)) '*' PN)  '+' (0.5 '*' ('PSCAL' VN GN
  80.    ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  81.  
  82.  CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 5 'RN' 0.91 'PN' 3.0
  83.          'UX' 0.35 'UY' 0.3 'UZ' 1.0 ;
  84.  
  85.  LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  86.  LISTP = 'MOTS' 'RN  ' 'UX  ' 'UY  ' 'UZ  ' 'PN  ' ;
  87.  
  88.  RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  89.                  LISTC LISTP
  90.                  RN VN PN GAMN CHPLI 'INRI'  ;
  91.  
  92. *
  93. **** Jacobian with respect to primitive variable
  94. *
  95.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOPRIM' MDOMINT MENTREE
  96.           LISTC LISTP
  97.           RN VN PN GAMN CHPLI 'INRI'  ;
  98.  
  99.  EPS=1.0D-6 ;
  100.  
  101. * RN
  102.  
  103.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  104.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  105.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  106.  
  107.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  108.                  LISTC LISTP
  109.                  (RN '+' DR) VN PN GAMN CHPLI 'INRI'  ;
  110.  
  111.  DRESDRN = ( RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  112.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  113.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  114.  
  115.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  116.  
  117.  'SI' (ERRO > 1.0D-3) ;
  118.      'ERREUR' 5 ;
  119.  'FINSI' ;
  120.  
  121.  
  122. * UX
  123.  
  124.  DUXSCAL = EPS '*'  ('EXTRAIRE' VN 'UX' PCEN) ;
  125.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  126.  
  127.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  128.                  LISTC LISTP
  129.                  RN (VN '+' DUX) PN GAMN CHPLI 'INRI'  ;
  130.  
  131.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  132.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  133.            'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  134.  
  135.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  136.  
  137.  'SI' (ERRO > 1.0D-3) ;
  138.      'ERREUR' 5 ;
  139.  'FINSI' ;
  140.  
  141. * UY
  142.  
  143.  DUYSCAL = EPS '*'  ('EXTRAIRE' VN 'UY' PCEN) ;
  144.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  145.  
  146.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  147.                  LISTC LISTP
  148.                  RN (VN '+' DUY) PN GAMN CHPLI 'INRI'  ;
  149.  
  150.  DRESDUYN = ( RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  151.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  152.            'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  153.  
  154.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  155.  
  156.  'SI' (ERRO > 1.0D-3) ;
  157.      'ERREUR' 5 ;
  158.  'FINSI' ;
  159.  
  160. * UZ
  161.  
  162.  DUZSCAL = EPS '*'  ('EXTRAIRE' VN 'UZ' PCEN) ;
  163.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  164.  
  165.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  166.                  LISTC LISTP
  167.                  RN (VN '+' DUZ) PN GAMN CHPLI 'INRI'  ;
  168.  
  169.  DRESDUZN = ( RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  170.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  171.            'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  172.  
  173.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  174.  
  175.  'SI' (ERRO > 1.0D-3) ;
  176.      'ERREUR' 5 ;
  177.  'FINSI' ;
  178.  
  179. * PN
  180.  
  181.  DPSCAL = EPS '*'  ('EXTRAIRE' PN 'SCAL' PCEN) ;
  182.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  183.  DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  184.  
  185.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  186.                  LISTC LISTP
  187.                  RN  VN (PN '+' DP) GAMN CHPLI 'INRI'  ;
  188.  
  189.  DRESDPN = ( RCHRES1 '-'  RCHRES) '/' DPSCAL ;
  190.  DRESDP  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  191.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  192.  
  193.  ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  194.  
  195.  'SI' (ERRO > 1.0D-3) ;
  196.      'ERREUR' 5 ;
  197.  'FINSI' ;
  198.  
  199. *
  200. **** Jacobian with respect to conservative variables
  201. *
  202.  
  203.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOCONS' MDOMINT MENTREE
  204.           LISTC LISTP
  205.           RN VN PN GAMN CHPLI 'INRI'  ;
  206.  
  207.  EPS=1.0D-6 ;
  208.  
  209. * RN
  210.  
  211.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  212.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  213.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  214.  
  215.  RN1 = RN '+' DR ;
  216.  GN1 = GN ;
  217.  RETN1 = RETN ;
  218.  
  219.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  220.  
  221.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  222.                  LISTC LISTP
  223.                  RN1 VN1 PN1 GAMN CHPLI 'INRI'  ;
  224.  
  225.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  226.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  227.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  228.  
  229.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  230.  
  231.  'SI' (ERRO > 1.0D-3) ;
  232.      'ERREUR' 5 ;
  233.  'FINSI' ;
  234.  
  235. * RUX
  236.  
  237.  DUXSCAL = EPS '*'  ('EXTRAIRE' GN 'UX' PCEN) ;
  238.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  239.  
  240.  RN1 = RN ;
  241.  GN1 = GN '+' DUX ;
  242.  RETN1 = RETN ;
  243.  
  244.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  245.  
  246.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  247.                  LISTC LISTP
  248.                  RN1 VN1 PN1 GAMN CHPLI 'INRI'  ;
  249.  
  250.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  251.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  252.            'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  253.  
  254.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  255.  
  256.  'SI' (ERRO > 1.0D-3) ;
  257.      'ERREUR' 5 ;
  258.  'FINSI' ;
  259.  
  260. * RUY
  261.  
  262.  DUYSCAL = EPS '*'  ('EXTRAIRE' GN 'UY' PCEN) ;
  263.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  264.  
  265.  RN1 = RN ;
  266.  GN1 = GN '+' DUY ;
  267.  RETN1 = RETN ;
  268.  
  269.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  270.  
  271.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  272.                  LISTC LISTP
  273.                  RN1 VN1 PN1 GAMN CHPLI 'INRI'  ;
  274.  
  275.  DRESDUYN = (RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  276.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  277.            'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  278.  
  279.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  280.  
  281.  'SI' (ERRO > 1.0D-3) ;
  282.      'ERREUR' 5 ;
  283.  'FINSI' ;
  284.  
  285.  
  286. * RUZ
  287.  
  288.  DUZSCAL = EPS '*'  ('EXTRAIRE' GN 'UZ' PCEN) ;
  289.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  290.  
  291.  RN1 = RN ;
  292.  GN1 = GN '+' DUZ ;
  293.  RETN1 = RETN ;
  294.  
  295.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  296.  
  297.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  298.                  LISTC LISTP
  299.                  RN1 VN1 PN1 GAMN CHPLI 'INRI'  ;
  300.  
  301.  DRESDUZN = (RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  302.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  303.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  304.  
  305.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  306.  
  307.  'SI' (ERRO > 1.0D-3) ;
  308.      'ERREUR' 5 ;
  309.  'FINSI' ;
  310.  
  311. * RETN
  312.  
  313.  DRSCAL = EPS '*'  ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  314.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  315.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  316.  
  317.  RN1 = RN ;
  318.  GN1 = GN ;
  319.  RETN1 = RETN '+' DR ;
  320.  
  321.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  322.  
  323.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  324.                  LISTC LISTP
  325.                  RN1 VN1 PN1 GAMN CHPLI 'INRI'  ;
  326.  
  327.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  328.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  329.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  330.  
  331.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  332.  
  333.  'SI' (ERRO > 1.0D-3) ;
  334.      'ERREUR' 5 ;
  335.  'FINSI' ;
  336.  
  337. **********************************************************************************
  338. **********************************************************************************
  339. ***************************  Riemann invariants, outlet **************************
  340. **********************************************************************************
  341. **********************************************************************************
  342.  
  343.  RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.15 ;
  344.  PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.011 ;
  345.  GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  346.  VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 0.19 'UX' 0.15
  347.     'UZ' -0.11 ;
  348.  
  349.  GN = RN '*' VN ;
  350.  RETN = (('INVERSE' (GAMN '-' 1)) '*' PN)  '+' (0.5 '*' ('PSCAL' VN GN
  351.    ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  352.  
  353.  CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 5 'RN' 0.91 'PN' 1.10
  354.          'UX' 0.35 'UY' 0.3 'UZ' 0.011 ;
  355.  
  356.  LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  357.  LISTP = 'MOTS' 'RN  ' 'UX  ' 'UY  ' 'UZ  ' 'PN  ' ;
  358.  
  359.  RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  360.                  LISTC LISTP
  361.                  RN VN PN GAMN CHPLI 'OUTRI'  ;
  362.  
  363. *
  364. **** Jacobian with respect to primitive variable
  365. *
  366.  
  367.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOPRIM' MDOMINT MENTREE
  368.           LISTC LISTP
  369.           RN VN PN GAMN CHPLI 'OUTRI'  ;
  370.  
  371.  EPS=1.0D-6 ;
  372.  
  373. * RN
  374.  
  375.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  376.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  377.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  378.  
  379.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  380.                  LISTC LISTP
  381.                  (RN '+' DR) VN PN GAMN CHPLI 'OUTRI'  ;
  382.  
  383.  DRESDRN = ( RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  384.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  385.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  386.  
  387.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  388.  
  389.  'SI' (ERRO > 1.0D-3) ;
  390.      'ERREUR' 5 ;
  391.  'FINSI' ;
  392.  
  393. * UX
  394.  
  395.  DUXSCAL = EPS '*'  ('EXTRAIRE' VN 'UX' PCEN) ;
  396.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  397.  
  398.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  399.                  LISTC LISTP
  400.                  RN (VN '+' DUX) PN GAMN CHPLI 'OUTRI'  ;
  401.  
  402.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  403.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  404.            'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  405.  
  406.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  407.  
  408.  'SI' (ERRO > 1.0D-3) ;
  409.      'ERREUR' 5 ;
  410.  'FINSI' ;
  411.  
  412. * UY
  413.  
  414.  DUYSCAL = EPS '*'  ('EXTRAIRE' VN 'UY' PCEN) ;
  415.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  416.  
  417.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  418.                  LISTC LISTP
  419.                  RN (VN '+' DUY) PN GAMN CHPLI 'OUTRI'  ;
  420.  
  421.  DRESDUYN = ( RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  422.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  423.            'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  424.  
  425.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  426.  
  427.  'SI' (ERRO > 1.0D-3) ;
  428.      'ERREUR' 5 ;
  429.  'FINSI' ;
  430.  
  431. * UZ
  432.  
  433.  DUZSCAL = EPS '*'  ('EXTRAIRE' VN 'UZ' PCEN) ;
  434.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  435.  
  436.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  437.                  LISTC LISTP
  438.                  RN (VN '+' DUZ) PN GAMN CHPLI 'OUTRI'  ;
  439.  
  440.  DRESDUZN = ( RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  441.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  442.            'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  443.  
  444.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  445.  
  446.  'SI' (ERRO > 1.0D-3) ;
  447.      'ERREUR' 5 ;
  448.  'FINSI' ;
  449.  
  450. * PN
  451.  
  452.  DPSCAL = EPS '*'  ('EXTRAIRE' PN 'SCAL' PCEN) ;
  453.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  454.  DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  455.  
  456.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  457.                  LISTC LISTP
  458.                  RN  VN (PN '+' DP) GAMN CHPLI 'OUTRI'  ;
  459.  
  460.  DRESDPN = ( RCHRES1 '-'  RCHRES) '/' DPSCAL ;
  461.  DRESDP  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  462.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  463.  
  464.  ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  465.  
  466.  'SI' (ERRO > 1.0D-3) ;
  467.      'ERREUR' 5 ;
  468.  'FINSI' ;
  469.  
  470. *
  471. **** Jacobian with respect to conservative variables
  472. *
  473.  
  474.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOCONS' MDOMINT MENTREE
  475.           LISTC LISTP
  476.           RN VN PN GAMN CHPLI 'OUTRI'  ;
  477.  
  478.  EPS=1.0D-6 ;
  479.  
  480. * RN
  481.  
  482.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  483.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  484.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  485.  
  486.  RN1 = RN '+' DR ;
  487.  GN1 = GN ;
  488.  RETN1 = RETN ;
  489.  
  490.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  491.  
  492.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  493.                  LISTC LISTP
  494.                  RN1 VN1 PN1 GAMN CHPLI 'OUTRI'  ;
  495.  
  496.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  497.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  498.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  499.  
  500.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  501.  
  502.  'SI' (ERRO > 1.0D-3) ;
  503.      'ERREUR' 5 ;
  504.  'FINSI' ;
  505.  
  506. * RUX
  507.  
  508.  DUXSCAL = EPS '*'  ('EXTRAIRE' GN 'UX' PCEN) ;
  509.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  510.  
  511.  RN1 = RN ;
  512.  GN1 = GN '+' DUX ;
  513.  RETN1 = RETN ;
  514.  
  515.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  516.  
  517.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  518.                  LISTC LISTP
  519.                  RN1 VN1 PN1 GAMN CHPLI 'OUTRI'  ;
  520.  
  521.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  522.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  523.            'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  524.  
  525.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  526.  
  527.  'SI' (ERRO > 1.0D-3) ;
  528.      'ERREUR' 5 ;
  529.  'FINSI' ;
  530.  
  531. * RUY
  532.  
  533.  DUYSCAL = EPS '*'  ('EXTRAIRE' GN 'UY' PCEN) ;
  534.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  535.  
  536.  RN1 = RN ;
  537.  GN1 = GN '+' DUY ;
  538.  RETN1 = RETN ;
  539.  
  540.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  541.  
  542.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  543.                  LISTC LISTP
  544.                  RN1 VN1 PN1 GAMN CHPLI 'OUTRI'  ;
  545.  
  546.  DRESDUYN = (RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  547.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  548.            'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  549.  
  550.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  551.  
  552.  'SI' (ERRO > 1.0D-3) ;
  553.      'ERREUR' 5 ;
  554.  'FINSI' ;
  555.  
  556.  
  557. * RUZ
  558.  
  559.  DUZSCAL = EPS '*'  ('EXTRAIRE' GN 'UZ' PCEN) ;
  560.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  561.  
  562.  RN1 = RN ;
  563.  GN1 = GN '+' DUZ ;
  564.  RETN1 = RETN ;
  565.  
  566.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  567.  
  568.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  569.                  LISTC LISTP
  570.                  RN1 VN1 PN1 GAMN CHPLI 'OUTRI'  ;
  571.  
  572.  DRESDUZN = (RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  573.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  574.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  575.  
  576.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  577.  
  578.  'SI' (ERRO > 1.0D-3) ;
  579.      'ERREUR' 5 ;
  580.  'FINSI' ;
  581.  
  582. * RETN
  583.  
  584.  DRSCAL = EPS '*'  ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  585.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  586.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  587.  
  588.  RN1 = RN ;
  589.  GN1 = GN ;
  590.  RETN1 = RETN '+' DR ;
  591.  
  592.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  593.  
  594.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  595.                  LISTC LISTP
  596.                  RN1 VN1 PN1 GAMN CHPLI 'OUTRI'  ;
  597.  
  598.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  599.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  600.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  601.  
  602.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  603.  
  604.  'SI' (ERRO > 1.0D-3) ;
  605.      'ERREUR' 5 ;
  606.  'FINSI' ;
  607.  
  608. **********************************************************************************
  609. **********************************************************************************
  610. ***************************  Supersonic inlet ************************************
  611. **********************************************************************************
  612. **********************************************************************************
  613.  
  614.  RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  615.  PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  616.  GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  617.  VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 15. 'UX' 1. 'UZ' 1.7;
  618.  
  619.  GN = RN '*' VN ;
  620.  RETN = (('INVERSE' (GAMN '-' 1)) '*' PN)  '+' (0.5 '*' ('PSCAL' VN GN
  621.    ('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
  622.  
  623.  CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 5 'RN' 1.11 'PN' 2.13
  624.          'UX' -11. 'UY' 2. 'UZ' 1.8 ;
  625.  
  626.  LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  627.  LISTP = 'MOTS' 'RN  ' 'UX  ' 'UY  ' 'UZ  ' 'PN  ' ;
  628.  
  629.  RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  630.                  LISTC LISTP
  631.                  RN VN PN GAMN CHPLI 'INSS'  ;
  632.  
  633.  ERRO = 'MAXIMUM' (RCHLIM '-' CHPLI) 'ABS' ;
  634.  
  635.  'SI' (ERRO > 1.0D-12) ;
  636.      'ERREUR' 5 ;
  637.  'FINSI' ;
  638.  
  639.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOCONS' MDOMINT MENTREE
  640.           LISTC LISTP
  641.           RN VN PN GAMN CHPLI 'INSS'  ;
  642.  
  643.  AA = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 5 'VAR1' 1.11 'VAR2' 1.13
  644.  'VAR3' 1.15 'VAR4' 1.17 'VAR5' 1.37 ;
  645.  
  646.  ERRO = 'MAXIMUM' ('KOPS' 'MULT' RJACO AA) 'ABS' ;
  647.  
  648.  'SI' (ERRO > 1.0D-12) ;
  649.      'ERREUR' 5 ;
  650.  'FINSI' ;
  651.  
  652.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOPRIM' MDOMINT MENTREE
  653.           LISTC LISTP
  654.           RN VN PN GAMN CHPLI 'INSS'  ;
  655.  
  656.  AA = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'RN' 1.11 'UX' 1.13
  657.  'UY' 1.15 'PN' 1.17 ;
  658.  
  659.  ERRO = 'MAXIMUM' ('KOPS' 'MULT' RJACO AA) 'ABS' ;
  660.  
  661.  'SI' (ERRO > 1.0D-12) ;
  662.      'ERREUR' 5 ;
  663.  'FINSI' ;
  664.  
  665.  
  666. **********************************************************************************
  667. **********************************************************************************
  668. ***************************  Supersonic outlet ***********************************
  669. **********************************************************************************
  670. **********************************************************************************
  671.  
  672.  RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  673.  PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  674.  GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  675.  VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' -15. 'UX' 1. 'UZ' -5.
  676.     ;
  677.  
  678.  GN = RN '*' VN ;
  679.  RETN = (('INVERSE' (GAMN '-' 1)) '*' PN)  '+' (0.5 '*' ('PSCAL' VN GN
  680.    ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  681.  
  682.  CHPLI CACCA = 'KOPS' MATRIK ;
  683.  
  684.  LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  685.  LISTP = 'MOTS' 'RN  ' 'UX  ' 'UY  ' 'UZ  ' 'PN  ' ;
  686.  
  687.  RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  688.                  LISTC LISTP
  689.                  RN VN PN GAMN CHPLI 'OUTSS'  ;
  690.  
  691. *
  692. **** Jacobian with respect to primitive variable
  693. *
  694.  
  695.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOPRIM' MDOMINT MENTREE
  696.           LISTC LISTP
  697.           RN VN PN GAMN CHPLI 'OUTSS'  ;
  698.  
  699.  EPS=1.0D-6 ;
  700.  
  701. * RN
  702.  
  703.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  704.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  705.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  706.  
  707.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  708.                  LISTC LISTP
  709.                  (RN '+' DR) VN PN GAMN CHPLI 'OUTSS'  ;
  710.  
  711.  DRESDRN = ( RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  712.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  713.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  714.  
  715.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  716.  
  717.  
  718.  'SI' (ERRO > 1.0D-3) ;
  719.      'ERREUR' 5 ;
  720.  'FINSI' ;
  721.  
  722. * UX
  723.  
  724.  DUXSCAL = EPS '*'  ('EXTRAIRE' VN 'UX' PCEN) ;
  725.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  726.  
  727.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  728.                  LISTC LISTP
  729.                  RN (VN '+' DUX) PN GAMN CHPLI 'OUTSS'  ;
  730.  
  731.  DRESDUXN = (RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  732.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  733.            'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  734.  
  735.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  736.  
  737.  'SI' (ERRO > 1.0D-3) ;
  738.      'ERREUR' 5 ;
  739.  'FINSI' ;
  740.  
  741. * UY
  742.  
  743.  DUYSCAL = EPS '*'  ('EXTRAIRE' VN 'UY' PCEN) ;
  744.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  745.  
  746.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  747.                  LISTC LISTP
  748.                  RN (VN '+' DUY) PN GAMN CHPLI 'OUTSS'  ;
  749.  
  750.  DRESDUYN = ( RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  751.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  752.            'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  753.  
  754.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  755.  
  756.  'SI' (ERRO > 1.0D-3) ;
  757.      'ERREUR' 5 ;
  758.  'FINSI' ;
  759.  
  760. * UZ
  761.  
  762.  DUZSCAL = EPS '*'  ('EXTRAIRE' VN 'UZ' PCEN) ;
  763.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  764.  
  765.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  766.                  LISTC LISTP
  767.                  RN (VN '+' DUZ) PN GAMN CHPLI 'OUTSS'  ;
  768.  
  769.  DRESDUZN = ( RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  770.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  771.            'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  772.  
  773.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  774.  
  775.  'SI' (ERRO > 1.0D-3) ;
  776.      'ERREUR' 5 ;
  777.  'FINSI' ;
  778.  
  779. * PN
  780.  
  781.  DPSCAL = EPS '*'  ('EXTRAIRE' PN 'SCAL' PCEN) ;
  782.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  783.  DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  784.  
  785.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  786.                  LISTC LISTP
  787.                  RN  VN (PN '+' DP) GAMN CHPLI 'OUTSS'  ;
  788.  
  789.  DRESDPN = ( RCHRES1 '-'  RCHRES) '/' DPSCAL ;
  790.  DRESDP  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  791.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  792.  
  793.  ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  794.  
  795.  'SI' (ERRO > 1.0D-3) ;
  796.      'ERREUR' 5 ;
  797.  'FINSI' ;
  798.  
  799. *
  800. **** Jacobian with respect to conservative variables
  801. *
  802.  
  803.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOCONS' MDOMINT MENTREE
  804.           LISTC LISTP
  805.           RN VN PN GAMN CHPLI 'OUTSS'  ;
  806.  
  807.  EPS=1.0D-6 ;
  808.  
  809. * RN
  810.  
  811.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  812.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  813.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  814.  
  815.  RN1 = RN '+' DR ;
  816.  GN1 = GN ;
  817.  RETN1 = RETN ;
  818.  
  819.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  820.  
  821.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  822.                  LISTC LISTP
  823.                  RN1 VN1 PN1 GAMN CHPLI 'OUTSS'  ;
  824.  
  825.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  826.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  827.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  828.  
  829.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  830.  
  831.  'SI' (ERRO > 1.0D-3) ;
  832.      'ERREUR' 5 ;
  833.  'FINSI' ;
  834.  
  835. * RUX
  836.  
  837.  DUXSCAL = EPS '*'  ('EXTRAIRE' GN 'UX' PCEN) ;
  838.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  839.  
  840.  RN1 = RN ;
  841.  GN1 = GN '+' DUX ;
  842.  RETN1 = RETN ;
  843.  
  844.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  845.  
  846.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  847.                  LISTC LISTP
  848.                  RN1 VN1 PN1 GAMN CHPLI 'OUTSS'  ;
  849.  
  850.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  851.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  852.            'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  853.  
  854.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  855.  
  856.  'SI' (ERRO > 1.0D-3) ;
  857.      'ERREUR' 5 ;
  858.  'FINSI' ;
  859.  
  860. * RUY
  861.  
  862.  DUYSCAL = EPS '*'  ('EXTRAIRE' GN 'UY' PCEN) ;
  863.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  864.  
  865.  RN1 = RN ;
  866.  GN1 = GN '+' DUY ;
  867.  RETN1 = RETN ;
  868.  
  869.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  870.  
  871.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  872.                  LISTC LISTP
  873.                  RN1 VN1 PN1 GAMN CHPLI 'OUTSS'  ;
  874.  
  875.  DRESDUYN = (RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  876.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  877.            'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  878.  
  879.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  880.  
  881.  'SI' (ERRO > 1.0D-3) ;
  882.      'ERREUR' 5 ;
  883.  'FINSI' ;
  884.  
  885.  
  886. * RUZ
  887.  
  888.  DUZSCAL = EPS '*'  ('EXTRAIRE' GN 'UZ' PCEN) ;
  889.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  890.  
  891.  RN1 = RN ;
  892.  GN1 = GN '+' DUZ ;
  893.  RETN1 = RETN ;
  894.  
  895.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  896.  
  897.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  898.                  LISTC LISTP
  899.                  RN1 VN1 PN1 GAMN CHPLI 'OUTSS'  ;
  900.  
  901.  DRESDUZN = (RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  902.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  903.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  904.  
  905.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  906.  
  907.  'SI' (ERRO > 1.0D-3) ;
  908.      'ERREUR' 5 ;
  909.  'FINSI' ;
  910.  
  911. * RETN
  912.  
  913.  DRSCAL = EPS '*'  ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  914.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  915.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  916.  
  917.  RN1 = RN ;
  918.  GN1 = GN ;
  919.  RETN1 = RETN '+' DR ;
  920.  
  921.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  922.  
  923.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  924.                  LISTC LISTP
  925.                  RN1 VN1 PN1 GAMN CHPLI 'OUTSS'  ;
  926.  
  927.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  928.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  929.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  930.  
  931.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  932.  
  933.  'SI' (ERRO > 1.0D-3) ;
  934.      'ERREUR' 5 ;
  935.  'FINSI' ;
  936.  
  937. **********************************************************************************
  938. **********************************************************************************
  939. ***************************  Pressure at the outlet ******************************
  940. **********************************************************************************
  941. **********************************************************************************
  942. *
  943. * We solve a Riemann probleme between the state in the center and
  944. * the same state but the pressure equal to the pressure at the outlet
  945. * We use AUSMPLUS
  946. *
  947. *
  948. ***** Calcul du residu dans le cas particulier PF = PC
  949. *
  950.  
  951.  RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  952.  PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  953.  GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  954.  VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' -15. 'UX' 1.
  955.       'UZ' -3.1 ;
  956.  
  957.  GN = RN '*' VN ;
  958.  RETN = (('INVERSE' (GAMN '-' 1)) '*' PN)  '+' (0.5 '*' ('PSCAL' VN GN
  959.    ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  960.  
  961.  CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 1 'PN' 2.1 ;
  962.  
  963.  LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  964.  LISTP = 'MOTS' 'RN  ' 'UX  ' 'UY  ' 'UZ  ' 'PN  ' ;
  965.  
  966.  RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  967.                  LISTC LISTP
  968.                  RN VN PN GAMN CHPLI 'OUTP'  ;
  969.  
  970. *
  971. **** Jacobian with respect to primitive variable
  972. *
  973.  
  974.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOPRIM' MDOMINT MENTREE
  975.           LISTC LISTP
  976.           RN VN PN GAMN CHPLI 'OUTP'  ;
  977.  
  978.  EPS=1.0D-6 ;
  979.  
  980. * RN
  981.  
  982.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  983.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  984.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  985.  
  986.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  987.                  LISTC LISTP
  988.                  (RN '+' DR) VN PN GAMN CHPLI 'OUTP'  ;
  989.  
  990.  DRESDRN = ( RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  991.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  992.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  993.  
  994.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  995.  
  996.  'SI' (ERRO > 1.0D-3) ;
  997.      'ERREUR' 5 ;
  998.  'FINSI' ;
  999.  
  1000.  
  1001. * UX
  1002.  
  1003.  DUXSCAL = EPS '*'  ('EXTRAIRE' VN 'UX' PCEN) ;
  1004.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1005.  
  1006.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1007.                  LISTC LISTP
  1008.                  RN (VN '+' DUX) PN GAMN CHPLI 'OUTP'  ;
  1009.  
  1010.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  1011.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1012.            'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1013.  
  1014.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1015.  
  1016.  'SI' (ERRO > 1.0D-3) ;
  1017.      'ERREUR' 5 ;
  1018.  'FINSI' ;
  1019.  
  1020. * UY
  1021.  
  1022.  DUYSCAL = EPS '*'  ('EXTRAIRE' VN 'UY' PCEN) ;
  1023.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1024.  
  1025.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1026.                  LISTC LISTP
  1027.                  RN (VN '+' DUY) PN GAMN CHPLI 'OUTP'  ;
  1028.  
  1029.  DRESDUYN = ( RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  1030.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1031.            'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  1032.  
  1033.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1034.  
  1035.  'SI' (ERRO > 1.0D-3) ;
  1036.      'ERREUR' 5 ;
  1037.  'FINSI' ;
  1038.  
  1039. * UZ
  1040.  
  1041.  DUZSCAL = EPS '*'  ('EXTRAIRE' VN 'UZ' PCEN) ;
  1042.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1043.  
  1044.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1045.                  LISTC LISTP
  1046.                  RN (VN '+' DUZ) PN GAMN CHPLI 'OUTP'  ;
  1047.  
  1048.  DRESDUZN = ( RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  1049.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1050.            'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  1051.  
  1052.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1053.  
  1054.  'SI' (ERRO > 1.0D-3) ;
  1055.      'ERREUR' 5 ;
  1056.  'FINSI' ;
  1057.  
  1058. * PN
  1059.  
  1060.  DPSCAL = EPS '*'  ('EXTRAIRE' PN 'SCAL' PCEN) ;
  1061.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1062.  DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  1063.  
  1064.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1065.                  LISTC LISTP
  1066.                  RN  VN (PN '+' DP) GAMN CHPLI 'OUTP'  ;
  1067.  
  1068.  DRESDPN = ( RCHRES1 '-'  RCHRES) '/' DPSCAL ;
  1069.  DRESDP  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1070.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  1071.  
  1072.  ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  1073.  
  1074.  'SI' (ERRO > 1.0D-3) ;
  1075.      'ERREUR' 5 ;
  1076.  'FINSI' ;
  1077.  
  1078. *
  1079. **** Jacobian with respect to conservative variables
  1080. *
  1081.  
  1082.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOCONS' MDOMINT MENTREE
  1083.           LISTC LISTP
  1084.           RN VN PN GAMN CHPLI 'OUTP'  ;
  1085.  
  1086.  EPS=1.0D-6 ;
  1087.  
  1088. * RN
  1089.  
  1090.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1091.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1092.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1093.  
  1094.  RN1 = RN '+' DR ;
  1095.  GN1 = GN ;
  1096.  RETN1 = RETN ;
  1097.  
  1098.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1099.  
  1100.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1101.                  LISTC LISTP
  1102.                  RN1 VN1 PN1 GAMN CHPLI 'OUTP'  ;
  1103.  
  1104.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  1105.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  1106.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1107.  
  1108.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1109.  
  1110.  'SI' (ERRO > 1.0D-3) ;
  1111.      'ERREUR' 5 ;
  1112.  'FINSI' ;
  1113.  
  1114. * RUX
  1115.  
  1116.  DUXSCAL = EPS '*'  ('EXTRAIRE' GN 'UX' PCEN) ;
  1117.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1118.  
  1119.  RN1 = RN ;
  1120.  GN1 = GN '+' DUX ;
  1121.  RETN1 = RETN ;
  1122.  
  1123.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1124.  
  1125.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1126.                  LISTC LISTP
  1127.                  RN1 VN1 PN1 GAMN CHPLI 'OUTP'  ;
  1128.  
  1129.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  1130.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1131.            'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1132.  
  1133.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1134.  
  1135.  'SI' (ERRO > 1.0D-3) ;
  1136.      'ERREUR' 5 ;
  1137.  'FINSI' ;
  1138.  
  1139. * RUY
  1140.  
  1141.  DUYSCAL = EPS '*'  ('EXTRAIRE' GN 'UY' PCEN) ;
  1142.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1143.  
  1144.  RN1 = RN ;
  1145.  GN1 = GN '+' DUY ;
  1146.  RETN1 = RETN ;
  1147.  
  1148.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1149.  
  1150.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1151.                  LISTC LISTP
  1152.                  RN1 VN1 PN1 GAMN CHPLI 'OUTP'  ;
  1153.  
  1154.  DRESDUYN = (RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  1155.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1156.            'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1157.  
  1158.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1159.  
  1160.  'SI' (ERRO > 1.0D-3) ;
  1161.      'ERREUR' 5 ;
  1162.  'FINSI' ;
  1163.  
  1164.  
  1165. * RUZ
  1166.  
  1167.  DUZSCAL = EPS '*'  ('EXTRAIRE' GN 'UZ' PCEN) ;
  1168.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1169.  
  1170.  RN1 = RN ;
  1171.  GN1 = GN '+' DUZ ;
  1172.  RETN1 = RETN ;
  1173.  
  1174.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1175.  
  1176.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1177.                  LISTC LISTP
  1178.                  RN1 VN1 PN1 GAMN CHPLI 'OUTP'  ;
  1179.  
  1180.  DRESDUZN = (RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  1181.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1182.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  1183.  
  1184.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1185.  
  1186.  'SI' (ERRO > 1.0D-3) ;
  1187.      'ERREUR' 5 ;
  1188.  'FINSI' ;
  1189.  
  1190. * RETN
  1191.  
  1192.  DRSCAL = EPS '*'  ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  1193.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1194.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1195.  
  1196.  RN1 = RN ;
  1197.  GN1 = GN ;
  1198.  RETN1 = RETN '+' DR ;
  1199.  
  1200.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1201.  
  1202.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1203.                  LISTC LISTP
  1204.                  RN1 VN1 PN1 GAMN CHPLI 'OUTP'  ;
  1205.  
  1206.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  1207.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1208.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  1209.  
  1210.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1211.  
  1212.  'SI' (ERRO > 1.0D-3) ;
  1213.      'ERREUR' 5 ;
  1214.  'FINSI' ;
  1215.  
  1216. **********************************************************************************
  1217. **********************************************************************************
  1218. ***************************  INSU ************************************************
  1219. **********************************************************************************
  1220. **********************************************************************************
  1221. *
  1222. * Total enthalpy, entropy, tangential velocity at the inlet
  1223. * Normal velocity inside
  1224. *
  1225. * We solve a Riemann probleme between the state at the interface thus
  1226. * reconstructed and the state at the center
  1227. *
  1228.  
  1229. *
  1230. **** Jacobian with respect to primitive variable
  1231. *
  1232.  
  1233.  RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  1234.  PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  1235.  GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  1236.  VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 1.13 'UX' 1.16
  1237.    'UZ' 1.11  ;
  1238.  GN = RN * VN ;
  1239.  RETN = (('INVERSE' (GAMN '-' 1)) '*' PN)  '+' (0.5 '*' ('PSCAL' VN GN
  1240.    ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  1241.  
  1242.  RHTN = RETN '+' PN ;
  1243.  HTN = RHTN '/' RN ;
  1244.  
  1245.  SN = PN '/' (RN '**' ('MAXI' GAMN)) ;
  1246.  
  1247.  CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 2 'HT' (('MAXIMUM' HTN) *
  1248.    1.11) 'S' (('MAXIMUM' SN) '*' 1.13) ;
  1249.  
  1250.  
  1251.  LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  1252.  LISTP = 'MOTS' 'RN  ' 'UX  ' 'UY  ' 'UZ  ' 'PN  ' ;
  1253.  
  1254.  RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1255.                  LISTC LISTP
  1256.                  RN VN PN GAMN CHPLI 'INSU'  ;
  1257.  
  1258.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOPRIM' MDOMINT MENTREE
  1259.           LISTC LISTP
  1260.           RN VN PN GAMN CHPLI 'INSU'  ;
  1261.  
  1262.  EPS=1.0D-6 ;
  1263.  
  1264. * RN
  1265.  
  1266.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1267.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1268.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1269.  
  1270.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1271.                  LISTC LISTP
  1272.                  (RN '+' DR) VN PN GAMN CHPLI 'INSU'  ;
  1273.  
  1274.  DRESDRN = ( RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  1275.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  1276.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1277.  
  1278.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1279.  
  1280.  'SI' (ERRO > 1.0D-3) ;
  1281.      'ERREUR' 5 ;
  1282.  'FINSI' ;
  1283.  
  1284.  
  1285. * UX
  1286.  
  1287.  DUXSCAL = EPS '*'  ('EXTRAIRE' VN 'UX' PCEN) ;
  1288.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1289.  
  1290.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1291.                  LISTC LISTP
  1292.                  RN (VN '+' DUX) PN GAMN CHPLI 'INSU'  ;
  1293.  
  1294.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  1295.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1296.            'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1297.  
  1298.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1299.  
  1300.  'SI' (ERRO > 1.0D-3) ;
  1301.      'ERREUR' 5 ;
  1302.  'FINSI' ;
  1303.  
  1304. * UY
  1305.  
  1306.  DUYSCAL = EPS '*'  ('EXTRAIRE' VN 'UY' PCEN) ;
  1307.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1308.  
  1309.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1310.                  LISTC LISTP
  1311.                  RN (VN '+' DUY) PN GAMN CHPLI 'INSU'  ;
  1312.  
  1313.  DRESDUYN = ( RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  1314.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1315.            'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  1316.  
  1317.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1318.  
  1319.  'SI' (ERRO > 1.0D-3) ;
  1320.      'ERREUR' 5 ;
  1321.  'FINSI' ;
  1322.  
  1323. * UZ
  1324.  
  1325.  DUZSCAL = EPS '*'  ('EXTRAIRE' VN 'UZ' PCEN) ;
  1326.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1327.  
  1328.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1329.                  LISTC LISTP
  1330.                  RN (VN '+' DUZ) PN GAMN CHPLI 'INSU'  ;
  1331.  
  1332.  DRESDUZN = ( RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  1333.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1334.            'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  1335.  
  1336.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1337.  
  1338.  'SI' (ERRO > 1.0D-3) ;
  1339.      'ERREUR' 5 ;
  1340.  'FINSI' ;
  1341.  
  1342. * PN
  1343.  
  1344.  DPSCAL = EPS '*'  ('EXTRAIRE' PN 'SCAL' PCEN) ;
  1345.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1346.  DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  1347.  
  1348.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1349.                  LISTC LISTP
  1350.                  RN  VN (PN '+' DP) GAMN CHPLI 'INSU'  ;
  1351.  
  1352.  DRESDPN = ( RCHRES1 '-'  RCHRES) '/' DPSCAL ;
  1353.  DRESDP  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1354.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  1355.  
  1356.  ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  1357.  
  1358.  'SI' (ERRO > 1.0D-3) ;
  1359.      'ERREUR' 5 ;
  1360.  'FINSI' ;
  1361.  
  1362. *
  1363. **** Jacobian with respect to conservative variables
  1364. *
  1365.  
  1366.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOCONS' MDOMINT MENTREE
  1367.           LISTC LISTP
  1368.           RN VN PN GAMN CHPLI 'INSU'  ;
  1369.  
  1370.  EPS=1.0D-6 ;
  1371.  
  1372. * RN
  1373.  
  1374.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1375.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1376.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1377.  
  1378.  RN1 = RN '+' DR ;
  1379.  GN1 = GN ;
  1380.  RETN1 = RETN ;
  1381.  
  1382.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1383.  
  1384.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1385.                  LISTC LISTP
  1386.                  RN1 VN1 PN1 GAMN CHPLI 'INSU'  ;
  1387.  
  1388.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  1389.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  1390.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1391.  
  1392.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1393.  
  1394.  'SI' (ERRO > 1.0D-3) ;
  1395.      'ERREUR' 5 ;
  1396.  'FINSI' ;
  1397.  
  1398. * RUX
  1399.  
  1400.  DUXSCAL = EPS '*'  ('EXTRAIRE' GN 'UX' PCEN) ;
  1401.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1402.  
  1403.  RN1 = RN ;
  1404.  GN1 = GN '+' DUX ;
  1405.  RETN1 = RETN ;
  1406.  
  1407.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1408.  
  1409.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1410.                  LISTC LISTP
  1411.                  RN1 VN1 PN1 GAMN CHPLI 'INSU'  ;
  1412.  
  1413.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  1414.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1415.            'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1416.  
  1417.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1418.  
  1419.  'SI' (ERRO > 1.0D-3) ;
  1420.      'ERREUR' 5 ;
  1421.  'FINSI' ;
  1422.  
  1423. * RUY
  1424.  
  1425.  DUYSCAL = EPS '*'  ('EXTRAIRE' GN 'UY' PCEN) ;
  1426.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1427.  
  1428.  RN1 = RN ;
  1429.  GN1 = GN '+' DUY ;
  1430.  RETN1 = RETN ;
  1431.  
  1432.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1433.  
  1434.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1435.                  LISTC LISTP
  1436.                  RN1 VN1 PN1 GAMN CHPLI 'INSU'  ;
  1437.  
  1438.  DRESDUYN = (RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  1439.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1440.            'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1441.  
  1442.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1443.  
  1444.  'SI' (ERRO > 1.0D-3) ;
  1445.      'ERREUR' 5 ;
  1446.  'FINSI' ;
  1447.  
  1448.  
  1449. * RUZ
  1450.  
  1451.  DUZSCAL = EPS '*'  ('EXTRAIRE' GN 'UZ' PCEN) ;
  1452.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1453.  
  1454.  RN1 = RN ;
  1455.  GN1 = GN '+' DUZ ;
  1456.  RETN1 = RETN ;
  1457.  
  1458.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1459.  
  1460.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1461.                  LISTC LISTP
  1462.                  RN1 VN1 PN1 GAMN CHPLI 'INSU'  ;
  1463.  
  1464.  DRESDUZN = (RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  1465.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1466.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  1467.  
  1468.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1469.  
  1470.  'SI' (ERRO > 1.0D-3) ;
  1471.      'ERREUR' 5 ;
  1472.  'FINSI' ;
  1473.  
  1474. * RETN
  1475.  
  1476.  DRSCAL = EPS '*'  ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  1477.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1478.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1479.  
  1480.  RN1 = RN ;
  1481.  GN1 = GN ;
  1482.  RETN1 = RETN '+' DR ;
  1483.  
  1484.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1485.  
  1486.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1487.                  LISTC LISTP
  1488.                  RN1 VN1 PN1 GAMN CHPLI 'INSU'  ;
  1489.  
  1490.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  1491.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1492.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  1493.  
  1494.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1495.  
  1496.  'SI' (ERRO > 1.0D-3) ;
  1497.      'ERREUR' 5 ;
  1498.  'FINSI' ;
  1499.  
  1500.  
  1501. **********************************************************************************
  1502. **********************************************************************************
  1503. ***************************  INJE ************************************************
  1504. **********************************************************************************
  1505. **********************************************************************************
  1506. *
  1507. * We impose RT and the normal momontum at the inlet.
  1508.  
  1509.  RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  1510.  PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  1511.  GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  1512.  VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 1.11 'UX' -1.3 'UZ'
  1513.    1.5 ;
  1514.  
  1515.  GN = RN '*' VN ;
  1516.  RETN = (('INVERSE' (GAMN '-' 1)) '*' PN)  '+' (0.5 '*' ('PSCAL' VN GN
  1517.    ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  1518.  
  1519.  MOME = 0.111 ;
  1520.  RT   = 4.1 ;
  1521.  
  1522.  CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 2 'RT' RT
  1523.     'MOME' MOME ;
  1524.  
  1525.  LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  1526.  LISTP = 'MOTS' 'RN  ' 'UX  ' 'UY  ' 'UZ  ' 'PN  ' ;
  1527.  
  1528.  RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1529.                  LISTC LISTP
  1530.                  RN VN PN GAMN CHPLI 'INJE'  ;
  1531.  
  1532. *
  1533. **** Jacobian with respect to primitive variable
  1534. *
  1535.  
  1536.  
  1537.  RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1538.                  LISTC LISTP
  1539.                  RN VN PN GAMN CHPLI 'INJE'  ;
  1540.  
  1541.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOPRIM' MDOMINT MENTREE
  1542.           LISTC LISTP
  1543.           RN VN PN GAMN CHPLI 'INJE'  ;
  1544.  
  1545.  EPS=1.0D-6 ;
  1546.  
  1547. * RN
  1548.  
  1549.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1550.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1551.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1552.  
  1553.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1554.                  LISTC LISTP
  1555.                  (RN '+' DR) VN PN GAMN CHPLI 'INJE'  ;
  1556.  
  1557.  DRESDRN = ( RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  1558.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  1559.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1560.  
  1561.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1562.  
  1563.  'SI' (ERRO > 1.0D-3) ;
  1564.      'ERREUR' 5 ;
  1565.  'FINSI' ;
  1566.  
  1567.  
  1568. * UX
  1569.  
  1570.  DUXSCAL = EPS '*'  ('EXTRAIRE' VN 'UX' PCEN) ;
  1571.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1572.  
  1573.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1574.                  LISTC LISTP
  1575.                  RN (VN '+' DUX) PN GAMN CHPLI 'INJE'  ;
  1576.  
  1577.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  1578.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1579.            'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1580.  
  1581.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1582.  
  1583.  'SI' (ERRO > 1.0D-3) ;
  1584.      'ERREUR' 5 ;
  1585.  'FINSI' ;
  1586.  
  1587. * UY
  1588.  
  1589.  DUYSCAL = EPS '*'  ('EXTRAIRE' VN 'UY' PCEN) ;
  1590.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1591.  
  1592.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1593.                  LISTC LISTP
  1594.                  RN (VN '+' DUY) PN GAMN CHPLI 'INJE'  ;
  1595.  
  1596.  DRESDUYN = ( RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  1597.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1598.            'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  1599.  
  1600.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1601.  
  1602.  'SI' (ERRO > 1.0D-3) ;
  1603.      'ERREUR' 5 ;
  1604.  'FINSI' ;
  1605.  
  1606. * UZ
  1607.  
  1608.  DUZSCAL = EPS '*'  ('EXTRAIRE' VN 'UZ' PCEN) ;
  1609.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1610.  
  1611.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1612.                  LISTC LISTP
  1613.                  RN (VN '+' DUZ) PN GAMN CHPLI 'INJE'  ;
  1614.  
  1615.  DRESDUZN = ( RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  1616.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1617.            'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  1618.  
  1619.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1620.  
  1621.  'SI' (ERRO > 1.0D-3) ;
  1622.      'ERREUR' 5 ;
  1623.  'FINSI' ;
  1624.  
  1625. * PN
  1626.  
  1627.  DPSCAL = EPS '*'  ('EXTRAIRE' PN 'SCAL' PCEN) ;
  1628.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1629.  DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  1630.  
  1631.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1632.                  LISTC LISTP
  1633.                  RN  VN (PN '+' DP) GAMN CHPLI 'INJE'  ;
  1634.  
  1635.  DRESDPN = ( RCHRES1 '-'  RCHRES) '/' DPSCAL ;
  1636.  DRESDP  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1637.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  1638.  
  1639.  ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  1640.  
  1641.  'SI' (ERRO > 1.0D-3) ;
  1642.      'ERREUR' 5 ;
  1643.  'FINSI' ;
  1644.  
  1645. *
  1646. **** Jacobian with respect to conservative variables
  1647. *
  1648.  
  1649.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOCONS' MDOMINT MENTREE
  1650.           LISTC LISTP
  1651.           RN VN PN GAMN CHPLI 'INJE'  ;
  1652.  
  1653.  EPS=1.0D-6 ;
  1654.  
  1655. * RN
  1656.  
  1657.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1658.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1659.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1660.  
  1661.  RN1 = RN '+' DR ;
  1662.  GN1 = GN ;
  1663.  RETN1 = RETN ;
  1664.  
  1665.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1666.  
  1667.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1668.                  LISTC LISTP
  1669.                  RN1 VN1 PN1 GAMN CHPLI 'INJE'  ;
  1670.  
  1671.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  1672.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  1673.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1674.  
  1675.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1676.  
  1677.  'SI' (ERRO > 1.0D-3) ;
  1678.      'ERREUR' 5 ;
  1679.  'FINSI' ;
  1680.  
  1681. * RUX
  1682.  
  1683.  DUXSCAL = EPS '*'  ('EXTRAIRE' GN 'UX' PCEN) ;
  1684.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1685.  
  1686.  RN1 = RN ;
  1687.  GN1 = GN '+' DUX ;
  1688.  RETN1 = RETN ;
  1689.  
  1690.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1691.  
  1692.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1693.                  LISTC LISTP
  1694.                  RN1 VN1 PN1 GAMN CHPLI 'INJE'  ;
  1695.  
  1696.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  1697.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1698.            'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1699.  
  1700.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1701.  
  1702.  'SI' (ERRO > 1.0D-3) ;
  1703.      'ERREUR' 5 ;
  1704.  'FINSI' ;
  1705.  
  1706. * RUY
  1707.  
  1708.  DUYSCAL = EPS '*'  ('EXTRAIRE' GN 'UY' PCEN) ;
  1709.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1710.  
  1711.  RN1 = RN ;
  1712.  GN1 = GN '+' DUY ;
  1713.  RETN1 = RETN ;
  1714.  
  1715.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1716.  
  1717.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1718.                  LISTC LISTP
  1719.                  RN1 VN1 PN1 GAMN CHPLI 'INJE'  ;
  1720.  
  1721.  DRESDUYN = (RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  1722.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1723.            'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1724.  
  1725.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1726.  
  1727.  'SI' (ERRO > 1.0D-3) ;
  1728.      'ERREUR' 5 ;
  1729.  'FINSI' ;
  1730.  
  1731.  
  1732. * RUZ
  1733.  
  1734.  DUZSCAL = EPS '*'  ('EXTRAIRE' GN 'UZ' PCEN) ;
  1735.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1736.  
  1737.  RN1 = RN ;
  1738.  GN1 = GN '+' DUZ ;
  1739.  RETN1 = RETN ;
  1740.  
  1741.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1742.  
  1743.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1744.                  LISTC LISTP
  1745.                  RN1 VN1 PN1 GAMN CHPLI 'INJE'  ;
  1746.  
  1747.  DRESDUZN = (RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  1748.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1749.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  1750.  
  1751.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1752.  
  1753.  'SI' (ERRO > 1.0D-3) ;
  1754.      'ERREUR' 5 ;
  1755.  'FINSI' ;
  1756.  
  1757. * RETN
  1758.  
  1759.  DRSCAL = EPS '*'  ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  1760.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1761.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1762.  
  1763.  RN1 = RN ;
  1764.  GN1 = GN ;
  1765.  RETN1 = RETN '+' DR ;
  1766.  
  1767.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1768.  
  1769.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1770.                  LISTC LISTP
  1771.                  RN1 VN1 PN1 GAMN CHPLI 'INJE'  ;
  1772.  
  1773.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  1774.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1775.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  1776.  
  1777.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1778.  
  1779.  'SI' (ERRO > 1.0D-3) ;
  1780.      'ERREUR' 5 ;
  1781.  'FINSI' ;
  1782.  
  1783.  
  1784. **********************************************************************************
  1785. **********************************************************************************
  1786. ***************************  INJELM **********************************************
  1787. **********************************************************************************
  1788. **********************************************************************************
  1789. *
  1790. * We impose RT and the normal momontum at the inlet.
  1791.  
  1792.  RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  1793.  PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  1794.  GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  1795.  VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 1.11 'UX' -1.3 'UZ'
  1796.    1.5 ;
  1797.  
  1798.  GN = RN '*' VN ;
  1799.  RETN = (('INVERSE' (GAMN '-' 1)) '*' PN)  '+' (0.5 '*' ('PSCAL' VN GN
  1800.    ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  1801.  
  1802.  MOME = 0.111 ;
  1803.  RT   = 4.1 ;
  1804.  
  1805.  CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 2 'RT' RT
  1806.     'MOME' MOME ;
  1807.  
  1808.  LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  1809.  LISTP = 'MOTS' 'RN  ' 'UX  ' 'UY  ' 'UZ  ' 'PN  ' ;
  1810.  
  1811.  RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1812.                  LISTC LISTP
  1813.                  RN VN PN GAMN CHPLI 'INJELM'  ;
  1814.  
  1815. *
  1816. **** Jacobian with respect to primitive variable
  1817. *
  1818.  
  1819.  
  1820.  RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1821.                  LISTC LISTP
  1822.                  RN VN PN GAMN CHPLI 'INJELM'  ;
  1823.  
  1824.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOPRIM' MDOMINT MENTREE
  1825.           LISTC LISTP
  1826.           RN VN PN GAMN CHPLI 'INJELM'  ;
  1827.  
  1828.  EPS=1.0D-6 ;
  1829.  
  1830. * RN
  1831.  
  1832.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1833.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1834.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1835.  
  1836.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1837.                  LISTC LISTP
  1838.                  (RN '+' DR) VN PN GAMN CHPLI 'INJELM'  ;
  1839.  
  1840.  DRESDRN = ( RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  1841.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  1842.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1843.  
  1844.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1845.  
  1846.  'SI' (ERRO > 1.0D-3) ;
  1847.      'ERREUR' 5 ;
  1848.  'FINSI' ;
  1849.  
  1850.  
  1851. * UX
  1852.  
  1853.  DUXSCAL = EPS '*'  ('EXTRAIRE' VN 'UX' PCEN) ;
  1854.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1855.  
  1856.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1857.                  LISTC LISTP
  1858.                  RN (VN '+' DUX) PN GAMN CHPLI 'INJELM'  ;
  1859.  
  1860.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  1861.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1862.            'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1863.  
  1864.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1865.  
  1866.  'SI' (ERRO > 1.0D-3) ;
  1867.      'ERREUR' 5 ;
  1868.  'FINSI' ;
  1869.  
  1870. * UY
  1871.  
  1872.  DUYSCAL = EPS '*'  ('EXTRAIRE' VN 'UY' PCEN) ;
  1873.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1874.  
  1875.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1876.                  LISTC LISTP
  1877.                  RN (VN '+' DUY) PN GAMN CHPLI 'INJELM'  ;
  1878.  
  1879.  DRESDUYN = ( RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  1880.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1881.            'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  1882.  
  1883.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1884.  
  1885.  'SI' (ERRO > 1.0D-3) ;
  1886.      'ERREUR' 5 ;
  1887.  'FINSI' ;
  1888.  
  1889. * UZ
  1890.  
  1891.  DUZSCAL = EPS '*'  ('EXTRAIRE' VN 'UZ' PCEN) ;
  1892.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1893.  
  1894.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1895.                  LISTC LISTP
  1896.                  RN (VN '+' DUZ) PN GAMN CHPLI 'INJELM'  ;
  1897.  
  1898.  DRESDUZN = ( RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  1899.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1900.            'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  1901.  
  1902.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1903.  
  1904.  'SI' (ERRO > 1.0D-3) ;
  1905.      'ERREUR' 5 ;
  1906.  'FINSI' ;
  1907.  
  1908. * PN
  1909.  
  1910.  DPSCAL = EPS '*'  ('EXTRAIRE' PN 'SCAL' PCEN) ;
  1911.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1912.  DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  1913.  
  1914.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1915.                  LISTC LISTP
  1916.                  RN  VN (PN '+' DP) GAMN CHPLI 'INJELM'  ;
  1917.  
  1918.  DRESDPN = ( RCHRES1 '-'  RCHRES) '/' DPSCAL ;
  1919.  DRESDP  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1920.            'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  1921.  
  1922.  ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  1923.  
  1924.  'SI' (ERRO > 1.0D-3) ;
  1925.      'ERREUR' 5 ;
  1926.  'FINSI' ;
  1927.  
  1928. *
  1929. **** Jacobian with respect to conservative variables
  1930. *
  1931.  
  1932.  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOCONS' MDOMINT MENTREE
  1933.           LISTC LISTP
  1934.           RN VN PN GAMN CHPLI 'INJELM'  ;
  1935.  
  1936.  EPS=1.0D-6 ;
  1937.  
  1938. * RN
  1939.  
  1940.  DRSCAL = EPS '*'  ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1941.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  1942.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1943.  
  1944.  RN1 = RN '+' DR ;
  1945.  GN1 = GN ;
  1946.  RETN1 = RETN ;
  1947.  
  1948.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1949.  
  1950.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1951.                  LISTC LISTP
  1952.                  RN1 VN1 PN1 GAMN CHPLI 'INJELM'  ;
  1953.  
  1954.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  1955.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  1956.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1957.  
  1958.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1959.  
  1960.  'SI' (ERRO > 1.0D-3) ;
  1961.      'ERREUR' 5 ;
  1962.  'FINSI' ;
  1963.  
  1964. * RUX
  1965.  
  1966.  DUXSCAL = EPS '*'  ('EXTRAIRE' GN 'UX' PCEN) ;
  1967.  DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1968.  
  1969.  RN1 = RN ;
  1970.  GN1 = GN '+' DUX ;
  1971.  RETN1 = RETN ;
  1972.  
  1973.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1974.  
  1975.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  1976.                  LISTC LISTP
  1977.                  RN1 VN1 PN1 GAMN CHPLI 'INJELM'  ;
  1978.  
  1979.  DRESDUXN = ( RCHRES1 '-'  RCHRES) '/' DUXSCAL ;
  1980.  DRESDUX  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1981.            'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1982.  
  1983.  ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1984.  
  1985.  'SI' (ERRO > 1.0D-3) ;
  1986.      'ERREUR' 5 ;
  1987.  'FINSI' ;
  1988.  
  1989. * RUY
  1990.  
  1991.  DUYSCAL = EPS '*'  ('EXTRAIRE' GN 'UY' PCEN) ;
  1992.  DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1993.  
  1994.  RN1 = RN ;
  1995.  GN1 = GN '+' DUY ;
  1996.  RETN1 = RETN ;
  1997.  
  1998.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1999.  
  2000.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  2001.                  LISTC LISTP
  2002.                  RN1 VN1 PN1 GAMN CHPLI 'INJELM'  ;
  2003.  
  2004.  DRESDUYN = (RCHRES1 '-'  RCHRES) '/' DUYSCAL ;
  2005.  DRESDUY  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  2006.            'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  2007.  
  2008.  ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  2009.  
  2010.  'SI' (ERRO > 1.0D-3) ;
  2011.      'ERREUR' 5 ;
  2012.  'FINSI' ;
  2013.  
  2014.  
  2015. * RUZ
  2016.  
  2017.  DUZSCAL = EPS '*'  ('EXTRAIRE' GN 'UZ' PCEN) ;
  2018.  DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  2019.  
  2020.  RN1 = RN ;
  2021.  GN1 = GN '+' DUZ ;
  2022.  RETN1 = RETN ;
  2023.  
  2024.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  2025.  
  2026.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  2027.                  LISTC LISTP
  2028.                  RN1 VN1 PN1 GAMN CHPLI 'INJELM'  ;
  2029.  
  2030.  DRESDUZN = (RCHRES1 '-'  RCHRES) '/' DUZSCAL ;
  2031.  DRESDUZ  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  2032.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  2033.  
  2034.  ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  2035.  
  2036.  'SI' (ERRO > 1.0D-3) ;
  2037.      'ERREUR' 5 ;
  2038.  'FINSI' ;
  2039.  
  2040. * RETN
  2041.  
  2042.  DRSCAL = EPS '*'  ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  2043.  ELE1 = 'MANUEL' 'POI1' PCEN ;
  2044.  DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  2045.  
  2046.  RN1 = RN ;
  2047.  GN1 = GN ;
  2048.  RETN1 = RETN '+' DR ;
  2049.  
  2050.  VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  2051.  
  2052.  RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'RESI' MDOMINT MENTREE
  2053.                  LISTC LISTP
  2054.                  RN1 VN1 PN1 GAMN CHPLI 'INJELM'  ;
  2055.  
  2056.  DRESDRN = (RCHRES1 '-'  RCHRES) '/' DRSCAL ;
  2057.  DRESDR  = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  2058.            'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  2059.  
  2060.  ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2061.  
  2062.  'SI' (ERRO > 1.0D-3) ;
  2063.      'ERREUR' 5 ;
  2064.  'FINSI' ;
  2065.  
  2066.  'FIN' ;
  2067.  
  2068.  
  2069.  
  2070.  
  2071.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales