Skip to content
Snippets Groups Projects
Commit e3217c01 authored by Andrew Boyd's avatar Andrew Boyd
Browse files

Merge branch 'master' of https://git.ecdf.ed.ac.uk/s2006749/Incompact3d

parents 5fba44f8 fc46de70
No related branches found
No related tags found
1 merge request!1Adding ellip_utils
Hide whitespace changes
Inline Side-by-side
Showing
with 2 additions and 0 deletions
src/BC-Cylinder.f90
+ 2
− 0
View file @ e3217c01
...@@ -84,6 +84,8 @@ contains ...@@ -84,6 +84,8 @@ contains
point=[xm, ym, zm] point=[xm, ym, zm]
call EllipsoidalRadius(point, ce, orientation, shape, r) call EllipsoidalRadius(point, ce, orientation, shape, r)
! r=sqrt_prec((xm-cexx)**two+(ym-ceyy)**two+(zm-cezz)**two) ! r=sqrt_prec((xm-cexx)**two+(ym-ceyy)**two+(zm-cezz)**two)
! r=sqrt_prec((xm-cexx)**two+(ym-ceyy)**two)
if (r-ra.gt.zeromach) then if (r-ra.gt.zeromach) then
! write(*,*) i, j, k ! write(*,*) i, j, k
cycle cycle
......
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