hyperresi_strict_2nd

Module

Description

$Id$

This module applies a sixth order hyperresistivity to the induction equation (following Brandenburg & Sarson 2002). This hyperresistivity ensures that the energy dissipation rate is positive define everywhere.

Spatial derivatives are accurate to second order.

** AUTOMATIC CPARAM.INC GENERATION ************************ Declare (for generation of cparam.inc) the number of f array variables and auxiliary variables added by this module

CPARAM logical, parameter :: lhyperresistivity_strict=.true.

MVAR CONTRIBUTION 0 MAUX CONTRIBUTION 3

Quick access

Routines:: hyperresistivity_strict(), register_hyperresi_strict()

Needed modules

cdata: $Id$
messages: $Id$

Variables

hyperresi_strict/del2_2nd [private]
hyperresi_strict/del2v_2nd [private]
hyperresi_strict/der2_2nd [private]

Subroutines and functions

subroutine hyperresi_strict/register_hyperresi_strict()

Set up indices for hyperresistivity auxiliary slots.

23-aug-07/anders: adapted from register_hypervisc_strict

Use :: farraymanager
Call to:: svn_id()

subroutine hyperresi_strict/hyperresistivity_strict(f)

Apply sixth order hyperresistivity with positive definite heating rate (see Brandenburg & Sarson 2002).

To avoid communicating ghost zones after each operator, we use derivatives that are second order in space.

23-aug-07/anders: adapted from hyperviscosity_strict

Parameters:: f (mx,my,mz,mfarray) [real]