read_precip_a.c

Go to the documentation of this file.

#include "prototypes_new.h"
 
int read_precip_a(char *fname,time_t tget,int i)
   
{
 extern struct dval dval;
 extern struct pdata pdata[10];
 extern struct station station[3000];
 extern int max_stations;
 FILE *fr;
 int j,k,ier,m,qual,len,len2;
 char kbuf[100],*p,*q,*r,buf[100],hb5[10];
 long t,u;
 float conv=.0174;
 double lat,lon;
 int number_found[5];
 char parm[10];
 int type;
 char pc,datbuf[50],parmbuf[50];
 int maxk,startk;
 struct tm *gm;
 
 fr=fopen(fname,"r");
 pdata[i].data_time=tget;
 
 gm=gmtime(&tget);
 
 pdata[i].ztime=(gm->tm_year+1900)*10000+(gm->tm_mon+1)*100+gm->tm_mday;
 
 printf("ztime %d\n", pdata[i].ztime);
 
 for(j=0;j<5;j++)
       number_found[j]=0;
 
 for(k=0;k<max_stations;k++) {
 
                  for(m=0;m<5;m++) {
 
                               pdata[i].stn[k].rain[m].data=-1;
                               pdata[i].stn[k].rain[m].qual=-1;
 
                 }
 
        }
 
if(fr==NULL) {
 
              for(k=0;k<5;k++)
                        pdata[i].used[k]=0;
 
              pdata[i].level=0;
            
              return(-1);
 
        }
 
 for(k=0;k<5;k++)
       pdata[i].used[k]=1;
 
 pdata[i].level=1;
 
/* initialize structure */
 
 for(;;) {
 
         bad:
         p=fgets(kbuf,100,fr);
         if(p==NULL)
                   break;
 
         if(kbuf[0]==':')
                        continue;
 
         if(kbuf[0]=='.' && kbuf[1]=='E')
                                  type=1;
 
         if(kbuf[0]=='.' && kbuf[1]=='A') 
                                  type=0;
 
         ier=sscanf(&kbuf[2],"%s %s %s",hb5,datbuf,parmbuf);
         if(ier==0)
                         continue;
 
         p=strchr(parmbuf,'/');
         if(p==NULL)
                 continue;
 
         pc=*(p+5);
 
         for(j=0;j<max_stations;j++) {
 
                  if(strcmp(hb5,station[j].hb5)==0 && 
                     pc==station[j].parm[4]) {  
 
                             break;
 
               }
 
        }
 
         if(j==max_stations) {
 
                      continue;
 
            }
 
         if(type==0) {
 
                     p=strchr(kbuf,'/');
                     if(p==NULL)
                             continue;
 
                     q=strchr(p,' ');
                     if(p==NULL)
                             continue;
                    
                   maxk=5;
                   startk=4;
 
           }
 
         else if(type==1) {
 
         q=kbuf;
         for(k=0;k<3;k++) {
 
                   p=strchr(q,'/');
                   if(p==NULL)
                             break;
 
                   q=p+1;
                  
           }
 
                   if(k!=3)
                          continue;  /* missing */
 
                   maxk=4;
                   startk=0;
 
            }
 
         for(k=startk;k<maxk;k++) {
 
                   pdata[i].stn[j].rain[k].qual=0;
 
                   if((p=strchr(q,'/'))==NULL &&
                      (p=strchr(q,'\n'))==NULL)
                                          goto bad;
 
                   *p=0;
 
                   strcpy(buf,q);
 
                   if((p=strchr(buf,'.'))==NULL) {
 
                              if((p=strchr(buf,'m'))!=NULL ||
                                 (p=strchr(buf,'M'))!=NULL) {
                                 
                                      pdata[i].stn[j].rain[k].data=-1;
                                      pdata[i].stn[j].rain[k].qual=-1;
 
                         }
 
                              else             {
 
                                             
                                               pdata[i].stn[j].rain[k].data=-2;
 
                         }
 
                }
 
                   else  {
 
                        number_found[k]++;
 
                        pdata[i].stn[j].rain[k].data=atof(buf);
 
                        r=p+3;
                        qual=8;
                        while(*r != 0) {
 
                          if(*r != ' '){
 
                             if(*r=='S')
                                qual=8;
 
                             else if(*r=='F')
                                qual=1;
 
                             else if(*r=='L')
                                qual=2;
 
                             else if(*r=='Q')
                                qual=3;
 
                             else if(*r=='F')
                                qual=4;
 
                             else if(*r=='V')
                                qual=8;
 
                             else if(*r=='E')
                                qual=5;
 
                             else if(*r=='T')
                                qual=6;
                                    
                             pdata[i].stn[j].rain[k].qual=qual;
 
                             break;
 
                   }
 
                                 else 
                                        r++;
 
              }
 
 
         }
 
          p=strchr(q,0);
          q=p+1;      
         
         }
 
       }
 
 
fclose(fr);
 
for(j=0;j<5;j++) {
 
     if(number_found[j] == 0)
              pdata[i].used[j]=0;
 
   }
 
 
return(1);
 
}