def radial_flare(script, flare_radius=None, start_radius=None, end_radius=None,<EOL>end_height=None):
<EOL>effective_radius = '<STR_LIT>'<EOL>r_func = '<STR_LIT>'<EOL>z_func = '<STR_LIT>'<EOL>r_func = r_func.replace('<STR_LIT>', str(effective_radius)).replace('<STR_LIT>', str(start_radius)).replace('<STR_LIT>', str(flare_radius))<EOL>z_func = z_func.replace('<STR_LIT>', str(effective_radius)).replace('<STR_LIT>', str(start_radius)).replace('<STR_LIT>', str(flare_radius))<EOL>function_cyl_co(script=script, r_func=r_func, z_func=z_func)<EOL>return None<EOL>
flare_radius must be >= z2 (height)
r2 max = flare_radius + r

r2 (num): radius of mesh at end of flare

+15 r= 8.8205
-15 r= 1.1795

z=10, 5 +/-15 - +/-15*0.74535599249992989880305788957709