我试图从这个源(第88页)重新创建图6.2中的图形。他们通过使用以下速度和温度的解(都是从求解基本的平行边板斜坡问题中得到的)得到这些图。< / p > < p > <跨类=“math-container”> $ v_x = \压裂{2(\ρg sinα(\))^ {n}} {n + 1} (H ^ {n + 1} - (H-z) ^ {n + 1}) $ < / span > < br / > <跨类=“math-container”> $ T = T_s + \压裂{问^{\补}_ {geo}} {\ kappa} (H-z) + \压裂{2啊^ {n + 3}(\ρg sinα(\))^ {n + 1}}{\卡帕(n + 2)}[1 - \压裂{z} {H} \压裂{1}{n + 3}(\压裂{H-z} {H}) ^ {n + 3}] $ < / span > < / p > < p >他们说为了创建图6.2,他们使用以下值。李> < / p > < ul > < <跨类=“math-container”> H = 100美元< / span > < /李> <李> <跨类=“math-container”> $ \α= 10 ^{\保监会}$ < / span > < /李> <李> <跨类= >“math-container T_s美元= -10 ^{\保监会}加元< / span > < /李> <李> <跨类=“math-container”>美元问^{\补}_ {geo} = 50 W m ^ {2} $ < / span > < /李> <李> <跨类=“math-container”> n = 3美元< / span > < /李> <李> <跨类=“math-container”> = 10美元^{-16}一个^ {1}Pa ^{3} = 10 ^{-16}一个^{1}公斤^ {3}m ^ {3} s ^ {6} $ < / span > < /李> <李> <跨类=“math-container”> $ \ρ= 910公斤m^{-3}$

我注意到，当我把这些值代入上述方程时，我没有得到与图中相同的结果。我看到过其他有类似图形的来源，例如图2显示$v_x$here，所以我知道他们所做的是正确的，我觉得我错过了一些东西。我仔细检查了单位，他们对$v_x$方程正确排列。对于$T$方程，我注意到我需要先将$T_s$转换为开尔文(即添加273)，然后再插入方程，然后我需要将$T$方程的最终结果转换为摄氏度(即减去273)。这仍然没有给我一个与图6.2中提供的匹配图。

当我把这些值代入方程时，为什么这些数字不匹配?我已经通过了一堆数学和编码，它们看起来不错，所以我觉得我的问题源于对冰川动力学的误解。 Is there some trick that I'm missing? Thanks!

Here is the original Figure 6.2

Here is what I got

Here is the Matlab code I used to get my results in case that helps.

H=100; alpha = 10; A = 10e-16; rho = 910; g = 9.81; kappa = 2.1; n = 3; q_geo = 50; Ts = -10; figure() z_grid1 = linspace(0,100,401); v_x = (2.*A.*(1/(n+1)).*(rho.*g.*sin(alpha)).^n)*(H^(n+1) -(H-z_grid1).^(n+1)); %equation 6.13 plot(v_x,z_grid1) title('v_x graph') xlabel('velocity [m/a]') ylabel('Height above the base z [m]') figure() z_grid2 = linspace(0,100,401); T = (Ts+273) + (q_geo./kappa).*(H-z_grid2)+(2.*A.*H.^(n+3)).* (((rho.*g.*sin(alpha))^(n+1))./ (kappa.*(n+2))).*(1-(z_grid2 ./H) - (1/(n+3)).*((H-z_grid2)./H).^(n+3)); %equation 6.28 plot(T-273,z_grid2) title('T graph') xlabel('T [Degrees Celcius]') ylabel('Height above the base z [m]') 

所以你要么保持你的代码，但转换你的10°角在弧度(更改为alpha = 0.174533)，或者你可以使用信德函数，它也计算正弦，但以角度在角度作为输入。$ a $: $10^{-16}$应该被注意1e-16。这修复了第一个图形(速度配置文件):

第二看起来仍然很奇怪…你在$q^{\perp}_{geo}$的定义中也有一个错误:单位是mW，所以它应该转换为q_geo=0.05;。这修复了方程的第一行:如果将其余部分放在注释中，它将很好地再现线性虚线，直到$z$ = 0 m的温度为-7.619°C。所以问题一定在第二部分，虽然我不知道在哪里