make_map.c

Go to the documentation of this file.

#include "prototypes.h"
 
void make_map(int num)
 
{
 
extern int pcp_in_use[200];
extern struct map map[1500];
extern struct pcp *pcp;
extern struct hrap_grid *hrap_grid;
 
char fbuf[100];
int x,y,ix,iy,l;
int minx,miny,totx,toty,uzpts,lzpts,mzpts,gzpts;
float lz,uz,mz,gz;
int ib;
 
if(pcp_in_use[num]==-1)
                           return;
 
strcpy(fbuf,"pcp");
 
read_file(fbuf,num,pcp);
 
minx=hrap_grid->hrap_minx;
miny=hrap_grid->hrap_miny;
totx=hrap_grid->maxi;
toty=hrap_grid->maxj;
 
for(ib=0;map[ib].hb5[0]!=0;ib++) {
 
        mz=0;uz=0;lz=0;gz=0;uzpts=0;lzpts=0;mzpts=0;gzpts=0;
 
        for(l=0;l<map[ib].hrap_points;l++) {
 
                x=map[ib].hrap_data[l].x;  
                y=map[ib].hrap_data[l].y;  
 
                ix=x-minx;
                iy=y-miny;
 
                if(ix < 0 || iy < 0 || ix > totx || iy > toty) {
 
                         continue;
 
               }
 
                if(map[ib].hrap_data[l].zone[3]==1) {
                            
                                  gz=gz+(float)pcp->value[ix][iy]/100;
                                  gzpts++;
 
                  }
 
                if(map[ib].hrap_data[l].zone[2]==1) {
                            
                                  uz=uz+(float)pcp->value[ix][iy]/100;
                                  uzpts++;
 
                  }
 
                if(map[ib].hrap_data[l].zone[1]==1) {
 
                                  mz=mz+(float)pcp->value[ix][iy]/100;
                                  mzpts++;
 
                       }
 
                if(map[ib].hrap_data[l].zone[0]==1){
 
                                  lz=lz+(float)pcp->value[ix][iy]/100;
                                  lzpts++;
 
                   }
                 
         
       
          }
 
         if(gzpts ==0)
                    gz=-1;
         else
                    gz=gz/(float)uzpts;
  
         if(uzpts ==0)
                    uz=-1;
         else
                    uz=uz/(float)uzpts;
  
         if(mzpts==0)
                    mz=-1;
         else
                    mz=mz/(float)mzpts;
         
         if(lzpts==0)
                   lz=-1;
 
         else
                   lz=lz/(float)lzpts;
         
              map[ib].uz[num]=uz;
              map[ib].lz[num]=lz;
              map[ib].mz[num]=mz;
              map[ib].gz[num]=gz;
 
              map[ib].maps_done[num]=1;
 
      }
 
return;
 
}