get_station_list_build.c

Go to the documentation of this file.

#include "prototypes_new.h"
 
void get_station_list(char *fname,int smonth,int emonth)
 
{
 
 extern struct dval dval;
 extern char station_list_custom_file[100];
 extern char station_climo_file[100];
 int i,ier,sflag,field,j,xadd,yadd;
 extern struct station station[3000];
 extern int max_stations;
 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 f[15];
 int ifound;
 int oldk;
 float elev;
 
 max_stations=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==3000)
              break;
 
    p=strchr(kbuf,'\n');
 
    if(p!=NULL)
             *p=0;
 
    j=0;
    while(kbuf[j] != 0) {
    
             if(kbuf[j]=='\t')
                  kbuf[j]=' ';
                  
             j++;             
 
      }
 
    station[i].isoh=calloc(24,sizeof(float));
    station[i].hb5=calloc(10,sizeof(char));
    station[i].name=calloc(50,sizeof(char));
    station[i].parm=calloc(10,sizeof(char));
    station[i].cparm=calloc(10,sizeof(char));
    station[i].index=calloc(30,sizeof(short int));
    station[i].zindex=calloc(5,sizeof(short int));
    
    ier=sscanf(kbuf,"%s %s %f %f %f %d ",station[i].hb5,
               station[i].parm,&station[i].lat,
               &station[i].lon,&elev,&station[i].tip);
 
    station[i].elev=(int)elev;
    
    if(strncasecmp("PPD",station[i].parm,3)!=0 &&
       strncasecmp("PPM",station[i].parm,3)!=0 )
                      continue;
                      
    for(m=0;m<i;m++) {
    
            if(strcmp(station[i].hb5,station[m].hb5)==0 &&
               station[i].parm[4]==station[m].parm[4]){
               
                                    printf("duplicate station %s\n",station[i].hb5);
                                    break;
                                    
                                    }
                                    
                      }
                      
    
               
    if(i!=m)
           continue;
 
    if(ier != 6)
               continue;
 
    lat=station[i].lat;
    lon=station[i].lon;
 
    t=dval.a * cos(lat*conv)/(1+sin(lat*conv))
      * cos((lon-dval.lo-90)*conv) + dval.xo +.5;
 
    station[i].x=t;
 
    u=dval.a * cos(lat*conv)/(1+sin(lat*conv))
      * sin((lon-dval.lo-90)*conv) + dval.yo + .5;
       
    station[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(station[i].name,&kbuf[j]);
                          break;
 
            }
 
                }
 
  i++;
 
  }
 
max_stations=i;
 
fclose(fp);
 
for(i=0;i<max_stations;i++){
         
                 for(m=0;m<24;m++)
                             station[i].isoh[m]=-99; 
                             
 
}
 
fp=NULL;
if(station_climo_file[0]!=0)
        fp=fopen(station_climo_file,"r");
       
if(fp==NULL) {
 
          printf("no climo file\n");
 
}
 
else {
 
     for(;;) {
     
         p=fgets(kbuf,200,fp);
         if(p==NULL)
                 break;
                 
         ier=sscanf(kbuf,"%s %s %f %f %f %f %f %f %f %f %f %f %f %f",
             hb5,parm,&f[0],&f[1],&f[2],&f[3],&f[4],&f[5],&f[6],&f[7],&f[8],
             &f[9],&f[10],&f[11]);
         
         if(ier != 14)
                    continue;
             
         for(i=0;i<max_stations;i++) {
           
                 if(strcmp(station[i].hb5,hb5)==0 &&
            station[i].parm[4]==parm[4]) {
                 
                              for(m=0;m<12;m++)
                                  station[i].isoh[m]=f[m];
                                  
                              strcpy(station[i].cparm,parm); 
                  
                              }
                              
                          }
 
                            
               }
           
     rewind(fp);
     
     for(;;) {
    
         p=fgets(kbuf,200,fp);
         if(p==NULL)
                 break;
                 
         ier=sscanf(kbuf,"%s %s %f %f %f %f %f %f %f %f %f %f %f %f",
             hb5,parm,&f[0],&f[1],&f[2],&f[3],&f[4],&f[5],&f[6],&f[7],&f[8],
             &f[9],&f[10],&f[11]);
         
         if(ier != 14)
                    continue;
           
         if(strcmp(parm,"PPMPBCM") != 0 &&
        strcmp(parm,"PPMRZCM") != 0)
                                continue;          
        
         for(i=0;i<max_stations;i++) {
           
                 if(strcmp(station[i].hb5,hb5)==0 &&
            (station[i].isoh[0] < 0 ||
            strcmp(station[i].cparm,"PPMPBCM")==0)) {
            
             if(strcmp(station[i].cparm,"PPMPBCM")==0)
                                  printf("overwrite\n");
                 
                              for(m=0;m<12;m++)
                                  station[i].isoh[m]=f[m];
                                  
                              strcpy(station[i].cparm,parm);
                     
                             
                              }
                              
                          }
 
                            
               }    
           
fclose(fp);        
           
 
}
 
}