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);