c c ========================================================== subroutine flgout1ac(q,mx,lb,ub,qo,mxo,lbo,ubo,aux, & lbr,ubr,shaper,meqn,maux,nc,t) c ========================================================== c implicit double precision(a-h,o-z) integer meqn, mx, mxo dimension q(meqn,mx), qo(mxo), aux(maux,mx) integer lb(1), ub(1), lbo(1), ubo(1), lbr(1), ubr(1), shaper(1), & mresult, stride, imin(1), imax(1), i, getindx, d c stride = (ub(1) - lb(1))/(mx-1) imin(1) = max(lb(1), lbr(1)) imax(1) = min(ub(1), ubr(1)) if (mod(imin(1)-lb(1),stride) .ne. 0) then imin(1) = imin(1) + stride - mod(imin(1)-lb(1),stride) endif imin(1) = getindx(imin(1), lb(1), stride) if (mod(imax(1)-lb(1),stride) .ne. 0) then imax(1) = imax(1) - mod(imax(1)-lb(1),stride) endif imax(1) = getindx(imax(1), lb(1), stride) c do 10 i = imin(1), imax(1) c # Density if (nc.eq.1) qo(i) = aux(1,i) c # Velocity u fluctations if (nc.eq.2) qo(i) = q(2,i) c # Pressure fluctations if (nc.eq.5) qo(i) = q(1,i) c # Speed of sound if (nc.eq.7) qo(i) = aux(2,i) 10 continue return end