From 1e3b7bc1ca27af1f8a4afece0086292ed90f2408 Mon Sep 17 00:00:00 2001 From: J Todd Date: Wed, 13 Jun 2018 16:59:00 +0100 Subject: [PATCH] Removed hardcoded undercut shape from glas.f90 --- src/glas.f90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/glas.f90 b/src/glas.f90 index a855203..b664ccb 100644 --- a/src/glas.f90 +++ b/src/glas.f90 @@ -217,8 +217,8 @@ SUBROUTINE Initializefcc(box,l,xo,ip,myid,maxx,maxy,maxz,minx,miny,minz,ntasks,S lc=4420.0+1.5e-04*(x-3300.0)**2+0.42*exp((x-3700.0)/2.0e+02) UCV=lc-1500.0*exp(-(x-3500.0)**2/50000.0) ! If (y.lt.lc-UC.or.y.gt.lc.or.z.gt.bint+3.0*SQRT(y-(lc-UC)).or.z.ge.WL-20.0) then - If (y.lt.lc-UC.or.z.gt.bint+3.0*SQRT(y-(lc-UC)).or.z.ge.WL-40.0) then - If (y.lt.UCV.or.(z.gt.bint+3.0*sqrt(y-UCV)).or.z.ge.WL-40.0) then + ! If (y.lt.lc-UC.or.z.gt.bint+3.0*SQRT(y-(lc-UC)).or.z.ge.WL-40.0) then + ! If (y.lt.UCV.or.(z.gt.bint+3.0*sqrt(y-UCV)).or.z.ge.WL-40.0) then ip=ip+1 xo(1,ip) = x0(1,k1) + Float(i-1)*b @@ -230,8 +230,8 @@ SUBROUTINE Initializefcc(box,l,xo,ip,myid,maxx,maxy,maxz,minx,miny,minz,ntasks,S if (xo(1,ip).lt.minx) minx=xo(1,ip) if (xo(2,ip).lt.miny) miny=xo(2,ip) if (xo(3,ip).lt.minz) minz=xo(3,ip) - EndIF - EndIF + ! EndIF + ! EndIF EndIF ! EndIF EndDo