c* Arrays dynamically allocated in this routine: c* c* real w (ia(_JT_nrows_)) c* real p (ia(_JT_nrows_)) c* real r (ia(_JT_nrows_)) c* real xold (ia(_JT_nrows_))