get_hrap_coord.c

Go to the documentation of this file.

#define _POSIX_SOURCE
#include "prototypes.h"
 
void get_hrap_coord(double maximum_latitude,double minimum_latitude,double center_longitude)
 
{
 
 extern struct grid *grid;
 extern char station_list_file[256];
 extern char grid_file[100];
 extern Widget scrollbar;
 extern Widget top_level;
 extern int posit;
 extern Widget s_text;
 int i,k,kk,maxi,maxj;
 struct HRAP hrap;
 struct HRAP irap;
 float cx,cy;
 extern struct dval dval;
 char dbuf[200],kbuf[200];
 extern struct hrap_grid *hrap_grid; 
 long r,s;
 extern struct pcp *pcp;
 extern struct pcp *spf;
 float conv=.0174,totlat;
 float dist,dist1,dist2,distance,value;
 float isohlat,isohlon;
 int jj,latindex,lonindex;
 int mini,minj;
 int slider_size,scrollm,ivalue,increment,page_increment;
 double pw;
 int l,ll;
 FILE *fp;
 char *p;
 float lat,lon;
 
 hrap_grid=(struct hrap_grid *) calloc(1,sizeof(struct hrap_grid));
 if(hrap_grid==NULL) {
 
                   printf("could not allocate hrap_grid space\n");
                   exit(1);
 
         }
         
pcp=(struct pcp *) calloc(1,sizeof(struct pcp));
 if(pcp==NULL) {
 
                   printf("could not allocate pcp space\n");
                   exit(1);
 
                  }          
 
 spf=(struct pcp *) calloc(1,sizeof(struct pcp));
 if(spf==NULL) {
 
                   printf("could not allocate pcp space\n");
                   exit(1);
 
                  }          
 
 grid=(struct grid *) calloc(1,sizeof(struct grid));
 if(grid==NULL) {
 
                   printf("could not allocate grid space\n");
                   exit(1);
 
                  }          
 
 
 maxi=-1;maxj=-1;mini=999999;minj=99999;
 
 strcpy(dbuf,"HRAP Coordinates\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);
 
 /* need longitude of 4 corners of display */
 
 totlat=maximum_latitude-minimum_latitude+2.0;
 
 hrap=LatLongToHrap((float)maximum_latitude,(float)center_longitude-5);
 
 cx=(int)hrap.x+1;
 cy=(int)hrap.y;
 
 if(cx > maxi)
           maxi=cx;
 
 if(cy > maxj)
           maxj=cy;
 
 if(cx < mini)
           mini=cx;
 
 if(cy < minj)
           minj=cy;
 
 hrap=LatLongToHrap((float)maximum_latitude,(float)center_longitude+5);
 
 cx=(int)hrap.x+1;
 cy=(int)hrap.y;
 
 if(cx > maxi)
           maxi=cx;
 
 if(cy > maxj)
           maxj=cy;
 
 if(cx < mini)
           mini=cx;
 
 if(cy < minj)
           minj=cy;
 
 hrap=LatLongToHrap((float)minimum_latitude,(float)center_longitude-5);        
 cx=(int)hrap.x+1;
 cy=(int)hrap.y;
 
 if(cx > maxi)
           maxi=cx;
 
 if(cy > maxj)
           maxj=cy;
 
 if(cx < mini)
           mini=cx;
 
 if(cy < minj)
           minj=cy;
 
 hrap=LatLongToHrap((float)minimum_latitude,(float)center_longitude+5);
 
 cx=(int)hrap.x+1;
 cy=(int)hrap.y;
 
 if(cx > maxi)
           maxi=cx;
 
 if(cy > maxj)
           maxj=cy;
 
 if(cx < mini)
           mini=cx;
 
 if(cy < minj)
           minj=cy;
 
/*allocate necessary space */
 
hrap_grid->hrap_minx=mini;
hrap_grid->hrap_miny=minj;
hrap_grid->maxi=maxi-mini;
hrap_grid->maxj=maxj-minj;
 
maxi=hrap_grid->maxi;
maxj=hrap_grid->maxj;
 
grid->value=(float **) calloc(maxi,sizeof(float *));
if(grid->value==NULL) {
 
                   printf("could not allocate grid space\n");
                   exit(1);
 
                  }   
 
pcp->value=(short int **) calloc(maxi,sizeof(short int *));
if(pcp->value==NULL) {
 
                   printf("could not allocate pcp space\n");
                   exit(1);
 
                  }          
 
spf->value=(short int **) calloc(maxi,sizeof(short int *));
if(spf->value==NULL) {
 
                   printf("could not allocate pcp space\n");
                   exit(1);
 
                  }          
    
                   
hrap_grid->coord=(struct coord**) calloc(maxi,sizeof(struct coord *));
 
if(hrap_grid->coord==NULL) {
 
                            printf("no memory for hrap array\n");
                            exit(1);
 
              }
 
for(i=0;i<maxi;i++) {
 
      hrap_grid->coord[i]=(struct coord *) calloc(maxj,sizeof(struct coord));
     
      
      if(hrap_grid->coord[i]==NULL) {
 
                printf("no memory for hrap array\n");
                exit(1);
 
          }
 
    }
 
 
 for(i=0;i<maxi;i++) {
 
          grid->value[i]=(float *) calloc(maxj,sizeof(float)); 
 
          if(grid->value[i]==NULL) {
 
                printf("no memory for grid array\n");
                exit(1);
 
          }
 
    }
 
 for(i=0;i<maxi;i++) {
 
          pcp->value[i]=(short int *) calloc(maxj,sizeof(short int)); 
 
          if(pcp->value[i]==NULL) {
 
                printf("no memory for pcp array\n");
                exit(1);
 
          }
 
    }
 
 for(i=0;i<maxi;i++) {
 
          spf->value[i]=(short int *) calloc(maxj,sizeof(short int)); 
 
          if(spf->value[i]==NULL) {
 
                printf("no memory for pcp array\n");
                exit(1);
 
          }
 
    }
 
 
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);
 
        r=dval.a * cos(hrap.y*conv)/(1+sin(hrap.y*conv))
          * cos((hrap.x-dval.lo-90)*conv) + dval.xo +.5;
   
        hrap_grid->coord[i][k].x=r;
  
        s=dval.a * cos(hrap.y*conv)/(1+sin(hrap.y*conv))
          * sin((hrap.x-dval.lo-90)*conv) + dval.yo + .5;
       
        hrap_grid->coord[i][k].y=s; 
 
   }
 
   }
 
return;
 
}