init_theta_grid_salpha Subroutine

  subroutine init_theta_grid_salpha(theta_grid_salpha_config_in)
    use theta_grid_params, only: init_theta_grid_params
    implicit none
    type(theta_grid_salpha_config_type), intent(in), optional :: theta_grid_salpha_config_in

    if (initialized) return
    initialized = .true.

    call init_theta_grid_params
    call read_parameters(theta_grid_salpha_config_in)
  end subroutine init_theta_grid_salpha