potential2D_polarization Module Procedure

module procedure potential2D_polarization module function potential2D_polarization(srcterm, SigP2, SigP3, SigH, Cm, v2, v3, Vminx2, Vmaxx2, Vminx3, Vmaxx3, dt, x, Phi0, perflag, it)

SOLVE IONOSPHERIC POTENTIAL EQUATION IN 2D USING MUMPS INCLUDES FULL OF POLARIZATION CURRENT, INCLUDING CONVECTIVE TERMS. VELOCITIES SHOULD BE TRIMMED (WITHOUT GHOST CELLS). THIS VERSION OF THE INTEGRATED POTENTIAL SOLVER OBVIATES ALL OTHERS SINCE A PURELY ELECTRSTATIC FORM CAN BE RECOVERED BY ZEROING OUT THE INERTIAL CAPACITANCE.

THIS FORM IS INTENDED TO WORK WITH CURVILINEAR MESHES. NOTE THAT THE FULL GRID VARIABLES (X%DX3ALL, ETC.) MUST BE USED HERE!!!

note the alt since we need to use dx2 as differential... Tricky bug/feature awkward way of handling this special case derivative which uses x3 as the differential to operate on a 2D array.

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(:,:):: 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
type(curvmesh), intent(in) :: x
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))


Calls

proc~~potential2d_polarization~~CallsGraph proc~potential2d_polarization potential2D_polarization interface~grad2d1_curv_alt grad2D1_curv_alt proc~potential2d_polarization->interface~grad2d1_curv_alt interface~grad2d3 grad2D3 proc~potential2d_polarization->interface~grad2d3 interface~elliptic2d_polarization elliptic2D_polarization proc~potential2d_polarization->interface~elliptic2d_polarization proc~grad2d1_curv_alt_23 grad2D1_curv_alt_23 interface~grad2d1_curv_alt->proc~grad2d1_curv_alt_23 proc~grad2d3_curv_23 grad2D3_curv_23 interface~grad2d3->proc~grad2d3_curv_23 proc~elliptic2d_polarization elliptic2D_polarization interface~elliptic2d_polarization->proc~elliptic2d_polarization mumps_exec mumps_exec proc~elliptic2d_polarization->mumps_exec proc~check_mumps_status check_mumps_status proc~elliptic2d_polarization->proc~check_mumps_status mumps_perm mumps_perm proc~elliptic2d_polarization->mumps_perm proc~quiet_mumps quiet_mumps proc~elliptic2d_polarization->proc~quiet_mumps

Called by

proc~~potential2d_polarization~~CalledByGraph proc~potential2d_polarization potential2D_polarization interface~potential2d_polarization potential2D_polarization interface~potential2d_polarization->proc~potential2d_polarization proc~potential_root_mpi_curv potential_root_mpi_curv proc~potential_root_mpi_curv->interface~potential2d_polarization interface~potential_root_mpi_curv potential_root_mpi_curv interface~potential_root_mpi_curv->proc~potential_root_mpi_curv interface~potential_root_mpi potential_root_mpi interface~potential_root_mpi->interface~potential_root_mpi_curv proc~electrodynamics_curv electrodynamics_curv proc~electrodynamics_curv->interface~potential_root_mpi interface~electrodynamics electrodynamics interface~electrodynamics->proc~electrodynamics_curv program~gemini3d Gemini3D program~gemini3d->interface~electrodynamics

Contents

None