stuff.c

Go to the documentation of this file.

 
int read_cqpf(time_t tim,int st,int fi)
 
{
 
struct tm *gmtim;
char points[100],kbuf[100],*p;
long r,s;
float conv=.0174;
float lat,lon,pcp;
int k,ier,i,numpoints;
struct abogus *cbogus;
FILE *fp;
int m;
 
wfobogus[0].owner=-1;
 
gmtim=gmtime(&tim);
 
for(m=0;m<20;m++){
 
   for(k=st;k<fi;k++) {
 
   for(i=0;i<100;i++) {
        wfobogus[m].bogus[k].inuse[i]=-1;
      }
}
 
  }
 
 
for(m=0;m<20;m++){
 
for(k=st;k<fi;k++) 
       wfobogus[m].bogus_flag[k]=-1;
 
 }
 
/* memory leak need to free bbogus area */
 
sprintf(points,"%s/fpointa.%s.%02d-%02d-%02d",scratch_file,arfc,
        gmtim->tm_mon+1,gmtim->tm_mday,gmtim->tm_year);
 
fp=fopen(points,"r");
 
if(fp==NULL)
          return(-1);
 
m=0;
for(;;) {
 
   p=fgets(kbuf,100,fp);
   
   if(p==NULL)
           break;
           
   p=strchr(kbuf,'\n');
   if(p!=NULL)
            *p=0;
 
   if(strcmp(kbuf,".END")==0)
                      continue;
 
   /* get owner here */
   
   wfobogus[m].owner=-1;
  
   for(k=0;k<20;k++) {  
 
              if(valid_list[k].wfo[0]==0)
                                     break;
 
              if(strcmp(valid_list[k].wfo,kbuf)==0) {
 
                            if(strcmp(wfo,kbuf)==0) 
                                                   continue;
 
                            wfobogus[m].owner=valid_list[k].alias;
                           
                            break;
 
            }
 
   }
 
   if(wfobogus[m].owner==-1)
                     continue;
   
  
 
   for(k=st;k<fi;k++) {
 
         wfobogus[m].bbogus[k]=NULL;
 
         numpoints=0;
 
         for(;;) {
 
                 p=fgets(kbuf,100,fp);
                 if(p==NULL)
                               return(1);
 
                 if(strncmp("END",kbuf,3)==0)
                                 break;
 
                 if(wfobogus[m].bbogus[k]==NULL){
 
                       wfobogus[m].bogus_flag[k]=0;
                       wfobogus[m].bbogus[k]=calloc(1,sizeof(struct abogus));
                       cbogus=wfobogus[m].bbogus[k];
 
             }
 
                 else  {
 
                 cbogus->next=calloc(1,sizeof(struct abogus));
                 cbogus=cbogus->next;
      
           }
 
                 numpoints++;
                 ier=sscanf(kbuf,"%f %f",&lat,&lon);
 
                 r=dval.a * cos(lat*conv)/(1+sin(lat*conv))
                       * cos((lon-dval.lo-90)*conv) + dval.xo +.5;
 
                 s=dval.a * cos(lat*conv)/(1+sin(lat*conv))
                       * sin((lon-dval.lo-90)*conv) + dval.yo + .5;
       
                 cbogus->x=r;
                 cbogus->y=s;
                 cbogus->next=NULL;
 
           }
 
               for(i=0;i<20;i++) {
 
                     p=fgets(kbuf,100,fp);
                     if(p==NULL)
                               return(1);
 
                     if(strncmp("END",kbuf,3)==0)
                                 break;
 
                     ier=sscanf(kbuf,"%f %f %f",&lat,&lon,&pcp);
 
                     r=dval.a * cos(lat*conv)/(1+sin(lat*conv))
                       * cos((lon-dval.lo-90)*conv) + dval.xo +.5;
 
                     s=dval.a * cos(lat*conv)/(1+sin(lat*conv))
                       * sin((lon-dval.lo-90)*conv) + dval.yo + .5;
       
                     wfobogus[m].bogus_flag[k]=1;
                     wfobogus[m].bogus[k].inuse[i]=1;
                     wfobogus[m].bogus[k].x[i]=r;
                     wfobogus[m].bogus[k].y[i]=s;
                     wfobogus[m].bogus[k].value[i]=pcp;
                     
           }
 
       if(numpoints==0)
               wfobogus[m].bogus_flag[k]=-1;
 
       }
 
m++;
 
}
 
wfobogus[m].owner=-1;
 
 
fclose(fp);
 
}