loc_icfl=max(abs(utrans(i,j,bi,bj)),
loc_icfl=max(loc_icfl,max(abs(vtrans(i,j,bi,bj)),
if (loc_icfl.gt.max_icfl) then
max_icfl = loc_icfl
_RL thick_bd, uflux, vflux, max_icfl, loc_icfl