ROAD DATA:
roadtype = Example_y,z-data_01
NX = number data points in x-direction (column)= 350
NY = number data points in y-direction (row)= 350
element read as (ix,iy) = (1,1), (1,2) ...
lattice constant (m) = 0.1000000E-03
scaling factor to give height in meter = 0.1000000E+01
computational time (minutes) = 0.1170E+00
FLAG PARAMETERS:
iflag1 (= 1 ---> no WINDOWING) = 3
(= 2 ---> WINDOWING with Bartlett)
(= 3 ---> WINDOWING with Welch)
iflag2 (= 1 ---> h - data (one single row, usually AFM data) = 2
(= 2 ---> y, h - data (two column))
(= 3 ---> x, y, h - data (three column))
iflag3 (= 1 ---> gives C(q)) = 1
(= 2 ---> gives CT(q) for TOP-profile)
(= 3 ---> gives CB(q) for BOTTOM-profile)
iflag4 (iflag4 = 1 (and iflag3 = 1) ---> wet surface) = 0
iflag5 (= 1 ---> gives G(r) = - ) = 0
(and G3(r) = - ^2)
iflag6 (= 1 ---> calculate power spectrum from 1D-data) = 0
(typically with NX=1, NY=N, N=number of point)
iswitch (= 1 ---> h(i,j)-->h(j,i); NX-->NY, NY-->NX) = 0
(useful if iflag6 = 1)
ido (= 1 or 2, fill out undefined data points hfill) = 2
ireduce (= 0 without anisotropy and topography files, = 1 with) = 1
itilt (= 0 remove nothing) = 4
(= 1 remove tilt of line(s) 1D; use only if iflag6=1)
(= 2 remove tilt of surface plane 2D; usually if iflag6=0)
(= 3 remove tilt and curvature of line(s) 1D; use only if iflag6=1)
(= 4 remove tilt and curvature of surface plane 2D; usually if iflag6=0)
Cqpoints (number of points for which C(q) is given) = 100
Phpoints (number of points for which P(h) and P(slope) is given) = 200
nreduce1 (factor reducing size of height data-set for plotting, use, e.,g., 8) = 2
PREPARE DATA: REMOVE UNDEFINED NUMBERS, REDUCE SLOPE:
ido = index telling what to do = 2
ido=0 --> do nothing
ido=1 --> fill undefined hfill
ido=2 --> fill undefined hfill, indicate large slope, fill undefined
undefined data point in input h-file must be denoted hfill =-0.1000000E+01
nit = number of iterations during filling= 1 1
radius (in units of lattice constant) for averaging = 0.1200000E+01
slopemax = maximum slope = 0.1000000E+02
fraction of total points which are undefined = 0.7906939E-01 0.7265306E-03
CALCULATED PARAMETERS 2D FFT:
radius of curvatures of fit surface (when itilt=4) (cm) = 0.1850799E+03 -0.1999966E+03
Ra = <|h-|> road roughness from measured data (m) = 0.2201408E-03
Rq/Ra (for Gaussian surface Rq = Ra*sqrt(pi/2)= 1.25..*Ra) = 0.1276153E+01
root mean square roughness Rq from measured data (m) = 0.2809335E-03
root mean square roughness, measured, after windowing (m) = 0.2557162E-03
root mean square road roughness from C(q) (m) = 0.2530915E-03
highest point (above average plane) in units of rms and in m = 0.3906125E+01 0.1097361E-02
lowest point (below average plane) in units of rms and in m =-0.4366060E+01 -0.1226572E-02
average slope <|grad h|> from measured data = 0.5789009E+00
maximum slope from measured data = 0.1361846E+02
rms of slope from measured data = 0.7613261E+00
fraction of surface area with slope above 10.0 = 0.8210113E-04
rms of slope from measured data removing slopes above 10.0 = 0.7528205E+00
rms of slope from C(q) = 0.7218490E+00
skewness /(rms)^3 measured data = 0.7838698E+00
for Gaussian surface, ensamble averaged skewness = 0
kurtosis (or sharpness) /(rms)^4 from measured data = 0.3350821E+01
for Gaussian surface, ensamble averaged kurtosis = 3
A_total/A_0 from measured data = 0.1206357E+01
A_total/A_0 from measured data removing slopes above 10.0 = 0.1205426E+01
A_total/A_0 from C(q) = 0.1218402E+01
TrippAverageLocal, Tripp and Peklenik gamma-parameter (related
to C(q)-anisotropy) = 0.1082001E+01 0.1543529E+01 0.1168605E+01
TrippAverageLocal, Tripp and Peklenik angle between gamma-ellipse and
x-axis (degree) =-0.1561626E+02 0.3867908E+02 0.1316322E+02