void write_fd(register FILE *f, int r0, int c0, int rn, int cn) {
register struct ent **pp;
int r, c;
(void) fprintf(f, "# This data file was generated by SCIM.\n");
(void) fprintf(f, "# You almost certainly shouldn't edit it.\n\n");
print_options(f);
for (c = 0; c < COLFORMATS; c++)
if (colformat[c])
(void) fprintf (f, "format %d = \"%s\"\n", c, colformat[c]);
for (c = c0; c <= cn; c++)
if (fwidth[c] != DEFWIDTH || precision[c] != DEFPREC || realfmt[c] != DEFREFMT)
(void) fprintf (f, "format %s %d %d %d\n", coltoa(c), fwidth[c], precision[c], realfmt[c]);
for (c = c0; c <= cn; c++)
if (col_hidden[c])
(void) fprintf(f, "hide %s\n", coltoa(c));
for (r = r0; r <= rn; r++)
if (row_hidden[r])
(void) fprintf(f, "hide %d\n", r);
//write_ranges(f);
write_marks(f);
if (mdir)
(void) fprintf(f, "mdir \"%s\"\n", mdir);
if (autorun)
(void) fprintf(f, "autorun \"%s\"\n", autorun);
for (c = 0; c < FKEYS; c++)
if (fkey[c]) (void) fprintf(f, "fkey %d = \"%s\"\n", c + 1, fkey[c]);
write_cells(f, r0, c0, rn, cn, r0, c0);
// write blocked cells
for (r = r0; r <= rn; r++) {
pp = ATBL(tbl, r, c0);
for (c = c0; c <= cn; c++, pp++)
if (*pp) {
if ((*pp)->flags & is_locked)
(void) fprintf(f, "lock %s%d\n", coltoa((*pp)->col), (*pp)->row);
/*if ((*pp)->nrow >= 0) {
(void) fprintf(f, "addnote %s ", v_name((*pp)->row, (*pp)->col));
(void) fprintf(f, "%s\n", r_name((*pp)->nrow, (*pp)->ncol, (*pp)->nlastrow, (*pp)->nlastcol));
}*/
}
}
/*
* Don't try to combine these into a single fprintf(). v_name() has
* a single buffer that is overwritten on each call, so the first part
* needs to be written to the file before making the second call.
*/
fprintf(f, "goto %s", v_name(currow, curcol));
//fprintf(f, " %s\n", v_name(strow, stcol));
fprintf(f, "\n");
}
void write_file0(const char* fn) {
FILE* fd = safe_fopen(fn, "w");
write_file_version(fd);
write_header(fd);
write_format(fd);
write_vertices(fd);
write_cells(fd);
write_cells_attributes(fd);
fclose(fd);
}
int main(int argc, char **argv)
{
struct GModule *module;
struct
{
struct Option *input;
struct Option *output;
} params;
G_gisinit(argv[0]);
module = G_define_module();
G_add_keyword(_("raster"));
G_add_keyword(_("hydrology"));
module->description =
_("Creates a topographic index raster map from an elevation raster map.");
params.input = G_define_standard_option(G_OPT_R_ELEV);
params.input->key = "input";
params.output = G_define_standard_option(G_OPT_R_OUTPUT);
params.output->description = _("Name for output topographic index raster map");
if (G_parser(argc, argv))
exit(EXIT_FAILURE);
/* Make sure that the current projection is not lat/long */
if ((G_projection() == PROJECTION_LL))
G_fatal_error(_("Lat/Long location is not supported by %s. Please reproject map first."),
G_program_name());
input = params.input->answer;
output = params.output->answer;
G_get_window(&window);
read_cells();
initialize();
calculate_atanb();
write_cells();
exit(EXIT_SUCCESS);
}
