cb - radial Grad B drift
  cf2bf=1.0     #turns on v2cb - perp. Grad B drift (nearly pol)
  cfqybf=1.0    #turns on vycb contrib to radial current
  cfq2bf=1.0    #turns on v2cb contrib to perp("2") current
  cfydd=0.0     #turns off divergence free diamagnetic current
  cf2dd=0.0     #turns off divergence free perp diagmatic current
  cfqybbo=0     #turn off Grad B current on boundary
  cfqydbo=1     #use full diagmagetic current on boundary to force j_r=0
  cfniybbo=1.   # use to avoid artificial source a