/com, add external poroidal B field to coil and air
B_Source = arg1		! Define the B, B_mag
allsel,all
csys,0
esel,s,type,,2
nsle,s
nsel,r,loc,Y,-57*INM,-56*INM
n_str=0
*get,n_count,node,0,count
*do,i,1,n_count,1
	  *get,n_next,node,n_str,NXTH
		*get,n_x,node,n_next,loc,x
		*get,n_z,node,n_next,loc,z
		D,n_next,AX,(n_z*B_Source/2)
		D,n_next,AZ,(-n_x*B_Source/2)
		D,n_next,AY,0.0
		n_str=n_next
*enddo
allsel,all
csys,0
esel,s,type,,2
nsle,s
nsel,r,loc,Y,56*INM,57*INM
!csys,12
!nsel,u,loc,y,-1e-5,1e-5
!nsel,u,loc,y,120-1e-5,120+1e-5
!csys,0
n_str=0
*get,n_count,node,0,count
*do,i,1,n_count,1
	  *get,n_next,node,n_str,NXTH
		*get,n_x,node,n_next,loc,x
		*get,n_z,node,n_next,loc,z
		D,n_next,AX,(n_z*B_Source/2)
		D,n_next,AZ,(-n_x*B_Source/2)
		D,n_next,AY,0.0
		n_str=n_next
*enddo