elliptic2d Submodule


Uses

  • module~~elliptic2d~~UsesGraph module~elliptic2d elliptic2d module~pdeelliptic PDEelliptic module~elliptic2d->module~pdeelliptic module~mumps_interface~2 mumps_interface module~pdeelliptic->module~mumps_interface~2 iso_fortran_env iso_fortran_env module~pdeelliptic->iso_fortran_env module~phys_consts phys_consts module~pdeelliptic->module~phys_consts module~mpimod mpimod module~pdeelliptic->module~mpimod module~mumps_rl mumps_rl module~mumps_interface~2->module~mumps_rl module~phys_consts->iso_fortran_env module~mpimod->iso_fortran_env module~mpimod->module~phys_consts mpi mpi module~mpimod->mpi

Contents


Module Procedures

module procedure elliptic2D_polarization module function elliptic2D_polarization(srcterm, SigP2, SigP3, SigH, gradSigH2, gradSigH3, Cm, v2, v3, Vminx2, Vmaxx2, Vminx3, Vmaxx3, dt, dx1, dx1i, dx2all, dx2iall, dx3all, dx3iall, Phi0, perflag, it)

interior+boundary-x3_adj-x2_adj. Note that are 3 sets of entries for each adjacent point note the different conductiances here ot be associated with derivatives in different directions BOTTOM GRID POINTS + CORNER TOP GRID POINTS + CORNER LEFT BOUNDARY RIGHT BOUNDARY INTERIOR LOCATION ix2-1,ix3-2 grid point ix2,ix3-2 grid point bit of intentional code duplication here and in the following sections to keep things organized in a way I can debug... ix2+1,ix3-2 grid point ix2-2,ix3-1 !ix2,ix3-1 grid point in ix2,ix3 equation add in polarziation terms that include previous time step potential at this grid point

Read more…

Arguments

Type IntentOptional AttributesName
real(kind=wp), intent(in), dimension(:,:):: srcterm

ZZZ - THESE WILL NEED TO BE MODIFIED CONDUCTIVITIES, AND WE'LL NEED THREE OF THEM

real(kind=wp), intent(in), dimension(:,:):: SigP2

ZZZ - THESE WILL NEED TO BE MODIFIED CONDUCTIVITIES, AND WE'LL NEED THREE OF THEM

real(kind=wp), intent(in), dimension(:,:):: SigP3

ZZZ - THESE WILL NEED TO BE MODIFIED CONDUCTIVITIES, AND WE'LL NEED THREE OF THEM

real(kind=wp), intent(in), dimension(:,:):: SigH

ZZZ - THESE WILL NEED TO BE MODIFIED CONDUCTIVITIES, AND WE'LL NEED THREE OF THEM

real(kind=wp), intent(in), dimension(:,:):: gradSigH2

ZZZ - THESE WILL NEED TO BE MODIFIED CONDUCTIVITIES, AND WE'LL NEED THREE OF THEM

real(kind=wp), intent(in), dimension(:,:):: gradSigH3

ZZZ - THESE WILL NEED TO BE MODIFIED CONDUCTIVITIES, AND WE'LL NEED THREE OF THEM

real(kind=wp), intent(in), dimension(:,:):: Cm

ZZZ - THESE WILL NEED TO BE MODIFIED CONDUCTIVITIES, AND WE'LL NEED THREE OF THEM

real(kind=wp), intent(in), dimension(:,:):: v2

ZZZ - THESE WILL NEED TO BE MODIFIED CONDUCTIVITIES, AND WE'LL NEED THREE OF THEM

real(kind=wp), intent(in), dimension(:,:):: v3

ZZZ - THESE WILL NEED TO BE MODIFIED CONDUCTIVITIES, AND WE'LL NEED THREE OF THEM

real(kind=wp), intent(in), dimension(:):: Vminx2
real(kind=wp), intent(in), dimension(:):: Vmaxx2
real(kind=wp), intent(in), dimension(:):: Vminx3
real(kind=wp), intent(in), dimension(:):: Vmaxx3
real(kind=wp), intent(in) :: dt
real(kind=wp), intent(in), dimension(0:):: dx1
real(kind=wp), intent(in), dimension(:):: dx1i
real(kind=wp), intent(in), dimension(0:):: dx2all
real(kind=wp), intent(in), dimension(:):: dx2iall
real(kind=wp), intent(in), dimension(0:):: dx3all
real(kind=wp), intent(in), dimension(:):: dx3iall
real(kind=wp), intent(in), dimension(:,:):: Phi0
logical, intent(in) :: perflag
integer, intent(in) :: it

Return Value real(kind=wp), dimension(size(SigP2,1),size(SigP2,2))

module procedure elliptic2D_polarization_periodic module function elliptic2D_polarization_periodic(srcterm, SigP, SigH, gradSigH2, gradSigH3, Cm, v2, v3, Vminx2, Vmaxx2, Vminx3, Vmaxx3, dt, dx1, dx1i, dx2all, dx2iall, dx3all, dx3iall, Phi0, perflag, it)

static model

Read more…

Arguments

Type IntentOptional AttributesName
real(kind=wp), intent(in), dimension(:,:):: srcterm
real(kind=wp), intent(in), dimension(:,:):: SigP
real(kind=wp), intent(in), dimension(:,:):: SigH
real(kind=wp), intent(in), dimension(:,:):: gradSigH2
real(kind=wp), intent(in), dimension(:,:):: gradSigH3
real(kind=wp), intent(in), dimension(:,:):: Cm
real(kind=wp), intent(in), dimension(:,:):: v2
real(kind=wp), intent(in), dimension(:,:):: v3
real(kind=wp), intent(in), dimension(:):: Vminx2
real(kind=wp), intent(in), dimension(:):: Vmaxx2
real(kind=wp), intent(in), dimension(:):: Vminx3
real(kind=wp), intent(in), dimension(:):: Vmaxx3
real(kind=wp), intent(in) :: dt
real(kind=wp), intent(in), dimension(0:):: dx1
real(kind=wp), intent(in), dimension(:):: dx1i
real(kind=wp), intent(in), dimension(0:):: dx2all
real(kind=wp), intent(in), dimension(:):: dx2iall
real(kind=wp), intent(in), dimension(0:):: dx3all
real(kind=wp), intent(in), dimension(:):: dx3iall
real(kind=wp), intent(in), dimension(:,:):: Phi0
logical, intent(in) :: perflag
integer, intent(in) :: it

Return Value real(kind=wp), dimension(size(SigP,1),size(SigP,2))

module procedure elliptic2D_cart module function elliptic2D_cart(srcterm, sig0, sigP, Vminx1, Vmaxx1, Vminx3, Vmaxx3, dx1, dx1i, dx3all, dx3iall, flagdirich, perflag, gridflag, it)

SOLVE IONOSPHERIC POTENTIAL EQUATION IN 2D USING MUMPS.

Read more…

Arguments

Type IntentOptional AttributesName
real(kind=wp), intent(in), dimension(:,:,:):: srcterm
real(kind=wp), intent(in), dimension(:,:,:):: sig0
real(kind=wp), intent(in), dimension(:,:,:):: sigP
real(kind=wp), intent(in), dimension(:,:):: Vminx1
real(kind=wp), intent(in), dimension(:,:):: Vmaxx1
real(kind=wp), intent(in), dimension(:,:):: Vminx3
real(kind=wp), intent(in), dimension(:,:):: Vmaxx3
real(kind=wp), intent(in), dimension(0:):: dx1
real(kind=wp), intent(in), dimension(:):: dx1i
real(kind=wp), intent(in), dimension(0:):: dx3all
real(kind=wp), intent(in), dimension(:):: dx3iall
integer, intent(in) :: flagdirich
logical, intent(in) :: perflag
integer, intent(in) :: gridflag
integer, intent(in) :: it

Return Value real(kind=wp), dimension(size(sig0,1),1,size(sig0,3))