%matplotlib inline
This script produces a default parameter file then uses this to set up a model run. The default parameter file reproduces the constant material properties case in Murphy Quinlan et al., 2021.
In order for this to run without installing the package, we put the directory on the PYTHONPATH. This should be removed if you have installed the package and are running this from your machine:
# Setup PYTHONPATH to allow pytesimal
# import without install of package
import setup_path
import pytesimal.load_plot_save
import pytesimal.quick_workflow
Define a filepath and filename for the parameter file:
folderpath = 'example_default/'
filename = 'example_parameters'
filepath = f'{folderpath}{filename}.txt'
Save a default parameters file to the filepath:
pytesimal.load_plot_save.make_default_param_file(filepath)
You can open this json
file in a text editor and change the default values.
For this example, we're just leaving the default values as they are and
loading it without editing, and starting a model run:
pytesimal.quick_workflow.workflow(filename, folderpath)
Just wait for a minute or two for your planetesimal to evolve!
# Once 400 millions years has passed and your planetesimal has cooled down,
# we can load the results in to analyse and plot:
filepath = 'example_default/example_parameters_results.npz'
(temperatures,
coretemp,
dT_by_dt,
dT_by_dt_core) = pytesimal.load_plot_save.read_datafile(filepath)
We can visualise the cooling history of the planetesimal:
# Specify a figure width and height
fig_w = 6
fig_h = 9
pytesimal.load_plot_save.two_in_one(fig_w,
fig_h,
temperatures,
coretemp,
dT_by_dt,
dT_by_dt_core)