#include <stdio.h>

extern FILE *inf, *opf, *stf;

stout(infile,sensor,outfile,stfile,minlon,maxlon,minlat,maxlat,nout,amin,amax,
      gmin,gmax,bmin,bmax,scale,offset,nlat,nlon,total,count,intcnt,intdat,
      nscan)

/*********************************/ 
/* declare the input parameters! */
/*********************************/ 
char *infile,*sensor,*outfile,*stfile;
int nlon,nlat,nscan,nout,intdat;
int bmin[],bmax[],total[],count[],intcnt[];
float minlon,maxlon,minlat,maxlat,scale[],offset[];
float amin[],amax[],gmin[],gmax[];
{
/*********************************/ 
/* declare the local variables!  */
/*********************************/ 
   int i,istat;
   float xs,ys;

/**************************************/ 
/* begin printing output status file! */
/**************************************/ 
   fprintf(stf,"Input file:  %s",infile);
   fprintf(stf,"\nSensor:      %s",sensor);
   fprintf(stf,"\nOutput file: %s",outfile);
   fprintf(stf,"\nStatus file: %s",stfile);
   fprintf(stf,"\n\nCharacteristics of the input data: "); 
   fprintf(stf,"\n\n           Data             ");
   fprintf(stf,"\n   Band    Minimum   Maximum");
   fprintf(stf,"\n   ----    -------   -------");
   for (i=0; i<nout; i++) {
     if (amin[i] >  1.0e10) amin[i] = -1.0;
     if (amax[i] < -1.0e10) amax[i] = -1.0;
     if (gmin[i] >  1.0e10) gmin[i] = -1.0;
     if (gmax[i] < -1.0e10) gmax[i] = -1.0;
     if (bmin[i] ==  9000000) bmin[i] = -1.0;
     if (bmax[i] == 0) bmax[i] = -1.0;
     fprintf(stf,"\n    %2d)   %8.2f  %8.2f",i+1,
             amin[i],amax[i]);
    
   }
   fprintf(stf,"\n\nCharacteristics of the output grid: ");
   if (minlon < maxlon)
     xs = ( maxlon - minlon ) / (float) nlon;
   else
     xs = ( 360.0 + maxlon - minlon ) / (float) nlon;
   ys = ( maxlat - minlat ) / (float) nlat;
   fprintf(stf,"\n\n   Longitude: %4d bins from %7.2f to %7.2f, width  = %7.4f",
     nlon,minlon,maxlon,xs);
   fprintf(stf,"\n   Latitude:  %4d bins from %7.2f to %7.2f, height = %7.4f",
     nlat,minlat,maxlat,ys);
   if (intdat == 1) 
      fprintf(stf,"\n   Interpolation -> ON");
   else
      fprintf(stf,"\n   Interpolation -> OFF");
   fprintf(stf,"\n");
   fprintf(stf,"\n         Number  Number   Number   Output grid       "); 
   fprintf(stf,"\n         of pts  of bins  of bins  data              "); 
   fprintf(stf,"\n   Band  grided  w/data   inter.   minimum  maximum  Scale  Offset");
   fprintf(stf,"\n   ----  ------  -------  -------  -------  -------  -----  ------");
   for (i=0; i<nout; i++) {
     fprintf(stf,"\n    %2d)%8d%9d%9d  %7.2f  %7.2f %7.2f %7.2f",
        i+1,total[i],count[i],intcnt[i],gmin[i],gmax[i],scale[i],offset[i]);
   }
   fprintf(stf,"");
   fprintf(stf,"\n\nTotal of %d scans processed!\n",nscan);

}