# Expand phase function in Legendre polynomials
for i in 4 8 16 32 64 128 3000; do pmom -l $i -n phase_function_500nm_reff10.dat > nmom$i.dat; done

# Calulate phase function from Legendre polynomials
for i in 4 8 16 32 64 128 3000; do phase -c -d -s 0.1 nmom$i.dat > phase$i.dat; done

# Calulate phase function from Legendre polynomials with delta scaling
for i in 4 8 16 32 64 128 3000; do phase -c -d -s 0.1 -f nmom$i.dat > phase$i\_ds.dat; done