get_zgage_grid.c

Go to the documentation of this file.

#define _POSIX_SOURCE
#include "prototypes_new.h"
 
void get_zgage_grid(double maximum_latitude,double minimum_latitude,
                    double center_longitude)
 
{
 
 extern char station_list_file[100];
 extern char hrap_zgage_file[100];
 extern struct station zstation[500];
 extern struct station station[3000];
 extern Widget scrollbar;
 extern Widget top_level;
 extern int posit;
 extern int max_zstations;
 extern int max_stations;
 extern Widget s_text;
 int i,k,kk,maxi,maxj,l,m,h;
 struct HRAP hrap;
 struct HRAP irap;
 float cx,cy;
 extern struct dval dval;
 char dbuf[200];
 extern struct hrap_grid *hrap_grid; 
 long r,s;
 float conv=.0174,totlat,sorted[30];
 float dist,dist1,dist2,distance,value;
 float isohlat,isohlon;
 int jj,latindex,lonindex;
 int slider_size,scrollm,ivalue,increment,page_increment;
 struct stat statbuf;
 double omaximum_latitude,ominimum_latitude,ocenter_longitude;
 int ier,newflag,mm,mer;
 FILE *fz;
 time_t ost_mtime;
 char *p,kbuf[200];
 struct tm *tm;
 
newflag=0;
 
if(hrap_zgage_file[0]==0)
          strcpy(hrap_zgage_file,"map_zqc.grid");
 
ier=stat(hrap_zgage_file,&statbuf);
 
/* get creation time of station_list file */
 
mer=stat(station_list_file,&statbuf);
 
/* grid file does not exist */
 
if(ier==-1)
          newflag=2;
 
else {
 
     fz=fopen(hrap_zgage_file,"r");
 
     p=fgets(dbuf,100,fz);
 
     ier=sscanf(dbuf,"%ld %lf %lf %lf",&ost_mtime,&omaximum_latitude,
                                &ominimum_latitude,&ocenter_longitude);
 
     tm=gmtime(&statbuf.st_mtime);
 
     if(ost_mtime != statbuf.st_mtime ||
        omaximum_latitude != maximum_latitude ||
        ominimum_latitude != minimum_latitude ||
        ocenter_longitude != center_longitude)
                      newflag=1;
 
     fclose(fz);
 
      }
 
if(newflag==0) {
 
           fz=fopen(hrap_zgage_file,"r");
           p=fgets(dbuf,100,fz);
 
         }
 
else if(newflag==1)
       fz=fopen(hrap_zgage_file,"w+");
 
else
       fz=fopen(hrap_zgage_file,"w");
 
if(newflag!=0) {
 
             strcpy(dbuf,"ZGage to grid recalculation\n");
             posit=posit+strlen(dbuf);
             XmTextInsert(s_text,posit,dbuf);
 
             XtVaGetValues(scrollbar,XmNmaximum,&scrollm,NULL);
 
             XmScrollBarGetValues(scrollbar,&ivalue,&slider_size,
                     &increment,&page_increment);
 
             XmScrollBarSetValues(scrollbar,scrollm-slider_size,
                      slider_size, increment,page_increment,True);
 
             XmUpdateDisplay(top_level);
 
             fprintf(fz,"%ld %f %f %f\n",statbuf.st_mtime,maximum_latitude,
                        minimum_latitude,center_longitude);
 
       }
    
for (i=0; i < hrap_grid->maxi; i++){     
 
for (k=0; k < hrap_grid->maxj ;k++){
 
        irap.x=hrap_grid->hrap_minx+i;
        irap.y=hrap_grid->hrap_miny+k;
        
        hrap=HrapToLatLong(irap);
        
        for(l=0;l<30;l++)
                    sorted[l]=9999999;
 
        if(newflag==0){
 
                     m=0;
                     fread(kbuf,sizeof(char),25,fz);
                     kbuf[25]=0;
                     for(mm=0;mm<5;mm++) {
 
                     hrap_grid->gage[i][k].zindex[mm]=atoi(&kbuf[m]);
 
                     m=m+5;
 
 
           }
 
 
                 }            
 
        else {
 
        for(m=0;m<max_zstations;m++) {
        
               dist1=hrap.y-zstation[m].lat;
               dist2=(hrap.x-zstation[m].lon)*cos((hrap.y+zstation[m].lat)/2*conv);
              
           dist=pow(dist1,2)+pow(dist2,2);
  
               for(l=0;l<5;l++) {
 
                       if(dist < sorted[l]) {
 
                               for(h=4; h > l; h--) {
 
                                     sorted[h]=sorted[h-1];
                                              
                                     hrap_grid->gage[i][k].zindex[h]=
                                     hrap_grid->gage[i][k].zindex[h-1];
   
                   }
 
                               sorted[l]=dist;
                                              
                               hrap_grid->gage[i][k].zindex[l]=m;
 
                               break;
 
                 }
                      
 
             }
 
         }
   
             for(l=0;l<5;l++)
                     fprintf(fz,"%04d ",hrap_grid->gage[i][k].zindex[l]);
 
      }
 
 
        }
 
   }
 
fclose(fz);
 
for(i=0;i<max_stations;i++) {
 
        for(l=0;l<5;l++)
                    sorted[l]=9999999;
 
        for(m=0;m<max_zstations;m++) {
        
               if(i==m)
                     continue; 
 
               dist1=station[i].lat-zstation[m].lat;
               dist2=station[i].lon-zstation[m].lon;
 
               dist=pow(dist1,2)+pow(dist2,2);
 
               for(l=0;l<5;l++) {
 
                       if(dist < sorted[l]) {
 
                               for(h=4; h > l; h--) {
 
                                     sorted[h]=sorted[h-1];
                                              
                                     station[i].zindex[h]=
                                     station[i].zindex[h-1];
   
                   }
 
                               sorted[l]=dist;
                                              
                               station[i].zindex[l]=m;
 
                               break;
 
                 }
                      
 
             }
 
         }
 
      }
 
return;
 
}