diff --git a/improve_osolve.f90 b/improve_osolve.f90
index 0cf227c81953a74c4abb15e47c77e5904ba1edfc..ca598e68307a139dd58320f0be421bc32090955e 100644
--- a/improve_osolve.f90
+++ b/improve_osolve.f90
@@ -107,11 +107,16 @@ shift=' '
       
 call show_time (total,step,inc,1,'compute ref. criterion on ov$')
 
-allocate (x(mpe),y(mpe),z(mpe),stat=threadinfo%err)            ; call heap (threadinfo,'x/y/z','improve_osolve',size(x)+size(y)+size(z),'dp',+1)
-allocate (vx(mpe),vy(mpe),vz(mpe),stat=threadinfo%err)         ; call heap (threadinfo,'vx/y/z','improve_osolve',size(vx)+size(vy)+size(vz),'dp',+1)
-allocate (dhdx(mpe),dhdy(mpe),dhdz(mpe),stat=threadinfo%err)   ; call heap (threadinfo,'dhdx/y/z','improve_osolve',size(dhdx)+size(dhdy)+size(dhdz),'dp',+1)
-allocate (crit(ov%nleaves),stat=threadinfo%err)                ; call heap (threadinfo,'crit','improve_osolve',size(crit),'dp',+1)
-allocate (crittemp(ov%nleaves),stat=threadinfo%err)            ; call heap (threadinfo,'crittemp','improve_osolve',size(crittemp),'dp',+1)
+allocate (x(mpe),y(mpe),z(mpe),stat=threadinfo%err)
+call heap (threadinfo,'x/y/z','improve_osolve',size(x)+size(y)+size(z),'dp',+1)
+allocate (vx(mpe),vy(mpe),vz(mpe),stat=threadinfo%err)
+call heap (threadinfo,'vx/y/z','improve_osolve',size(vx)+size(vy)+size(vz),'dp',+1)
+allocate (dhdx(mpe),dhdy(mpe),dhdz(mpe),stat=threadinfo%err)
+call heap (threadinfo,'dhdx/y/z','improve_osolve',size(dhdx)+size(dhdy)+size(dhdz),'dp',+1)
+allocate (crit(ov%nleaves),stat=threadinfo%err)
+call heap (threadinfo,'crit','improve_osolve',size(crit),'dp',+1)
+allocate (crittemp(ov%nleaves),stat=threadinfo%err)
+call heap (threadinfo,'crittemp','improve_osolve',size(crittemp),'dp',+1)
 
 crit=0.d0
 crittemp=0.d0