new_read_qpf_grids.c

Go to the documentation of this file.

 
#include "prototypes.h"
 
int read_qpf_grids(int k,char *fname)
 
{
extern int pcp_in_use[200];
FILE *fp;
extern struct hrap_grid *hrap_grid;
int i,j,h;
char kbuf[10000];
extern struct pcp *pcp;
int minhrapi,minhrapj,maxhrapi,maxhrapj;
int ghrapi,ghrapj;
int ii,jj;
int ier;
char *p;
int gmini,gminj,gmaxi,gmaxj;
 
fp=fopen(fname,"r");
 
 
if(fp==NULL) {
 
             printf("could not open %s\n",fname);
             
             return(-1);            
 
       }
 
for(i=0;i<hrap_grid->maxi;i++) {
 
for(j=0;j<hrap_grid->maxj;j++) {
 
          pcp->value[i][j]=0;
 
    }
 
}
 
minhrapi=hrap_grid->hrap_minx;
minhrapj=hrap_grid->hrap_miny;
maxhrapi=hrap_grid->hrap_minx + hrap_grid->maxi;
maxhrapj=hrap_grid->hrap_miny + hrap_grid->maxj;
 
p=fgets(kbuf,100,fp);
 
ier=sscanf(kbuf,"%d %d %d %d %d\n",&gmini,&gminj,&gmaxi,&gmaxj,&iflag);
 
if(ier==4)
        iflag=0;
 
for(i=0;i<gmaxi;i++) {
 
     if(iflag==0)
     fread(kbuf,sizeof(char),gmaxj*4,fp);
 
     else
     fread(kbuf,sizeof(char),gmaxj*6,fp);
 
     /* get hrap coord of gridded data */
 
     ghrapi=gmini+i;
 
     ii=ghrapi-minhrapi;
 
     if(ghrapi >= minhrapi && ghrapi < maxhrapi) {  
 
     h=0;
 
     for(j=0;j<gmaxj;j++) {
 
            ghrapj=gminj+j;
 
            jj=ghrapj-minhrapj;
 
 
 
            if(ghrapj >= minhrapj && ghrapj < maxhrapj)  {
 
                        pcp->value[ii][jj]=atoi(&kbuf[h]);;
                        
 
                }
 
            h=h+4;
 
 
      }
 
   }
 
   }
 
fclose(fp);
 
/*copy to internal file */
 
write_file("pcp",k,pcp);
 
pcp_in_use[k]=1;
 
return(1);
 
}