18
if(NSubA[NNit]=0) and (VarRoA[NNit]=l) then
RoM:=sqr(0.5*(sqrt(RoH)+sqrt(RoL)));
NRQ:=1;
QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,
QC,VCV,VC,VM,RO,KAPPA,MU,HS);
C(BET,RE,DPH,KSH,KK, ERCJERR);
if IERR = 0 then begin DP(EDPA,DPM,EDPS,EDPR,EDP);
P(EPA,PMG,EPS, EPR, EP);
T(ETA, TMG, ETS, ETR, ET);
ROS(EROSA, RIM, EROSS, EROSR, EROS);
EPSI(BET.DPM,PMG,EEPS);
EPH(RIM,TMG,PMG,ROSC, ERO,EKAP,EVIS);
if (NSubA[NNIT]=0 and (VarRoA[NNit]=l) then begin
ROW(EROWA, EROWS, EROWR, EROW); ERO := EROWS
end
else EROWR := 0.0;
NRQ:=0;
TETDP:=0.0;
ifEDP <> 0.0 then begin
DPMD:=0.001*EDP*DPM; DPM1:=DPM+DPMD; DPM2:=DPM-DPMD;
QCalc(NRQ,TMG,PMG,DPMl,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,
QC1,VCV1,VC1,VM1,RO,KAPPA,MU,HS);
QCalc(NRQ,TMG,PMG,DPM2,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,
QC2,VCV2,VC2,VM2,RO.KAPPA,MU,HS);
TETDP:=DPM*(QC1-QC2)/(2.0*DPMD)/QC;
end:
DDD:=0.001*EDD*DD20: DD01:=DD20+DDD; DD02:=DD20-DDD;
QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD01,DT20,RIM,BET,RE,KSH,KK,ROSC,
QC1,VCV1,VC1,VM1,RO,KAPPA.MU,HS):
QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD02,DT20,RIM,BET,R£,KSH,KK,ROSC,
QC2,VCV2,VC2,VM2,RO,KAPPA,MU,HS);
TETDD:=DD20*(QC1-QC2)/(2.0*DDD)/QC;
DDT:=0.001*EDT*DT20; DT01:=DT20+DDT; DT02:=DT20-DDT;
QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT01,RIM,BET,RE,KSH,KK,ROSC,
QC 1,VCV1,VC1,VM1,RO, KAPPA, MU,HS);
QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT02,RIM,BET,RE,KSH,KK,ROSC,
QC2,VCV2,VC2,VM2,RO,KAPPA,MU,HS);
TETDT:=DT20*(QC1-QC2)/(2.0*DDT)/QC;
DKAP:=0.001*EKAP* KAPPA; KAP1:=KAPPA+DKAP; KAP2:=KAPPA-DKAP;
QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,
QC1,VCV1,VC1,VM1,RO,KAP1,MU,HS);
QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,
QC2,VCV2,VC2,VM2,RO,KAP2,MU,HS);
TETKAP:=KAPPA*(QC1-QC2)/(2.0*DKAP)/QC;
DMU:=0.001*EVIS*MU; MU1:=MU+DMU; MU2:=MU-DMU;
QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,
QC1,VCV1,VC1,VM1,RO,KAPPA,MU1,HS);
QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,
QC2,VCV2,VC2,VM2,RO,KAPPA,MU2,HS);
TETMU:=MU*(QC1-QC2)/(2.0*DMU)/QC;
TETC:=1.0; TETEPS:=1.0; TETRO:=0.5;
if (NSubA[NNit]=0) and (VarRoA[NNit]=l) then begin