Commit 6f63e1bb authored by sbiri's avatar sbiri
Browse files

changed temp to tmp in line 947-950 not to confuse temporary var. for temperature

parent 843eb946
......@@ -944,10 +944,10 @@ def get_L(L, lat, usr, tsr, qsr, hin, Ta, sst, qair, qsea, wind, monob, zo,
psim_calc(hin[1]/monob, meth)))
Rb = g*dthv*hin[1]/(tv*uz*uz)
if (L == "tsrv"):
temp = (g*kappa*tsrv /
tmp = (g*kappa*tsrv /
np.maximum(np.power(usr, 2)*Ta*(1+0.6077*qair), 1e-9))
temp = np.minimum(np.abs(temp), 200)*np.sign(temp)
monob = 1/np.copy(temp)
tmp = np.minimum(np.abs(tmp), 200)*np.sign(tmp)
monob = 1/np.copy(tmp)
elif (L == "Rb"):
zol = (Rb*(np.power(np.log((hin[1]+zo)/zo)-psim_calc((hin[1]+zo) /
monob, meth) +
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment