diff --git a/src/BC-Ellipsoid.f90 b/src/BC-Ellipsoid.f90
index c8d72d1026fc20747529c05662f7c7576e6fe7bc..0abe7501cfbffcbf3e5af079bd7730732429f328 100644
--- a/src/BC-Ellipsoid.f90
+++ b/src/BC-Ellipsoid.f90
@@ -84,7 +84,7 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
         if (ym /= ym) then
             write(*,*) "ym = ", ym, " should be ", ((real(j-1,mytype))*dy), ", as j = ", j,  " (or maybe it should be ", ((real(j-1,mytype))*dy)/real(nraf,mytype), ")"
             if (j.lt.(nyraf)) then 
-                write(*,*) "yp(j-1) = ", yp(j-1), " yp(j+1) = ", yp(j+1)
+                write(*,*) "yp(j-1) = ", yp(j-1), " yp(j+1,..) = ", yp(j+1:)
             endif
         endif
 
@@ -99,12 +99,12 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
                     !     write(*,*) point,position(i_body,:),orientation(i_body,:),shape(i_body,:)
                     ! endif
                     call EllipsoidalRadius(point,position(i_body,:),orientation(i_body,:),shape(i_body,:),r)
-                    if (r /= r) then 
-                        write(*,*) "Point, position, orientation, shape for body", i_body, " = "
-                        write(*,*) point,position(i_body,:),orientation(i_body,:),shape(i_body,:)
-                        write(*,*) "R calculated = ", r
-                        write(*,*) "Timestep = ", itime
-                    endif
+                    ! if (r /= r) then 
+                    !     write(*,*) "Point, position, orientation, shape for body", i_body, " = "
+                    !     write(*,*) point,position(i_body,:),orientation(i_body,:),shape(i_body,:)
+                    !     write(*,*) "R calculated = ", r
+                    !     write(*,*) "Timestep = ", itime
+                    ! endif
 
                     is_inside = (r-ra(i_body)).lt.zeromach
                     ! if (is_inside) then 
diff --git a/src/genepsi3d.f90 b/src/genepsi3d.f90
index 51b0092a38bb1feb90d18b908e7df369647c974f..6bd7f4b83cea3dbcbf43d29b9daf2b17ca336059 100644
--- a/src/genepsi3d.f90
+++ b/src/genepsi3d.f90
@@ -276,7 +276,7 @@ contains
     real(mytype),dimension(ysize(1),nyraf,ysize(3))    :: yepsi
     real(mytype),dimension(zsize(1),zsize(2),nzraf)    :: zepsi
     real(mytype),dimension(ny)                         :: yp
-    real(mytype),dimension(nyraf)                      :: ypraf
+    real(mytype),dimension(nyraf+20)                      :: ypraf
     real(mytype)                     :: dxraf,dyraf,dzraf
     integer                          :: i,j,k
     integer                          :: ii,jj,kk