[Har] =  Hardat_to_Normal_and_Skew_dev(3, 16)

stuct_dim = size(Har);


for i = 1:stuct_dim(2)
    norm1(i) = Har(i).norm(1);
    skew1(i) = Har(i).skew(1);
     norm2(i) = Har(i).norm(2);
     skew2(i) = Har(i).skew(2);
     main(i) =  Har(i).Imag(1);
     x_field_trim(i)  =  Har(i).Imag(2);
     y_field_trim(i)  =  Har(i).Imag(3);
end

x_field = norm1 - norm1(1);
y_field = skew1 - skew1(1);



