#include "IWInclude.h"

main()
{
   int nxyz[3],mode,count;
   char label[80];
   short lilarray[128*128];
   int i,j,k;
   float dmin,dmax,dmean,secmean;

   nxyz[0] = 128;
   nxyz[1] = 128;
   nxyz[2] = 20;

   mode = 1;
   sprintf(label, "demo");
   IMOpen(1, "/pop/diggs/junk.dat", "new");  /* open a new file */
   IMCrHdr(1, nxyz, nxyz, mode, label, 1);   /* create header for file */
   IMAlCon(1,0);
   dmin = 10000;
   dmax = -10000;
   dmean = 0;
   IMWrHdr(1, label, 0, dmin, dmax, dmean);
   for(i=0;i<nxyz[2];i++)
   {
      count = 0;
      secmean = 0;
      for(j=0;j<nxyz[1];j++)
      {
         for(k=0;k<nxyz[0];k++)
         {
           lilarray[count] = (i * 10) + k;
           if(lilarray[count] < dmin)
              dmin = lilarray[count];
           if(lilarray[count] > dmax)
              dmax = lilarray[count];
           secmean += lilarray[count]; 
           count++;
         }
      }
      secmean = secmean / (nxyz[0] * nxyz[1]);
      dmean += secmean;
      IMWrSec(1, lilarray);                      /* write the sections one by one */
   }
   dmean = dmean / (nxyz[2]);
   IMWrHdr(1, label, 1, dmin, dmax, dmean);
   IMClose(1);
}