read_zlevel_b.c

Go to the documentation of this file.

#include "prototypes_new.h"
 
int read_zlevel_b(char *fname,time_t tget,int i)
   
{
 extern struct dval dval;
 extern struct zdata zdata[10];
 extern struct station zstation[500];
 extern int max_zstations;
 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;
 int uflag[4];
 
 fr=fopen(fname,"r");
 zdata[i].data_time=tget;
 
 gm=gmtime(&tget);
 
 zdata[i].ztime=(gm->tm_year+1900)*10000+(gm->tm_mon+1)*100+gm->tm_mday;
 
 for(j=0;j<4;j++) {
 
       number_found[j]=0;
       uflag[j]=0;
       
       }
 
 for(k=0;k<max_zstations;k++) {
 
                  for(m=0;m<5;m++) {
 
                               zdata[i].stn[k].zlevel2[m].data=-1;
                               zdata[i].stn[k].zlevel2[m].qual=-1;
 
                 }
 
        }
 
if(fr==NULL) {
 
               for(k=0;k<max_zstations;k++) {
 
                  for(m=0;m<4;m++) {
 
                               zdata[i].stn[k].zlevel2[m].data=
                               zdata[i].stn[k].zlevel1[m].data;
                               
                               zdata[i].stn[k].zlevel2[m].qual=
                               zdata[i].stn[k].zlevel1[m].qual;
                               
                               if(zdata[i].stn[k].zlevel2[m].data>=0)
                                         number_found[m]++;
                               
 
                 }
 
        }
        
     for(m=0;m<4;m++) {
 
              if(number_found[m] == 0)
                        zdata[i].used[m]=0;
            
              zdata[i].level[m]=1;
 
                   }
                   
    
     return(-1);
 
        }
 
 for(k=0;k<5;k++) {
 
       zdata[i].used[k]=1;
       zdata[i].level[k]=2;
       
       }
 
/* initialize structure */
 
 for(;;) {
 
         bad:
         p=fgets(kbuf,100,fr);
         if(p==NULL)
                   break;
 
         if(kbuf[0]==':')
                        continue;
 
         ier=sscanf(&kbuf[3],"%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_zstations;j++) {
 
                  if(strcmp(hb5,zstation[j].hb5)==0 && 
                     pc==zstation[j].parm[4]) {  
 
                             break;
 
               }
 
        }
 
         if(j==max_zstations) {
 
                      continue;
 
            }
 
         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++) {
 
                   zdata[i].stn[j].zlevel2[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) {
 
                              uflag[k]=1;
                      
                     
 
                              }
 
                   else  {
 
                        number_found[k]++;
 
                        zdata[i].stn[j].zlevel2[k].data=atof(buf);
                        
                        r=p+2;
                        qual=8;
                        while(*r != 0) {
 
                          if(*r != ' '){
 
                             if(*r=='S')
                                qual=8;
 
                             else if(*r=='E')
                                qual=5;
 
                             else if(*r=='W')
                                qual=2;
 
                             else if(*r=='F')
                                qual=1;
 
                             zdata[i].stn[j].zlevel2[k].qual=qual;
                             
                        if((qual == 8) &&
                              zdata[i].stn[j].zlevel2[k].data !=
                              zdata[i].stn[j].zlevel1[k].data) {
                  
                                  uflag[k]=1;
                  
                       
                                                   
                              }                     
 
                             break;
 
              }
 
                          else 
                                        r++;
 
                  }
 
 
         }
 
          p=strchr(q,0);
          q=p+1;      
         
         }
 
       }
 
 
fclose(fr);
 
for(m=0;m<4;m++) {
 
              if(uflag[m]==1) {
 
                   number_found[m]=0;
           
                   for(k=0;k<max_zstations;k++) {
 
                               zdata[i].stn[k].zlevel2[m].data=
                               zdata[i].stn[k].zlevel1[m].data;
                               
                               zdata[i].stn[k].zlevel2[m].qual=
                               zdata[i].stn[k].zlevel1[m].qual;
                               
                               if(zdata[i].stn[k].zlevel2[m].data>=0)
                                         number_found[m]++;
                               
 
                 }
                    
                 zdata[i].level[m]=1;
        
 
 }
              
 
}
 
for(j=0;j<4;j++) {
 
     if(number_found[j] == 0) {
     
              zdata[i].used[j]=0;
              zdata[i].level[j]=1;
 
              }
    
   }
 
return(1);
 
}