|
12 | 12 | #include "AuxiliaryState.h" |
13 | 13 | #include "DataTypes.h" |
14 | 14 | #include "HorzMesh.h" |
| 15 | +#include "HorzOperators.h" |
15 | 16 | #include "OceanState.h" |
16 | 17 | #include "Tracers.h" |
17 | 18 | #include "VertCoord.h" |
18 | | -#include "HorzOperators.h" |
19 | 19 |
|
20 | 20 | namespace OMEGA { |
21 | 21 |
|
@@ -64,13 +64,15 @@ TracerHighOrderHorzAdvOnCell::TracerHighOrderHorzAdvOnCell(const HorzMesh *Mesh) |
64 | 64 | : HorzontalMesh(Mesh), |
65 | 65 | NAdvCellsForEdge("NumberOfCellsContribToAdvectionAtEdge", |
66 | 66 | Mesh->NEdgesOwned), |
67 | | - AdvCellsForEdge("IndexOfCellsContributingToAdvection", |
68 | | - Mesh->NEdgesOwned, Mesh->MaxEdges2 + 2), |
| 67 | + AdvCellsForEdge("IndexOfCellsContributingToAdvection", Mesh->NEdgesOwned, |
| 68 | + Mesh->MaxEdges2 + 2), |
69 | 69 | AdvMaskHighOrder("MaskForHighOrderAdvectionTerms", Mesh->NEdgesAll), |
70 | | - AdvCoefs("CommonAdvectionCoefficients", Mesh->MaxEdges2+2, Mesh->NEdgesAll), |
71 | | - AdvCoefs3rd("CommonAdvectionCoeffsForHighOrder", Mesh->MaxEdges2+2, |
| 70 | + AdvCoefs("CommonAdvectionCoefficients", Mesh->MaxEdges2 + 2, |
| 71 | + Mesh->NEdgesAll), |
| 72 | + AdvCoefs3rd("CommonAdvectionCoeffsForHighOrder", Mesh->MaxEdges2 + 2, |
72 | 73 | Mesh->NEdgesAll), |
73 | | - HighOrderFlxHorz("HigherOrderHorizontalFlux",Tracers::getNumTracers(), Mesh->NEdgesAll,Mesh->NVertLevels/VecLength), |
| 74 | + HighOrderFlxHorz("HigherOrderHorizontalFlux", Tracers::getNumTracers(), |
| 75 | + Mesh->NEdgesAll, Mesh->NVertLayers / VecLength), |
74 | 76 | NEdgesOnCell(Mesh->NEdgesOnCell), EdgesOnCell(Mesh->EdgesOnCell), |
75 | 77 | CellsOnEdge(Mesh->CellsOnEdge), EdgeSignOnCell(Mesh->EdgeSignOnCell), |
76 | 78 | DvEdge(Mesh->DvEdge), AreaCell(Mesh->AreaCell) { |
@@ -98,7 +100,7 @@ void TracerHighOrderHorzAdvOnCell::init() { |
98 | 100 | // Allocate Kokkos arrays in member data |
99 | 101 |
|
100 | 102 | SecondDerivativeOnCell secondDerivativeOnCell(Mesh); |
101 | | - Array3DReal DerivTwo("DerivTwo", MaxEdges2+2, 2, NEdgesAll); |
| 103 | + Array3DReal DerivTwo("DerivTwo", MaxEdges2 + 2, 2, NEdgesAll); |
102 | 104 | parallelFor( |
103 | 105 | {NCellsOwned}, |
104 | 106 | KOKKOS_LAMBDA(int ICell) { secondDerivativeOnCell(DerivTwo, ICell); }); |
|
0 commit comments