diff --git a/raft/raft_model.py b/raft/raft_model.py
index 97dd07e..a66878f 100644
--- a/raft/raft_model.py
+++ b/raft/raft_model.py
@@ -847,7 +847,7 @@ def step_func_equil(X, args, Y, oths, Ytarget, err, tol_, iter, maxIter):
         
         for i, fowt in enumerate(self.fowtList):
             print(f"Found mean offets of FOWT {i+1} with surge = {fowt.Xi0[0]: .2f} m,  sway  = {fowt.Xi0[1]: .2f},  and heave = {fowt.Xi0[2]: .2f} m")
-            print(f"                                 roll  = {fowt.Xi0[3]*180/np.pi: .2f} deg, pitch = {fowt.Xi0[4]*180/np.pi: .2f} deg, and yaw   = {fowt.Xi0[5]*180/np.pi: .2f} deg")
+            print(f"                                 roll  = {fowt.Xi0[3]*180/np.pi: .2f} deg, pitch = {fowt.Xi0[4]*180/np.pi: .2f}, and yaw   = {fowt.Xi0[5]*180/np.pi: .2f} deg")
 
         
         #dsolvePlot(info) # plot solver convergence trajectories
diff --git a/raft/raft_rotor.py b/raft/raft_rotor.py
index 8a97238..56743c1 100644
--- a/raft/raft_rotor.py
+++ b/raft/raft_rotor.py
@@ -636,7 +636,7 @@ def calcHydroConstants(self, dgamma=0, rho=1025, g=9.81):
         return A_hydro, I_hydro
 
 
-    def calcCavitation(self, case, azimuth=0, clearance_margin=1.0, Patm=101325, Pvap=2300, error_on_cavitation=False):
+    def calcCavitation(self, case, azimuth=0, clearance_margin=1.0, Patm=101325, Pvap=2500, error_on_cavitation=False):
         ''' Method to calculate the cavitation number of the rotor
         (wind speed (m/s), rotor speed (RPM), pitch angle (deg), azimuth (deg))
 
@@ -838,7 +838,7 @@ def calcAero(self, case, current=False, display=0):
         # Set up vectors in axis frame. Assuming CCBlade forces (but not 
         # moments) are relative to the rotor axis
         forces_axis = np.array([loads["T"][0], loads["Y"][0], loads["Z" ][0]])
-        moments_axis = np.array([loads["Q"][0], loads["My"][0], loads["Mz"][0]]) 
+        moments_axis = np.array([loads["My"][0], loads["Q"][0], loads["Mz"][0]])        
         
         # Rotate forces and moments to be relative to global orientation (but still wrt hub)
         self.f0[:3] = np.matmul(self.R_q, forces_axis)