大西洋中脊横截剖面轮廓形状如上右图红色曲线所示。对于某处剖面,以中脊对称轴向上方向为 y 轴、以海平面为 x 轴建立直角坐标系,则轮廓线上某点 A 的横坐标 x 表示 A 离对称轴的距离(单位为千米),而纵坐标 y 表示 A 点的海拔(单位为千米,一般为负数),x 和 y 存在如下关系:
其中,∣x∣ 表示求 x 的绝对值,D 为对称轴处离海面的距离(单位为 km)、C 为给定的常量、u 为海底减薄速率(单位为 km/Myr),假设 D=2.5km、C=0.11、u=50km/Myr,本关任务是绘制这种情况下的剖面轮廓线。
import matplotlib.pyplot as plt
import math
########## Begin ##########
D,C,U=2.5,0.11,50
X=[]
Y=[]
d,delta=1200,200
for x in range(-d,d+1,delta):
X.append(x)
y=-(D+C*(10*abs(x)/U)**0.5)
Y.append(y)
########## End ##########
#以下代码用于打印结果
for i in range(len(X)):
print('%5s\t%.2f' % (int(X[i]),Y[i]))
#以下代码用于保存图片
plt.plot(x,y,'r-')
plt.axis([-d,d,-4.2,0])
plt.savefig( 'src/step7/student/pic.png' )
plt.close()