ElectronConfinementTimeFirstApprox.pl

#!/bin/perl

do "parameters.pl";



$MeanUloop=`cat MeanUloop`;
$MeanPlasmaCurrent=`cat MeanPlasmaCurrent`;
$StateEqElectronDensity=`cat StateEqElectronDensity`;
$ElectronTemperature=`cat ElectronTemperature`;

#Ref: PhD Brotankova eq. 3.27 p. 27 (electron density modification)
open RESULT, ">ElectronConfinementTimeFirstApprox";
printf RESULT "%4.1f",3/8*$PlasmaVolume*($ElectronTemperature*$eV*$StateEqElectronDensity)/($MeanPlasmaCurrent*$MeanUloop);
close RESULT;