Pro g3bCut, gfile
; ~tdet/HHMS/SSBC/g3blook.pro
; //saturn/~tdet/HHMS/SSBC/g3bstribs.pro previously
; //europa/~tdet/HHMS/SSurf/RnD/stribs.pro previously
; //europa/~tdet/SSurf/MapStats/stribs.pro previously
; Date: 1999 12 22
; make scatter plots of params at inner and outer limits of grid
; Note: B is in units of deci-gauss (dg), 1 dg = 10^-5 T = 10^4 nT
; Example: g3bCut, 'g3bss.1952_005.wso'
; Example: g3bCut, 'g3bss.1951_030.wso'

COMMON BLOC0, Ro0,Vr0,Vt0,Vp0,Br0,Bt0,Bp0,Tp0
COMMON BLOC1, Ro, Vr, Vt, Vp, Br, Bt, Bp, Tp
COMMON BLOC2, Ltime, drx, dtx, dpx
COMMON BLOC3, Km, Jm, Im
COMMON BLOCX, RC

readg3b, gfile

Rs = 6.95E5                ; ( km )
Rad = 21.5 + drx*Findgen(Im)
Plot, Rad, RC, psym=-4
Stop, 'move on?'

The = 42.5 + dtx*Findgen(Jm)
Phi = 2.5 + dpx*Findgen(Km)
AU = 90    ; index of Rad corresponding to 1 AU (= 215 Rs)
;AU = 0    ; index of Rad corresponding whatever
Fk = 0
Lk = Km-1
Fj = 0
Lj = Jm-1

Nz = lonarr(Im)
For i=0,Im-1 Do Begin
  zot = where(Tp(*,*,I) le 0., Nzot)
  Nz(i) = Nzot
EndFor

!P.Charsize = 1.4
PlotLoop:
  Plot, The,Vp(Fk+1,*,AU), Title=gfile $
      , xTitle='Theta', xRange=[30,150], xTicks=8,xMinor=3, /xStyle $
      , yTitle='Vp 1AU,Fk & Lk', yRange=[-10,10] $
      , psym=6,symsize=.4
  oPlot, The,Vp(Lk-1,*,AU), psym=7,symsize=.4
  oPlot, [90,90],[-15,15], linestyle=1

  If (!d.name eq 'X') Then Stop, 'Vr?'
  Plot, The,Vr(Fk,*,AU), Title=gfile $
      , xTitle='Theta', xRange=[30,150], xTicks=8,xMinor=3, /xStyle $
      , yTitle='Vr 1AU,Lk', yRange=[200,800] $
      , psym=6,symsize=.4
  oPlot, The,Vr(Lk,*,AU), psym=7,symsize=.4
  oPlot, [90,90],[200,800], linestyle=1

  If (!d.name eq 'X') Then Stop, 'Vt?'
  Plot, The,Vt(Fk,*,AU), Title=gfile $
      , xTitle='Theta', xRange=[30,150], xTicks=8,xMinor=3, /xStyle $
      , yTitle='Vt 1AU,Lk', yRange=[-10,10] $
      , psym=6,symsize=.4
  oPlot, The,Vt(Lk,*,AU), psym=7,symsize=.4
  oPlot, [90,90],[-20,20], linestyle=1
  oPlot, [0,180],[0,0], linestyle=1

  If (!d.name eq 'X') Then Stop, 'Ro?'
  Plot_io, The,Ro(Fk,*,AU), Title=gfile $
         , xTitle='Theta', xRange=[30,150], xTicks=8,xMinor=3, /xStyle $
	 , yTitle='Ro (1AU)', yRange=[1E-12,1E-8] $
         , psym=6,symsize=.2
  oPlot, The,Ro(Lk,*,AU), psym=7,symsize=.2
  oPlot, [90,90],[1E-13,1E-8], linestyle=1

  If (!d.name eq 'X') Then Stop, 'Br?'
  Plot, The,1E4*Br(Fk,*,AU), Title=gfile ,psym=6,symsize=.2 $
      ,xTitle='Theta', xRange=[30,150], xTicks=8,xMinor=3, /xStyle $
      ,yTitle='Br (nT,1AU)', yRange=[-5,5], /yStyle
  oPlot, The,1E4*Br(Lk,*,AU), psym=7,symsize=.2
  oPlot, [90,90],[-15,15], linestyle=1
  oPlot, [0,180],[0,0], linestyle=1

  If (!d.name eq 'X') Then Stop, 'Bt?'
  Plot, The,1E4*Bt(Fk,*,AU), Title=gfile $
      , xTitle='Theta', xRange=[30,150], xTicks=8,xMinor=3, /xStyle $
      , yTitle='Bt (nT,1AU)', yRange=[-.3,.3], /yStyle $
      , psym=6,symsize=.2
  oPlot, The,1E4*Bt(Lk,*,AU), psym=7,symsize=.2
  oPlot, [90,90],[-15,15], linestyle=1
  oPlot, [0,180],[0,0], linestyle=1

  If (!d.name eq 'X') Then Stop, 'Bp?'
  Plot, The,1E4*Bp(Fk,*,AU), Title=gfile $
      , xTitle='Theta', xRange=[30,150], xTicks=8,xMinor=3, /xStyle $
      , yTitle='Bp (nT,1AU)', yRange=[-5,5], /yStyle $
      , psym=6,symsize=.2
  oPlot, The,1E4*Bp(Lk,*,AU), psym=7,symsize=.2
  oPlot, [90,90],[-15,15], linestyle=1
  oPlot, [0,180],[0,0], linestyle=1

  If (!d.name eq 'X') Then Stop, 'Bs?'
  Plot, The,1E4*(Br(Fk,*,AU)-Bp(Fk,*,AU)), Title=gfile $
      , xTitle='Theta', xRange=[30,150], xTicks=8,xMinor=3, /xStyle $
      , yTitle='Bs (nT,1AU)', yRange=[-10,10], /yStyle $
      , psym=6,symsize=.2
  oPlot, The,1E4*(Br(Lk,*,AU)-Bp(Lk,*,AU)), psym=7,symsize=.2
  oPlot, [90,90],[-15,15], linestyle=1
  oPlot, [0,180],[0,0], linestyle=1

  Stop, 'Replay?'
GoTo, PlotLoop

End