pcOpacP = 0.5
if (Quiet):
LaunchNowin()
else:
Launch()
OpenDatabase(fIn)
md = GetMetaData(fIn) #Metadata
#Import definitions
pyv.lfmExprsEB()
#Start plotting
pyv.setAtts() #Some defaults
Nphi = len(PhiCs)
radScl = np.pi / 180.0
#Create points
for k in range(Nphi):
PhiC = PhiCs[k]
LatC = LatCs[k]
fOut = "Slc.P%d.png" % np.int(PhiC)
#Generate seed points
Cp = np.cos(PhiC * radScl)
Sp = np.sin(PhiC * radScl)
ThCp = 90 - LatC #Critical latitudes in theta
w3d.windowValid = 0
else:
#Doing H+
w3d.viewNormal = (0.973617, 0.0921289, 0.208763)
w3d.focus = (-1.025, 0, 0.00223994)
w3d.viewUp = (-0.209201, -0.00498575, 0.97786)
w3d.viewAngle = 30
w3d.parallelScale = 25.6503
w3d.nearPlane = -51.3007
w3d.farPlane = 51.3007
w3d.imagePan = (-0.0645008, 0.0175247)
w3d.imageZoom = 3.7975
w3d.perspective = 1
w3d.eyeAngle = 2
w3d.centerOfRotationSet = 0
w3d.centerOfRotation = (-1.025, 0, 0.00223994)
w3d.axis3DScaleFlag = 0
w3d.axis3DScales = (1, 1, 1)
w3d.shear = (0, 0, 1)
w3d.windowValid = 1
SetView3D(w3d)
DrawPlots()
#pyv.cleanLegends(plXs,plYs,plTits)
pyv.setAtts()
SaveWindow()
if (not Quiet):
OpenGUI()
v.imageZoom = 1.0
SetView3D(v)
light = LightAttributes()
light.type = 1 # object
light.direction = (0, -0.75, -1)
light.brightness = 0.3
SetLight(1, light)
pyv.genPolarCoords2(dbs[2], 2, 10, 45, contours=True)
speciesStr = "Species: " + species
pyv.genTit(t**S=speciesStr, Pos=(0.7, 0.86), height=0.02)
pyv.genTit(t**S="Energies: All", Pos=(0.7, 0.89), height=0.02)
tit = pyv.genTit(t**S=title, Pos=(0.25, 0.955))
pyv.cleanLegends(plXs, plYs, plTits, plHt=0.015)
pyv.setAtts(pWidth=1000, VidDir=tmpvid)
DrawPlots()
md = GetMetaData(dbs[2])
dt = md.times[1] - md.times[0]
pyv.doTimeLoop(Ninit=0,
Nfin=None,
T0=T0,
dt=dt,
Save=True,
tLabPos=(0.35, 0.05))
outVid = "vid_r_phi_mag_" + varNames[vi] + "_all.mp4"
pyv.makeVid(outVid=outVid, vidDir=tmpvid, overwrite=True)