get_NEXRAD_coord.c

Go to the documentation of this file.

#include "prototypes.h"
 
void get_NEXRAD_coord()
 
{
 extern Widget scrollbar;
 extern Widget top_level;
 extern Display *display;
 extern int posit;
 extern Widget s_text;
 int i,k,kk;
 struct HRAP hrap;
 struct HRAP irap;
 float cx,cy,lat,lon;
 extern struct nexcoord nexcoord[2];
 FILE *fp;
 char *p,kbuf[120],dbuf[80];
 int slider_size,scrollm,ivalue,increment,page_increment;
 
 fp=fopen("NEXRAD_radars","r");
 
 if(fp==NULL) {
 
              printf("could not open NEXRAD_radars\n");
              exit(1);
 
              }
 
 kk=0;
 for(;;) {
 
     p=fgets(kbuf,80,fp);
     if(p==NULL) {
 
                 fclose(fp);
                 return;
 
           }
 
     p=pars_line(kbuf,"radar=",dbuf);
 
     if(p==NULL)
               continue;
                   
     strcpy(nexcoord[kk].name,dbuf); 
 
     p=pars_line(kbuf,"latitude=",dbuf);
 
     if(p==NULL)
               continue;
 
     lat=atof(p);
 
     p=pars_line(kbuf,"longitude=",dbuf);
 
     if(p==NULL)
               continue;
 
     strcpy(dbuf,"Coordinates for ");
     strcat(dbuf,nexcoord[kk].name);
     strcat(dbuf,"\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);
 
     XFlush(display);
     XmUpdateDisplay(top_level);
 
     lon=atof(p);
 
     if(lon > 0)
              lon=-lon; 
    
     hrap=LatLongToHrap(lat,-lon);
 
     cx=(int)hrap.x+1;
     cy=(int)hrap.y;
 
     for (i=0; i < 131; i++)
     {     
     for (k=0; k < 131 ;k++)
     {
      
        irap.x=cx + (float)(k - 65);
        irap.y=cy + (float)(65 -  i);
        hrap=HrapToLatLong(irap);
 
        get_NEXRAD_pixels(i,k,hrap.y,hrap.x,kk);
 
    }
 
   }
 
}
return;
 
}
void get_NEXRAD_pixels(int i, int k, float lat, float lon, int kk)
 
  {
 
     extern struct dval dval;
     long r,s;
     float conv=.0174;
     extern struct nexcoord nexcoord[2];
     
     r=dval.a * cos(lat*conv)/(1+sin(lat*conv))
       * cos((lon-dval.lo-90)*conv) + dval.xo +.5;
 
     nexcoord[kk].boxes[130-i][k].x=r;
 
     s=dval.a * cos(lat*conv)/(1+sin(lat*conv))
       * sin((lon-dval.lo-90)*conv) + dval.yo + .5;
       
     nexcoord[kk].boxes[130-i][k].y=s;
 
     return;
 
   }