build_hrap.c

Go to the documentation of this file.

#include "prototypes.h"
 
struct top *topo;  
struct h_grid *h_grid;
double pgon[5000][2];
 
 
void main(int argc, char **argv)
 
{
 
FILE *fp,*fr,*fh,*fb,*fn;
double max_lat,min_lat,cen_lon,lat,lon,point[2],area2,area3,felev,area4;
double totlat,totlon;
int i,numverts,j,ier,ii,jj,lower,mid,upper,numpoints,lpoints,mpoints,upoints,good,zones,glacial,gpoints,icount;
char *p,kbuf[100],ibuf[100],dbuf[100],hb5[100];
long fpos,fpos1;
 
topo=NULL;
h_grid=NULL;
 
p=pars_argv(argv,"max_lat=");
 
if(p==NULL) {
 
            printf("no max_lat specified\n");
            exit(1);
 
      }
 
max_lat=atof(p);
 
p=pars_argv(argv,"min_lat=");
 
if(p==NULL) {
 
            printf("no min_lat specified\n");
            exit(1);
 
      }
 
 
min_lat=atof(p);
 
p=pars_argv(argv,"cen_lon=");
 
if(p==NULL) {
 
            printf("no cen_lon specified\n");
            exit(1);
 
      }
 
 
cen_lon=atof(p);
 
p=pars_argv(argv,"index_file=");
 
if(p==NULL) {
 
            printf("no index file specified\n");
            exit(1);
 
      }
 
 
fr=fopen(p,"r");
 
if(fr==NULL) {
 
            printf("could not open index file %s\n",p);
 
            exit(1);
 
      }
 
 
p=pars_argv(argv,"topo=");
 
if(p==NULL){
 
           printf("no topo file specified\n");
 
           exit(1);
 
     }
 
get_topo(p);
 
p=pars_argv(argv,"hrap=");
 
if(p==NULL){
 
           printf("no hrap file specified\n");
 
           exit(1);
 
     }
 
fh=fopen(p,"w");
 
p=pars_argv(argv,"basin=");
 
if(p==NULL){
 
           printf("no basin file specified\n");
 
           exit(1);
 
     }
 
fb=fopen(p,"w");
 
p=pars_argv(argv,"basin_name=");
 
if(p==NULL){
 
           printf("no basin_name file specified\n");
 
           exit(1);
 
     }
 
fn=fopen(p,"w");
 
get_coord(max_lat,min_lat,cen_lon);
 
for(;;) {
 
        icount=0;
        p=fgets(ibuf,100,fr);
 
        if(p==NULL)
                exit(1);
 
        p=pars_line(ibuf,"i=",hb5);
 
       
 
        if(p==NULL){
 
             printf("no input file specified\n");
 
             continue;
 
     }
 
         p=strchr(hb5,'\n');
         if(p!=NULL)
                *p=0;
 
        fp=fopen(hb5,"r");
 
        if(fp==NULL) {
 
              printf("could not open input file %s\n",p);
 
              continue;
 
      }
 
        p=pars_line(ibuf,"upper=",dbuf);
 
        if(p==NULL){
 
           printf("no upper area specified - assume none exists\n");
           area3=2000000;
 
     }
 
         else
           area3=atof(p);
 
        p=pars_line(ibuf,"middle=",dbuf);
 
        if(p==NULL){
 
           printf("no middle area specified - assume none exists\n");
           area2=2000000;
 
     }
 
         else
           area2=atof(p);
 
        if(area2 < 30000 && area3 > 30000) {
 
                           printf("middle area with no upper area\n");
 
                           exit(1);
 
             }
 
        p=pars_line(ibuf,"glacial=",dbuf);
 
        if(p==NULL){
 
           printf("no glacial specified - assume none exists\n");
           area4=2000000;
 
     }
 
         else
           area2=atof(p);
 
         fpos=ftell(fh);
         
         fprintf(fh,"                    \n");
 
         zones=1;
 
         if(area2 < 30000)
                      zones++;
 
         if(area3 < 30000)
                      zones++;
 
         if(area4 < 30000)
                      zones++;
 
         p=fgets(kbuf,100,fp);
         p=strchr(kbuf,'\n');
         *p=0;
 
         strcat(kbuf,"  ");
         strcat(kbuf,hb5);
         strcat(kbuf,"\n");
 
         sscanf(kbuf,"%d",&numverts);
         fputs(kbuf,fb);
 
         sprintf(kbuf,"%d   %s   0   0  1  -center\n",icount,hb5);         
         icount++;
         fputs(kbuf,fn);        
 
         i=0;
         for(i=0;i<numverts;i++) {
 
             p=fgets(kbuf,100,fp);
 
             if(p==NULL)
                     break;
 
             fputs(kbuf,fb);
             ier=sscanf(kbuf,"%lf %lf",&lon,&lat);
 
             if(ier != 2)
                continue;
 
             pgon[i][0]=-lon;
             pgon[i][1]=lat;
 
             }
 
         fclose(fp);
 
         printf("number of points in region %d\n",numverts);
         printf("area2=%f area3=%f\n",area2,area3);
 
         numpoints=0;lpoints=0;mpoints=0;upoints=0;gpoints=0;
         totlat=0;totlon=0;
 
         for(i=0;i<h_grid->maxi-1;i++)
         {
         for(j=0;j<h_grid->maxj-1;j++)
         {
 
           good=0;
           for(ii=i;ii<=i+1;ii++){
           for(jj=j;jj<=j+1;jj++){
 
               point[0]=h_grid->latlon[ii][jj].x;
               point[1]=h_grid->latlon[ii][jj].y;
 
               if(macmartintest(pgon,numverts,point)!=0)
                    good++;
 
 
         }
     }
   
          if(good==0)
                   continue;
 
           numpoints++;
 
           printf("point %f %f\n",point[0],point[1]);
           totlon=totlon+point[0];
           totlat=totlat+point[1];
 
           lower=0;
           mid=0;
           upper=0;
           glacial=0;
  
           for(ii=i;ii<=i+1;ii++){
           for(jj=j;jj<=j+1;jj++){
 
             felev=(float)h_grid->elev[ii][jj] * 32.1;
 
             if(zones==1)
                          lower++;
 
             else if(zones==2) {
 
                        if(felev < area3)
                                  lower++;
 
                        else
                                  upper++;
 
              }
 
             else if(zones==3) {
 
                        if(felev < area2)
                                  lower++;
 
                        else if(felev < area3)
                                  mid++;
 
                        else
                                  upper++;
 
              }
 
 
 
             else if(zones==4) {
 
                        if(felev < area2)
                                  lower++;
 
                        else if(felev < area3)
                                  mid++;
 
                        else if(felev < area4)
                                  upper++;
 
                        else
                                  glacial++;
 
              }
 
 
            }
     }
 
        fprintf(fh,"%d %d ",h_grid->box[i][j].x,h_grid->box[i][j].y);
 
        if(lower > 0) {
 
           fprintf(fh,"1 ");
           lpoints++;
 
            } 
 
        if(mid > 0) {
 
           fprintf(fh,"2 ");
           mpoints++;
 
     }
 
        if(upper > 0) {
 
           fprintf(fh,"3 ");
 
           upoints++;
 
     }
 
        if(glacial > 0) {
 
           fprintf(fh,"4 ");
 
           upoints++;
 
     }
 
        fprintf(fh,"\n");
 
             
    }
 
       }
 
      if(numpoints > 0) {
 
      totlat=totlat/(float)numpoints;
      totlon=totlon/(float)numpoints;
 
    }
 
      fprintf(fn,"%f   %f\n",-totlon,totlat);
 
      printf("%f %f\n",totlat,totlon);
 
      fpos1=ftell(fh);
      fseek(fh,fpos,SEEK_SET);
 
      fprintf(fh,"%d %d %s",numpoints,zones,hb5);
 
      printf("number of hrap boxes inside region %d\n",numpoints);
      printf("lower %d middle %d upper %d\n",lpoints,mpoints,upoints);
 
      fseek(fh,fpos1,SEEK_SET);
      icount=0;
 
 
      }
 
fclose(fh);
fclose(fb);
 
}