more removals

common/ModBasicMath.F90

@@ -191,7 +191,6 @@ subroutine QuadFit_1D(x, f, a0, a1, a2)
     lhs = 0.
     rhs = 0.
     do i = 1, size(x)
-      print *, "i = ", i
       xi = x(i)
       xi2 = xi*xi
       xi3 = xi*xi2
@@ -206,13 +205,7 @@ subroutine QuadFit_1D(x, f, a0, a1, a2)
 
       lhs(3, 3) = lhs(3, 3) + xi4
 
-      print *, "lhs"
-      do j = 1, 3
-        print *, lhs(j, :)
-      end do
-
       rhs = rhs + (/1._WP, xi, xi2/)*f(i)
-      print *, "rhs = ", rhs
     end do ! i
 
     lhs(2, 1) = lhs(1, 2)

examples/carotid_web/initcond.F90

@@ -12,7 +12,7 @@ program randomized_cell_gen
   integer, parameter :: ranseed = 112
 
   !initial condition setup parameters
-  real(WP), parameter :: hematocrit = 0.16
+  real(WP), parameter :: hematocrit = 0.2
   real(WP) :: tuber = 4
   real(WP), parameter :: tubelen = 30
 
@@ -34,14 +34,10 @@ program randomized_cell_gen
   nodeNum = nodeNum + 1
 
   !calculate number of cells for the defined hematocrit, assuming all blood cells are healthy RBCs for volume
-  !hematocrit = 4 * nrbc / (tube_radius^2 * tube_length)
-  print *, "hematocrit", hematocrit
+  !hematocrit = 4 * nrbc / (3 * tube_radius^2 * tube_length)
   nrbcMax = ((3*(tubelen*tuber**2*hematocrit))/4)
-  ! nrbcMax = 500
 
-  if (rootWorld) then
-    write (*, *) "Num RBCs in simulation is ", nrbcMax
-  end if
+  if (rootWorld) write (*, *) "Num RBCs in simulation is ", nrbcMax
 
   !set other initialization params
   vBkg(1:2) = 0.; vBkg(3) = 8.

examples/randomized_case/initcond.F90

@@ -34,7 +34,7 @@ program randomized_cell_gen
   nodeNum = nodeNum + 1
 
   !calculate number of cells for the defined hematocrit, assuming all blood cells are healthy RBCs for volume
-  !hematocrit = 4 * nrbc / (tube_radius^2 * tube_length)
+  !hematocrit = 4 * nrbc / (3 * tube_radius^2 * tube_length)
   nrbcMax = ((3*(tubelen*tuber**2*hematocrit))/4)
 
   if (rootWorld) write (*, *) "Num RBCs in simulation is ", nrbcMax