Date: Sun, 15 Jul 2001 18:26:51 -0600
From: "Roger S. Miller"
X-Mailer: Mozilla 4.76 [en] (X11; U; Linux 2.0.27 i586)
X-Accept-Language: en
To: Markus Neteler
Subject: r.drain -- new version
Markus,
I worked on the original r.drain for a while and wasn't really able to
understand what it was doing, much less where the problem was. I
figured that I might be able to spend the rest of the weekend figuring
it out and in the end I would still have a program with unusual behavior
in flat spots. I decided instead to finish the version that I started
based on r.fill.dir.
My version of r.drain (which I called r.d2) is in the attached zip file.
r.d2 reproduces most of the behavior of r.drain. It does not use the
"null_value" option and it won't read input from a sites file. It uses
all of the flags defined for r.drain. My first draft of the model
didn't correctly handle the ew,ns and diagonal resolutions (as
r.fill.dir does not). I believe that I fixed that problem. There are
also some other minor behaviors that differ a little from r.drain.
r.d2 runs pretty quickly on a small, clean dem, but it can take quite a
while to run on a large cell map - particularly one with very many
areas where the flow directions are ambiguous. I didn't write anything
into the program to warn the user about the execution time.
I tested it on the fcell map that you sent me when this issue first came
up and on a lat-long cell map. It worked well in both cases. The
version of r.fill.dir that I started with was the version I had archived
here, which was not updated for all the changes needed to make
r.fill.dir run on all of our platforms. I built in the changes that I
could remember, but there may be a few more that come up when you try to
compile it on the Cray.
I hope this meets your needs.
Roger