diff --git a/compute_divergence.f90 b/compute_divergence.f90
index 4fe9affe593befc32b3382deca7ff073ede1c0b5..af55c12bebe5a8ba5c046e2e772ad31af3012bff 100644
--- a/compute_divergence.f90
+++ b/compute_divergence.f90
@@ -62,11 +62,16 @@ nint=8
 
 mpe=params%mpe
 
-allocate (divergence(osolve%nleaves),stat=err)           ; if (err.ne.0) call stop_run ('Error alloc divergence in compute_divergence$')
-allocate (divergencetemp(osolve%nleaves),stat=err)       ; if (err.ne.0) call stop_run ('Error alloc divergencetemp in compute_divergence$')
-allocate (x(mpe),y(mpe),z(mpe),stat=err)          ; if (err.ne.0) call stop_run ('Error alloc x,y,z in compute_divergence$')
-allocate (vx(mpe),vy(mpe),vz(mpe),stat=err)       ; if (err.ne.0) call stop_run ('Error alloc vx,vy,vz in compute_divergence$')
-allocate (dhdx(mpe),dhdy(mpe),dhdz(mpe),stat=err) ; if (err.ne.0) call stop_run ('Error alloc dhdx,dhdy,dhdz in compute_divergence$')
+allocate (divergence(osolve%nleaves),stat=err)
+if (err.ne.0) call stop_run ('Error alloc divergence in compute_divergence$')
+allocate (divergencetemp(osolve%nleaves),stat=err)
+if (err.ne.0) call stop_run ('Error alloc divergencetemp in compute_divergence$')
+allocate (x(mpe),y(mpe),z(mpe),stat=err)
+if (err.ne.0) call stop_run ('Error alloc x,y,z in compute_divergence$')
+allocate (vx(mpe),vy(mpe),vz(mpe),stat=err)
+if (err.ne.0) call stop_run ('Error alloc vx,vy,vz in compute_divergence$')
+allocate (dhdx(mpe),dhdy(mpe),dhdz(mpe),stat=err)
+if (err.ne.0) call stop_run ('Error alloc dhdx,dhdy,dhdz in compute_divergence$')
 
 
 divergence=0.d0