shared_data.f90

Go to the documentation of this file.

!*******************************************************************************
!
!  Note separating the Doxygen comment block here so detailed decription is
!  found in the Module not the file.
!
!*******************************************************************************
      MODULE shared_data
      USE stel_kinds
      USE stel_constants
 
      IMPLICIT NONE
 
!*******************************************************************************
!  shared_data module parameters
!*******************************************************************************
!  Solver parameters
      INTEGER, PARAMETER   :: ngmres_steps = 100
      INTEGER, PARAMETER   :: prediag = 1
      INTEGER, PARAMETER   :: preblock = 2
      REAL (dp), PARAMETER :: fsq_res = 1.e-16_dp
      REAL (dp), PARAMETER :: levm_ped = 1.e-10_dp
      REAL (dp), PARAMETER :: mu_ped = 1.e-8_dp
      INTEGER, PARAMETER   :: gmres_peak = 2
      INTEGER, PARAMETER   :: gmres_no_peak = 1
 
!  IO Parameters
      INTEGER, PARAMETER   :: unit_out = 336
 
!  Force calulation flags.
      LOGICAL, PARAMETER   :: l_pedge = .true.
      LOGICAL, PARAMETER   :: l_natural = .true.
 
!*******************************************************************************
!  shared_data module variables
!*******************************************************************************
!  Size variables.
      INTEGER                                     :: neqs
      INTEGER                                     :: ndims
      INTEGER                                     :: niter
      INTEGER                                     :: mblk_size
      INTEGER                                     :: nsp
 
!  Solver control variables.
      INTEGER                                     :: nprecon
      INTEGER                                     :: nprecon_type
      INTEGER                                     :: ngmres_type = gmres_peak
      INTEGER                                     :: iortho = 3
      INTEGER                                     :: hesspass_test = -1
      INTEGER                                     :: in_hess_nfunct
      INTEGER                                     :: out_hess_nfunct
      REAL (dp)                                   :: mupar_test
      REAL (dp)                                   :: fsq_total
      REAL (dp)                                   :: fsq_total1
 
!  Solver work variables.
      REAL(dp), DIMENSION(:), ALLOCATABLE         :: xc
      REAL(dp), DIMENSION(:), ALLOCATABLE, TARGET :: gc
      REAL(dp), DIMENSION(:), ALLOCATABLE, TARGET :: gc_sup
      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE     :: buv_res
      REAL(dp), DIMENSION(:), ALLOCATABLE         :: xc0
      REAL(dp), DIMENSION(:), ALLOCATABLE         :: gc0
      REAL(dp), DIMENSION(:,:,:,:), ALLOCATABLE   :: col_scale
 
      REAL (dp)                                   :: fsq_gmres
      REAL (dp)                                   :: fsq_lin
      REAL (dp)                                   :: etak_tol
      REAL (dp)                                   :: levm_scale = 1
      REAL (dp)                                   :: wtotal
      REAL (dp)                                   :: wtotal0
      REAL (dp)                                   :: delta_t
      REAL (dp)                                   :: fsqvs
      REAL (dp)                                   :: fsqvu
      REAL (dp)                                   :: fsqvv
      REAL (dp)                                   :: ste(4)
      REAL (dp)                                   :: bs0(12)
      REAL (dp)                                   :: bu0(12)
      REAL (dp)                                   :: bsbu_ratio_s
      REAL (dp)                                   :: jsju_ratio_s
      REAL (dp)                                   :: bsbu_ratio_a
      REAL (dp)                                   :: jsju_ratio_a
      REAL (dp)                                   :: scale_s
      REAL (dp)                                   :: scale_u
 
!  Cordinate basis.
      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE     :: r1_i
      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE     :: z1_i
      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE     :: ru_i
      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE     :: zu_i
      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE     :: rv_i
      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE     :: zv_i
 
!  Shared quantities.
      REAL (dp)                                   :: jsupvdota
      REAL (dp), DIMENSION(:), ALLOCATABLE        :: torflux
      REAL (dp), DIMENSION(:), ALLOCATABLE        :: polflux
      REAL (dp), DIMENSION(:,:,:), POINTER        :: asubsmnsf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: asubsmncf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: asubumncf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: asubumnsf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: asubvmncf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: asubvmnsf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: fsupsmnsf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: fsupsmncf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: fsupumncf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: fsupumnsf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: fsupvmncf
      REAL (dp), DIMENSION(:,:,:), POINTER        :: fsupvmnsf
      
!  Flags to control evolving origin and edge.
      LOGICAL :: l_push_edge
      LOGICAL :: l_push_s
      LOGICAL :: l_push_u
      LOGICAL :: l_push_v
      LOGICAL :: l_linearize
      LOGICAL :: l_conjgrad
      LOGICAL :: l_getwmhd
      LOGICAL :: l_getfsq
      LOGICAL :: l_applyprecon
      LOGICAL :: l_printoriginforces = .false.
      LOGICAL :: l_init_state
      LOGICAL :: l_update_state = .false.
      LOGICAL :: l_par_state
      LOGICAL :: lcolscale
      LOGICAL :: lasym = .false.
      LOGICAL :: lrecon = .false.
      LOGICAL :: lverbose = .true.
      LOGICAL :: lequi1 = .true.
 
      REAL (dp) :: siesta_curtor = 0.0
 
      END MODULE shared_data