get_valid_gages.c

Go to the documentation of this file.

#define _POSIX_SOURCE
 
#include "prototypes_new.h"
 
void get_valid_gages(double maximum_latitude,double minimum_latitude,
                     double center_longitude)
 
{
 
 extern char station_list_file[100];
 extern char forecast_station_list_file[100];
 extern  char hrap_gage_file[100];
 extern int wfo_orig;
 extern int max_fstations;
 extern struct fstation fstation[500];
 extern struct station station[3000];
 int i,k,kk,maxi,maxj,l,m,h;
 struct HRAP hrap;
 struct HRAP irap;
 char dbuf[200];
 extern struct hrap_grid *hrap_grid; 
 float conv=.0174,totlat,sorted[20];
 float dist,dist1,dist2,distance,value;
 float isohlat,isohlon;
 int jj,latindex,lonindex;
 int mini,minj;
 int mm,hh;
 struct stat statbuf;
 int newflag;
 time_t ltime1,ltime2;
 time_t oltime1,oltime2;
 int ier,mer;
 FILE *fz;
 char *p;
 double omaximum_latitude,ominimum_latitude,ocenter_longitude;
 char kbuf[200];
 
 mini=hrap_grid->hrap_minx;
 minj=hrap_grid->hrap_miny;
 
 newflag=0;
 
 if(hrap_gage_file[0]==0)
                     newflag=1;
 
 else {
 
     ier=stat(hrap_gage_file,&statbuf);
     if(ier==-1)
                     newflag=1;
                     
 }
 
 
 mer=stat(station_list_file,&statbuf);
 
 if(mer==-1)
          ltime1=-1;
          
 else
          ltime1=statbuf.st_mtime;
          
 mer=stat(forecast_station_list_file,&statbuf);
 
 if(mer==-1)
          ltime2=-1;
          
 else
          ltime2=statbuf.st_mtime;
 
 
 if(newflag==0) {
 
    fz=fopen(hrap_gage_file,"r");
 
    p=fgets(dbuf,100,fz);
 
    ier=sscanf(dbuf,"%ld %ld %lf %lf %lf",&oltime1,&oltime2,&omaximum_latitude,
                                &ominimum_latitude,&ocenter_longitude);
     
 
    if(oltime1 != ltime1 || oltime2 != ltime2 ||
        omaximum_latitude != maximum_latitude ||
        ominimum_latitude != minimum_latitude ||
        ocenter_longitude != center_longitude)
                      newflag=1;
 
    fclose(fz);
 
 
    }
 
if(newflag==0) {
 
           fz=fopen(hrap_gage_file,"r");
           p=fgets(dbuf,100,fz);                       
 
         }
 
else if(newflag==1) {
 
       fz=fopen(hrap_gage_file,"w+");
 
       fprintf(fz,"%ld %ld %f %f %f\n",ltime1,ltime2,maximum_latitude,
                        minimum_latitude,center_longitude);
 
       }
 
 
 
 for (i=0; i < hrap_grid->maxi; i++)
     {     
 for (k=0; k < hrap_grid->maxj ;k++)
     {
 
        irap.x=mini+i;
        irap.y=minj+k;
        
        hrap=HrapToLatLong(irap);
        /* we have lat long of point we need seasonal isohyet */
 
        for(l=0;l<5;l++)
                    sorted[l]=9999999;
                    
        if(newflag==0){
 
                     m=0;
                     fread(kbuf,sizeof(char),50,fz);
                     kbuf[50]=0;
                     for(mm=0;mm<10;mm++) {
 
                     hrap_grid->gage[i][k].index[mm]=atoi(&kbuf[m]);
 
                     m=m+5;
 
 
           }
 
 
                 }   
                         
        else {
 
        for(mm=0;mm<max_fstations;mm++) {
 
              if(fstation[mm].owner != wfo_orig)
                             continue;
                  
              if(fstation[mm].qpf==-1)
                        continue;
        
               m=fstation[mm].pindex;
 
               dist1=hrap.y-station[m].lat;
               dist2=hrap.x-station[m].lon;
 
               dist=pow(dist1,2)+pow(dist2,2);
 
               for(l=0;l<5;l++) {
 
                       if(dist < sorted[l]) {
 
                               for(hh=4; hh > l; hh--) {
 
                                     sorted[hh]=sorted[hh-1];
                                             
                                              
                                     hrap_grid->gage[i][k].index[hh]=
                                     hrap_grid->gage[i][k].index[hh-1];
   
                   }
 
                               sorted[l]=dist;
                                              
                               hrap_grid->gage[i][k].index[l]=m;
 
                               break;
 
                 }
                      
 
             }
             
 
       }
    
        for(l=5;l<10;l++) 
                    sorted[l]=9999999;
                    
     
        for(mm=0;mm<max_fstations;mm++) {
          
              if(fstation[mm].qpf==-1)
                        continue;
        
               m=fstation[mm].pindex;
 
               dist1=hrap.y-station[m].lat;
               dist2=hrap.x-station[m].lon;
 
               dist=pow(dist1,2)+pow(dist2,2);
 
               for(l=0;l<5;l++) {
 
                       if(dist < sorted[l+5]) {
 
                               for(hh=9; hh > l+5; hh--) {
 
                                     sorted[hh]=sorted[hh-1];
                                              
                                     hrap_grid->gage[i][k].index[hh]=
                                     hrap_grid->gage[i][k].index[hh-1];
                    
 
                   }
 
                               sorted[l+5]=dist;
                                              
                               hrap_grid->gage[i][k].index[l+5]=m;
 
                               break;
 
                 }
                      
 
             }
             
 
    }
    
        for(l=0;l<10;l++)
                     fprintf(fz,"%04d ",hrap_grid->gage[i][k].index[l]);
 
    }
       
     }
 
}
 
fclose(fz);
 
return;
}