ax2.grid(True) ax2.set_yticks(0.5 + np.arange(nlayers)) labels = [] for i in range(int(winddata['Nz'][0])): labels.append("%s" % winddata['z'][i][0]) ax2.set_yticklabels(labels) ax1.set_ylim([0, 15]) ax1.scatter(ut, tau) if outfile: f.set_figheight(9) f.set_figwidth(6.5) savefig(outfile) if imageshow: show(block=True) if __name__ == '__main__': indt = True winddata = readMASS(sys.argv[1], 'WL', rtype='dict', indatetime=indt) outfile = None dimmfile = None if len(sys.argv) == 3: outfile = sys.argv[2] windplot(winddata, outfile, indatetime=indt, imageshow=True)
ax1.set_ylim([0, 4]) ax1.errorbar(mut, fseeing, yerr=e_fseeing, fmt=None, label="MASS") if dimmfile: ax1.scatter(dimm_ut, dimmdict["seeing"], alpha=0.3, label="DIMM", color="red") ax1.legend(loc=0) if outfile: f.set_figheight(9) f.set_figwidth(6.5) savefig(outfile) if imageshow: show(block=True) if __name__ == "__main__": indt = True cn2data = readMASS(sys.argv[1], "TX", rtype="dict", indatetime=indt) massdata = readMASS(sys.argv[1], "A", rtype="dict", indatetime=indt) outfile = None dimmfile = None if len(sys.argv) == 3: outfile = sys.argv[2] if len(sys.argv) == 4: outfile = sys.argv[3] dimmfile = sys.argv[2] cn2plot_image(massdata, cn2data, dimmfile, outfile, indatetime=indt, imageshow=True)