get_zstation_list.c

Go to the documentation of this file.

#include "prototypes_new.h"
 
void get_zstation_list(char *fname)
 
{
 
 extern char zstation_list_custom_file[100];
 extern struct dval dval;
 int i,ier,sflag,field,j,xadd,yadd;
 extern struct station zstation[750];
 extern int max_zstations;
 FILE *fp;
 char *p,kbuf[200],hb5[10],parm[10];
 double dist1,dist2,dist,sorted[30];
 int m,l,h,k;
 int kk;
 float lat,lon;
 float conv=.0174;
 long t,u;
 float elev;
 
 max_zstations=0;
 
 fp=fopen(fname,"r");
 
 if(fp==NULL) {
 
              printf("could not open %s\n",fname);
              exit(1);
 
              }
 
 i=0;
 for(;;) {
 
    p=fgets(kbuf,120,fp);
    if(p==NULL)
              break;
 
    if(i==750)
              break;
          
    p=strchr(kbuf,'\n');
 
    if(p!=NULL)
             *p=0;
 
    j=0;
    while(kbuf[j] != 0) {
    
             if(kbuf[j]=='\t')
                  kbuf[j]=' ';
                  
             j++;             
 
      }
 
    zstation[i].hb5=calloc(10,sizeof(char));
    zstation[i].name=calloc(50,sizeof(char));
    zstation[i].parm=calloc(10,sizeof(char));
 
    ier=sscanf(kbuf,"%s %s %f %f %f %d ",zstation[i].hb5,
               zstation[i].parm,&zstation[i].lat,
               &zstation[i].lon,&elev,&zstation[i].tip);
           
    if(strncasecmp("HZI",zstation[i].parm,3)!=0)     
                      continue;
        
     
    for(m=0;m<i;m++) {
    
            if(strcmp(zstation[i].hb5,zstation[m].hb5)==0 &&
               zstation[i].parm[4]==zstation[m].parm[4]){
               
                                    printf("duplicate zstation %s\n",zstation[i].hb5);
                                    break;
                                    
                                    }
                                    
                      }          
          
 
    if(ier != 6)
               continue;
 
    lat=zstation[i].lat;
    lon=zstation[i].lon;
 
    t=dval.a * cos(lat*conv)/(1+sin(lat*conv))
      * cos((lon-dval.lo-90)*conv) + dval.xo +.5;
 
    zstation[i].x=t;
 
    u=dval.a * cos(lat*conv)/(1+sin(lat*conv))
      * sin((lon-dval.lo-90)*conv) + dval.yo + .5;
       
    zstation[i].y=u;
 
    sflag=1;
    field=0;
    for(j=0;j<strlen(kbuf);j++) {
 
                 if(kbuf[j] !=' ' && sflag==1) {
                                    sflag=0; 
                                    field++;
 
                  }
      
                 else if(kbuf[j]==' ')
                          sflag=1;
 
                 if(field==7) {
 
                          kbuf[j+49]=0;
 
                          strcpy(zstation[i].name,&kbuf[j]);
                          break;
 
            }
 
                }
 
  i++;
 
  }
 
max_zstations=i;
 
fclose(fp);
 
strcpy(kbuf,fname);
strcat(kbuf,".zcustom");
strcpy(zstation_list_custom_file,kbuf);
 
fp=fopen(kbuf,"r");
 
if(fp!=NULL) {
 
     for(;;) {
 
           p=fgets(kbuf,80,fp);
           if(p==NULL)
                   break;
 
           p=strchr(kbuf,'\n');
           *p=0;
 
           ier=sscanf(kbuf,"%s %s %d %d",hb5,parm,&xadd,&yadd);
 
           for(j=0;j<max_zstations;j++) {
 
                   if(strcmp(hb5,zstation[j].hb5)==0) {
                     
                      if(strncmp(parm,zstation[j].parm,3)==0 &&
                              parm[4]==zstation[j].parm[4]) {
 
                              zstation[j].xadd=xadd;
                              zstation[j].yadd=yadd;
                              break;
 
                }
 
         }
 
     }
     }
     
   fclose(fp); 
 
   }        
 
return;
 
}