hourly_rain.c

Go to the documentation of this file.

/*  possible error returns
 
    -1 no report after begin time
    -2 no report after end time
    -3 time gap between begin time and first report after begin time too large
    -4 data not reliable
 
    r->rn    amount of rain
    r->tdif    time gap of end time and last report
    
                                                                 */
 
#include "prototypes.h"
 
hourly_rain(struct stn_values *stn_values,int inum,
           time_t tstart,time_t tstop,struct rainfall *r)
 
{
 
int j,eflag;
int xpctd;
 
r->rn=0.0;
r->rel=0;
r->obs=0;
r->tskips=0;
r->tdif=0L;
 
eflag=-1;
    
xpctd=(tstart-tstop)/3600;
 
if(xpctd==0)
        xpctd=1;
          
for(j=0;j<inum;j++) {
 
         /* observation time before start time */
 
                if(stn_values[j].clock > tstart) 
                                      continue;
 
                if(stn_values[j].clock <= tstop) {
 
                                      /*
                                      printf("pph %d %d %f\n",xpctd,r->obs,                                      
                                      (float)r->obs/(float)xpctd);                
                                      */
                 
                                      if((float)r->obs/(float)xpctd > .75)
                                      return(1);
                                      
                                      else
                                      return(-1);
 
                                                   }
 
                if((r->obs==0 || j==0) && (tstart - stn_values[j].clock) > r->maxtdif){         
 
                                                  r->tdif=tstart-stn_values[j].clock;
                                                  return(-3);
 
                                                                            }
 
                r->obs++;
                eflag=-2;
 
                if(stn_values[j].value < 0 || stn_values[j].value >= r->maxgap) {
 
                                         r->rel++;
                                         if(r->rel >= 2)
                                                    eflag=-4;
 
                                                             }
 
                else
                                         r->rn=r->rn+stn_values[j].value;
 
                                 }
 
return(eflag);
 
}