# collect samples for pdf estimation from ic1 case
# sampling region (xs, xe) hard-coded for ic1 case

[u, v, w, y]

# field_name (arg_name=arg (, pair)* ) = expression

c = 1.0 - y

ddxc (export=true, xs=176, xe=400) = ddx(c)

fsd (export=true, xs=176, xe=400) = sqrt (sqr(ddx(c)) + sqr(ddy(c)) + sqr(ddz(c)))

avg { c, fsd, ddxc }