diff --git a/code/termspecs/termspec.vnsd_fluctuation b/code/termspecs/termspec.vnsd_fluctuation index d3e3fca..36598af 100644 --- a/code/termspecs/termspec.vnsd_fluctuation +++ b/code/termspecs/termspec.vnsd_fluctuation @@ -12,7 +12,10 @@ nx = - ddx(c) / fsd ny = - ddy(c) / fsd nz = - ddz(c) / fsd -vn = u * nx + v * ny + w * nz +unx = u * nx +vny = v * ny +wnz = w * nz +vn = unx + vny + wnz divn = ddx(nx) + ddy(ny) + ddz(nz) @@ -35,7 +38,7 @@ gradt_vnsd1 = ddx(vnsd) - nx * gradn_vnsd divn2_dndn = divn2 - gradn_gradn -sign_divn2_dndn = $dsign(divn2_dndn, divn) +sign_divn2_dndn = $dsign(1.0, divn) * divn2_dndn vnsd_divn = vnsd * divn @@ -55,7 +58,10 @@ vnsd_sign_divn2_dndn = vnsd * sign_divn2_dndn avg { c, u, v, w, fsd, lapc, ddx_c, ddy_c, ddz_c, d2dx_c, d2dy_c, d2dz_c } -avg fsd { u, v, w, nx, ny, nz, absk, divn, sd, gradt_vnsd1, +avg fsd { u, v, w, nx, ny, nz, absk, divn, +unx, vny, wnz, +vn, sd, +gradt_vnsd1, vnsd, vnsd_divn, vnsd_divn_cor,