! ! Local parameters ! ! maxstep : Over how many 'shells' do you want to calculate the SF's? ! max_moment : How many moments do you want to calculate? ! integer, parameter :: maxmom=6,maxstep=2,minstep=0