read_t_a.c

Go to the documentation of this file.

#include "prototypes_new.h"
 
int read_t_a(char *fname,time_t tget,int i)
   
{
 extern struct dval dval;
 extern struct tdata tdata[10];
 extern struct station tstation[1000];
 extern int max_tstations;
 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[6];
 char parm[10];
 int type;
 char pc,datbuf[50],parmbuf[50];
 int maxk,startk;
 struct tm *gm;
 float fval;
 int kread;
 int offset;
 
 fr=fopen(fname,"r");
 
 tdata[i].data_time=tget;
 
 gm=gmtime(&tget);
 
 tdata[i].ztime=(gm->tm_year+1900)*10000+(gm->tm_mon+1)*100+gm->tm_mday;
 
 for(j=0;j<6;j++)
       number_found[j]=0;
 
 for(k=0;k<max_tstations;k++) {
 
                  for(m=0;m<6;m++) {
 
                               tdata[i].stn[k].tlevel1[m].data=-99;
                               tdata[i].stn[k].tlevel1[m].qual=-99;
 
                 }
 
        }
 
if(fr==NULL) {
 
              for(k=0;k<6;k++) {
          
                        tdata[i].used[k]=0;
                        tdata[i].level[k]=0;
            
            }
            
              return(-1);
 
        }
 
 for(k=0;k<6;k++) {
 
       tdata[i].used[k]=1;
       tdata[i].level[k]=1;
 
        }
 
/* initialize structure */
 
 for(;;) {
 
         bad:
         p=fgets(kbuf,100,fr);
         if(p==NULL)
                   break;
 
         if(kbuf[0]==':')
                        continue;
 
         ier=sscanf(&kbuf[2],"%s %s %s",hb5,datbuf,parmbuf);
         if(ier==0)
                         continue;         
 
         p=strchr(parmbuf,'/');
         if(p==NULL)
                 continue;
 
         if(*(p+6)=='X') {
             
           offset=7;    
                   kread=1;
                   startk=4;
           maxk=5;
           
           }
           
     else if(*(p+6)=='N') {
     
               offset=7;
               kread=1;
               startk=5;
           maxk=6;
           
           }
           
         else      {
 
                   offset=0;
                   kread=3;  
               startk=0;
           maxk=4;
     
     
     }         
                 
         pc=*(p+5);
 
         for(j=0;j<max_tstations;j++) {
 
                  if(strcmp(hb5,tstation[j].hb5)==0 && 
                     pc==tstation[j].parm[4]) {  
 
                             break;
 
               }
 
        }
 
         if(j==max_tstations) {
 
                      continue;
 
            }
        
         q=kbuf;
         for(k=0;k<kread;k++) {
 
                   p=strchr(q,'/');
                   if(p==NULL)
                             break;
 
                   q=p+1;
                  
           }
 
         if(k!=kread)
                          continue;  /* missing */
 
         for(k=startk;k<maxk;k++) {
 
                   tdata[i].stn[j].tlevel1[k].qual=0;
 
                   if((p=strchr(q,'/'))==NULL &&
                      (p=strchr(q,'\n'))==NULL)
                                          goto bad;
 
                   *p=0;
 
                   strcpy(buf,q+offset);
           
                   if((p=strchr(buf,'m'))!=NULL ||
                      (p=strchr(buf,'M'))!=NULL) {
                                 
                                      tdata[i].stn[j].tlevel1[k].data=-99;
                                      tdata[i].stn[j].tlevel1[k].qual=-99;
 
                         }
                         
               else {
                             
           p=buf;
 
                   number_found[k]++;
 
                 
                   tdata[i].stn[j].tlevel1[k].data=atoi(buf);           
            
                   r=p;
                   qual=8;
           
                   while(*r != 0) {
 
                if(isalpha(*r) != 0) {
 
                             if(*r=='S' || *r=='A')
                                qual=8;
 
                             else if(*r=='F')
                                qual=1;
 
                             else if(*r=='V')
                                qual=8;
 
                             else if(*r=='E')
                                qual=5;
                
                 tdata[i].stn[j].tlevel1[k].qual=qual;  
 
 
 
                             break;
 
              }
 
                     else 
                                        r++;
 
            }
 
               }
        
          p=strchr(q,0);
          q=p+1;      
         
        }
 
       }
 
 
fclose(fr);
 
for(j=0;j<6;j++) {
 
     if(number_found[j] == 0)
              tdata[i].used[j]=0;
 
   }
 
 
return(1);
 
}