Hello all, I have produced an animation of a fireball using the comet3 function (see code and question below)
kt = 'Please enter the yield in kilotonnes you wish to test: ';
kt = input(kt);
Fireball = 'Your fireball volume for your chosen yield: ';
if (kt<=1);
fireball_radius1 = 60;
elseif(kt<=2>=1);
fireball_radius2 = 80;
elseif(kt<=3>=2);
fireball_radius3 = 90;
end
z = linspace(-1,1,5000);
r = sqrt(1-z.^2);
t = linspace(0,150*pi,5000);
x = r.*cos(t);
y = r.*sin(t);
comet3(x,y,z)
I want the fireball radius to be 60, 80 or 90 dependant on which section the yield the user inputs falls into. How do I re-call the radius values? So for example if the user inputs a value of 1.5Kt the radius of the fireball needs to be 80Km. How is this recalled to plot that value. Thanks
