#include <stdio.h>

write_grid(opf,opf2,grid,pop,dscale,doffset,nlat,nlon,nout,gmin,gmax,bmin,bmax,totpts,pop_flag)

FILE *opf, *opf2;
float **grid[];
short int **pop[];
float dscale[],doffset[],gmin[],gmax[];
int nlat,nlon,nout,bmin[],bmax[],*totpts,pop_flag;
{

   int i,j,k,ival4,istat;
   short int ival;

   for (k=0; k<nout; k++) {  
     for (i=nlat-1; i>=0; i--) {
       for (j=0; j<nlon; j++) {
         if (grid[k][i][j] > gmax[k]) gmax[k] = grid[k][i][j];
         if ((grid[k][i][j] < gmin[k]) && (grid[k][i][j] >= 0.0))
              gmin[k] = grid[k][i][j];
         if (grid[k][i][j] == -1.0)
           ival4 = -1;
         else if (grid[k][i][j] == -3.0)
           ival4 = -3;
         else if (grid[k][i][j] == -4.0)
           ival4 = -4;
         else if (grid[k][i][j] == -5.0)
           ival4 = -5;
         else if (grid[k][i][j] >= 0.0) {
           if (grid[k][i][j] > 655.35) grid[k][i][j] = 655.35;
           ival4 = (int)((grid[k][i][j] + doffset[k])*dscale[k] + 0.5);
           if (ival4 < bmin[k]) bmin[k] = ival4;
           if (ival4 > bmax[k]) bmax[k] = ival4;
         }
         else
           ival4 = -1;
         ival = ival4;
         istat = fwrite(&ival,2,1,opf);
         if (pop_flag == 1) istat = fwrite(&pop[k][i][j],2,1,opf2);
         *totpts = *totpts + 1;
       }
     }
   }
   return;
}