get_forecast_station_list.c

Go to the documentation of this file.

#include "prototypes_new.h"
 
void get_forecast_station_list(char *fname)
 
{
 
 extern struct fdata fdata[60];
 extern int isom;
 extern struct tag tag[20];
 extern struct station station[3000];
 extern int max_stations;
 extern struct fstation fstation[500];
 extern int max_fstations;
 FILE *fp;
 char *p,kbuf[150];
 int i,ier;
 int m,l,h,k,j;
 int kk,ii;
 char dbuf[10],tbuf[10],wfo[10];
 extern char *txbuf[3];
 int n;
 
 max_fstations=0;
 
 fp=fopen(fname,"r");
 
 if(fp==NULL) {
 
              printf("could not open forecast station_list %s\n",fname);
              exit(1);
 
              }
 
 i=0;
 
 for(ii=0;ii<500;ii++) {
 
    p=fgets(kbuf,100,fp);
    if(p==NULL)
              break;
      
    p=strchr(kbuf,'\n');
 
    if(p!=NULL)
             *p=0;
 
 
    *(txbuf[0])=0;
    *(txbuf[1])=0;
    *(txbuf[2])=0;
 
    ier=sscanf(kbuf,"%s %s %s %s %s",
               fstation[i].hb5,wfo,txbuf[0],txbuf[1],txbuf[2]);
    
    if(ier < 1)
               continue;
 
    j=0;
    while(wfo[j] != 0) {
 
              wfo[j]=tolower(wfo[j]);
              j++;
 
            }
 
    for(j=0;j<max_stations;j++) {
 
                 if(strcmp(fstation[i].hb5,station[j].hb5)==0 ) {  
                                            
                     fstation[i].owner=-1;
                           
                     for(m=0;m<20;m++) {
 
                          if(strcmp(wfo,tag[m].wfo)==0) {
 
                          fstation[i].owner=m;
                          break;
 
                        }
      
            }
            
                     fstation[i].qpf=-1;
                     fstation[i].frz=-1;
                     fstation[i].txn=-1;
                     
                     fstation[i].index=j;
                     
                     for(n=0;n<60;n++) {
                     
                         fdata[n].fstn[j].rain=calloc(5,sizeof(float));                    
                         fdata[n].fstn[j].frz=calloc(5,sizeof(float));                     
                         fdata[n].fstn[j].tx=calloc(5,sizeof(int));                     
                         fdata[n].fstn[j].tn=calloc(5,sizeof(int));
                         
                     
                       for(k=0;k<5;k++) {
                       
                         fdata[n].fstn[j].rain[k]=-9999;
                         fdata[n].fstn[j].frz[k]=-9999;
                         fdata[n].fstn[j].tx[k]=-9999;
                         fdata[n].fstn[j].tn[k]=-9999;
                         
                                        }
 
 
                                  }
                                                           
                     if(ier==5 || ier==2) {
 
                                  fstation[i].qpf=1;
                                  fstation[i].frz=1;
                                  fstation[i].txn=1;
          
                }
  
                     else {
                     
                                for(m=0;m<3;m++) {
                        
                                  if(*(txbuf[m])==0)
                                                continue;
 
                                       if(strcmp(txbuf[m],"HZ")==0)
                                                  fstation[i].frz=1;
 
                                       else if(strcmp(txbuf[m],"PP")==0)
                                                  fstation[i].qpf=1;
                                   
                                       else if(strcmp(txbuf[m],"TA")==0)
                                                  fstation[i].txn=1;
                                                  
                                        
                                        
                  }
                  
              }    
 
                     i++;
                     break;
 
                           }
 
                }
 
  }
 
max_fstations=i;
fclose(fp);
 
return;
 
}