diff --git a/src/navier-stokes.c b/src/navier-stokes.c
index e071e39..453e2c2 100644
--- a/src/navier-stokes.c
+++ b/src/navier-stokes.c
@@ -49,7 +49,7 @@ int uk(
   }
 
   // iterate
-  for(t=starting_time;t<starting_time+nsteps;t++){
+  for(t=starting_time;nsteps==0 || t<starting_time+nsteps;t++){
     ns_step(u, K1, K2, N1, N2, nu, delta, L, g, fft1, fft2, ifft, tmp1, tmp2, tmp3, irreversible);
     
     if(t%print_freq==0){
@@ -125,7 +125,7 @@ int eea(
   uint64_t first_box = print_freq - (starting_time % print_freq);
 
   // iterate
-  for(t=starting_time;t<starting_time+nsteps;t++){
+  for(t=starting_time;nsteps==0 || t<starting_time+nsteps;t++){
     ns_step(u, K1, K2, N1, N2, nu, delta, L, g, fft1, fft2, ifft, tmp1, tmp2, tmp3, irreversible);
 
     // convolution term
@@ -230,7 +230,7 @@ int quiet(
   copy_u(u, u0, K1, K2);
 
   // iterate
-  for(t=0;t<nsteps;t++){
+  for(t=0;nsteps==0 || t<nsteps;t++){
     ns_step(u, K1, K2, N1, N2, nu, delta, L, g, fft1, fft2, ifft, tmp1, tmp2, tmp3, irreversible);
   }