BUGFIX: ind_r should be treated the same for 1,3 and 4

75a4dedb · Sabine · d56bfae2 · 75a4dedb
Commit 75a4dedb authored Feb 09, 2017 by Sabine
--- a/src/releaseparticles.f90
+++ b/src/releaseparticles.f90
@@ -335,7 +335,7 @@ subroutine releaseparticles(itime)

  !Af ind_rel is defined in readcommand.f

-            if (ind_rel .eq. 1) then
+            if ((ind_rel .eq. 1).or.(ind_rel .eq. 3).or.(ind_rel .eq. 4)) then

  ! Interpolate the air density
  !****************************