function indxrjunk (nrho, ntheta, nchi, irho, itheta, ichi) IMPLICIT NONE INTEGER :: indxrjunk, nrho, ntheta, nchi, irho, itheta, ichi indxrjunk = (irho-1)*ntheta*nchi + (itheta - 1) * nchi + ichi RETURN END function indxrjunk