/* fixnode.c -- traverse the node file and fix the elevation of all the new
 *              interpolated points.
 *
 * Written by Curtis Olson, started November 1997.
 *
 * Copyright (C) 1997  Curtis L. Olson  - curt@infoplane.com
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
 *
 * $Id$
 * (Log is kept at end of this file)
 */


#include <stdio.h>
#include <string.h>

#include "fixnode.h"
#include "../Dem2node/mesh.h"
#include "triload.h"


/* load the node information */
void fixnodes(char *basename, struct MESH *m) {
    char file[256];
    FILE *fd;
    int i;

    printf("Fixing up node elevations\n");

    /* we could just fix the new nodes as the first "for" statement
     * would do, but that leads to funny results (I need to figure out
     * why.)  So, let's try fixing all of them */

    /* for ( i = origcount + 1; i <= nodecount; i++ ) { */
    for ( i = 1; i <= nodecount; i++ ) {
	/* printf("Current: %d %.2f %.2f %.2f\n", i, nodes[i][0],
	       nodes[i][1], nodes[i][2]); */

	nodes[i][2] = mesh_altitude(m, nodes[i][0], nodes[i][1]);

	/* printf("Fixed: %d %.2f %.2f %.2f\n", i, nodes[i][0],
	       nodes[i][1], nodes[i][2]); */
    }

    strcpy(file, basename);
    strcat(file, ".1.node");

    printf("Overwriting original node file: %s\n", file);

    fd = fopen(file, "w");

    fprintf(fd, "%d 2 1 0\n", nodecount);

    for ( i = 1; i <= nodecount; i++ ) {
	fprintf(fd, "%d %.2f %.2f %.2f 0\n", i, nodes[i][0],
	       nodes[i][1], nodes[i][2]);
    }

    fclose(fd);
}


/* $Log$
/* Revision 1.2  1997/12/02 13:12:07  curt
/* Updated to fix every node.
/*
 * Revision 1.1  1997/11/27 00:17:33  curt
 * Initial revision.
 *
 */