Télécharger rovap.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : rovap.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ***************** CAS TEST : rovap.dgibi **************************
  5. *
  6. * Test de l'opérateur VARI ROVAP(P,T)
  7. * Les données sont un FLOTTANT, un LISTREEL ou un CHPO
  8. *
  9. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  10. 'OPTI' 'ECHO' 0 ;
  11. *
  12. *------------------------------------------------------> FLOTTANT
  13. *
  14. X1 = 1.0133D5 ;
  15. X2 = 373.15D0 ;
  16. X3 = VARI ROVAP X1 X2 ;
  17. XSTO1 = X3 ;
  18. *
  19. X1 = 5.021D5 ;
  20. X2 = 425.15D0 ;
  21. X3 = VARI ROVAP X1 X2 ;
  22. XSTO2 = X3 ;
  23. *
  24. X1 = 10.027D5 ;
  25. X2 = 453.15D0 ;
  26. X3 = VARI ROVAP x1 x2 ;
  27. XSTO3 = X3 ;
  28. *------------------------------------------------------> LISTREEL
  29. Y1 = 'PROG' 1.0133d5 5.021d5 10.027d5 ;
  30. Y2 = 'PROG' 373.15D0 425.15D0 453.15D0 ;
  31. Y3 = VARI ROVAP Y1 Y2 ;
  32. YSTO1 = EXTR Y3 1 ;
  33. YSTO2 = EXTR Y3 2 ;
  34. YSTO3 = EXTR Y3 3 ;
  35. *------------------------------------------------------> CHPO
  36. P1 = 1.0133 373.15 ;
  37. P2 = 5.021 425.15 ;
  38. P3 = 10.027 453.15 ;
  39. P1P3 = P1 'DROI' 1 P2 'DROI' 1 P3 ;
  40. Z1 = 1.D5 * ('COOR' 1 P1P3) ;
  41. Z2 = 1. * ('COOR' 2 P1P3) ;
  42. Z3 = VARI ROVAP Z1 Z2 ;
  43. ZSTO1 = 'EXTR' Z3 'SCAL' P1 ;
  44. ZSTO2 = 'EXTR' Z3 'SCAL' P2 ;
  45. ZSTO3 = 'EXTR' Z3 'SCAL' P3 ;
  46. *-------------------------------------------> Controle
  47. CTRL1 = XSTO1 + YSTO1 + ZSTO1 / 3.d0 - XSTO1 / XSTO1 ;
  48. CTRL2 = XSTO2 + YSTO2 + ZSTO2 / 3.d0 - XSTO2 / XSTO2 ;
  49. CTRL3 = XSTO3 + YSTO3 + ZSTO3 / 3.d0 - XSTO3 / XSTO3 ;
  50. CTRL4 = CTRL1 + CTRL2 + CTRL3 ;
  51. XREF1 = 0.597728631 ;
  52. XREF2 = 2.679528403 ;
  53. XREF3 = 5.15995872 ;
  54. CTRL5 = XSTO1 + YSTO1 + ZSTO1 / 3.d0 - XREF1 / XREF1 ;
  55. CTRL6 = XSTO2 + YSTO2 + ZSTO2 / 3.d0 - XREF2 / XREF2 ;
  56. CTRL7 = XSTO3 + YSTO3 + ZSTO3 / 3.d0 - XREF3 / XREF3 ;
  57. CTRL8 = CTRL5 + CTRL6 + CTRL7 ;
  58. *-------------------------------------------> Affichage
  59. 'MESS' ' ' ;
  60. 'MESS' ' ' ;
  61. 'MESS' 'Pt 1 VARI ROVAP ----->' XSTO1 YSTO1 ZSTO1 XREF1 ;
  62. 'MESS' 'Pt 2 VARI ROVAP ----->' XSTO2 YSTO2 ZSTO2 XREF2 ;
  63. 'MESS' 'Pt 3 VARI ROVAP ----->' XSTO3 YSTO3 ZSTO3 XREF3 ;
  64. 'MESS' ' ' ;
  65. 'MESS' 'Comparaison calculs rel.----->' CTRL4 ;
  66. 'MESS' 'Erreur rel. VDI ----->' CTRL8 ;
  67. 'MESS' ' ' ;
  68. 'MESS' ' ' ;
  69. 'MESS' ' ' ;
  70. 'MESS' ' ' ;
  71. 'MESS' ' ' ;
  72. *-------------------------------------------> Compte-rendu et sortie
  73. * EPS8 : Incertitude relative sur les data vdi
  74. EPS4 = 1.E-14 ;
  75. EPS8 = 6.E-4 ;
  76. CTRL4 = ABS CTRL4 ;
  77. CTRL8 = ABS CTRL8 ;
  78. LOG4 = CTRL4 > EPS4 ;
  79. LOG8 = CTRL8 > EPS8 ;
  80. L0 = LOG4 'OU' LOG8 ;
  81. 'SI' L0 ;
  82. 'ERREUR' 5 ;
  83. 'SINON' ;
  84. 'ERREUR' 0 ;
  85. 'FINSI' ;
  86. 'FIN' ;
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  

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