make_mat.c

Go to the documentation of this file.

#include "prototypes_new.h"
 
void make_mat(int num,int mnum)
 
{
 
extern int render_all;
extern wfo_orig;
extern int wfo_in_use[20];
extern int pcp_in_use[500];
extern struct map map[1500];
extern struct pcp *pcp;
extern struct hrap_grid *hrap_grid;
int h;
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++) {
 
        if(render_all == 0) {
        
        /* skip if you are not the owner and there is no value */  
        /* should not affect daily_qc */
        /* may affect verify */
        /* should fix specify */
        /* auto_specify ??? */
       
        if(map[ib].owner!=wfo_orig &&                            
           map[ib].tuz[num] < 0 && 
           map[ib].tlz[num] < 0 &&
           map[ib].tmz[num] < 0 &&
           map[ib].tgz[num] < 0)
                             continue;       
 
        }
 
        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 ||
                   pcp->value[ix][iy] < -500) {
 
                         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].tuz[mnum]=uz;
         map[ib].tlz[mnum]=lz;
         map[ib].tmz[mnum]=mz;
         map[ib].tgz[mnum]=gz;
 
         map[ib].tmaps_done[mnum]=1;
     
      }
 
return;
 
}