get_rain.c

Go to the documentation of this file.

#include "prototypes.h"
 
float get_rain(struct stn_values *stn_values,char *pc,int inum,time_t etime,time_t btime,struct rainfall *r)
 
{
 
float val;
int ier;
struct tm *tm;
int isave;
int i;
int zdate;
 
if(pc[2]=='D') {
 
               ier=-1;
               tm=gmtime(&etime);
               
               zdate=(tm->tm_year+1900)*10000+(tm->tm_mon+1)*100+tm->tm_mday;
               
               for(i=0;i<inum;i++) {
              
               if(stn_values[i].zdate == zdate &&
                   stn_values[i].ztime/10000 > 11 && 
                    stn_values[i].ztime/10000 < 14 ) {
                    
                       ier=1;
                       r->rn=stn_values[i].value;
                       
                       /* printf("found daily rain %d %f\n",stn_values[i].zdate,r->rn); */
 
                      break; 
                       }    
                       
                       }                        
                  
               }
 
if(pc[2]=='H') {
 
              r->maxtdif=1L*3600L;
              r->maxgap=3.0;
              ier=hourly_rain(stn_values,inum,etime,btime,r);
 
               }
 
else if(strncmp(&pc[3],"RG",2)==0) {
 
               r->maxtdif=7L*3600L;
               r->firstdif=1L*3600L;
               r->maxgap=3.0;
               ier=domsat_rain(stn_values,inum,etime,btime,r);
 
 
                                 }
 
else if(strncmp(&pc[3],"RP",2)==0) {
 
                r->maxtdif=7L*3600L;
                r->firstdif=1L*3600L;
                r->maxgap=3.0;
                ier=domsat_rain(stn_values,inum,etime,btime,r);
 
                 }
 
else if(strncmp(&pc[3],"RR",2)==0) {
                   
                r->maxtdif=36L*3600L;
                r->maxgap=2.0;
                ier=alert_rain(stn_values,inum,etime,btime,r);
 
                                   }
 
else if(strncmp(&pc[3],"RM",2)==0) {
                  
                   
                r->maxtdif=6L*3600L;
                r->maxgap=5.0;
                r->firstdif=1L*3600L;
                
                if((etime-btime)/3600 == 24) {
                
               
                            etime=etime-4*3600;
                            btime=btime-4*3600;
                            r->maxtdif=6L*3600L;
                            r->maxgap=5.0;
                            r->firstdif=1L*3600L;
                
                
                                              }
                
                ier=snotel_rain(stn_values,inum,etime,btime,r); 
 
                
                if(r->rn > 0 && (etime-btime)/3600 == 24) {
                
                                isave=(r->rn*10);
                                r->rn=(float)isave/10;                
                 
 
                               }
 
     }
 
r->ier=ier;
 
if(ier != 1)
         val=-99;
 
else      
         val=r->rn;
 
 
if(val >= 0 && val < .005)
              val=0;
              
return(val);
 
}