check_consistency.c

Go to the documentation of this file.

#include "prototypes_new.h" 
 
void check_consistency(int j)
 
{
 
int k,m;
float rtotal;
extern int max_stations;
extern struct pdata pdata[10];
extern struct station station[3000];
int i,ii;
 
for(k=0;k<max_stations;k++) {
 
        pdata[j].stn[k].tcons=1;
     
    if(pdata[j].stn[k].frain[4].data < 0 ||
       pdata[j].stn[k].frain[4].qual==1)
                                     continue;
    
        rtotal=0;
        for(m=0;m<4;m++) {
     
                if(pdata[j].stn[k].frain[m].data >= 0)
                            rtotal=rtotal+pdata[j].stn[k].frain[m].data;
        
             }
          
        if(fabs(rtotal-pdata[j].stn[k].frain[4].data) > .01) 
                    pdata[j].stn[k].tcons=-1;
            
            
  }
 
for(k=0;k<max_stations;k++) {
 
        station[k].lat;
        station[k].lon;
    
        for(m=0;m<5;m++) {
 
            if(pdata[j].stn[k].frain[m].qual==1 ||
           pdata[j].stn[k].frain[m].data < 0)
                               continue;
 
             pdata[j].stn[k].scons[m]=1;
         
         for(ii=0;ii<30;ii++){
   
                 i=station[k].index[ii];
 
                 if(station[i].lat == station[k].lat &&
                    station[i].lon == station[k].lon &&
            pdata[j].stn[i].frain[m].qual != 1 &&
            pdata[j].stn[i].frain[m].data >=0 &&
            pdata[j].stn[i].frain[m].data !=
            pdata[j].stn[k].frain[m].data) {
            
               pdata[j].stn[k].scons[m]=-1;
               
               }
                 
            }
            
            
        }    
            
  }
           
                
}