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.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
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 |