make_fmap.c

Go to the documentation of this file.

#include "prototypes.h"
 
void make_fmap(int k)
 
{
 
extern struct map map[500];
extern struct pcp *qpf;
extern struct hrap_grid *hrap_grid;
int m,ib,l;
int x,y,ix,iy;
int minx,miny,totx,toty,uzpts,lzpts,mzpts,gzpts;
float lz,uz,mz,gz;
time_t tim;
struct tm *gmtim; 
 
tim=time(NULL);
 
gmtim=gmtime(&tim);
     
minx=hrap_grid->hrap_minx;
miny=hrap_grid->hrap_miny;
totx=hrap_grid->maxi;
toty=hrap_grid->maxj;
 
read_file("qpf",k,qpf); 
      
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)qpf->value[ix][iy]/100;
                                  gzpts++;
 
                  }
 
                if(map[ib].hrap_data[l].zone[2]==1) {
                            
                                  uz=uz+(float)qpf->value[ix][iy]/100;
                                  uzpts++;
 
                  }
 
                if(map[ib].hrap_data[l].zone[1]==1) {
 
                                  mz=mz+(float)qpf->value[ix][iy]/100;
                                  mzpts++;
 
                       }
 
                 if(map[ib].hrap_data[l].zone[0]==1) {
 
                                  lz=lz+(float)qpf->value[ix][iy]/100;
                                  lzpts++;
 
                  
                 
          }
       
          }
 
         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;
 
         if(gzpts==0)
                    gz=-1;
         else
                    gz=gz/(float)gzpts;
 
 
         map[ib].gz[k]=gz;
         map[ib].uz[k]=uz;
         map[ib].mz[k]=mz;
         map[ib].lz[k]=lz;
         map[ib].maps_done[k]=1;
 
       }
 
}