result750=loadtxt('result750nm.dat')

figure(1)
contourf(result750[:,2].reshape(30,25))
colorbar()
xlabel('effective radius [um]')
ylabel('cloud optical thickness')
title('radiance [mW / (m^2 nm sr)] @ 750 nm')