mod_resistivity Module

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



Contents


Variables

TypeVisibilityAttributesNameInitial
type(settings_t), private, pointer:: settings=> null()
type(background_t), private, pointer:: background=> null()

Derived Types

type, public :: resistivity_t

Components

TypeVisibilityAttributesNameInitial
procedure(realdp), public, pointer, nopass:: eta
procedure(realdp), public, pointer, nopass:: detadT
procedure(realdp), public, pointer, nopass:: detadr

Type-Bound Procedures

procedure, public :: delete

Functions

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

Arguments

TypeIntentOptionalAttributesName
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

TypeIntentOptionalAttributesName
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

TypeIntentOptionalAttributesName
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

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)


Subroutines

private subroutine delete(this)

Arguments

TypeIntentOptionalAttributesName
class(resistivity_t), intent(inout) :: this