Corrected the sign in turn with sign for vertical stability current

matlab/IMAS/data_coils.m

@@ -12,7 +12,7 @@ Z = [ 5.4640; 3.2780; 1.0920; -1.0720; -3.2580; -5.4440;...
 dR = [0.7400; 0.7400;0.7400;0.7400;0.7400;0.7400;0.9590;0.5801;0.6963;0.6382;0.8125;1.5590;0.132;0.132];
 dZ = [2.0930; 2.093; 2.093; 2.093; 2.093; 2.093; 0.9841; 0.7146; 0.9538; 0.9538; 0.9538; 1.1075; 0.132;0.132];
 
-N = [554;554;554;554;554;554;248.6;115.2;185.9;169.9;216.8;459.4; 4; 4];
+N = [554;554;554;554;554;554;248.6;115.2;185.9;169.9;216.8;459.4; 4; -4];
 
 N = round(N);
 

matlab/IMAS/gen_filament.m

@@ -82,7 +82,7 @@ for ii = 1:numel(AC) % loop all the coils given vvstructure
     
     % Try to use a number of filaments that respect the shape of the parallelogram
     ratio = sqrt((AD*AD'))/sqrt((AB*AB'));
-    num1 = ceil(sqrt(num(ii)*ratio));
+    num1 = ceil(sqrt(abs(num(ii))*ratio));
     num2 = ceil(num1/ratio);
     
     pointx = zeros(num1*num2, 1);
@@ -102,7 +102,7 @@ for ii = 1:numel(AC) % loop all the coils given vvstructure
     % but will not necessarily generate a number of windings = number of
     % turns. The following normalization garanties that the total current
     % flowing in the cross section of coil Iw*nw = Ia*nturns respects the input data.
-    Twc(end+1:end+num1*num2,end+1) = 1*sign(num(ii))/(num1*num2)*num(ii);
+    Twc(end+1:end+num1*num2,end+1) = 1*abs(num(ii)/(num1*num2))*sign(num(ii));
 end
 
 % Reshape filaments(windings) into unique vector