Notes from meeting. These cover questions I had while using Butterpy, including what variables correspond with those listed in the paper, deprecated keywords, etc:
regions.py
activity_rate (regions.py) and Activity Level (paper) are the same
diffrot_shear is expressed as a fraction of the rotation frequency (unitless, 0 is no sheer, 1 means the eq rotates twice as fast as the pole, negative values allowed) the same as latitudinal rot shear (in paper)
alpha_med paramterizes spot contrast and spot area
Zach has always kept it constant so far... probably varies with stellar mass irl but this is probably hard to get good numbers for
decay_timescale in spots is as a fraction of the rotation period (unitless)
this is Spot Lifetime from the paper
decay_time in regions, units of day, part of a hacky solution and only matters at the highest activity levels (doesn't matter for activity_level <7)
keep it constant
Behavior: above activity level 7, you get saturation of spots and then it'll stop emerging new spots for a while
spots.py
Note: Need parameter descriptions in the docstring for this class.
dur: this doesn't get used anymore, ignore. Used to be used so that they could define a "burn in time" to get to equilibrium.
tau_evolve: relic, need to delete
threshold: don't have to worry about this either
Best Practices:
Run your regions for longer than you need, then take a radom snippet from somewhere in the middle. "Burn in" is typically less than a year, usually .1-.5 depending on activity rate and cycle length.
Notes from meeting. These cover questions I had while using Butterpy, including what variables correspond with those listed in the paper, deprecated keywords, etc:
regions.py
spots.py Note: Need parameter descriptions in the docstring for this class.
Best Practices: