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