def _scale_prof(step, rprof, rad=None):
rbot, rtop = misc.get_rbounds(step)<EOL>if rbot == <NUM_LIT:0>:  <EOL><INDENT>return rprof<EOL><DEDENT>if rad is None:<EOL><INDENT>rad = step.rprof['<STR_LIT:r>'].values + rbot<EOL><DEDENT>return rprof * (<NUM_LIT:2> * rad / (rtop + rbot))**<NUM_LIT:2><EOL>
Scale profile to take sphericity into account