psc_bnd_fields_impl: fix scalar t

Author: JamesMcClung

src/libpsc/psc_bnd_fields/psc_bnd_fields_impl.hxx

@@ -557,8 +557,10 @@ struct BndFields_ : BndFieldsBase
       real_t s = 0.0;
       real_t p = 0.0;
       if (radiation) {
-        Real3 x3_s = (Real3(edge_idx) + Real3::unit(d1) * 0.5) * grid.domain.dx;
-        Real3 x3_p = (Real3(edge_idx) + Real3::unit(d2) * 0.5) * grid.domain.dx;
+        Real3 x3_s =
+          (Real3(edge_idx) + Real3::unit(d1) * real_t(0.5)) * grid.domain.dx;
+        Real3 x3_p =
+          (Real3(edge_idx) + Real3::unit(d2) * real_t(0.5)) * grid.domain.dx;
         s = radiation->pulse_s_lower(grid.time(), d0, p, x3_s);
         p = radiation->pulse_p_lower(grid.time(), d0, p, x3_p);
       }
@@ -606,8 +608,10 @@ struct BndFields_ : BndFieldsBase
       real_t s = 0.0;
       real_t p = 0.0;
       if (radiation) {
-        Real3 x3_s = (Real3(edge_idx) + Real3::unit(d1) * 0.5) * grid.domain.dx;
-        Real3 x3_p = (Real3(edge_idx) + Real3::unit(d2) * 0.5) * grid.domain.dx;
+        Real3 x3_s =
+          (Real3(edge_idx) + Real3::unit(d1) * real_t(0.5)) * grid.domain.dx;
+        Real3 x3_p =
+          (Real3(edge_idx) + Real3::unit(d2) * real_t(0.5)) * grid.domain.dx;
         s = radiation->pulse_s_upper(grid.time(), d0, p, x3_s);
         p = radiation->pulse_p_upper(grid.time(), d0, p, x3_p);
       }