plot two or more optical spectrum

[mbjavan]

please help me to plot  optical spectrums, calculated and saved in some *.nc files at a single figure.
I need a script 

[zh]

Please check here:
http://www.quantumwise.com/documents/tutorials/latest/SiliconOptical/index.html/chap.Si.Band.html#sect1.running.si

[marmotte]

Dear ATK team,

I'm trying to get the absorption for my system. I've already calculated the dielectric constants and it is OK. My problem is when I apply the script of : http://www.quantumwise.com/documents/tutorials/latest/SiliconOptical/index.html/chap.Si.Band.html#sect1.running.si, I have RuntimeWarning: divide by zero encountered in divide.

I guess I have to change values in this script and this values must be extracted from optical spectrum. Could you please tell me what is the meaning

of the values in red, I mean what the corresponding values in the dielectric constant spectrum

ax.axis([180,1000,2.2,6.4])
ax.set_ylabel(r"$n$", size=16)
ax.tick_params(axis='x', labelbottom=False, labeltop=True)
ax = pylab.subplot(212)
ax.plot(l,alpha,'r')
ax.axis([180,1000,0,0.24])


Thx

[Anders Blom]

Where is the warning? Probably from the line dividing by energy. If so, it's not so critical, it just means you have a zero energy.

180 and 1000 are wavelengths, computed from the energies on the line commented as "# Calculate the wavelength" and 2.2 and 6.4 and the other ones are needed because axis() take a 4-element list (xmin, xmax, ymin, ymax).

[marmotte]

Dear Anders,

Actually it is critical because I can't have the absorption ! how we can avoid the problem. all the optical properties are extracted from dielectric constant and this later is correct. I mean how to avoid this problem.

Best regards,

[marmotte]

Dear Ander,

Now it's working, please avoid my last post.

best