mod_resistivity Module

Module containing resistivity-related routines, calculates and sets the resistivity values based on the equilibrium configuration.



Variables

Type Visibility Attributes Name Initial
type(settings_t), private, pointer :: settings => null()
type(background_t), private, pointer :: background => null()

Derived Types

type, public ::  resistivity_t

Components

Type Visibility Attributes Name Initial
procedure(real(dp)), public, pointer, nopass :: eta
procedure(real(dp)), public, pointer, nopass :: detadT
procedure(real(dp)), public, pointer, nopass :: detadr

Type-Bound Procedures

procedure, public :: delete

Functions

public function new_resistivity(settings_tgt, background_tgt) result(resistivity)

Arguments

Type IntentOptional Attributes Name
type(settings_t), intent(in), target :: settings_tgt
type(background_t), intent(in), target :: background_tgt

Return Value type(resistivity_t)

private function spitzer_eta(x)

Default profile for the resistivity $\eta$. Returns one of the following: - The Spitzer resistivity based on the equilibrium temperature profile - a fixed resistivity value - zero (if resistivity is disabled) Values are normalised on return.

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)

private function spitzer_detadT(x)

Default profile for the derivative of the resistivity $\eta$ with respect to the temperature $T$. Returns one of the following: - zero (if resistivity is disabled or for fixed resistivity) - derivative of Spitzer resistivity with respect to T Values are normalised on return.

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)

private function spitzer_detadr(x)

Default profile for the derivative of the resistivity $\eta$ with respect to position. Returns zero unless a dropoff profile was chosen.

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)


Subroutines

private subroutine delete(this)

Arguments

Type IntentOptional Attributes Name
class(resistivity_t), intent(inout) :: this