mod_physicaldata.t

Go to the documentation of this file.

module mod_physicaldata
  implicit none
  save
 
  type state
     !> Variables, normally cell center conservative values
     double precision, dimension(:^D&,:), allocatable :: w
     !> Variables, cell face values
     double precision, dimension(:^D&,:), allocatable :: ws
     !> Variables, cell edge values
     double precision, dimension(:^D&,:), allocatable :: we
     !> Variables, cell corner values
     double precision, dimension(:^D&,:), allocatable :: wc
     !> extra variables do not need ghost cell and equation flux
     double precision, dimension(:^D&,:), pointer :: wextra=>null()
     !> debug scratch: arbitrary intermediate fields captured via debug_store,
     !> flushed to a standalone .dat by save_wdebug (see mod_input_output)
     double precision, dimension(:^D&,:), allocatable :: wdebug
     !> Time-independent magnetic field at cell center and cell interface
     double precision, dimension(:^D&,:,:), pointer :: b0=>null()
     !> Time-independent electric current density at cell center
     double precision, dimension(:^D&,:), pointer :: j0=>null()
     !> Time-independent equi vars (B0 is not into this array)
     double precision, dimension(:^D&,:,:), pointer :: equi_vars=>null()
     !> Cell-center positions
     double precision, dimension(:^D&,:), pointer :: x=>null()
     !> Cell sizes in coordinate units
     double precision, dimension(:^D&,:), pointer :: dx=>null()
     !> Cell local timesteps
     double precision, dimension(:^D&), pointer :: dt=>null()
     !> Cell sizes at cell center in length unit
     double precision, dimension(:^D&,:), pointer :: ds=>null()
     !> Cell sizes at cell face in length unit
     double precision, dimension(:^D&,:), pointer :: dsc=>null()
     !> Volumes of a cell
     double precision, dimension(:^D&), pointer :: dvolume=>null()
     !> Areas of cell-center surfaces
     double precision, dimension(:^D&,:), pointer :: surface=>null()
     !> Areas of cell-face surfaces
     double precision, dimension(:^D&,:), pointer :: surfacec=>null()
     !> special values for a block
     double precision, dimension(:), pointer :: special_values=>null()
     !> ID of a grid block
     integer :: igrid=-1
     !> index range of block array in cell centers
     integer :: ixg^l
     !> index range of block array in cell faces
     integer :: ixgs^l
     !> level of AMR
     integer :: level
     !> If it face a physical boundary
     logical, dimension(:), pointer :: is_physical_boundary(:) =>null()
  end type state
 
{^nooned
  type state_sub
     !> Variables, normally center
     double precision, dimension(:^DE&,:), allocatable :: w
     !> Variables for the cornerpositions on the slice 
     double precision, dimension(:^DE&,:), allocatable :: wc
     !> Variables, normally center, one level coarser representative
     double precision, dimension(:^DE&,:), allocatable :: wcoarse
     !> Cell-center positions
     double precision, dimension(:^DE&,:), allocatable :: x
     !> Corner positions on the slice
     double precision, dimension(:^DE&,:), allocatable :: xc
     !> Cell-center positions, one level coarser representative
     double precision, dimension(:^DE&,:), allocatable :: xcoarse
     !> Cell sizes
     double precision, dimension(:^DE&,:), allocatable :: dx
     !> Cell sizes, one level coarser
     double precision, dimension(:^D&,:), allocatable :: dxcoarse
     !> Cell sizes in length unit
     double precision, dimension(:^D&,:), allocatable :: ds
     !> Volumes of a cell
     double precision, dimension(:^DE&), allocatable :: dvolume
     !> Volumes of a cell, one level coarser representative
     double precision, dimension(:^DE&), allocatable :: dvolumecoarse
     !> Areas of cell-center surfaces 
     double precision, dimension(:^DE&,:), allocatable :: surface
     !> Areas of cell-face surfaces
     double precision, dimension(:^DE&,:), allocatable :: surfacec
     !> ID of a grid block
     integer :: igrid=-1
  end type state_sub
}
{^ifoned
  type state_sub
     !> Variables, normally center
     double precision, dimension(:), allocatable :: w
     !> Variables for the cornerpositions on the slice 
     double precision, dimension(:), allocatable :: wc
     !> Variables, normally center, one level coarser representative
     double precision, dimension(:), allocatable :: wcoarse
     !> Cell-center positions
     double precision, dimension(:), allocatable :: x
     !> Corner positions on the slice
     double precision, dimension(:), allocatable :: xc
     !> Cell-center positions, one level coarser representative
     double precision, dimension(:), allocatable :: xcoarse
     !> ID of a grid block
     integer :: igrid=-1
  end type state_sub
}
  type grid_field
     !> Variables new state
     double precision, dimension(:^D&,:), allocatable :: w
     !> Variables old state
     double precision, dimension(:^D&,:), allocatable :: wold
  end type grid_field
  !> buffer for pole boundary
  type(state) :: pole_buf
 
  !> --- debug field dump (flushed by save_wdebug in mod_input_output) ---
  !> when .true., intermediate fields captured in ps(:)%wdebug are written out
  logical :: wdebug_on = .false.
  !> number of debug slots allocated in ps(:)%wdebug
  integer :: n_wdebug = 0
 
  !> array of physical states for all blocks on my processor
  type(state), dimension(:), allocatable, target :: ps
  !> array of physical states, temp 1 for multi-step time integrator
  type(state), dimension(:), allocatable, target :: ps1
  !> array of physical states, temp 2 for multi-step time integrator
  type(state), dimension(:), allocatable, target :: ps2
  !> array of physical states, temp 3 for multi-step time integrator
  type(state), dimension(:), allocatable, target :: ps3
  !> array of physical states, temp 4 for multi-step time integrator
  type(state), dimension(:), allocatable, target :: ps4
  !> array of physical blocks, one level coarser representative
  type(state), dimension(:), allocatable, target :: psc
 
  !> array of physical blocks in reduced dimension
  type(state_sub), dimension(:), allocatable, target :: ps_sub
 
{^ifoned
  double precision, dimension(:), allocatable :: collapseddata
}
{^nooned
  double precision, dimension(:^DE&,:), allocatable :: collapseddata
}
  !> array of physical blocks of meshed fields for particles
  type(grid_field), dimension(:), allocatable, target :: gridvars
 
  !> velocities store for constrained transport
  type ct_velocity
    double precision, dimension(:^D&,:), allocatable :: vnorm,cbarmin,cbarmax
    double precision, dimension(:^D&,:,:), allocatable :: vbarc,vbarlc,vbarrc
  end type ct_velocity
 
end module mod_physicaldata