Commit 0b1dcca7 authored by Markus Gaug's avatar Markus Gaug
Browse files

don't know what is different

parent befd6651
Pipeline #7373 canceled with stages
......@@ -51,7 +51,7 @@ z0_guess = 3.3 # z0 in beta inversion (referene point with assumed mol.
# initialize the axis ranges
x_min = 0 # in km
x_max = 6.3 # in km
y_min = 11.8
y_min = 12.
y_max = 22.5
# molecular fit ranges
......@@ -76,8 +76,10 @@ min_ranges_beta = [ 0.5 ]
max_ranges_beta = [ 4. ]
sigma_min = -3e-4 *1e3
sigma_max = 0.003*1e3
sigma_max = 0.0038*1e3
lr_min = -5.
lr_max = 110.
# standard channel configuration
labels = [ 'BRL 355 nm (elastic)', 'BRL 355 nm (photon counting)',
......@@ -386,7 +388,7 @@ def show_file(filename):
alpha_err = np.sqrt(aod_err**2/(x_mol2 - x_mol1)**2 + aod_fitted**2/(x_mol2 - x_mol1)**2 /12)
print ('AOD fitted: ', aod_fitted, ' ans2: ', ans2[0], ' ans1: ', ans1[0], 'alpha_fitted:', alpha_fitted, ' alpha_err:', alpha_err)
ax[1].errorbar(alpha_fitted, 0.5*(x_mol1+x_mol2)+offset, (x_mol2-x_mol1)/np.sqrt(12), alpha_err, label=r'$\alpha_{aer}$'+' ({:.0f} nm)'.format(channel.wavelength_nm))
ax[1].errorbar(alpha_fitted, 0.5*(x_mol1+x_mol2)+offset, (x_mol2-x_mol1)/np.sqrt(12), alpha_err, label=r'$\alpha_{aer}$'+' ({:.0f} nm)'.format(channel.wavelength_nm), color=colors[channel_id])
if (channel_id == raman_id+1):
print ('b_widths: len: ', len(b_widths[ids]), ' cont: ', b_widths[ids], ' x: len: ',len(x_height), ' cont: ', x_height)
......@@ -473,7 +475,7 @@ def show_file(filename):
#print ('alphaaer_filt', alphaaer_filt)
#ax2.plot(beta_aerosol,x_raman_filt/1000.*costheta,label=r'$\beta_{back}$ (355 nm)'+'z0={:.2f}'.format(z0), color=colors[channel_id])
ax2.plot(beta_aerosol_uv,x_raman_filt/1000.*costheta,label=r'$\beta_{back}$ (355 nm)', color=colors[channel_id])
ax2.plot(beta_aerosol_uv,x_raman_filt/1000.*costheta,label=r'$\beta_{back}$ (355 nm)', color='maroon') #colors[channel_id])
ids_lr = np.where((x_raman_filt/1000.*costheta > 1.4) & (x_raman_filt/1000.*costheta < 1.55))
mean_uv = np.mean(alphaaer_filt[ids_lr]/beta_aerosol_uv[ids_lr])
......@@ -557,17 +559,17 @@ def show_file(filename):
ids_magic = np.where(height_magic > min_magic)
if (vertical_axis):
ax[0].plot(sr_magic[ids_magic]+np.log(7.5/48.),height_magic[ids_magic],label=u"MAGIC, 21:09 UTC, 0.4\u00B0", color='gray', linestyle='dashed')
ax[0].plot(sr_magic[ids_magic]+np.log(7.5/48.),height_magic[ids_magic],label=u"MAGIC, 21:09 UTC, 0.4\u00B0\n (preliminary)", color='gray', linestyle='dashed')
else:
ax[0].plot(height_magic[ids_magic],sr_magic[ids_magic]+np.log(7.5/48.),label=u"MAGIC, 21:09 UTC, 0.4\u00B0", color='gray', linestyle='dashed')
height_smagic = np.array(magic_stab.x/1000.)
sigma_magic = np.array(magic_stab.sigma)*1000.
ids_smagic = np.where(height_smagic > min_magic)
ax[1].plot(sigma_magic[ids_smagic],height_smagic[ids_smagic], label=r"$\alpha_{aer}$ (MAGIC)",color='gray', linestyle='dashed')
ax[1].plot(sigma_magic[ids_smagic],height_smagic[ids_smagic], label=r"$\alpha_{aer}$ (MAGIC, preliminary)",color='gray', linestyle='dashed')
#ax[0].legend(loc='lower right')
ax[0].legend(loc='lower left', bbox_to_anchor=(0.465,0.03))
ax[0].legend(loc='lower left', bbox_to_anchor=(0.465,0.026))
#ax[1].legend(loc='lower right')
ax[1].legend(loc='lower left', bbox_to_anchor=(0.465,0.03))
ax2.legend(loc='upper right', bbox_to_anchor=(0.95,0.95))
......@@ -576,8 +578,12 @@ def show_file(filename):
ax[2].legend(loc='lower left', bbox_to_anchor=(0.465,0.03))
ax[0].grid(axis='y')
ax[1].grid(axis='y')
#ax[2].grid(axis='x')
ax[2].grid(axis='y')
ax[1].set_xlim(sigma_min, sigma_max)
ax[2].set_xlim(lr_min, lr_max)
ax2.set_xlim(sigma_min*1e-2*9, sigma_max*1e-2*9)
ax3.set_xlim(lr_min*0.01, lr_max*0.01)
if (vertical_axis):
ax[0].set_ylabel('Altitude above LIDAR (km)')
ax[1].set_ylabel('Altitude above LIDAR (km)')
......@@ -590,7 +596,6 @@ def show_file(filename):
ax[2].set_xlabel(r"LIDAR Ratio (sr)")
ax2.set_xlabel(r"Backscatter coefficient $\beta$ (km$^{-1}$ sr$^{-1}$)")
ax3.set_xlabel(r"$\AA$ngstr$\"o$m exponent")
ax3.set_xlim(0.,1.)
else:
ax[0].set_ylabel(r'Logarithm range-corrected signal (log(p.e. $\cdot m) )$')
plt.savefig('BRL_RamanLIDAR_{:s}_{:s}.pdf'.format(ddate,dtime))
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment