diff --git a/src/forces.f90 b/src/forces.f90
index 961f26074b02bcb2989c77785fc3c5643467d86b..44915cf6161cb50ded0bebc219edd4d9b96dc501 100644
--- a/src/forces.f90
+++ b/src/forces.f90
@@ -1803,9 +1803,9 @@ contains
       mom2 = sum(tunsty(:)) + sum(tconvy(:)) + (-2.0*tconv2_sign+1.0)*sum(tconvy2(:))
       mom3 = sum(tunstz(:)) + sum(tconvz(:)) + (-2.0*tconv2_sign+1.0)*sum(tconvz2(:))
  
-      dra1 = (sum(tdiffx) + sum(tdiffx2) + tp1 - mom1)
-      dra2 = (sum(tdiffy) + sum(tdiffy2) + tp2 - mom2)
-      dra3 = (sum(tdiffz) + sum(tdiffz2) + tp3 - mom3)
+      dra1 = -(sum(tdiffx) + sum(tdiffx2) + tp1 - mom1)*(1.0-2.0*torq_flip)
+      dra2 = -(sum(tdiffy) + sum(tdiffy2) + tp2 - mom2)*(1.0-2.0*torq_flip)
+      dra3 = -(sum(tdiffz) + sum(tdiffz2) + tp3 - mom3)*(1.0-2.0*torq_flip)
       
       ! do k=1,zsize(3)
 
diff --git a/src/module_param.f90 b/src/module_param.f90
index ff3167d7a310893c83ddf0d842ffe34fafbcd057..1342d389aacb65a502ea7e249650864aa5bfb28e 100644
--- a/src/module_param.f90
+++ b/src/module_param.f90
@@ -622,6 +622,6 @@ module ibm_param
   real(mytype) :: position_1(3),linearVelocity_1(3),orientation_1(4),angularVelocity_1(4)
   real(mytype) :: chord,thickness,omega, tconv2_sign, shear_velocity
   integer :: inana ! Analytical BC as Input
-  integer :: imove, nozdrift, force_csv, bodies_fixed, cube_flag, torques_flag,orientations_free, shear_flow_ybc, shear_flow_zbc,torq_debug
+  integer :: imove, nozdrift, force_csv, bodies_fixed, cube_flag, torques_flag,orientations_free, shear_flow_ybc, shear_flow_zbc,torq_debug, torq_flip
 end module ibm_param
 !############################################################################
diff --git a/src/parameters.f90 b/src/parameters.f90
index e5d1090f5d8229cc91c8327d2a519847a76df376..a7818016fa52c59515ff5a90f604e798ce07f279 100644
--- a/src/parameters.f90
+++ b/src/parameters.f90
@@ -63,7 +63,7 @@ subroutine parameter(input_i3d)
   NAMELIST /ibmstuff/ cex,cey,cez,shx,shy,shz,oriw,orii,orij,orik,lvx,lvy,lvz,avx,avy,avz,ra, &
       nobjmax,nraf,nvol,iforces, cvl_scalar, npif, izap, ianal, imove, thickness, chord, omega , &
       ubcx,ubcy,ubcz,rads,rho_s, c_air, grav_x,grav_y,grav_z, nozdrift, force_csv, bodies_fixed, cube_flag, tconv2_sign, &
-      torques_flag, orientations_free, shear_flow_ybc, shear_flow_zbc, shear_velocity, torq_debug
+      torques_flag, orientations_free, shear_flow_ybc, shear_flow_zbc, shear_velocity, torq_debug, torq_flip
   NAMELIST /ForceCVs/ xld, xrd, yld, yud, zld, zrd
   NAMELIST /LMN/ dens1, dens2, prandtl, ilmn_bound, ivarcoeff, ilmn_solve_temp, &
        massfrac, mol_weight, imultispecies, primary_species, &