Skip to content
Snippets Groups Projects
Commit a5ce915e authored by s2006749's avatar s2006749
Browse files

more debugging outputs

parent 6de9114d
No related branches found
No related tags found
1 merge request!1Adding ellip_utils
Hide whitespace changes
Inline Side-by-side
Showing
with 5 additions and 0 deletions
src/BC-Ellipsoid.f90
+ 5
− 0
View file @ a5ce915e
...@@ -55,6 +55,7 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp) ...@@ -55,6 +55,7 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
! Intitialise epsi ! Intitialise epsi
epsi(:,:,:)=zero epsi(:,:,:)=zero
! Update center of moving ellipsoid ! Update center of moving ellipsoid
! if (t.ne.0.) then ! if (t.ne.0.) then
...@@ -79,6 +80,10 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp) ...@@ -79,6 +80,10 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
do j=nyi,nyf do j=nyi,nyf
! ym=(real(j-1,mytype))*dy ! ym=(real(j-1,mytype))*dy
ym=yp(j) ym=yp(j)
if (ym /= ym) then
write(*,*) "ym = ", ym, " should be ", ((real(j-1,mytype))*dy), ", as j = ", jj
endif
do i=nxi,nxf do i=nxi,nxf
xm=real(i-1,mytype)*dx xm=real(i-1,mytype)*dx
point=[xm, ym, zm] point=[xm, ym, zm]
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment