//  -*- Mode: C++;  -*-
//  File: ravel.cc
//  Author: Dino Bellugi (dino@Moray.Berkeley.EDU)
//  Copyright (C) Dino Bellugi, 1996
// *~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
// * FUNCTION:	Main ravel program
// *		
// * CLASSES:	Ravel_Class
// * 
// * RELATED PACKAGES: Grid, Ravel_Class
// *
// * HISTORY:
// * Created:	Thu Oct 24 12:48:27 1996 (dino)
// *~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

#include <stdlib.h>
#include <fstream.h>

#include "grid.h"
#include "ravel_class.h"

int main (int argc, char *argv[]) 
{
#if (STREAMS)
	if (argc != 5)
#else 
	if (argc != 4)
#endif
	{
		// error - print usage
#if (STREAMS)
		cerr << "Usage: " << argv[0] << " <elevation grid> <streams grid> <critical slope (%)> <ravel constant>\n";
#else 
		cerr << "Usage: " << argv[0] << " <elevation grid> <critical slope (%)> <ravel constant>\n";
#endif
		exit(1);
	}

	Grid z;
#if (MESSAGES)
	cerr << "Reading elevation grid ...\n";
#endif
	z.readGrid(argv[1]);

#if (STREAMS)
#if (MESSAGES)
	cerr << "Reading channels grid ...\n";
#endif
	c.readGrid(argv[2]);

	// check grids
	if (z.xdim != c.xdim || z.ydim != c.ydim)
	{
#if (MESSAGES)
		cerr << "Grids have different dimensions ...\n";
#endif
		exit(1);
	}

	Ravel_Class totalRavel(z, c, atof(argv[3]), atof(argv[4]));
#else
	Ravel_Class totalRavel(z, atof(argv[2]), atof(argv[3]));
#endif

#if (MESSAGES)
	cerr << "Writing ravel grid ...\n";
#endif
	totalRavel.coutGrid();
}