SUBROUTINE  gradconj(  u,  b,  p,  q,  r,  it_max,  prec,  itgc,  rnorm_k,  rnudt2  )
      USE numerics
      IMPLICIT  NONE
      REAL(rp),  DIMENSION(sx-1:ex+1,sy-1:ey+1,sz-1:ez+1),  INTENT(INOUT)  ::  u
      REAL(rp),  DIMENSION(sx-1:ex+1,sy-1:ey+1,sz-1:ez+1),  INTENT(INOUT)  ::  b
      REAL(rp),  DIMENSION(sx-1:ex+1,sy-1:ey+1,sz-1:ez+1),  INTENT(INOUT)  ::  p
      REAL(rp),  DIMENSION(sx-1:ex+1,sy-1:ey+1,sz-1:ez+1),  INTENT(INOUT)  ::  q
      REAL(rp),  DIMENSION(sx-1:ex+1,sy-1:ey+1,sz-1:ez+1),  INTENT(INOUT)  ::  r
      REAL(rp),  INTENT(IN)   ::  prec
      REAL(rp),  INTENT(OUT)  ::  rnorm_k
      INTEGER,   INTENT(IN)   ::  it_max
      INTEGER,   INTENT(OUT)  ::  itgc
      REAL(rp),  INTENT(IN)   ::  rnudt2
      END  SUBROUTINE  gradconj