コード例 #1
0
ファイル: windplot.py プロジェクト: saltastro/masspipe
    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)
コード例 #2
0
ファイル: cn2plot.py プロジェクト: saltastro/masspipe
    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)