qbasic代码

rho1=(rg*tsat(i )/g) ^-1 : pmean=(p(i+1)+p(i))/2
bigx(i)=2*muu(i) /(pmean*di*SQR( rho1* g))
mtot=ABS(mtot+ delmdot(i))
qsum=qsum-delq(i)
IF ABS(mtot-mdot(1)) <= (mdot(1)*.0001#) OR mtot > mdot(1) THEN PRINT" All mass removed. " :GOTO endroutine
IF tsat(i)<=tsink THEN PRINT" T at or below sink temperature ",tsat(i):GOTO endroutine
IF biqx(i)<slip THEN GOTO nextelem
IF blgx(i) >= slip AND bigx(i) < molec THEN GOTO 100
IF bigx(i) >= molec AND iflag <2 THEN PRINT"MOLECULAR FLOW ",biqx(i): iflag=2!
GOTO nextelem
100 : IF iflag< 1 THEN PRINT"SLIP FLOW ",bigx(i) : PRINT" @ z-",z(i)
iflag=1!
v2=mach(i+1) ^2*k*rg*tsat(i+1)
f=(dout/2#/rho(i+1)/v2)*(p(i)-p(i+1))/delz
PRINT" f @ i = "; f,i
nextelem:
nn=i
ncount=0
NEXT i
nn=i-1
endroutine:
PRINT" DONE AND DONE AT i =",i:PRINT"Qsonic = ", q(1)," Qtot = ', qsum
INPUT" hit RETURN ";DUMMY