/*----------------------------------------------------------------------*/
void
assess_init(assess_t *a)
{
integrate_init(&a->integ_aux, SIZE_EPS,
double_exponential);
dickman_init(&a->dickman_aux);
}
int initialize_grid()
{
set_bvals_init(&level0_Grid, &level0_Domain);
set_bvals(&level0_Grid, 0);
if(!is_dt_set_by_script) {
new_dt(&level0_Grid);
}
lr_states_init(level0_Grid.Nx1,level0_Grid.Nx2,level0_Grid.Nx3);
Integrate = integrate_init(level0_Grid.Nx1,level0_Grid.Nx2,level0_Grid.Nx3);
#ifdef SELF_GRAVITY
SelfGrav = selfg_init(&level0_Grid, &level0_Domain);
if(is_restart == 0) {
(*SelfGrav)(&level0_Grid, &level0_Domain);
}
/* Only bvals for Phi set when last argument of set_bvals = 1 */
set_bvals(&level0_Grid, 1);
#endif
return 0;
}
void att_ctrl_pitch_init(float dt ,float K1, float K2, float lambda){
pitch_dt = dt;
pitch_K1 = K1;
pitch_K2 = K2;
pitch_lambda = lambda;
derivative_init(&pitch_deriv, dt);
integrate_init (&pitch_inte, dt);
}
