processing
Computations of phyvars that are not output by StagYY.
Time series are returned along with the time at which the variables are evaluated. Radial profiles are returned along with the radial positions at which the variables are evaluated.
dtime
Compute time increment.
Compute dt as a function of time.
Parameters:
Name | Type | Description | Default |
---|---|---|---|
sdat |
StagyyData
|
a |
required |
Returns:
Type | Description |
---|---|
Tseries
|
dt and time arrays. |
dt_dt
Compute derivative of temperature.
Compute dT/dt as a function of time using an explicit Euler scheme.
Parameters:
Name | Type | Description | Default |
---|---|---|---|
sdat |
StagyyData
|
a |
required |
Returns:
Type | Description |
---|---|
Tseries
|
derivative of temperature and time arrays. |
ebalance
Compute energy balance.
Compute Nu_t - Nu_b + V*dT/dt as a function of time using an explicit Euler scheme. This should be zero if energy is conserved.
Parameters:
Name | Type | Description | Default |
---|---|---|---|
sdat |
StagyyData
|
a |
required |
Returns:
Type | Description |
---|---|
Tseries
|
energy balance and time arrays. |
mobility
Compute plates mobility.
Compute the ratio vsurf / vrms.
Parameters:
Name | Type | Description | Default |
---|---|---|---|
sdat |
StagyyData
|
a |
required |
Returns:
Type | Description |
---|---|
Tseries
|
mobility and time arrays. |
delta_r
diff_prof
diffs_prof
advts_prof
advds_prof
advas_prof
energy_prof
advth
init_c_overturn
c_overturned
Compute theoretical overturned concentration.
This compute the resulting composition profile if fractional crystallization of a SMO is assumed and then a purely radial overturn happens.
Parameters:
Name | Type | Description | Default |
---|---|---|---|
step |
Step
|
a |
required |
Returns:
Type | Description |
---|---|
Rprof
|
the composition and radius. |