distBot)             call set( Fwall, i, 0.0 + E0*( ((node_val(ll,i)/kappa)*( LLL ))**0 ))                end do     case ("Burchard00")         do i=0,nNodes             distTop = max(0.0,node_val(top,i))             distBot = max(0.0,node_val(bottom,i))             LLL = 0.0 / 