check paper test cases

DiffractionGratings/grating3D_parallel_Mmatrix.sh

@@ -4,20 +4,22 @@ export OMP_NUM_THREADS=1
 export NPROC=96
 
 ### choose configuration
-### theta loop study for the pillar case
+### theta loop study for the "U" case
+### config 1
 export flag_angle_study="theta"
 export loop_angle_max=50
 export fixed_angle=90
-export GRATING_CASE="pillar"
+export GRATING_CASE="U"
 
-### phi loop study for the half_ellipsoid case
+### phi loop study for the "U" case
+### config 2
 # export flag_angle_study="phi"
 # export loop_angle_max=360
 # export fixed_angle=50
-# export GRATING_CASE="half_ellipsoid"
+# export GRATING_CASE="U"
 
-export nb_angle=20
-export nb_lam=100
+export nb_angle=40
+export nb_lam=50
 export lambda_min=400
 export lambda_max=1200
 export FLAG_TOTAL=0

DiffractionGratings/grating3D_postplot_Mmatrix.py

@@ -164,7 +164,7 @@ ax=axes[2,1];zplot = ax.pcolormesh(L,P,Qscat_pin.real);add_colorbar(zplot);ax.ti
 plt.savefig('BALANCE_new_code.jpg')
 
 fig, axes = plt.subplots(4, 4, subplot_kw=dict(projection='polar') ,figsize=(12,12))
-flag_lam = False
+flag_lam = True
 rlabel = r"$\hbar \nu$" if not flag_lam else r"$\lambda_0$"
 anglelabel = r"$\varphi_0$" if flag_angle_study=='phi' else r"$\theta_0$"
 which_r = tab_hnu if not flag_lam else tab_lam
@@ -173,8 +173,8 @@ for i in range(4):
     for j in range(4):
         ax=axes[i,j]
         if i!=0 or j!=0:
-            sm=ax.contourf(tab_angle*np.pi/180,which_r,M[i,j]/M[0,0],cmap=plt.cm.bwr,levels=30,vmin=-1,vmax=1)
-            # sm=ax.contourf(tab_angle*np.pi/180,which_r,M[i,j]/M[0,0],cmap=plt.cm.bwr,levels=30)
+            # sm=ax.contourf(tab_angle*np.pi/180,which_r,M[i,j]/M[0,0],cmap=plt.cm.bwr,levels=30,vmin=-1,vmax=1)
+            sm=ax.contourf(tab_angle*np.pi/180,which_r,M[i,j]/M[0,0],cmap=plt.cm.bwr,levels=30)
             ax.text(0. , 1. , r"$M_{%d%d}$"%(i+1,j+1) , fontsize=16 , transform=ax.transAxes)
             ax.set_xticks([])
             ax.set_yticks([])