QXPLR(3)

SYNOPSIS
SUBROUTINE QXPLR (LUN, KPRINT, IPASS)

NAME
qxplr -

LIBRARY
SLATEC

KEYWORDS
QUICK CHECK

AUTHOR
(UNKNOWN)

DESCRIPTION
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * F I S H P A K * * * * * * A PACKAGE OF FORTRAN SUBPROGRAMS FOR THE SOLUTION OF * * * * SEPARABLE ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS * * * * (VERSION 3 , JUNE 1979) * * * * BY * * * * JOHN ADAMS, PAUL SWARZTRAUBER AND ROLAND SWEET * * * * OF * * * * THE NATIONAL CENTER FOR ATMOSPHERIC RESEARCH * * * * BOULDER, COLORADO (80307) U.S.A. * * * * WHICH IS SPONSORED BY * * * * THE NATIONAL SCIENCE FOUNDATION * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * PROGRAM TO ILLUSTRATE THE USE OF SUBROUTINE HWSPLR TO SOLVE THE EQUATION (1/R)(D/DR)(R*(DU/DR)) + (1/R**2)(D/DTHETA)(DU/DTHETA) = 16*R**2 ON THE QUARTER-DISK 0 .LT. R .LT. 1, 0 .LT. THETA .LT. PI/2 WITH WITH THE BOUNDARY CONDITIONS U(1,THETA) = 1 - COS(4*THETA), 0 .LE. THETA .LE. 1 AND (DU/DTHETA)(R,0) = (DU/DTHETA)(R,PI/2) = 0, 0 .LE. R .LE. 1. (NOTE THAT THE SOLUTION U IS UNSPECIFIED AT R = 0.) THE R-INTERVAL WILL BE DIVIDED INTO 50 PANELS AND THE THETA-INTERVAL WILL BE DIVIDED INTO 48 PANELS.

ROUTINES CALLED
HWSPLR, PIMACH

REVISION HISTORY (YYMMDD)
800103 DATE WRITTEN
890718 Changed computation of PI to use PIMACH. (WRB)
890911 Removed unnecessary intrinsics. (WRB)
890911 REVISION DATE from Version 3.2
891214 Prologue converted to Version 4.0 format. (BAB)
901010 Added PASS/FAIL message and cleaned up FORMATs. (RWC)