[GRASSLIST:648] Re: running r.topmodel

[Thomas Adams]

Wouter,

That's excellent; the choice of going to R is pretty interesting. Does 
this mean the Glue algorithms will be ported there as well?

Regards,
Tom


Wouter Buytaert wrote:
>
> Hi Jaime and the rest of the list,
>
> Just to say that I'm now one of the topmodel developers at Lancaster. 
> We are doing a major cleanup of the code and intend to release some 
> new versions soon, although in R rather than in GRASS (for several 
> reasons...) incorporating Monte Carlo and other stuff.
>
> BTW, there is a bug in the r.topmodel infiltration module which Huidae 
> Cho fixed some weeks ago but probably didn't make it to the official 
> version yet.
>
> Anyway, if you have more questions about topmodel, feel free to ask 
> via the list or directly to me.
>
> cheers
> wouter
>
> On Thu, 30 Mar 2006, Thomas Adams wrote:
>
>> Jaime,
>>
>> I believe I can help you, although:
>>
>> (1) I have never run TOPMODEL in GRASS (just the Fortran standalone 
>> version)
>>
>> (2) the link to Lancaster Univ. is old, in fact I'm not sure Kieth Bevan
>> Lancaster Univ. even make the original source available now, from 
>> what I can
>> see at: http://www.es.lancs.ac.uk/hfdg/freeware/hfdg_freeware_top.htm
>>
>> My guesses are these, since they conform to the GRASS r.topmodel 
>> manual page
>> and the original files for the standalone version of TOPMODEL:
>>
>> parameters=3Dstring
>> (i) TOPMODEL Parameters file, where the standalone file looked like 
>> this:
>>
>> SUBCATCHMENT 1
>> 0.032 5.0 50. 3600.0 3600.0 0.05 0.0000328 0.002 0 1.0 0.02 0.1
>>
>> I can not remember what's what; I'll have to dig this up, plus I'll 
>> have to
>> look at the r.topmodel source code more =97 I believe the data are 
>> read in this
>> order from what I can see in the source=85
>>
>> input=3Dstring
>> (i) Rainfall and potential evapotranspiration data file {probably a 
>> space
>> delimited file with 2 columns}
>> output=3Dstring
>> (o) Output file {looks like flow is cu meters/timestep}
>> Qobs=3Dstring
>> (i) OPTIONAL Observed flow file {single column file}
>>
>> My recollection is that the inputs are in terms of meters/time-step 
>> (flow
>> would be cu. meters/delta-t) =97 but I MUST CONFIRM THIS!
>>
>> The original input file in the standalone version combined the GRASS 
>> input
>> Qobs files and looked like this:
>>
>> 950 1.0 <=3D=3D=3D=3D=3D=3D Number of timesteps & Number of hours per 
>> time step (I
>> believe)
>> .0000000 .0000619 .0000328
>> .0000000 .0000647 .0000325
>> .0010000 .0000619 .0000330
>> .0010000 .0000560 .0000333
>> .0010000 .0000457 .0000348
>> .0005000 .0000323 .0000356
>> .0005000 .0000156 .0000356
>> .0000000 .0000000 .0000356
>> .0005000 .0000000 .0000351
>>
>> Notice that the values are very small, which makes me believe that I 
>> think my
>> guess about the units of the input are correct; this could be 
>> different in
>> the GRASS r.topmodel code, of course. The columns are, respectively,
>> precipitation, PE, and ObsQ.
>>
>> At least this can get you going. I'll try to get more on this in the 
> next few
>> days. With a little help of some of the more experienced GRASS module=20
>> contributors, I don't mind taking this on.
>>
>> Regards,
>> Tom
>>
>> Jaime Carrera wrote:
>>>  Hi,
>>>
>>>  Has any one used the r.topmodel module? If so, I have the following
>>>  questions:
>>>  1) What are the instructions that are required in the "topmodel 
>>> paramet=
> ers
>>>  file"?
>>>  2) How should the topographic index statistics file be written?
>>>
>>>  Cheers,
>>>
>>>  Jaime
>
>


-- 
Thomas E Adams
National Weather Service
Ohio River Forecast Center
1901 South State Route 134
Wilmington, OH 45177

EMAIL:	thomas.adams at noaa.gov

VOICE:	937-383-0528
FAX:	937-383-0033