get_hrap_mask.c

Go to the documentation of this file.

 
#include "prototypes_new.h"
 
int get_hrap_mask(char *fname)
 
{
 
extern struct tag tag[20];
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;
int iflag;
 
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++) {
 
          hrap_grid->owner[i][j]=-1;
          
 
    }
 
}
 
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);
 
for(i=0;i<20;i++)
     tag[i].wfo[0]=0;
     
ier=sscanf(kbuf,"%d %d %d %d %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s\n",
           &gmini,&gminj,&gmaxi,&gmaxj,
           tag[0].wfo,tag[1].wfo,tag[2].wfo,tag[3].wfo,
           tag[4].wfo,tag[5].wfo,tag[6].wfo,tag[7].wfo,
           tag[8].wfo,tag[9].wfo,tag[10].wfo,tag[11].wfo,
           tag[12].wfo,tag[13].wfo,tag[14].wfo,tag[15].wfo,
           tag[16].wfo,tag[17].wfo,tag[18].wfo,tag[19].wfo);
      
for(i=0;i<gmaxi;i++) {
 
     fread(kbuf,sizeof(char),gmaxj*6,fp);
     
     kbuf[gmaxj*6]=0; 
     /* 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)  {
 
                           hrap_grid->owner[ii][jj]=atoi(&kbuf[h]);
                        
 
                }
            
            h=h+6;
 
      }
 
   }
 
   }
 
fclose(fp);
 
return(1);
 
}