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
