@@ -367,8 +367,8 @@ Field3D Div_par_K_Grad_par(const Field3D& kY, const Field3D& f, CELL_LOC outloc)
367367 + Div_par (kY , outloc) * Grad_par (f, outloc);
368368}
369369
370- Field3D Div_par_K_Grad_par_mod (const Field3D& Kin, const Field3D& fin,
371- Field3D& flow_ylow, bool bndry_flux) {
370+ Field3D Div_par_K_Grad_par_mod (const Field3D& Kin, const Field3D& fin, Field3D& flow_ylow,
371+ bool bndry_flux) {
372372 TRACE (" FV::Div_par_K_Grad_par_mod" );
373373
374374 ASSERT2 (Kin.getLocation () == fin.getLocation ());
@@ -396,18 +396,22 @@ Field3D Div_par_K_Grad_par_mod(const Field3D& Kin, const Field3D& fin,
396396 const auto iym = i.ym ();
397397
398398 // Upper cell edge
399- const BoutReal c_up = 0.5 * (Kin[i] + K_up[iyp]); // K at the upper boundary
400- const BoutReal J_up = 0.5 * (coord->J [i] + coord->J .yup ()[iyp]); // Jacobian at boundary
399+ const BoutReal c_up = 0.5 * (Kin[i] + K_up[iyp]); // K at the upper boundary
400+ const BoutReal J_up =
401+ 0.5 * (coord->J [i] + coord->J .yup ()[iyp]); // Jacobian at boundary
401402 const BoutReal g_22_up = 0.5 * (coord->g_22 [i] + coord->g_22 .yup ()[iyp]);
402- const BoutReal gradient_up = 2 . * (f_up[iyp] - fin[i]) / (coord->dy [i] + coord->dy .yup ()[iyp]);
403+ const BoutReal gradient_up =
404+ 2 . * (f_up[iyp] - fin[i]) / (coord->dy [i] + coord->dy .yup ()[iyp]);
403405
404406 const BoutReal flux_up = c_up * J_up * gradient_up / g_22_up;
405407
406408 // Lower cell edge
407- const BoutReal c_down = 0.5 * (Kin[i] + K_down[iym]); // K at the lower boundary
408- const BoutReal J_down = 0.5 * (coord->J [i] + coord->J .ydown ()[iym]); // Jacobian at boundary
409+ const BoutReal c_down = 0.5 * (Kin[i] + K_down[iym]); // K at the lower boundary
410+ const BoutReal J_down =
411+ 0.5 * (coord->J [i] + coord->J .ydown ()[iym]); // Jacobian at boundary
409412 const BoutReal g_22_down = 0.5 * (coord->g_22 [i] + coord->g_22 .ydown ()[iym]);
410- const BoutReal gradient_down = 2 . * (fin[i] - f_down[iym]) / (coord->dy [i] + coord->dy .ydown ()[iym]);
413+ const BoutReal gradient_down =
414+ 2 . * (fin[i] - f_down[iym]) / (coord->dy [i] + coord->dy .ydown ()[iym]);
411415
412416 const BoutReal flux_down = c_down * J_down * gradient_down / g_22_down;
413417
@@ -468,7 +472,6 @@ Field3D Div_par_K_Grad_par_mod(const Field3D& Kin, const Field3D& fin,
468472 return result;
469473}
470474
471-
472475/* ******************************************************************************
473476* Delp2
474477* perpendicular Laplacian operator
0 commit comments