[Anuga-user] Extraction of runup from ANUGA sww file

[Muhammad I. <muh...@uq...>]

Hi,

I have been using ANUGA for quite a decent time on wave flume modelling. However, I am struggling to extract the runup time series (swash motion) of my simulation. The only way that I can think of  is to extract the data from the sww file and find the location of min depth (threshold depth, thr) at each time (Please refer to my matlab script below). This approach is convenient for regular waves but it is computationally expensive for random waves (t = 10mins). Note that depthc{1} = depth(t,x), i.e. 24000x750 array (0:0.025:600s and 0:0.01:7.5m).

Thr = 0.01 m corresponds to the elevation of my runup wires from the bed.

for i=1:nt;
            d=depthc{k}(i,:);
            d(d<0) = 0
            a=find(d<thr,1,'first');
            if isempty(a)==0; %%to omit the data with zero value.
                index(i)=a;
            else index(i)=nx;

            end
            d(index(i)+1:nx) = NaN
            dfilt(i,:) = d
        end
        xsh=xc(index);

I am not sure whether there is an output_export.py file in ANUGA example folder that may solve this problem faster than my matlab script. At the moment, it took me more than 24hours to get the output on my WORK PC. I hope you can give comments on this matter.



Regards,
Muhammad

Muhammad Shazril Idris Ibrahim
Doctor of Philosophy (Ph.D.) Student - School of Civil Engineering
Room 219 | Building 78
School of Civil Engineering | The University of Queensland | QLD 4072
Australia
M: +614 50220276 | E: Muh...@uq...