from matplotlib.font_manager import FontProperties
font = FontProperties(fname='/WORK/grmc_dxj_1/wangnan/data/font/arial.ttf')

fig =plt.figure(num=1,figsize=(18,5))   ###圖??陸?fig
plt.subplot(121)
plt.fill_between(range(len(p_var3)),p_var3_90,p_var3_10,facecolor='lightcoral',alpha=0.3)
plt.plot(range(len(p_var3)),p_var3,color='indianred',linestyle='-',linewidth=2.,label='PRD (cut_3day)')
plt.fill_between(range(len(p_var3)),p_var1_90,p_var1_10,facecolor='cornflowerblue',alpha=0.3)
plt.plot(p_var1,color='royalblue',linewidth=2.,label='PRD (cut_1day)')
plt.title('PRD',fontsize=26)
plt.text(0.,-18,'3-day advance control',c='indianred',fontsize=24)
plt.text(0.,-20.5,'1-day advance control',c='royalblue',fontsize=24)
plt.yticks( fontsize = 24)  #xy軸刻度數(shù)字
plt.xticks( fontsize = 24)
plt.ylabel(r'$\mathrm{ \Delta O_{3} \  (ppb)}$',fontsize=26)
plt.axhline(y=0,ls=":",c="k")
plt.xticks([0,4,8,12],['08:00','12:00','16:00','20:00'])
plt.ylim([-22,5])
plt.yticks([-15,0],[-15,0])

plt.subplot(122)
plt.fill_between(range(len(y_var3)),y_var3_90,y_var3_10,facecolor='lightcoral',alpha=0.3)
plt.plot(range(len(y_var3)),y_var3,color='indianred',linestyle='-',linewidth=2.)
plt.fill_between(range(len(y_var3)),y_var1_90,y_var1_10,facecolor='cornflowerblue',alpha=0.3)
plt.plot(y_var1,color='royalblue',linewidth=2.)
plt.ylim([-40,8])
plt.title('YRD',fontsize=26)

plt.yticks( size = 24)  #xy軸刻度數(shù)字
plt.xticks( size = 24)
plt.xticks([0,4,8,12],['08:00','12:00','16:00','20:00'])
plt.yticks([-30,-15,0],[-30,-15,0])
plt.ylabel(r'$\mathrm{ \Delta O_{3} \  (ppb)}$',fontsize=26)

plt.axhline(y=0,ls=":",c="k")

plt.savefig('./FigS_NOx_reduction_pannel.png',bbox_inches='tight' , dpi=400)
plt.show()

image.png