get_topo.c

Go to the documentation of this file.

#include "prototypes.h"
 
void get_topo(char *fname)
 
{
 
 int i,j,ier,maxj,maxi,max_value,k;
 float lat,lon;
 extern struct top *topo;
 FILE *fp;
 char *p,kbuf[10000];
 float max_lat,max_lon,delta_lat,delta_lon,total_lat,total_lon;
 
 topo=(struct top *) calloc(1,sizeof(struct top));
 if(topo==NULL) {
 
                   printf("could not allocate topo space\n");
                   exit(1);
 
                  }          
 
 max_value=0;
 
 fp=fopen(fname,"r");
 
 if(fp==NULL) {
 
              printf("could not open %s\n",fname);
              exit(1);
 
              }
 
 p=fgets(kbuf,80,fp);
 
 ier=sscanf(kbuf,"%f %f %f %f %f %f",&max_lat,&max_lon,&total_lat,&total_lon,&delta_lat,&delta_lon);
 delta_lat=delta_lat/60.;
 delta_lon=delta_lon/60.;
 
 maxi=total_lon/delta_lon+1;
 maxj=total_lat/delta_lat+1;
 
 topo->maxi=maxi;
 topo->maxj=maxj;
 topo->max_lat=max_lat;
 topo->max_lon=max_lon;
 topo->total_lat=total_lat;
 topo->total_lon=total_lon;
 topo->delta_lat=delta_lat;
 topo->delta_lon=delta_lon;
 
 topo->value=(short int**) calloc(maxi,sizeof(short int *));
 
 if( topo->value==NULL) {
 
                            printf("no memory for topo array\n");
                            exit(1);
 
              }
 
 for(i=0;i<maxi;i++) {
 
      topo->value[i]=(short int *) calloc(maxj,sizeof(short int));
 
      if( topo->value[i]==NULL) {
 
                printf("no memory for topo array\n");
                exit(1);
 
          }
 
    }
 
 for(j=0;j<maxj;j++) {
 
 fread(kbuf,sizeof(char),maxi*4,fp);
 
 k=0;
 
 for(i=0;i<maxi;i++) {
         
         topo->value[i][j]=atoi(&kbuf[k]);
 
         k=k+4;         
 
    }
 
   }
 
fclose(fp);
printf("return from topo\n");
return;
 
}