diff --git a/src/module_param.f90 b/src/module_param.f90
index 9ce0040be596773985784d8f91b1686a4d56407b..7c63620b77f3cf5b36a4c0bf4b00b63bed10661b 100644
--- a/src/module_param.f90
+++ b/src/module_param.f90
@@ -671,7 +671,7 @@ module ibm_param
   real(mytype) :: chord,thickness,omega, tconv2_sign, shear_velocity
   real(mytype) :: ce(30),sh(30),ori(40), lv(30), av(30)
   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, ztorq_only, nbody, inviscid_output
+  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, nbody, inviscid_output, div_visu_flag
 end module ibm_param
 !############################################################################
 !############################################################################
diff --git a/src/parameters.f90 b/src/parameters.f90
index 79b1920a83c1554a11e9cfe22555d2a816fdf8c0..7aac0dda6174a40e2cf2292ff208561e3e76d888 100644
--- a/src/parameters.f90
+++ b/src/parameters.f90
@@ -73,7 +73,7 @@ subroutine parameter(input_i3d)
       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, ztorq_only, nbody, &
-      inviscid_output
+      inviscid_output, div_visu_flag
   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, &
@@ -772,6 +772,7 @@ subroutine parameter_defaults()
   force_csv=0
   nbody=1
   inviscid_output=0
+  div_visu_flag=0
   ra(:) = 1.0
 
   !! Gravity field
diff --git a/src/visu.f90 b/src/visu.f90
index ed0b4ca920d3867bf7fea543366b215de3f75e40..1f46058e71bdd20dcd9471c04b255e698940d098 100644
--- a/src/visu.f90
+++ b/src/visu.f90
@@ -175,7 +175,7 @@ contains
     use decomp_2d_io, only : decomp_2d_start_io
 
     use param, only : nrhotime, ilmn, iscalar, ioutput, irestart
-    use ibm_param, only : inviscid_output
+    use ibm_param, only : inviscid_output,div_visu_flag
     use variables, only : sx, cifip6, cisip6, ciwip6, cifx6, cisx6, ciwx6
     use variables, only : sy, cifip6y, cisip6y, ciwip6y, cify6, cisy6, ciwy6
     use variables, only : sz, cifip6z, cisip6z, ciwip6z, cifz6, cisz6, ciwz6
@@ -267,26 +267,28 @@ contains
     ! Write pressure
       call write_field(ta1, ".", "pp", num, .true., flush=.true.)
     endif
-    !  ! Interpolate div_visu_var
-    ! !WORK Z-PENCILS
-    ! call interzpv(ppi3,div_visu_var(:,:,:,1),dip3,sz,cifip6z,cisip6z,ciwip6z,cifz6,cisz6,ciwz6,&
-    ! (ph3%zen(1)-ph3%zst(1)+1),(ph3%zen(2)-ph3%zst(2)+1),nzmsize,zsize(3),1)
-    ! !WORK Y-PENCILS
-    ! call transpose_z_to_y(ppi3,pp2,ph3) !nxm nym nz
-    ! call interypv(ppi2,pp2,dip2,sy,cifip6y,cisip6y,ciwip6y,cify6,cisy6,ciwy6,&
-    !         (ph3%yen(1)-ph3%yst(1)+1),nymsize,ysize(2),ysize(3),1)
-    ! !WORK X-PENCILS
-    ! call transpose_y_to_x(ppi2,pp1,ph2) !nxm ny nz
-    ! call interxpv(ta1,pp1,di1,sx,cifip6,cisip6,ciwip6,cifx6,cisx6,ciwx6,&
-    !         nxmsize,xsize(1),xsize(2),xsize(3),1)
-
-
-    ! ! Rescale div_visu_var
-    ! call rescale_pressure(ta1)
-
-    ! ! Write pressure
-    ! call write_field(ta1, ".", "div_visu", num, .true., flush=.true.)
 
+    if (div_visu_flag.eq.1) then 
+     ! Interpolate div_visu_var
+    !WORK Z-PENCILS
+    call interzpv(ppi3,div_visu_var(:,:,:,1),dip3,sz,cifip6z,cisip6z,ciwip6z,cifz6,cisz6,ciwz6,&
+    (ph3%zen(1)-ph3%zst(1)+1),(ph3%zen(2)-ph3%zst(2)+1),nzmsize,zsize(3),1)
+    !WORK Y-PENCILS
+    call transpose_z_to_y(ppi3,pp2,ph3) !nxm nym nz
+    call interypv(ppi2,pp2,dip2,sy,cifip6y,cisip6y,ciwip6y,cify6,cisy6,ciwy6,&
+            (ph3%yen(1)-ph3%yst(1)+1),nymsize,ysize(2),ysize(3),1)
+    !WORK X-PENCILS
+    call transpose_y_to_x(ppi2,pp1,ph2) !nxm ny nz
+    call interxpv(ta1,pp1,di1,sx,cifip6,cisip6,ciwip6,cifx6,cisx6,ciwx6,&
+            nxmsize,xsize(1),xsize(2),xsize(3),1)
+
+
+    ! Rescale div_visu_var
+    call rescale_pressure(ta1)
+
+    ! Write pressure
+    call write_field(ta1, ".", "div_visu", num, .true., flush=.true.)
+    endif
     ! LMN - write density
     if (ilmn) call write_field(rho1(:,:,:,1), ".", "rho", num)