热弹流润滑的温度计算问题

charaderrrrr

黄平老师的热弹流点接触程序运行后，发现温度输出都为0遂找温度T的相关参数倒推：子程序TRA3中的D（I），B（I,K）到TCAL中的A（I,K），R(K)一直为1,  A2和A5在第一次循环后趋于零...
发现这中间值都在第一次循环后归零或者趋于零，并且循环中温度各梯度没有变化
TRA3中储存节点温度的D输出一直是不对的，感觉是整个TRA3的问题
但对照了一下能量方程的撰写，除了用的是旧的语法格式，其他好像没什么问题

所以想请教一下，后面循环置零的原因是什么 和 温度梯度内没有变化的原因，以及相应的修改建议
并且在其他程序不变的情况下，添加输出语句看中间值时，会发生输出误差变化甚至NaN的情况，也想请问一下原因