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