DEM error in header file PreviousNext
CAESAR discussion board > Model Users >
Message/Author
Next message Pierluigi De Rosa (Pierluigi)  posted on Wednesday, August 30, 2006 - 05:42 pm
Dear all CAESAR Users,

I'm starting now to use CAESAR to model river processes and I have tried the sample files to begin to work with it.
It semms work fine, so I want to try CAESAR on my sample area and I created a DEM and rotate it in L to R direction.
below are the first 6 line of my DEM:

ncols 37
nrows 46
xllcorner 1780.835274
yllcorner 2128.812173
cellsize 2.850302
NODATA_value -9999
-9999 -9999 -9999 ....
-9999 178.1002 175.9218 173.9347 ....

but when I load the data CAESAR tell me that there is an eroor in header file, then I press OK and It will crash and give me this error:

"Indice oltre i limiti della matrice"
or in english:
"Index out matrix limit"

Where is my DEM wrong??

Thanks to all and my sincere compliments to Prof. Tom Coulthard for the dear job carried out with CAESAR

Pierluigi
Next message Tom Coulthard (Tom)  posted on Thursday, August 31, 2006 - 10:35 am
Hello Pierluigi,
Thanks for the post -
At first glance this may be because the cell size is not an integer.. when CAESAR loads up the data file, it assumes that the DEM spacing is in whole numbers..
Theres nothing wrong with having floating point cell spacing - just that I've never come across it before, and all the variables within the model that use this are integers!
Try changing the header so the cellsize is 3 and see if that works..
Its a quick fix, and if that does work, I would recommend re-sampling the DEM to 3m..
Hope this helps,
Tom
Next message Tom Coulthard (Tom)  posted on Thursday, August 31, 2006 - 10:38 am
Just tried to re-create that error, and that appears to be the problem, let me know if it still does that,
Tom
Next message Pierluigi De Rosa (Pierluigi)  posted on Friday, September 01, 2006 - 02:09 pm
Dear Tom,

I tried to resample a DEM to an integer value, but I have still the same problem. I tried also other ascii DEM. Here (http://www.gfosservices.it/dem_rotate.txt) there is a mine ascii dem I have rotated and clipped using DEMeditor.exe
When I load it into CAESAR seems work all fine but when I press "Load data" I have this error
"Index out matrix boundary"

Can you try to use my DEM? what's the problem in my DEM file?

Thanks

Pierluigi
Next message Tom Coulthard (Tom)  posted on Friday, September 01, 2006 - 04:12 pm
Hello Pierluigi,
On the Grid tab before you press the load data button, do you have the 'override header file' box checked? If so uncheck it,
Tom
Next message Tom Coulthard (Tom)  posted on Friday, September 01, 2006 - 04:18 pm
Hello again,
Just tried downloading your file and starting it in CAESAR and seemed to work fine...? Is it running OK for you now?
Tom
Next message Pierluigi De Rosa (Pierluigi)  posted on Monday, September 04, 2006 - 11:44 am
Dear Tom,
I still have problem...
when I press "Load data" I have an error and CAESAR will crash!!
The error is "Index out matrix limit" (you know what does it means??);
if I press Ignore CAESAR seems to be in stand by (It tells me "Initialising")

I use the default setting and Catchment mode.
These are my setting:
DEM data file:dem_rotate.txt
Grain data file:null
Bedrock data file: null
rainfall data file: hourly_rain_data.dat
tracer run? UNCHECKED
memory limit 1
init# of scans 20
override header file UNCHECKED

all other setting are leave by default

It seems you OK?

Maybe is some problem in my machine??

Thanks

Pierluigi
Next message Tom Coulthard (Tom)  posted on Tuesday, September 05, 2006 - 03:35 am
Hi Pierluigi,
It seems to work for me?? I'm in Singapore at the moment, so my replies will be out of sync, but I'll try and email you some files to see if that helps,
Tom
Next message Tom Coulthard (Tom)  posted on Tuesday, September 05, 2006 - 03:48 am
Hi Pierluigi,
Problem sorted..
Its because the time limit (max run duration) defaults to 1000 hours, and the rainfall data set is longer than that - so the program tries to load in too much data and goes over the extent of the hourly rainfall limit array (which is defined to 1000 hours). If you change the max run duration to greater than the length of the rainfall record - (make it 1000000) then it works fine.
I;ve copied a working config file int he post below... and emailed you a zip file of a working version. Looks fine -
Tom
Next message Tom Coulthard (Tom)  posted on Tuesday, September 05, 2006 - 03:49 am
<?xml>
- <parms>
- <general-parms>
<headeroverride>false</headeroverride>
<x-coordinate>593</x-coordinate>
<y-coordinate>358</y-coordinate>
<initscans>20</initscans>
<maxerodelimit>0.02</maxerodelimit>
<cellsize>5</cellsize>
<memorylimit>1</memorylimit>
<minq>0.01</minq>
<creeprate>0.0025</creeprate>
<lateralerosionrate>0.0000002</lateralerosionrate>
<maxiter>100000</maxiter>
<runstarttime>0</runstarttime>
<maxrunduration>1000000</maxrunduration>
<slopefailurethreshold>45</slopefailurethreshold>
<wssmoothingradius>2</wssmoothingradius>
<mvalue>0.005</mvalue>
<growgrasstime>5</growgrasstime>
<initialq>60</initialq>
<wssmoothing>false</wssmoothing>
<grass-sediment>false</grass-sediment>
<flowdistribution>2</flowdistribution>
<mintimestep>0</mintimestep>
<evaporation>0.0</evaporation>
<vegcritshear>180.0</vegcritshear>
<bedslope>true</bedslope>
<wsslope>false</wsslope>
<catchment_mode>true</catchment_mode>
<reach_mode>false</reach_mode>
<lat1>true</lat1>
<lat2>false</lat2>
<lat3>false</lat3>
<cross_stream_grad>2</cross_stream_grad>
- <saveoptions>
<option>elevations</option>
<checked>true</checked>
</saveoptions>
- <saveoptions>
<option>elevdiff</option>
<checked>true</checked>
</saveoptions>
- <saveoptions>
<option>grainsize</option>
<checked>true</checked>
</saveoptions>
- <saveoptions>
<option>total tracer g/s</option>
<checked>true</checked>
</saveoptions>
- <saveoptions>
<option>tracer layer 1</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>tracer layer 2</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>tracer layer 3</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>tracer layer 4</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>tracer layer 5</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>tracer layer 6</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>tracer layer 7</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>tracer layer 8</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>tracer layer 9</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>water depth</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>d50</option>
<checked>true</checked>
</saveoptions>
- <saveoptions>
<option>flow velocity</option>
<checked>false</checked>
</saveoptions>
- <saveoptions>
<option>soil saturation</option>
<checked>false</checked>
</saveoptions>
</general-parms>
- <grain-size>
<gs>0.0005</gs>
<gp>0.144</gp>
<ss>true</ss>
<fv>0.066</fv>
</grain-size>
- <grain-size>
<gs>0.001</gs>
<gp>0.022</gp>
<ss>true</ss>
<fv>0.109</fv>
</grain-size>
- <grain-size>
<gs>0.002</gs>
<gp>0.019</gp>
<ss>false</ss>
<fv>0.164</fv>
</grain-size>
- <grain-size>
<gs>0.004</gs>
<gp>0.029</gp>
<ss>false</ss>
<fv>0.237</fv>
</grain-size>
- <grain-size>
<gs>0.008</gs>
<gp>0.068</gp>
<ss>false</ss>
<fv>0.338</fv>
</grain-size>
- <grain-size>
<gs>0.016</gs>
<gp>0.146</gp>
<ss>false</ss>
<fv>0.479</fv>
</grain-size>
- <grain-size>
<gs>0.032</gs>
<gp>0.220</gp>
<ss>false</ss>
<fv>0.678</fv>
</grain-size>
- <grain-size>
<gs>0.064</gs>
<gp>0.231</gp>
<ss>false</ss>
<fv>0.959</fv>
</grain-size>
- <grain-size>
<gs>0.128</gs>
<gp>0.121</gp>
<ss>false</ss>
<fv>1.357</fv>
</grain-size>
- <file-parms>
<inputtimestep>1440</inputtimestep>
<saveinterval>1000</saveinterval>
<savetologfileinterval>60</savetologfileinterval>
<tracerrun>false</tracerrun>
<uniquefilecheck>false</uniquefilecheck>
- <filenames>
<desc>DEM Data File</desc>
<name>dem_rotate.txt</name>
</filenames>
- <filenames>
<desc>Grain Data File</desc>
<name>null</name>
</filenames>
- <filenames>
<desc>Bedrock Data File</desc>
<name>null</name>
</filenames>
- <filenames>
<desc>Rain Data File</desc>
<name>hourly_rain_data.dat</name>
</filenames>
- <filenames>
<desc>Tracer File</desc>
<name>null</name>
</filenames>
- <filenames>
<desc>Tracer Sed Vol File</desc>
<name>null</name>
</filenames>
- <filenames>
<desc>Tracer Grain Size Data File</desc>
<name>null</name>
</filenames>
- <sources>
<input>false</input>
<x>
<y>
<filename>
</sources>
- <sources>
<input>false</input>
<x>
<y>
<filename>
</sources>
- <sources>
<input>false</input>
<x>
<y>
<filename>
</sources>
- <sources>
<input>false</input>
<x>
<y>
<filename>
</sources>
</file-parms>
- <description>
<s>
</description>
- <outputfile-parms>
<generateavifile>false</generateavifile>
<avifile>out.avi</avifile>
<avifreq>1000</avifreq>
<generatetimeseriesfile>false</generatetimeseriesfile>
<timeseriesfile>catchment.dat</timeseriesfile>
<timeseriesfreq>60</timeseriesfreq>
<generateiterationsfile>false</generateiterationsfile>
<iterationsfile>iterout.dat</iterationsfile>
</outputfile-parms>
- <display>
<top>86</top>
<left>175</left>
<width>875</width>
<height>640</height>
</display>
- <lateral>
<oldlat>true</oldlat>
<newlat>false</newlat>
</lateral>
</parms>
Next message Tom Coulthard (Tom)  posted on Tuesday, September 05, 2006 - 04:39 am
PS - when I run your catchment, you get loads of landslides occuring - you might want to change the slope threshold (above which landslides happen) to 60 or 70 degrees for now (just to get it working..).
Good luck!
Tom
Next message Pierluigi De Rosa (Pierluigi)  posted on Tuesday, September 05, 2006 - 09:01 am
Dear Tom, Dear All CAESAR users,

I tried a greater value of "max run duration" parameter and all works fine!!

Thanks to the precious suggestions of Prof. Tom.

I hope this can help other CAESAR user.

Pierluigi
Back to top
Add Your Message Here
Post:
Bold text Italics Underline Create a hyperlink Insert a clipart image

Username: Posting Information:
This is a private posting area. Only registered users and moderators may post messages here.
Password:
Options: Enable HTML code in message
Automatically activate URLs in message
Action: