Pro g3bShell, gfile
; ~tdet/HHMS/igmvw/g3bshell.pro           2000-03-23
; ~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
; make plots of parameters on spherical shell(s)
; Note: B is in units of deci-gauss (dg), 1 dg = 10^-5 T = 10^4 nT
; Example: g3bShell, 'g3b.lnk'
; Example: g3bShell, '../igmvw/grids/g3b.005.004'

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

Rgas = 8.32E-3
Rs = 6.95E5                ; ( km )
Rad = 21.5 + drx*Findgen(Im)
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

Pb = (Br^2 + Bt^2 + Bp^2)/(8.*!pi)
Pp = Ro*Tp*Rgas
Pt = Pb + Pp
Print, max(Pb(*,*,0)), max(Pp(*,*,0)), max(Pt(*,*,0)), max(Tp(*,*,0))
Print, avg(Pb(*,*,0)), avg(Pp(*,*,0)), avg(Pt(*,*,0)), avg(Tp(*,*,0))
Print, min(Pb(*,*,0)), min(Pp(*,*,0)), min(Pt(*,*,0)), min(Tp(*,*,0))

Stop, 'Tp shell by shell?'
Window, /Free, Title='Shell Shuffle', xSize=10*Km, ySize=12*Jm

Qp = Pt
Buf = bytarr(10*Km,12*Jm,Im)
For I=0,Im-1 Do Begin
  minQp = min(Qp(*,*,I))
  maxQp = max(Qp(*,*,I))
  Print, I, minQp, maxQp
  Img = BytScl(Qp(*,*,I),min=minQp,max=maxQp)
  Buf(*,*,I) = ReBin(Img,10*Km,12*Jm)
  TV, Buf(*,*,I)
  xyOuts, .95,.98, strcompress(string(I+1)), color=0
EndFor

PlotLoop:
For I=0,Im-1 Do Begin
  TV, Buf(*,*,I)
  xyOuts, .96,.96, strcompress(string(I+1)), color=0
  If(I lt 3) Then wait, 1
  If(I lt 6) Then wait, .5
EndFor

  Stop, 'Replay?'
GoTo, PlotLoop

End