isphion=1
  b0=5.         #=1 for normal direction B field
  rsigpl=1.e-8  #anomalous cross field conductivity
  cfjhf=1.      #turn on heat flow from current (fqp)
  cfjve=1.      #makes vex=vix-cfjve*fqx
  jhswitch=1    #Joule Heating switch
  newbcl=1      #Sheath boundary condition (bcee, i) from current equation
  newbcr=1
  isfdiax=1     #Factor to turn on diamagnetic contribution to sheath
  cfyef=1.0     #EXB drift in y direction
  cf2ef=1.0     #EXB drift in 2 direction
  cfybf=1.0     #turns on vy