@@ -39,24 +39,24 @@ struct TestSetupPlane {
3939
4040 ErrorMeasures ExpectedNormalStressErrors = {0.0033910709836867704 ,
4141 0.0039954090464502795 };
42- KOKKOS_FUNCTION Real windStressX (Real X, Real Y) const {
42+ KOKKOS_FUNCTION Real sfcStressX (Real X, Real Y) const {
4343 return std::cos (TwoPi * X / Lx) * std::sin (TwoPi * Y / Ly);
4444 }
4545
46- KOKKOS_FUNCTION Real windStressY (Real X, Real Y) const {
46+ KOKKOS_FUNCTION Real sfcStressY (Real X, Real Y) const {
4747 return std::sin (TwoPi * X / Lx) * std::cos (TwoPi * Y / Ly);
4848 }
4949};
5050
5151struct TestSetupSphere {
5252 ErrorMeasures ExpectedNormalStressErrors = {0.0038588958862868362 ,
5353 0.003813760171030077 };
54- KOKKOS_FUNCTION Real windStressX (Real Lon, Real Lat) const {
54+ KOKKOS_FUNCTION Real sfcStressX (Real Lon, Real Lat) const {
5555 return -4 * std::sin (Lon) * std::cos (Lon) * std::pow (std::cos (Lat), 3 ) *
5656 std::sin (Lat);
5757 }
5858
59- KOKKOS_FUNCTION Real windStressY (Real Lon, Real Lat) const {
59+ KOKKOS_FUNCTION Real sfcStressY (Real Lon, Real Lat) const {
6060 return -std::pow (std::sin (Lon), 2 ) * std::pow (std::cos (Lat), 3 );
6161 }
6262};
@@ -83,8 +83,8 @@ int testSfcStressForcingVars(Real RTol) {
8383 Mesh->NEdgesOwned );
8484 Err += setVectorEdge (
8585 KOKKOS_LAMBDA (Real (&VecField)[2 ], Real X, Real Y) {
86- VecField[0 ] = Setup.windStressX (X, Y);
87- VecField[1 ] = Setup.windStressY (X, Y);
86+ VecField[0 ] = Setup.sfcStressX (X, Y);
87+ VecField[1 ] = Setup.sfcStressY (X, Y);
8888 },
8989 ExactNormalStressEdge, EdgeComponent::Normal, Geom, Mesh,
9090 ExchangeHalos::No);
@@ -94,11 +94,11 @@ int testSfcStressForcingVars(Real RTol) {
9494
9595 // Set inputs
9696 Err += setScalar (
97- KOKKOS_LAMBDA (Real X, Real Y) { return Setup.windStressX (X, Y); },
97+ KOKKOS_LAMBDA (Real X, Real Y) { return Setup.sfcStressX (X, Y); },
9898 SfcStressForcing.ZonalStressCell , Geom, Mesh, OnCell);
9999
100100 Err += setScalar (
101- KOKKOS_LAMBDA (Real X, Real Y) { return Setup.windStressY (X, Y); },
101+ KOKKOS_LAMBDA (Real X, Real Y) { return Setup.sfcStressY (X, Y); },
102102 SfcStressForcing.MeridStressCell , Geom, Mesh, OnCell);
103103
104104 // Compute numerical result
0 commit comments