(testparticles; #55) remove ifdefs

danieljprice · Jul 7, 2023 · 26cd79e · 26cd79e
1 parent 5077b45
commit 26cd79e
Showing 1 changed file with 8 additions and 12 deletions.
diff --git a/src/setup/setup_testparticles.F90 b/src/setup/setup_testparticles.F90
@@ -41,12 +41,8 @@ subroutine setpart(id,npart,npartoftype,xyzh,massoftype,vxyzu,polyk,gamma,hfact,
  use options,        only:iexternalforce,alpha,alphamax,alphau,beta,nfulldump
  use units,          only:set_units
  use physcon,        only:solarm
-#ifdef GR
- use externalforces, only:iext_gr
- use metric,         only:a
-#else
  use externalforces, only:iext_star
-#endif
+ use metric,         only:a
  use eos,            only:ieos
  use physcon,        only:pi
  use prompting,      only:prompt
@@ -70,7 +66,11 @@ subroutine setpart(id,npart,npartoftype,xyzh,massoftype,vxyzu,polyk,gamma,hfact,
  ! general parameters
  !
  time  = 0.
- gamma = 5./3. ! GR cannot have gamma=1
+ if (gr) then
+    gamma = 5./3. ! GR cannot have gamma=1
+ else
+    gamma = 1.
+ endif
  polyk = 0.
  npart = 10
  ieos  = 11
@@ -192,12 +192,8 @@ subroutine setpart(id,npart,npartoftype,xyzh,massoftype,vxyzu,polyk,gamma,hfact,
     vxyzu(1:3,i) = rtan*vcirc
  enddo
 
-#ifdef GR
- iexternalforce = iext_gr
- a              = spin
-#else
- iexternalforce = iext_star
-#endif
+ a = spin
+ if (.not.gr) iexternalforce = iext_star
 
 end subroutine setpart