def run_test_1d():
extra_model_kwargs = {
'rho_0': 0.1,
'onesided_flag': True,
'chi': 0.0,
'origin_flag': True,
}
model_kwargs = dict(default_model_1d_kwargs, **extra_model_kwargs)
m = model.Model1D(**model_kwargs)
t_output_every = 20.0
t_upto = 100.0
output_dir = 'test_1d'
force_resume = None
run_model(t_output_every, output_dir, m, force_resume, t_upto=t_upto)
def run_cannock_1d():
extra_model_kwargs = {
'rho_0': 1.0,
'onesided_flag': True,
'chi': 1.5,
'origin_flag': True,
}
model_kwargs = dict(default_model_1d_kwargs, **extra_model_kwargs)
m = model.Model1D(**model_kwargs)
t_output_every = 50.0
t_upto = 1e4
output_dir = '/Users/ewj/Desktop/cannock/agent_data/{}'.format(m)
force_resume = None
run_model(t_output_every, output_dir, m, force_resume, t_upto=t_upto)
def run_test_2d():
trap_n = 1
extra_model_kwargs = {
'rho_0': 2e-4,
'onesided_flag': True,
'chi': 0.0,
'walls': walls.Traps(L, dx, trap_n, trap_d, trap_w, trap_s),
'origin_flag': True,
}
model_kwargs = dict(default_model_2d_kwargs, **extra_model_kwargs)
m = model.Model2D(**model_kwargs)
t_output_every = 20.0
t_upto = 100.0
output_dir = 'test_2d'
force_resume = None
run_model(t_output_every, output_dir, m, force_resume, t_upto=t_upto)
def run_trap_no_chi():
trap_n = 1
extra_model_kwargs = {
'rho_0': 1e-1,
'onesided_flag': False,
'chi': 0.0,
'walls': walls.Traps(L, dx, trap_n, trap_d, trap_w, trap_s),
'origin_flag': False,
'p_0': 1.0,
'c_source': 0.0,
}
model_kwargs = dict(default_model_2d_kwargs, **extra_model_kwargs)
m = model.Model2D(**model_kwargs)
t_output_every = 100.0
t_upto = 4e4
output_dir = None
force_resume = None
run_model(t_output_every, output_dir, m, force_resume, t_upto=t_upto)
def run_field():
rho_0 = 0.002
c_delta_0 = 0.1
c_delta = c_delta_0 / rho_0
extra_model_kwargs = {
'rho_0': rho_0,
'spatial_flag': True,
'periodic_flag': True,
'pore_flag': True,
'origin_flags': np.array([True, False]),
'L': np.array([1000.0, 300.0]),
'pore_pf': 0.1,
'c_field_flag': True,
'c_dx': 50.0,
'c_D': 10.0,
'c_delta': c_delta,
'c_0': 1.0,
'tumble_flag': True,
'tumble_chemo_flag': True,
'chi': 100.0,
'temporal_chemo_flag': False,
}
model_kwargs = dict(default_model_kwargs, **extra_model_kwargs)
model = Model(**model_kwargs)
t_output_every = 1.0
t_upto = 50.0
output_dir = None
force_resume = None
run_utils.run_model(t_output_every, output_dir, m=model,
force_resume=force_resume, t_upto=t_upto)
def run_spatial():
extra_model_kwargs = {
'spatial_flag': True,
'periodic_flag': True,
'pore_flag': True,
'pore_pf': 0.0707,
'tumble_flag': True,
'tumble_chemo_flag': True,
'temporal_chemo_flag': True,
'onesided_flag': True,
'chi': 0.5,
}
model_kwargs = dict(default_model_kwargs, **extra_model_kwargs)
model = Model(**model_kwargs)
t_output_every = 50.0
t_upto = 100.0
output_dir = None
force_resume = None
run_utils.run_model(t_output_every, output_dir, m=model,
force_resume=force_resume, t_upto=t_upto)