Esempio n. 1
0
def setup_options(options):
    if options is None:
        parser= get_options()
        options, args= parser.parse_args([])
    #Set up to default fit
    options.potential= 'mpdiskplhalofixbulgeflat'
    options.aAmethod= 'staeckelg'
    options.dfeh=0.1
    options.dafe=0.05
    options.ngl= 20
    options.singles= True
    return options
Esempio n. 2
0
                        r'$[\alpha/\mathrm{Fe}] = %.3f$' % (afes[bin]),
                        size=16.,top_right=True)

    left, bottom, width, height= 0.1, 0.1, 0.8, 0.25
    ax2= pyplot.axes([left,bottom,width,height])
    fig= pyplot.gcf()
    fig.sca(ax2)
    bovy_plot.bovy_plot(rs,derivProps[:,2],'k-',lw=2.,
                        xrange=[4.,10.],
                        overplot=True)
    indx= numpy.argmin(numpy.fabs(derivProps[:,2]))
    bovy_plot.bovy_plot([4.,10.],[0.,0.],'-',color='0.5',overplot=True)
    bovy_plot.bovy_plot([rs[indx],rs[indx]],[derivProps[indx,2],1000.],
                        'k--',overplot=True)
    thisax= pyplot.gca()
    pyplot.ylabel(r'$\mathrm{Correlation\ between}$'+'\n'+r'$\!\!\!\!\!\!\!\!R_d\ \&\ \Sigma_{1.1}(R)$')
    pyplot.xlabel(r'$R\ (\mathrm{kpc})$')   
    pyplot.xlim(4.,10.)
    pyplot.ylim(-.5,.5)
    bovy_plot._add_ticks()
    pyplot.sca(axTop)
    bovy_plot.bovy_plot([rs[indx],rs[indx]],[0.01,sigs[3,indx]],
                        'k--',overplot=True)   
    bovy_plot.bovy_end_print(options.outfilename)
        

if __name__ == '__main__':
    parser= get_options()
    options,args= parser.parse_args()
    illustrateBestR(options,args)
    #First do all dens comparisons
    if not options.justvel:
        options.type= 'z'
        plotDensComparisonDFMulti(options,args)
        #options.type= 'R'
        #plotDensComparisonDFMulti(options,args)
    options.type= 'vt'
    plotVelComparisonDFMulti(options,args)
    options.type= 'vz'
    plotVelComparisonDFMulti(options,args)
    options.type= 'vr'
    plotVelComparisonDFMulti(options,args)
    return None

if __name__ == '__main__':
    (options,args)= get_options().parse_args()
    try:
        os.mkdir(os.path.join('..',args[0][0:-1]))
    except OSError:
        pass
    args[0]= os.path.join('..',args[0][0:-1],args[0])
    if options.allgroups:
        options.group= 'aenhanced'
        plotAllComparisonDFMulti(options,args)
        options.group= 'apoor'
        plotAllComparisonDFMulti(options,args)
        options.group= 'apoorfpoor'
        plotAllComparisonDFMulti(options,args)
        options.group= 'aintermediate'
        plotAllComparisonDFMulti(options,args)
    else: