;+
;
;Name: DIPOLE
;
; Purpose: Calculate a dipole field with hydrostatic density
;          differentiated by open vs. closed field
;          at a location rpB, thetapB, phipB (true spherical coordinates)
;
;;Keyword Inputs:
;
;  BO: magnitude of dipole at photosphere
;                       DEFAULT 10 Gauss
;
;  THETOPEN: latitude at photosphere (north and south)
;       above which field is open
;  HYDRO: how to handle the plasma in the closed region
;        The options are:
;
;             HYDRO=1; exponential isothermal hydrostatic equilibrium (DEFAULT)
;             HYDRO=2; radial hydrostatic equilibrium
;             HYDRO=3; Vasquez et al 2003 hydrostatic atmosphere
;
;  DENSPROF
;             For HYDRO=1, DENSPROF = density at coronal base in CGS units, scaled by 5d8
;             For HYDRO=2, DENSPROF can be array [A,B,C,D,E,F] (units in cgs)
;                       or multiplier of array [densprof*A,B,densprof*C,D,densprof*E,F]
;                               (note for widget it has to be multiplier, so, scalar)
;                    dens = A*r^-B + C*r^-D + E*r^-F
;                  and inside DIPOLE program A,B,C,D,E,F defaults to WSM streamer values
;            For HYDRO=3, DENSPROF can be array [A1,A2,A3,A4,A5,aa,bb,alpha,beta]
;                       or multiplier of array as above
;                     dens=A1*exp(A2/r)*r^-2*(1+A3/r+A4/r^2+A5/r^3)
;                     temp=TE*(aa+1)/(aa+bb*r^alpha + (1-bb)*r^-beta)
;                  and inside DIPOLE program A1,A2,A3,A4,A5,aa,bb,alpha,beta
;                       defaults to Vasquez streamer values
;               Note- default=1 for all
;
;  TE
;             isothermal temperature value
;                       DEFAULT 1.5D6
;
;  ODENSPROF, OT0 - sets plasma for open field
;               follows same form as DENSPROF - but HYDRO=1 ODENSPROF is CGS units scaled by 1d8
;                   and inside DIPOLE A,B,C,D,E,F defaults to WSM coronal hole values
;                   or A1,A2,A3,A4,A5,aa,bb,alpha,beta  to Vasquez 2003 values
;
;
;  VELIMPOSE - impose a velocity of constant magnitude VELIMPOSE directed along the field
;               in open regions.
;                       UNITS KM/SEC
;                       DEFAULT 0.d0