plt.grid()
plt.xticks(fontsize=ff)
plt.yticks(fontsize=ff)
plt.show()'''
from pcc import get_my_cmap2
xbins = 10**np.arange(-1, 2.08, 0.08)
ybins = 10**np.arange(-1, 2.08, 0.08)
counts, _, _ = np.histogram2d(A[maske_p], B[maske_p], bins=(xbins, ybins))
fig, ax = plt.subplots(figsize=(10, 10))
plt.pcolormesh(xbins, ybins, counts.T, cmap=get_my_cmap2(), vmin=0.1)
cb = plt.colorbar()
cb.set_label('number of samples', fontsize=ff)
ax.set_xscale('log')
ax.set_yscale('log')
plt.xlabel('DPR Rainrate in mm/h', fontsize=ff)
plt.ylabel('Radolan Rainrate in mm/h', fontsize=ff)
plt.xticks(fontsize=ff)
plt.yticks(fontsize=ff)
plt.title('Corr: ' + str(round(corr, 3)) + r'$\pm$' + str(round(eror, 3)) +
'\nBias: ' + str(round(bias, 3)),
fontsize=20)
plt.grid()
plt.xticks(fontsize=ff)
plt.yticks(fontsize=ff)
plt.title('Solid Corr: ' + str(round(r_corr2,3)) + r'$\pm$'+ str(round(r_eror2,3)), fontsize=20)
plt.grid()
#plt.savefig('/automount/ags/velibor/plot/IRS/NSsolidWS_RR.png')
#plt.close()
plt.show()
'''
from pcc import get_my_cmap2
fig = plt.figure(figsize=(8, 8))
ax2 = fig.add_subplot(111, aspect='equal')
plt.hist2d(C[maske_CD], D[maske_CD], bins=60, cmap=get_my_cmap2(), vmin=0.1)
cbar = plt.colorbar(shrink=0.7)
cbar.set_label('number of samples', fontsize=ff2)
cbar.ax.tick_params(labelsize=ff)
plt.title('Solid Phase Corr: ' + str(round(z_corr2, 3)) + r'$\pm$' +
str(round(z_eror2, 3)),
fontsize=ff2)
plt.xlim(15, 70)
plt.ylim(15, 70)
cx, cy = np.arange(0, 80, 1), np.arange(0, 80, 1)
plt.plot(cx, cy, color='black')
plt.xlabel('DPR Reflectivity in dBZ', fontsize=ff2)
plt.ylabel('Radolan Reflectivity in dBZ', fontsize=ff2)
plt.grid()
plt.xticks(fontsize=ff)
plt.yticks(fontsize=ff)
plt.title('Solid Corr: ' + str(round(r_corr2,3)) + r'$\pm$'+ str(round(r_eror2,3)), fontsize=20)
plt.grid()
#plt.savefig('/automount/ags/velibor/plot/IRS/NSsolidWS_RR.png')
#plt.close()
plt.show()
'''
from pcc import get_my_cmap2
fig =plt.figure(figsize=(8,8))
ax2 = fig.add_subplot(111, aspect='equal')
plt.hist2d(C[maske_CD],D[maske_CD], bins=60, cmap=get_my_cmap2(),vmin=0.1)
cbar = plt.colorbar(shrink=0.7)
cbar.set_label('number of samples', fontsize=ff2)
cbar.ax.tick_params(labelsize=ff)
plt.title('Solid Phase Corr: ' + str(round(z_corr2,3)) + r'$\pm$'+ str(round(z_eror2,3)), fontsize=ff2)
plt.xlim(15,70)
plt.ylim(15,70)
cx,cy = np.arange(0,80,1),np.arange(0,80,1)
plt.plot(cx,cy, color='black')
plt.xlabel('DPR Reflectivity in dBZ',fontsize=ff2)
plt.ylabel('Radolan Reflectivity in dBZ',fontsize=ff2)
plt.grid()
plt.xticks(fontsize=ff)
plt.yticks(fontsize=ff)
ax1.scatter(c,d,color='black',alpha=0.3, label='solid: '+ corcor(c,d))
ax1.legend()
plt.xlabel('dpr reflectivity in dbz')
plt.ylabel('radolan reflectivity in dbz')
plt.grid()
plt.show()'''
from pcc import get_my_cmap2
xbins = 10**np.arange(-1, 2.08, 0.08)
ybins = 10**np.arange(-1, 2.08, 0.08)
counts, _, _ = np.histogram2d(c[maske_cd],d[maske_cd], bins=(xbins, ybins))
fig, ax = plt.subplots(figsize=(10,10))
plt.pcolormesh(xbins, ybins, counts.T, cmap=get_my_cmap2(),vmin=0.1)
cb = plt.colorbar()
cb.set_label('number of samples', fontsize=ff)
ax.set_xscale('log')
ax.set_yscale('log')
plt.xlabel('DPR Rainrate in mm/h',fontsize=ff)
plt.ylabel('Radolan Rainrate in mm/h',fontsize=ff)
plt.xticks(fontsize=ff)
plt.yticks(fontsize=ff)
plt.title('Corr: ' + str(round(r_corr2,3)) + r'$\pm$'+ str(round(r_eror2,3))+
'\nBias: '+ str(round(r_bias2,3)), fontsize=20)
plt.grid()
plt.show()
plt.xticks(fontsize=ff)
plt.yticks(fontsize=ff)
plt.title('Winter Corr: ' + str(round(corr,3)) + r'$\pm$'+ str(round(eror,3)), fontsize=20)
plt.grid()
plt.show()
#plt.savefig('/automount/ags/velibor/plot/IRS/NSwinter_RR.png')
#plt.close()
'''
fig =plt.figure(figsize=(8,8))
ax2 = fig.add_subplot(111, aspect='equal')
plt.hist2d(C[maske_z],D[maske_z], bins=60, cmap=get_my_cmap2(),vmin=0.1, vmax=400)
cbar = plt.colorbar(shrink=0.7)
cbar.set_label('number of samples', fontsize=ff2)
#cbar.set_ticks(labelsize=20)
cbar.ax.tick_params(labelsize=ff)
cx,cy = np.arange(0,80,1),np.arange(0,80,1)
plt.plot(cx,cy, color='black')
plt.title('Winter Corr: ' + str(round(corr2,3)) + r'$\pm$'+ str(round(eror2,3)), fontsize=ff2)
plt.xlim(15,70)
plt.ylim(15,70)
plt.xlabel('DPR Reflectivity in dBZ',fontsize=ff2)
plt.ylabel('Radolan Reflectivity in dBZ',fontsize=ff2)