get_station_list_new.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) {
 
          for(i=0;i<max_stations;i++) {
         
          
          for(k=0;k<12;k++) {
 
            ier=is_good(k,smonth,emonth);
 
            if(ier == -1)
                   continue;
 
             station[i].isoh[k]=(get_lisohyet(station[i].lat,station[i].lon,k)/25.4);
 
             if(station[i].isoh[k] < 0)
                      printf("No PRISM climatology for %s\n",station[i].hb5);
                      
             strcpy(station[i].cparm,"PPMPBZZ");
 
       }
 
       }
 
}
 
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;
       
         if(parm[0]!='P')
                  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);        
           
 
}
 
/* finally sum up for seasonal totals */
 
for(i=0;i<max_stations;i++){
 
oldk=-1;
 
for(kk=0;kk<12;kk++) {
 
 if(kk < 3)
          k=9+kk;
          
 else
          k=kk-3;
 
 if(oldk==-1) {
 
      station[i].isoh[k+12]=station[i].isoh[k];
      
    
          }
    
     else {
    
       station[i].isoh[k+12]=station[i].isoh[k] + station[i].isoh[oldk];
 
      
    }
    
    oldk=k+12;
    
    }
 
 
}
 
strcpy(kbuf,fname);
strcat(kbuf,".custom");
strcpy(station_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_stations;j++) {
 
                   if(strcmp(hb5,station[j].hb5)==0) {
                     
                      if(strncmp(parm,station[j].parm,3)==0 &&
                              parm[4]==station[j].parm[4]) {
 
                              station[j].xadd=xadd;
                              station[j].yadd=yadd;
                              break;
 
                }
 
         }
 
     }
     }
     
 
fclose(fp);          
 
   }
 
for(i=0;i<max_stations;i++) {
 
        for(l=0;l<30;l++)
                    sorted[l]=9999999;
 
        for(m=0;m<max_stations;m++) {
        
               if(i==m)
                     continue; 
 
               dist1=station[i].lat-station[m].lat;
               dist2=(station[i].lon-station[m].lon)*
                  cos((station[i].lat+station[m].lat)/2*conv);
 
               dist=pow(dist1,2)+pow(dist2,2);
           
               for(l=0;l<30;l++) {
 
                       if(dist < sorted[l]) {
 
                               for(h=29; h > l; h--) {
 
                                     sorted[h]=sorted[h-1];
                                              
                                     station[i].index[h]=
                                     station[i].index[h-1];
   
                   }
 
                               sorted[l]=dist;
                                              
                               station[i].index[l]=m;
 
                               break;
 
                 }
                      
 
             }
 
         }
 
      }
 
return;
 
}