diff --git a/src/module_param.f90 b/src/module_param.f90
index 1342d389aacb65a502ea7e249650864aa5bfb28e..90386e2a12b049b79468101dacedca679e67157b 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, torq_flip
+  integer :: imove, nozdrift, force_csv, bodies_fixed, cube_flag, torques_flag,orientations_free, shear_flow_ybc, shear_flow_zbc,torq_debug, torq_flip, ztorq_only
 end module ibm_param
 !############################################################################
diff --git a/src/parameters.f90 b/src/parameters.f90
index a7818016fa52c59515ff5a90f604e798ce07f279..bdee324691f53953a21b2c2953eba8dc4f443b15 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, torq_flip
+      torques_flag, orientations_free, shear_flow_ybc, shear_flow_zbc, shear_velocity, torq_debug, torq_flip, ztorq_only
   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, &
diff --git a/src/xcompact3d.f90 b/src/xcompact3d.f90
index 8634f356f71cd3f768ee8cc63e101189b2b106ac..ffec99fb0b2cd7ae46938a1551297678ee374e51 100644
--- a/src/xcompact3d.f90
+++ b/src/xcompact3d.f90
@@ -149,6 +149,9 @@ program xcompact3d
         endif 
         if (orientations_free.eq.1) then 
          torque = [xtorq,ytorq,ztorq]
+         if (ztorq_only.eq.1) then
+            torque = [0.0_mytype, 0.0_mytype, ztorq]
+         endif
         else 
          torque(:) = zero
         endif