The line:
nm=l-1
will produce an out of bound array element access

Since l varies from k to 0, we have the line:
if (fabs(w[nm])+anorm == anorm) break;
accessing wm[k-1] to wm[-1]