 |
V3FIT
|
2 SUBROUTINE add_resistive_e
9 USE descriptor_mod,
ONLY: iam
11 USE siesta_displacement,
ONLY: update_upperv
18 INTEGER,
PARAMETER :: nits = 10
19 REAL(dp),
PARAMETER :: zero = 0
20 REAL(dp) :: save_eta, fmhd, eps
22 LOGICAL,
PARAMETER :: lcurrent_only=.true.
25 INTEGER,
SAVE :: nCheck = 0
49 CALL update_state (.false., zero, zero)
54 IF (ncheck .EQ. 0)
CALL checkforces(
xc,
gc)
66 IF (.NOT.
lverbose .AND. j.EQ.6) cycle
67 WRITE (j,
'(/,a,i3)')
' UPDATING RESISTIVE E-FIELD: ITERATIONS='
68 WRITE (j,
'(a,1p2e12.3)')
69 ' JUMP IN FSQ DUE TO RESISTIVE DIFFUSION: ', fmhd,
fsq_total1
80 END SUBROUTINE add_resistive_e
integer nprecon
Preconditioner flag.
logical lverbose
Use verbose screen output.
logical l_getfsq
Compute |F|^2.
Module contained subroutines and functions updating MHD forces and Wmhd.
real(dp), dimension(:), allocatable xc
1D array of Fourier mode displacement components.
This file contains subroutines for aupdating from t to t + delta_t the magnetic field and pressure as...
logical l_init_state
Store initial field/pressure state.
This file contains subroutines for updating half-grid magnetic fields.
subroutine update_bfield(l_add_res)
Update contravariant componets of the magnetic field.
subroutine, public init_state(lcurrent_only, lpar_in)
Initialize equilibrium state.
real(dp) eta_factor
Resistivity value.
integer, parameter unit_out
File output io unit.
real(dp) ftol
Force tolarance.
real(dp), dimension(:), allocatable, target gc
1D Array of Fourier mode MHD force components
logical l_update_state
Update the ste array.
real(dp) fsq_total1
|F|^2 WITHOUT column scaling.
This file contains all the variables and maximum sizes of the inputs for a SIESTA namelist input file...
Initializes unperturbed siesta fields and pressure in real space.
This file contains variables and parameters used by many modules in SIESTA.