FixNode/fixnode.c

   1 /* fixnode.c -- traverse the node file and fix the elevation of all the new
   2  *              interpolated points.
   3  *
   4  * Written by Curtis Olson, started November 1997.
   5  *
   6  * Copyright (C) 1997  Curtis L. Olson  - curt@infoplane.com
   7  *
   8  * This program is free software; you can redistribute it and/or modify
   9  * it under the terms of the GNU General Public License as published by
  10  * the Free Software Foundation; either version 2 of the License, or
  11  * (at your option) any later version.
  12  *
  13  * This program is distributed in the hope that it will be useful,
  14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16  * GNU General Public License for more details.
  17  *
  18  * You should have received a copy of the GNU General Public License
  19  * along with this program; if not, write to the Free Software
  20  * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  21  *
  22  * $Id$
  23  * (Log is kept at end of this file)
  24  */
  25
  26
  27 #include <stdio.h>
  28 #include <string.h>
  29
  30 #include "fixnode.h"
  31 #include "../Dem2node/mesh.h"
  32 #include "triload.h"
  33
  34
  35 /* load the node information */
  36 void fixnodes(char *basename, struct MESH *m) {
  37     char file[256];
  38     FILE *fd;
  39     int i;
  40
  41     printf("Fixing up node elevations\n");
  42
  43     /* we could just fix the new nodes as the first "for" statement
  44      * would do, but that leads to funny results (I need to figure out
  45      * why.)  So, let's try fixing all of them */
  46
  47     /* for ( i = origcount + 1; i <= nodecount; i++ ) { */
  48     for ( i = 1; i <= nodecount; i++ ) {
  49         /* printf("Current: %d %.2f %.2f %.2f\n", i, nodes[i][0],
  50                nodes[i][1], nodes[i][2]); */
  51
  52         nodes[i][2] = mesh_altitude(m, nodes[i][0], nodes[i][1]);
  53
  54         /* printf("Fixed: %d %.2f %.2f %.2f\n", i, nodes[i][0],
  55                nodes[i][1], nodes[i][2]); */
  56     }
  57
  58     strcpy(file, basename);
  59     strcat(file, ".1.node");
  60
  61     printf("Overwriting original node file: %s\n", file);
  62
  63     fd = fopen(file, "w");
  64
  65     fprintf(fd, "%d 2 1 0\n", nodecount);
  66
  67     for ( i = 1; i <= nodecount; i++ ) {
  68         fprintf(fd, "%d %.2f %.2f %.2f 0\n", i, nodes[i][0],
  69                nodes[i][1], nodes[i][2]);
  70     }
  71
  72     fclose(fd);
  73 }
  74
  75
  76 /* $Log$
  77 /* Revision 1.2  1997/12/02 13:12:07  curt
  78 /* Updated to fix every node.
  79 /*
  80  * Revision 1.1  1997/11/27 00:17:33  curt
  81  * Initial revision.
  82  *
  83  */