ENH: Updated tutorial solver wrt updated cloud functionality

tutorials/lagrangian/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam/UEqn.H

@@ -4,8 +4,8 @@
       + fvm::div(phi, U)
       + turbulence->divDevRhoReff(U)
      ==
-        thermoCloud1.SU()
-      + kinematicCloud1.SU()
+        thermoCloud1.SU(U)
+      + kinematicCloud1.SU(U)
       + rho.dimensionedInternalField()*g
     );
 

tutorials/lagrangian/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam/hsEqn.H

@@ -6,7 +6,7 @@
       - fvm::laplacian(turbulence->alphaEff(), hs)
      ==
         DpDt
-      + thermoCloud1.Sh()
+      + thermoCloud1.Sh(hs)
     );
 
     hsEqn.relax();