| Message/Author |
|
|
|
|
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 |
|
|
|
|
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 |
|
|
|
|
Just tried to re-create that error, and that appears to be the problem, let me know if it still does that,
Tom |
|
|
|
|
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 |
|
|
|
|
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 |
|
|
|
|
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 |
|
|
|
|
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 |
|
|
|
|
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 |
|
|
|
|
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 |
|
|
|
|
<?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> |
|
|
|
|
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 |
|
|
|
|
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 |
