在解决表面磁性水平z轴旋转问题时，可以使用各种编程语言和算法来实现。下面是一个使用Python编程语言的示例代码，其中使用了numpy库和matplotlib库来进行计算和可视化。

import numpy as np
import matplotlib.pyplot as plt

# 创建一个表示磁场的向量
magnetic_field = np.array([0, 0, 1])

# 创建一个表示表面的点的坐标数组
points = np.array([[1, 0, 0], [0, 1, 0], [-1, 0, 0], [0, -1, 0]])

# 创建一个表示旋转角度的数组
angles = np.linspace(0, 2*np.pi, 100)

# 初始化一个空的结果数组
results = np.zeros((len(angles), len(points)))

# 对于每个角度和每个点，计算旋转后的磁场向量在z轴上的投影
for i, angle in enumerate(angles):
    rotation_matrix = np.array([[np.cos(angle), -np.sin(angle), 0],
                                [np.sin(angle), np.cos(angle), 0],
                                [0, 0, 1]])
    rotated_field = np.dot(rotation_matrix, magnetic_field)
    results[i] = np.dot(rotated_field, points.T)

# 可视化结果
plt.plot(angles, results[:, 0], label='Point 1')
plt.plot(angles, results[:, 1], label='Point 2')
plt.plot(angles, results[:, 2], label='Point 3')
plt.plot(angles, results[:, 3], label='Point 4')
plt.xlabel('Rotation Angle (radians)')
plt.ylabel('Z-axis Projection')
plt.legend()
plt.show()

这段代码首先创建了一个表示磁场的向量magnetic_field和一个表示表面的点的坐标数组points。然后，通过使用旋转矩阵对磁场向量进行旋转，并计算旋转后的磁场向量在z轴上的投影。最后，通过使用matplotlib库将结果可视化出来。

请注意，这只是一个示例代码，具体的解决方法可能会因问题的具体要求而有所不同。