_RL xkar
& xkar, stefan,
WRITE(iUnit,*) 'BlkF: xkar =',xkar
& Tf0kel, Rgas, xkar, stefan,
xkar = 0.4 _d 0
czol = zref*xkar*gravity
rdn=xkar/(log(zref/zice))
rd = rdn/(1. _d 0 + rdn*(aln-psimh)/xkar)
rh = rhn/(1. _d 0 + rhn*(aln-psixh)/xkar)
usn = ws/(one + rdn/xkar*(zwln-psimh) )
rd = rdn/(1. _d 0 + rdn*(zwln-psimh)/xkar)
rh = rhn/(1. _d 0 + rhn*(ztln-psixh)/xkar)
re = ren/(1. _d 0 + ren*(ztln-psixh)/xkar)