subroutine osav(j)
c
c=======================================================================
c                                                                    ===
c  This routine stores the current  J slab data into the volume      ===
c  data and put it in common block VOLDAT. These volume data is      ===
c  for Shapiro filtering  and other volume data operations that      ===
c  may involve horizontal or vertical operators.                     ===
c                                                                    ===
c=======================================================================
c
c-----------------------------------------------------------------------
c  Define global data.
c-----------------------------------------------------------------------
c
#include <cdefs.h>
#include <param.h>
#include <onedim.h>
#include <workspa.h>
#include <voldat.h>
#include <options.h>
c
c-----------------------------------------------------------------------
c  Define local data.
c-----------------------------------------------------------------------
c
      integer i,ip,j,k,m
c
c-----------------------------------------------------------------------
c  Begin executable code.
c-----------------------------------------------------------------------
c
      do 100 k=1,km
      do 100 i=1,imt
        ip=(j-1)*imt+i
        xu(ip,k)=ua(i,k)
        xv(ip,k)=va(i,k)
 100  continue
      do 110 i=1,imt
        ip=(j-1)*imt+i
        bkeep(ip,1)=bcon(i,1)
        bkeep(ip,2)=bcon(i,2)
 110  continue
      do 120 m=1,nt
      do 120 k=1,km
      do 120 i=1,imt
        ip=(j-1)*imt+i
        xt(ip,k,m)=ta(i,k,m)
 120  continue
      return
      end