|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
calcwi.f里面的50行
ca(kk) = max(0.,rcl*qr(k+2)/ri(kk))
我认为错了,应该为 ca(kk) = max(0.,rcl*qr(k+2)/wm(kk)*1000)
ff = ener*dfa(nn)*tcl**dft(nn)*exp(freq*dfe(nn)/tcl)*tbc
rr = ener*dba(nn)*tcl**dbt(nn)*exp(freq*dbe(nn)/tcl)*tbc
do k=1,nst
kk = kmap(k)
ff = ff*ca(kk)**coff(kk,nn)
bb = bb*ca(kk)**cofb(kk,nn)
上面的rr应该为bb
还有reactck.f里面的subroutine fun的变量gamdt没有赋值就使用了
|
|