get_forecast_station_list_new.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 struct station zstation[750];
 extern struct station tstation[1000];
 extern int max_stations;
 extern int max_zstations;
 extern int max_tstations;
 extern struct fstation fstation[1000];
 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<1000;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++;
 
            }
 
 
    fstation[i].qpf=-1;
    fstation[i].frz=-1;
    fstation[i].txn=-1;
    fstation[i].pindex=-1;
    fstation[i].zindex=-1;
    fstation[i].tindex=-1;
 
                     
    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;
                                                  
                                        }
                                        
                  }
                  
    }
 
    if(fstation[i].qpf == 1) {
 
    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].pindex=j;
                  
                     for(n=0;n<60;n++) {
                     
                         fdata[n].fstn[j].rain=calloc(5,sizeof(float));                    
                         
                     
                     for(k=0;k<5;k++)
                         fdata[n].fstn[j].rain[k]=-9999;
                         
                                  }
                  
                break;        
                  
            }
            
            
        }         
                     
         } 
 
    if(fstation[i].frz == 1) {
 
    for(j=0;j<max_zstations;j++) {
 
                 if(strcmp(fstation[i].hb5,zstation[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].zindex=j;
                     
 
             
                     for(n=0;n<60;n++) {
                            
             
                
                         fdata[n].fstn[j].frz=calloc(5,sizeof(float));                     
                         
                     for(k=0;k<5;k++)         
                         fdata[n].fstn[j].frz[k]=-9999;
                                   
                                  }
                  
            break;    
                  
            }
            
                        
        }         
              
                     
         } 
 
 
    if(fstation[i].txn == 1) {
 
 
 
    for(j=0;j<max_tstations;j++) {
 
                 if(strcmp(fstation[i].hb5,tstation[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].tindex=j;
                     
                     for(n=0;n<60;n++) {
                                         
                         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].tx[k]=-9999;
                         fdata[n].fstn[j].tn[k]=-9999;
                         
                                        }
 
 
                                  }
                  
            break;    
                  
            }
            
                        
        }         
              
                     
         } 
     
  if(fstation[i].tindex == -1)
                fstation[i].txn=-1;
     
  if(fstation[i].pindex == -1)
                fstation[i].qpf=-1;
         
  if(fstation[i].zindex == -1)
                fstation[i].frz=-1;
     
       
  i++;
                    
 
  }
 
max_fstations=i;
fclose(fp);
 
 
return;
 
}