第一循环：calculateh.f

第一次运行循环时，四面体网格仍然是未细化的网格。在 calculateh.f 中执行以下操作：

• 确定每条边的长度并将其存储在字段 d(*) 中
• 对于每个节点，确定其所属边的平均长度。将此值乘以用户指定的极限值，并除以用户定义准则的实际值。将结果存储在字段 h(*) 中。例如，如果用户选择了应力准则，极限值为 50，节点中的实际应力值为 200，则该节点中 0.1 的平均边长将导致该节点中的值为 0.025，即局部所需的边长为 0.025。注意，细化信息存储在 filab(87) 中，这是一个长度为 87 的字符字符串。在位置 1:2 存储"RM"，应用细化准则的字段存储在位置 3:6，极限值存储在位置 7:26（f20.0），要细化的单元集存储在位置 27:87。
• 确定整个网格中 d(*) 的最小值并将其存储在 $d_{min}$ 中。