parser.add_argument('--fontsize', dest = "fontsize", type = int, default = 16, help = 'fontsize of plot annotations (default: 16)') parser.add_argument('--linewidth', dest = "linewidth", type = int, default = 3, help = 'fontsize of plot annotations (default: 3)') parser.add_argument('--dpi', dest = "dpi", type = int, default = 100, help = 'dpi of plot annotations (default: 100)') return parser ############################################################################### ### Parse Arguments ### args = new_argument_parser().parse_args() ### Get Fargo Parameters ### p = Parameters() num_rad = p.ny; num_theta = p.nx r_min = p.ymin; r_max = p.ymax surface_density_zero = p.sigma0 taper_time = p.masstaper viscosity = p.nu dt = p.ninterm * p.dt fargo_par = util.get_pickled_parameters() jupiter_mass = 1e-3
parser.add_argument('--fontsize', dest = "fontsize", type = int, default = 16, help = 'fontsize of plot annotations (default: 16)') parser.add_argument('--linewidth', dest = "linewidth", type = int, default = 3, help = 'fontsize of plot annotations (default: 3)') parser.add_argument('--dpi', dest = "dpi", type = int, default = 100, help = 'dpi of plot annotations (default: 100)') return parser ############################################################################### ### Parse Arguments ### args = new_argument_parser().parse_args() ### Get Fargo Parameters ### p = Parameters() num_rad = p.ny; num_theta = p.nx r_min = p.ymin; r_max = p.ymax surface_density_zero = p.sigma0 planet_mass = 1.0 taper_time = p.masstaper viscosity = p.nu dt = p.ninterm * p.dt fargo_par = util.get_pickled_parameters() jupiter_mass = 1e-3
parser.add_argument('--dpi', dest="dpi", type=int, default=100, help='dpi of plot annotations (default: 100)') return parser ############################################################################### ### Parse Arguments ### args = new_argument_parser().parse_args() ### Get Fargo Parameters ### p = Parameters() num_rad = p.ny num_theta = p.nx r_min = p.ymin r_max = p.ymax surface_density_zero = p.sigma0 taper_time = p.masstaper scale_height = p.aspectratio viscosity = p.nu dt = p.ninterm * p.dt
parser.add_argument('--dpi', dest="dpi", type=int, default=100, help='dpi of plot annotations (default: 100)') return parser ############################################################################### ### Parse Arguments ### args = new_argument_parser().parse_args() ### Get Fargo Parameters ### p = Parameters() num_rad = p.ny num_theta = p.nx r_min = p.ymin r_max = p.ymax surface_density_zero = p.sigma0 dust_surface_density_zero = p.sigma0 * p.epsilon taper_time = p.masstaper viscosity = p.nu scale_height = p.aspectratio dt = p.ninterm * p.dt
return parser ############################################################################### ### Parse Arguments ### args = new_argument_parser().parse_args() directories = master_directories[args.choice] accretion_rates = master_accretion_rates[args.choice] start_times = master_start_times[args.choice] end_times = master_end_times[args.choice] frame_ranges = master_frame_ranges[args.choice] ### Get Fargo Parameters ### p = Parameters(directory="../" + directories[0]) num_rad = p.ny num_theta = p.nx r_min = p.ymin r_max = p.ymax surface_density_zero = p.sigma0 taper_time = p.masstaper viscosity = p.nu dt = p.ninterm * p.dt fargo_par = util.get_pickled_parameters(directory="../" + directories[0])
from pylab import * import sys sys.path.insert(0, '../../utils/python') from reader import Fields from advanced import Parameters n = 1 outputdir = "../../outputs/fargo_dusty" fluids = ["gas", "dust1", "dust2", "dust3"] field = "dens" fig = figure(figsize=(20, 5)) p = Parameters(outputdir) for i, fluid in enumerate(fluids): data = Fields("../../outputs/fargo_dusty", fluid, n).get_field(field).reshape(p.ny, p.nx) ax = fig.add_subplot(1, len(fluids), i + 1) ax.imshow(log10(data), origin='lower', aspect='auto') show()
help='save density for num_grains grain sizes (default: False)') parser.add_argument('-o', dest="make_opacities", action='store_true', default=False, help='make opacities (default: False)') return parser ### Parse Arguments ### args = new_argument_parser().parse_args() ### Get FARGO Parameters ### p = Parameters() fargo_par = util.get_pickled_parameters() num_rad = p.ny num_theta = p.nx r_min = p.ymin r_max = p.ymax surface_density_zero = p.sigma0 dust_surface_density_zero = p.sigma0 * p.epsilon planet_mass = 1.0 taper_time = p.masstaper scale_height = p.aspectratio viscosity = p.nu
parser.add_argument('--dpi', dest="dpi", type=int, default=100, help='dpi of plot annotations (default: 100)') return parser ############################################################################### ### Parse Arguments ### args = new_argument_parser().parse_args() ### Get Fargo Parameters ### p = Parameters() num_rad = p.ny num_theta = p.nx r_min = p.ymin r_max = p.ymax surface_density_zero = p.sigma0 taper_time = p.masstaper scale_height = p.aspectratio viscosity = p.nu dt = p.ninterm * p.dt
parser.add_argument('--dpi', dest="dpi", type=int, default=100, help='dpi of plot annotations (default: 100)') return parser ############################################################################### ### Parse Arguments ### args = new_argument_parser().parse_args() ### Get Fargo Parameters ### p = Parameters() num_rad = p.ny num_theta = p.nx r_min = p.ymin r_max = p.ymax surface_density_zero = p.sigma0 taper_time = p.masstaper scale_height = p.aspectratio viscosity = p.nu dt = p.ninterm * p.dt