import numpy as np import matplotlib.pyplot as plt # 设置模拟参数 D = 0.12# 热扩散率 L = 20 # 模拟深度 T = 365 * 10 # 模拟时间 h = 0.1 # 栅格大小 dt = 0.01 # 时间步长 N = int(T / dt) # 时间步数 M = int(L / h) + 1 # 深度格点数 K = int(T / 90*dt) # 画图次数 # 初始化温度分布 T0 = np.ones(M) * 10 T0[0] = 10 + 12 # 表面温度 T0[-1] = 11 # 深度为20米处温度 # 初始化温度矩阵 T = np.zeros((M, N)) T[:, 0] = T0 # 进行数值求解 for n in range(1, N): for i in range(1, M - 1): T[i, n] = T[i, n-1] + D * dt / h**2 * (T[i+1, n-1] - 2*T[i, n-1] + T[i-1, n-1]) # 边界条件 if not np.all(T[:,n]==0): T[0, n] =10 + 12 * np.sin(2 * np.pi * n * dt / T) T[-1, n] =11 else: break # 每隔90天画一次图 if n % int(T / 90*dt) == 0: plt.plot(T[:, n], np.linspace(0, L, M), label=f'{n*dt/365:.0f} year') plt.legend() # 显示温度轮廓图 plt.xlabel('Temperature (°C)') plt.ylabel('Depth (m)') plt.title('Temperature Profile') plt.gca().invert_yaxis() plt.show()代码错误，如何改正

v options
are turned off.

INVOCATION
A shell is a program that provides the command line (i.e., the all-
text user interface) that lets you interact with your operating sys‐
tem. The 代码中的错误是在计算画图次数 `K` 时出现了错误，应该将 `int(T / shell waits for user input and then executes the commands
entered on the command line. The command line 90*dt)` 改为 `K = int(T / 90 / dt)`。另外，在画图时应该先创建一个 typically includes the
name of the program being run, followed by any settings or options
needed图像对象，然后将所有的曲线添加到该对象中，最后显示图像。修改后的代码如下 to make the program work correctly, followed by the arguments to
be passed to the program.

Bash can be started：

import numpy as np
import matplotlib.pyplot as plt

# 设置模拟参数
D = 0.12  #  in  several  ways.  The most common way is to
       start it as an interactive shell.  Interactive means that bash  热扩散率
L = 20  # 模拟深度
T = 365 * 10  # 模拟时间
 is  run
       with  its standard input connected to the terminal and that it displays
       the shell prompt to let the user knowh = 0.1  # 栅格大小
dt = 0.01  # 时间步长
N = int(T / it is ready to read a command.

       Bash can also be started non-interactively by running a shell script.
       The  script  dt)  # 时间步数
M = int(L / h) + 1  # 深度格点数
K = int will  be  run  in  much  the  same way as a normal shell
       command, but the shell will exit when the(T / 90 / dt)  # 画图次数

# 初始化温度分布
T0 = np.ones(M) script  completes.   This  is
       how  you  would  typically  run  a  shell script that sets environment
       * 10
T0[0] = 10 + 12  # 表面温度
T0[-1] =  variables, for example.

       When Bash starts, it reads and executes the commands in  these  files,
       in  this order:

11  # 深度为20米处温度

# 初始化温度矩阵
T = np.zeros((M, N))
T       ·  /etc/profile

       ·  ~/.bash_profile

       ·  ~/.bash_login

       ·  ~/.profile

      [:, 0] = T0

# 进行数值求解
for n in range(1, N):
    for i in The  /etc/profile  file  is  read  once, when Bash is started as a login
       shell.  The other files range(1, M - 1):
        T[i, n] = T[i, n-1] + D * dt / h are read when Bash is started as  an  interac‐
       tive  shell.  If any of the files exist but cannot**2 * (T[i+1, n-1] - 2*T[i, n-1] + T[i-1, be read, Bash reports
       an error.  This can happen, for example, if one of the files is not  a
       n-1])

    # 边界条件
    if not np.all(T[:,n]==0):
        T[0, n] = regular  file  but a symbolic link that points to a file that does not
       exist.

       When a shell script is run, 10 + 12 * np.sin(2 * np.pi * n * dt / T)
        T[-1, n] =  Bash reads and executes the commands in the
       script.   When  the  script  is  done, Bash exits.  If a11
    else:
        break

    # 每隔90天画一次图
    if n % K == 0:
        command in the
       script fails for any reason, Bash exits with an error code.

       When Bash is run non-interactively to run plt.plot(T[:, n], np.linspace(0, L, M), label=f'{n*dt/365:.0f} year')

# a shell script, the shell  is
       not  started  as  a  login  shell by default.  This means that the 显示温度轮廓图
plt.xlabel('Temperature (°C)')
plt.ylabel('Depth (m)')
plt.title('Temperature Profile')
plt shell
       does not read the /etc/profile or ~/.profile files.  However, it  does
       read  the  ~/.bash.gca().invert_yaxis()
plt.legend()
plt.show()

运行后将得到正确的温度轮廓图。