Removes drift dependency on toroidal current dir

Upon closer inspection/thought, the toroidal current induces a poloidal field. Reversing the direction of the toroidal current reverses the direction of the poloidal B-field. The impact of the drift flows is, then ExB = (Erad + Epol)x(Btor+Bpol) = Erad x Btor + Epol x Btor + Erad x Bpol + Epol x Bpol (-> 0). Thus, the only impact of reversing Itor on the drifts is a reversal of the toroidal shear flow. However, UEDGE assumes toroidicity, so this contribution is not considered.

bbb/odesetup.m

@@ -1113,12 +1113,12 @@ call syld96(matt,matp,cion,cizb,crmb)
 
 c...  rescale magnetic field quantities with b0
       call s2scal (nx+2, ny+2, abs(b0/b0old), btot, 1, nx+2)
-      call s2scal (nx+2, ny+2, sign(1.,upbparadir*b0/b0old), rbfbt, 1, nx+2)
-      call s2scal (nx+2, ny+2, upbparadir*b0old/b0, rbfbt2, 1, nx+2)
+      call s2scal (nx+2, ny+2, sign(1.,b0/b0old), rbfbt, 1, nx+2)
+      call s2scal (nx+2, ny+2, b0old/b0, rbfbt2, 1, nx+2)
       call s2scal (nx+2, ny+2, abs(b0/b0old), rbpol, 1, nx+2)
-      call s2scal (nx+2, ny+2, upbparadir*b0old/b0, curvrby, 1, nx+2)
+      call s2scal (nx+2, ny+2, b0old/b0, curvrby, 1, nx+2)
       call s2scal (nx+2, ny+2, abs(b0old/b0), curvrb2, 1, nx+2)
-      call s2scal (nx+2, ny+2, upbparadir*b0old/b0, gradby, 1, nx+2)
+      call s2scal (nx+2, ny+2, b0old/b0, gradby, 1, nx+2)
       call s2scal (nx+2, ny+2, abs(b0old/b0), gradb2, 1, nx+2)
          # v2cb sign change comes from rbfbt in forming uu
       b0old = b0