Thanks Dave.
Another way to find the variable in the code (the history field doesn't always have the same name as the variable in the code) is to find the history call for Z0MV (by grep or other means):
call hist_addfld1d (fname='Z0MV', units='m', &
avgflag='A', long_name='roughness length over vegetation, momentum', &
ptr_patch=this%z0mv_patch, default='inactive')
So, the variable in the code is z0mv_patch and you can find instances of this using grep or other means.
The z0mv is calculated first in biogeophys/CanopyTemperatureMod.F90:
z0m(p) = z0mr(patch%itype(p)) * htop(p)
z0mv(p) = z0m(p)
and then as Dave noted modified for canopy density in biogeophys/CanopyFluxesMod.F90:
z0mv(p) = exp(egvf * log(z0mv(p)) + (1._r8 - egvf) * log(z0mg(c)))