potential2d Submodule

  1. potential2d.f90
  2. potential2d

Uses

  • module~~potential2d~~UsesGraph module~potential2d potential2d module~calculus calculus module~potential2d->module~calculus module~grid grid module~potential2d->module~grid module~potential_mumps potential_mumps module~potential2d->module~potential_mumps module~pdeelliptic PDEelliptic module~potential2d->module~pdeelliptic module~phys_consts phys_consts module~calculus->module~phys_consts module~mesh mesh module~calculus->module~mesh module~reader reader module~grid->module~reader iso_fortran_env iso_fortran_env module~grid->iso_fortran_env module~grid->module~phys_consts module~grid->module~mesh module~mpimod mpimod module~grid->module~mpimod module~potential_mumps->module~calculus module~potential_mumps->module~pdeelliptic module~potential_mumps->iso_fortran_env module~potential_mumps->module~phys_consts module~interpolation interpolation module~potential_mumps->module~interpolation module~potential_mumps->module~mesh module~potential_mumps->module~mpimod module~pdeelliptic->iso_fortran_env module~pdeelliptic->module~phys_consts module~mumps_interface~2 mumps_interface module~pdeelliptic->module~mumps_interface~2 module~pdeelliptic->module~mpimod module~reader->iso_fortran_env module~reader->module~phys_consts module~phys_consts->iso_fortran_env module~mumps_rl mumps_rl module~mumps_interface~2->module~mumps_rl module~interpolation->module~phys_consts module~mesh->module~phys_consts module~mpimod->iso_fortran_env module~mpimod->module~phys_consts mpi mpi module~mpimod->mpi

Contents

Module Procedures


Module Procedures

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.

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(:,:):: 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))

module procedure potential2D_polarization_periodic module function potential2D_polarization_periodic(srcterm, SigP, 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.

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(:,:):: 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
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(SigP,1),size(SigP,2))

module procedure potential2D_fieldresolved module function potential2D_fieldresolved(srcterm, sig0, sigP, Vminx1, Vmaxx1, Vminx3, Vmaxx3, x, flagdirich, perflag, it)

SOLVE IONOSPHERIC POTENTIAL EQUATION IN 2D USING MUMPS ASSUME THAT WE ARE RESOLVING THE POTENTIAL ALONG THE FIELD LINE AND THAT IT VARIES IN X1 AND X3 (X2 IS NOMINALL JUST ONE ELEMENT. LEFT AND RIGHT BOUNDARIES (IN X3) ARE ASSUMED TO USE DIRICHLET BOUNARY CONDITIONS, WHILE THE (ALTITUDE) TOP CAN BE NEUMANN OR DIRICHLET. BOTTOM (ALTITUDE) IS ALWAYS ASSUMED TO BE DIRICHLET.

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
type(curvmesh), intent(in) :: x
integer, intent(in) :: flagdirich
logical, intent(in) :: perflag
integer, intent(in) :: it

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