surf=0;


while(surf<TheSystem.LDE.NumberOfSurfaces)


    if(TheSystem.LDE.GetSurfaceAt(surf).Comment=="TEST")


       surc=surf;


       break


    end


    surf=surf+1;


end

    max_field = -1E9;


    for i = 1:TheSystem.SystemData.Fields.NumberOfFields


       X_field= TheSystem.SystemData.Fields.GetField(i).X;


       Y_field= TheSystem.SystemData.Fields.GetField(i).Y;


       radial_field= sqrt(X_field^2+Y_field^2);


       if (radial_field> max_field) 


           max_field = radial_field; max_field_num=i;


       end


    end


    fprintf('%d    %f\n',max_field_num,max_field);