convergence order 2 geo - postplot

DiffractionGratings/grating3D_postplot.py

@@ -48,7 +48,10 @@ if myDir[6:]=='solarcell':
 elif myDir[6:]=='conv':
     print('cf pdf')
     pl.figure()
-    pl.plot(np.loadtxt(myDir+'/temp-Q_scat.txt',ndmin=2)[:,1])
+    data_abs = np.loadtxt(myDir+'/temp-Q_scat.txt',ndmin=2)[:,1]
+    pl.plot(data_abs[0:int(len(data_abs)/2)],label='linear elements')
+    pl.plot(data_abs[int(len(data_abs)/2):] ,label='curved elements')
+    pl.legend()
     pl.savefig('fig_convergence_absorption.pdf')
 else:
     print('===> Computed from diffraction efficiencies')