from pylab import *

i=5
l=4
mie_r1=loadtxt('/home/users/lainer/Strahlung/Mie/tmp_mie/mie1.out')
semilogx(mie_r1[:,0], mie_r1[:,i], label='reff=1nm') # change 1 to desired variable


title('Ice particles - variation of asymmetry parameter g')
xlabel('wavelength [nm]')
ylabel('g')
legend(loc=l)


mie_r5=loadtxt('/home/users/lainer/Strahlung/Mie/tmp_mie/mie5.out')
semilogx(mie_r5[:,0], mie_r5[:,i], label='reff=5nm') # change 1 to desired variable
#ylim(0,1.1)
legend(loc=l)


mie_r10=loadtxt('/home/users/lainer/Strahlung/Mie/tmp_mie/mie10.out')
semilogx(mie_r10[:,0], mie_r10[:,i], label='reff=10nm') # change 1 to desired variable
legend(loc=l)



mie_r15=loadtxt('/home/users/lainer/Strahlung/Mie/tmp_mie/mie15.out')
semilogx(mie_r15[:,0], mie_r15[:,i], label='reff=15nm') # change 1 to desired variable
legend(loc=l)



mie_r20=loadtxt('/home/users/lainer/Strahlung/Mie/tmp_mie/mie20.out')
semilogx(mie_r20[:,0], mie_r20[:,i], label='reff=20nm') # change 1 to desired variable
legend(loc=l)


savefig('mie_ssa.png')
"mie.py" 41L, 1031C