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 #include <unistd.h>
  30
  31 #include "fixnode.h"
  32 #include "../Dem2node/mesh.h"
  33 #include "triload.h"
  34
  35
  36 /* load the node information */
  37 void fixnodes(char *filename, struct MESH *m) {
  38     char toname[256];
  39     FILE *fd;
  40     int i;
  41
  42     printf("Fixing up node elevations\n");
  43
  44     /* we could just fix the new nodes as the first "for" statement
  45      * would do, but that leads to funny results (I need to figure out
  46      * why.)  So, let's try fixing all of them */
  47
  48     /* for ( i = origcount + 1; i <= nodecount; i++ ) { */
  49     for ( i = 1; i <= nodecount; i++ ) {
  50         /* printf("Current: %d %.2f %.2f %.2f\n", i, nodes[i][0],
  51                nodes[i][1], nodes[i][2]); */
  52
  53         nodes[i][2] = mesh_altitude(m, nodes[i][0], nodes[i][1]);
  54
  55         /* printf("Fixed: %d %.2f %.2f %.2f\n", i, nodes[i][0],
  56                nodes[i][1], nodes[i][2]); */
  57     }
  58
  59
  60     sprintf(toname, "%s.orig", filename);
  61     printf("Moving %s to %s\n", filename, toname);
  62     rename(filename, toname);
  63
  64     printf("Saving new node file: %s\n", filename);
  65
  66     fd = fopen(filename, "w");
  67
  68     fprintf(fd, "%d 2 1 0\n", nodecount);
  69
  70     for ( i = 1; i <= nodecount; i++ ) {
  71         fprintf(fd, "%d %.2f %.2f %.2f 0\n", i, nodes[i][0],
  72                nodes[i][1], nodes[i][2]);
  73     }
  74
  75     fclose(fd);
  76 }
  77
  78
  79 /* $Log$
  80 /* Revision 1.3  1998/01/09 23:03:08  curt
  81 /* Restructured to split 1deg x 1deg dem's into 64 subsections.
  82 /*
  83  * Revision 1.2  1997/12/02 13:12:07  curt
  84  * Updated to fix every node.
  85  *
  86  * Revision 1.1  1997/11/27 00:17:33  curt
  87  * Initial revision.
  88  *
  89  */