Pro g3bSlice, gfile
; ~tdet/HHMS/igmvw/g3bslice.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: g3bSlice, 'grids/g3b.000.014'

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)
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

Window, /Free, Title='Meridional Slice', xSize=600, ySize=600

R2Ro = Ro
For I=0,Im-1 Do R2Ro(*,*,I) = Rad(I)^2*Ro(*,*,I)*Vr(*,*,I)^2
;RTp = Tp
;For I=0,Im-1 Do RTp(*,*,I) = Rad(I)*Tp(*,*,I)


Qp = Reform(R2Ro(36,*,*),Jm,Im)
  minQp = min(Qp)
  maxQp = max(Qp)
  Print, minQp, maxQp

PlotLoop:
  Polar_TV, Qp,The,Rad, pos=1
  Stop, 'Replay?'
GoTo, PlotLoop

End