get_basin_data.c

Go to the documentation of this file.

#include "prototypes.h"
 
void get_basin_data(char *basin_file,char *hrap_file)
 
{
 
extern struct map map[500];
extern struct dval dval;
FILE *fr,*fp;
char *p,ibuf[100];
int ib,ier,l,numpts,ip,x,y,ip2,ip3,ip4,zones,i,maxib,dum,mm;
double lat,lon;
long r,s;
float conv=.0174;
 
if(basin_file[0]==0 ) {
 
     printf("no basin_file specified\n");
     exit(1);
 
   }
 
if(hrap_file[0]==0 ) {
 
     printf("no basin_file specified\n");
     exit(1);
 
   }
 
fr=fopen(basin_file,"r");
 
if(fr==NULL) {
 
            printf("could not open basin_file %s\n",basin_file);
 
            exit(1);
 
      }
 
fp=fopen(hrap_file,"r");
 
if(fp==NULL) {
 
            printf("could not open hrap_file %s\n",hrap_file);
 
            exit(1);
 
      }
 
ib=0;
 
for(;;) {
 
        p=fgets(ibuf,80,fr);
 
        if(p==NULL)
                  break;
 
        p=strchr(ibuf,'\n');
        if(p!=NULL)
                  *p=0;
 
        ier=sscanf(ibuf,"%d %d %s",&numpts,&dum,map[ib].hb5);
 
        map[ib].basin_points=numpts;
 
        map[ib].basin=calloc(numpts,sizeof(XPoint));
 
        for(l=0;l<numpts;l++) {
 
               p=fgets(ibuf,80,fr);
 
               if(p==NULL)
                        break;
 
               ier=sscanf(ibuf,"%lf %lf",&lon,&lat);
 
               lon=-lon;
 
               r=dval.a * cos(lat*conv)/(1+sin(lat*conv))
                       * cos((lon-dval.lo-90)*conv) + dval.xo +.5;
 
               s=dval.a * cos(lat*conv)/(1+sin(lat*conv))
                       * sin((lon-dval.lo-90)*conv) + dval.yo + .5;
 
               map[ib].basin[l].x=r;
               map[ib].basin[l].y=s;
 
         }
 
      p=fgets(ibuf,80,fr);
 
      ib++;
 
      }
 
map[ib].hb5[0]=0;     
maxib=ib;
 
for(ib=0;ib<maxib;ib++) {
  
        p=fgets(ibuf,80,fp);
 
        if(p==NULL)
                 break;
 
        ier=sscanf(ibuf,"%d %d",&numpts,&zones);
 
        map[ib].hrap_points=numpts;
 
        map[ib].hrap_data=calloc(numpts,sizeof(struct hrap_data));
 
        for(mm=0;mm<4;mm++) 
                     map[ib].zones[mm]=-1;
 
        for(l=0;l<numpts;l++) {
 
              for(mm=0;mm<4;mm++) 
                     map[ib].hrap_data[l].zone[mm]=-1;
 
              p=fgets(ibuf,100,fp);
 
              if(p==NULL)
                          break;
 
              ier=sscanf(ibuf,"%d %d %d %d %d %d\n",&x,&y,&ip,&ip2,&ip3,&ip4);
 
              if(ip < 0 || ip > 4) {
 
                     printf("hrap error\n");
                     exit(1);
 
           }
              
              map[ib].hrap_data[l].x=x;
              map[ib].hrap_data[l].y=y;
              map[ib].hrap_data[l].zone[ip-1]=1;
              map[ib].zones[0]=1;              
 
              if(ier>=4) {
 
                  if(ip2 < 0 || ip2 > 4) {
 
                     printf("hrap error\n");
                     exit(1);
 
           }
              
                    map[ib].hrap_data[l].zone[ip2-1]=1;
                    map[ib].zones[ip2-1]=1;
 
          }
 
 
              if(ier>=5) {
 
              if(ip3 < 0 || ip3 > 4) {
 
                     printf("hrap error\n");
                     exit(1);
 
           }
 
                    map[ib].hrap_data[l].zone[ip3-1]=1;
                    map[ib].zones[ip3-1]=1;
 
          }
 
              if(ier>=6) {
 
              if(ip4 < 0 || ip4 > 4) {
 
                     printf("hrap error\n");
                     exit(1);
 
           }
              
                    map[ib].hrap_data[l].zone[ip4-1]=1;
                    map[ib].zones[ip4-1]=1;
 
          }
    
         }
 
 
         for(i=0;i<51;i++)
                 map[ib].maps_done[i]=-1;       
 
 
      }
 
fclose(fr);
fclose(fp);
 
return;
 
}