Commit f813e57b authored by Lore, Jeremy's avatar Lore, Jeremy
Browse files

Bugfix for RK routines for dL integration that would not work for more than one system of equations

parent ad0373fb

matlab/bfield_library_jdl/fieldline_following/dl/fl_derivs_dl_g3d.m

+1 −1
Original line number Diff line number Diff line
@@ -5,7 +5,7 @@ end 
N = length(RPZ);



% Axisym part

[Bout,ierr_Bas] = bfield_geq_bicub(g,RPZ(1:2:N-2),RPZ(3:2:N),nowarn);

[Bout,ierr_Bas] = bfield_geq_bicub(g,RPZ(1:3:N-2),RPZ(3:3:N),nowarn);

if ierr_Bas == 1

    if ~nowarn

        warning('AS bfield error in fl_derivs_phi_g3d')

matlab/bfield_library_jdl/fieldline_following/dl/fl_derivs_dl_gfile.m

+1 −1
Original line number Diff line number Diff line
@@ -5,7 +5,7 @@ end 
N = length(RPZ);



% Axisym part

[Bout,ierr_Bas] = bfield_geq_bicub(bfield.g,RPZ(1:2:N-2),RPZ(3:2:N),nowarn);

[Bout,ierr_Bas] = bfield_geq_bicub(bfield.g,RPZ(1:3:N-2),RPZ(3:3:N),nowarn);

if ierr_Bas == 1

    if ~nowarn

        warning('AS bfield error in fl_derivs_dl_gfile')