Java转Python气象计算代码结果不一致，求助排查循环与参数问题

问题背景 #

现有一段可正常输出预期结果的气象计算Java代码，转换为Python后运行结果与原代码不符，主要表现为rv2等数值变化异常，与Java结果差异明显。

核心错误点分析 #

1. 参数P取值错误：Java中P=100000，Python误写为1000000，直接影响分子扩散系数Dp的计算，导致后续所有依赖该参数的结果出现系统性偏差。
2. 循环逻辑不一致：Java使用for(i=0;i<N;i++)遍历，i从0到437（共438次循环）；Python原代码的while循环先执行i +=1，导致i从1开始，且第一个循环结束后i=438，第二个循环因i<N不成立直接跳过，同时时间t的计算完全错位。
3. 变量引用错误：第二个循环计算RR2时，误用了第一个循环的rl变量，应该使用当前循环的rl2，导致rv2计算逻辑错误。
4. 输出格式不匹配：Java用格式化输出控制数值精度和格式，Python原代码直接打印，导致结果格式不一致，不利于对比验证。

修正后的Python代码 #

import math

N = 438
tmin = 0
dt1 = 47.25
dt2 = 10.22
E0 = 610.78  
S = 1
r0l = 2e-6
r0v = 4e-6
nopt1 = 1.1e+3
nopt2 = 1.1e+4
T = 268.15 
b = 0.62
nv = 1e+9
pl = 917
pv = 1000
Rp = 461
P = 100000  # 修正参数取值，与Java一致

# 初始计算部分，与Java逻辑对齐
Ev = E0 * math.pow(10, ((7.63 * (T - 273.15)) / (T - 31.25)))
El = E0 * math.pow(10, ((9.5 * (T - 273.15)) / (T - 7.65)))
Dp = (-1.89725e-5 + 1.5e-7 * T) * (101326 / P)
a1 = math.pow(r0v, 3) * ((pv * nv) / (pl * nopt1)) + math.pow(r0l, 3)
a2 = math.pow(r0v, 3) * ((pv * nv) / (pl * nopt2)) + math.pow(r0l, 3)
r1 = math.pow(a1, 1.0/3.0)
r2 = math.pow(a2, 1.0/3.0)
s1 = (math.pow(r1, 2) - math.pow(r0l, 2)) * pl * Rp * T
s2 = (math.pow(r2, 2) - math.pow(r0l, 2)) * pl * Rp * T
y = 2 * Dp * Ev * (1 - (El / Ev))
t1 = s1 / y
t2 = s2 / y

# 第一个循环：对应Java的nopt1遍历逻辑
with open("outfile1.txt", "w") as out1:
    for i in range(N):
        t = tmin + dt1 * i
        rl = math.pow(((r0l*r0l) + ((2*Dp*Ev)/(pl*Rp*T))*(S-(El/Ev))*t), 0.5)
        RR1 = (math.pow(r0v, 3)*pv*nv + pl*nopt1*(math.pow(r0l, 3)-math.pow(rl, 3))) / (pv*nv)
        rv1 = math.pow(RR1, 1.0/3.0)
        L1 = b / (nv*math.pow(rv1, 2) + nopt1*math.pow(rl, 2))
        # 与Java一致的格式化输出
        output_line = "%.2e\t%.2e\t%.2f\t%.2f\n" % (rv1, rl, t, L1)
        print(output_line.strip())
        out1.write(output_line)

# 第二个循环：对应Java的nopt2遍历逻辑
with open("outfile2.txt", "w") as out2:
    for i in range(N):
        t = tmin + dt2 * i
        rl2 = math.pow(((r0l*r0l) + ((2*Dp*Ev)/(pl*Rp*T))*(S-(El/Ev))*t), 0.5)
        RR2 = (math.pow(r0v, 3)*pv*nv + pl*nopt2*(math.pow(r0l, 3)-math.pow(rl2, 3))) / (pv*nv)
        rv2 = math.pow(RR2, 1.0/3.0)
        L2 = b / (nv*math.pow(rv2, 2) + nopt2*math.pow(rl2, 2))
        # 与Java一致的格式化输出
        output_line = "%.2e\t%.2e\t%.2f\t%.2f\n" % (rv2, rl2, t, L2)
        print(output_line.strip())
        out2.write(output_line)

修正说明 #

1. 修正P的取值为100000，确保Dp计算与Java一致。
2. 使用for i in range(N)替代while循环，完全对齐Java的循环逻辑，i从0到437，共438次迭代。
3. 第二个循环中修正变量引用，用rl2参与RR2的计算，确保逻辑正确。
4. 使用Python的格式化字符串，与Java的printf输出格式保持一致，便于结果对比。
5. 使用with语句管理文件，自动关闭文件，比原Java的手动关闭更简洁安全。

内容的提问来源于stack exchange，提问作者Dova