* fichier : rovapsat.dgibi ************************************************************************ * Section : Mathematiques Fonctions ************************************************************************ ***************** CAS TEST : rovapsat.dgibi ************************** * * Test de l'opérateur VARI ROVAP(P,TSAT(P)) * Les données sont un FLOTTANT, un LISTREEL ou un CHPO * 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ; 'OPTI' 'ECHO' 0 ; * *------------------------------------------------------> FLOTTANT * X1 = 1.0133D5 ; X3 = VARI ROVAP X1 ; XSTO1 = X3 ; * X1 = 5.021D5 ; X3 = VARI ROVAP X1 ; XSTO2 = X3 ; * X1 = 10.027D5 ; X3 = VARI ROVAP X1 ; XSTO3 = X3 ; *------------------------------------------------------> LISTREEL Y1 = 'PROG' 1.0133d5 5.021d5 10.027d5 ; Y3 = VARI ROVAP Y1 ; YSTO1 = EXTR Y3 1 ; YSTO2 = EXTR Y3 2 ; YSTO3 = EXTR Y3 3 ; *------------------------------------------------------> CHPO P1 = 1.0133 0. ; P2 = 5.021 0. ; P3 = 10.027 0. ; P1P3 = P1 'DROI' 1 P2 'DROI' 1 P3 ; Z1 = 1.D5 * ('COOR' 1 P1P3) ; Z3 = VARI ROVAP Z1 ; ZSTO1 = 'EXTR' Z3 'SCAL' P1 ; ZSTO2 = 'EXTR' Z3 'SCAL' P2 ; ZSTO3 = 'EXTR' Z3 'SCAL' P3 ; *-------------------------------------------> Controle CTRL1 = XSTO1 + YSTO1 + ZSTO1 / 3.d0 - XSTO1 / XSTO1 ; CTRL2 = XSTO2 + YSTO2 + ZSTO2 / 3.d0 - XSTO2 / XSTO2 ; CTRL3 = XSTO3 + YSTO3 + ZSTO3 / 3.d0 - XSTO3 / XSTO3 ; CTRL4 = CTRL1 + CTRL2 + CTRL3 ; XREF1 = 0.597728631 ; XREF2 = 2.679528403 ; XREF3 = 5.15995872 ; CTRL5 = XSTO1 + YSTO1 + ZSTO1 / 3.d0 - XREF1 / XREF1 ; CTRL6 = XSTO2 + YSTO2 + ZSTO2 / 3.d0 - XREF2 / XREF2 ; CTRL7 = XSTO3 + YSTO3 + ZSTO3 / 3.d0 - XREF3 / XREF3 ; CTRL8 = CTRL5 + CTRL6 + CTRL7 ; *-------------------------------------------> Affichage 'MESS' ' ' ; 'MESS' ' ' ; 'MESS' 'Pt 1 VARI ROVAPsat -->' XSTO1 YSTO1 ZSTO1 XREF1 ; 'MESS' 'Pt 2 VARI ROVAPsat -->' XSTO2 YSTO2 ZSTO2 XREF2 ; 'MESS' 'Pt 3 VARI ROVAPsat -->' XSTO3 YSTO3 ZSTO3 XREF3 ; 'MESS' ' ' ; 'MESS' 'Comparaison calculs rel.----->' CTRL4 ; 'MESS' 'Erreur rel. VDI ----->' CTRL8 ; 'MESS' ' ' ; 'MESS' ' ' ; 'MESS' ' ' ; 'MESS' ' ' ; 'MESS' ' ' ; *-------------------------------------------> Compte-rendu et sortie * EPS8 : Incertitude relative sur les data vdi EPS4 = 1.E-14 ; EPS8 = 6.E-4 ; CTRL4 = ABS CTRL4 ; CTRL8 = ABS CTRL8 ; LOG4 = CTRL4 > EPS4 ; LOG8 = CTRL8 > EPS8 ; L0 = LOG4 'OU' LOG8 ; 'SI' L0 ; 'ERREUR' 5 ; 'SINON' ; 'ERREUR' 0 ; 'FINSI' ; 'FIN' ;