*SYNTAXE_C3 : This means this input deck is written in C3 syntaxe
read result 21 ;
qbreak = chrono break_side qtotgas 0.001 1.d6 ;
tbreak = chrono break_side gastemp 0.001 1.d6 ;
tbreak2 = chrono break_side gastemp 0.00 1.d6 ;
Pbreak = chrono break_side pressure 0.00 1.d6 ;
vbreak = chrono break_side gasv 0.001 1.d6 ;
pland  = chrono land_side  pressure 0.001 1.d6 ;
mass_in = chrono  nspipe gasmass ;
diame = 1.153 ;
grav = 9.81 ;
liq_dens = 1000.0 ;
Pi = getpi ;
Section = Pi * 0.25 *  diame * diame ;
mass_flux = qbreak div Section ;
dens = mass_flux div vbreak ;
tmp1 = dens mult -1. ;
tmp2 = liq_dens add tmp1 ;
tmp3 = tmp2 div dens ;
tmp4 = tmp3 div diame ;
tmp5 = tmp4 div grav ;
tmp6 = vbreak div tmp5  ;
froude = vbreak mult tmp6 ;
FRVSQ = froude versus qbreak ;
FRVSV = froude versus vbreak ;
PT = Pbreak versus tbreak2 ; 
gnuplot qbreak tbreak vbreak pland mass_in froude FRVSQ FRVSV PT ;