;************************************************************************
;+
;Name: gibbaglowprams
;
; Purpose: Calculate the density, temperature, magnetic field of the
;          Gibson-Bagenal-Low model at a location rpB, thetapB, phipB
;       (true spherical coordinates)
;
;Keyword Inputs:
;
;         PHYSICAL PROPERTIES
;
;               APAR - (also affects MORPHOLOGY) just as in the Bogdan and Low model,
;                       this allows an expansion of the field and introduces
;                       currents. apar is the radial shift transformation used to make a
;                       tear-drop shape (r -> r + apar)
;                       NOTE APAR not equal to zero can introduce negative densities
;                       DEFAULT 0.0
;
;               BPAR - cusp height
;                       DEFAULT 2.78
;
;  HYDRO: how to handle the plasma in the closed region
;        The options are:
;
;             HYDRO=2; radial hydrostatic equilibrium
;             HYDRO=3; Vasquez et al 2003 hydrostatic atmosphere
;
;  DENSPROF
;             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 FOR_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 FOR_DIPOLE program A1,A2,A3,A4,A5,aa,bb,alpha,beta
;                       defaults to Vasquez streamer values
;               Note- default=1 for all
;
;  TE
;             For HYDRO=3, temperature parameter 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 FOR_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
;               GAMMA1, GAMMA3 - Dipole/Octopole terms (units Gauss)
;                       DEFAULT GBL:  -10. -1.2
;                         (note, original paper had 1.2 for the octopole
;                           so this is not an exact reproduction
;                           for now we are going to insist on same sign G1, G3 for GBL normalization
;                           because this introduces additional null at equator instead of pole)
;                       FOR TESTING ONLY -- uncomment further down where gamma1,gamma3,bpar set
;                         also in gibbaglow.pro where NORM1, NORM3 are set
;                         DEFAULT Low 1986: 45*b*b, -4
;                         (note, original paper had 3*b*b,4 - we believe this was a typo)
;                         (also note, to reproduce in normalization of GBL, use gamma1=45, gamma3=360)
;
;               VSCALE - scale for radial velocity profile along the open field (varies with height)
;                           UNITS KM/SEC
;                           DEFAULT 0.d0
;