Assign and use t_wm: temporal filtering for temperature
Allow y0 != z_os (thermal surface roughness length)
y0 != z_os
Why does u_star -> NaN? Only for statistics: compute u_star assuming log-law
add ffy component and rotation
map2d_{xm1,ym1}