Fix minor precision mismatch

Author: Mark Zhang

src/common/include/2dHardcodedIC.fpp

@@ -135,8 +135,8 @@
 
     case (207) ! Kelvin Helmholtz Instability
         sigma = 0.05_wp / sqrt(2.0_wp)
-        gauss1 = exp(- (y_cc(j) - 0.75_wp) ** 2 / (2 * sigma ** 2))
-        gauss2 = exp(- (y_cc(j) - 0.25_wp) ** 2 / (2 * sigma ** 2))
+        gauss1 = exp(- (y_cc(j) - 0.75_wp) ** 2 / (2.0_wp * sigma ** 2))
+        gauss2 = exp(- (y_cc(j) - 0.25_wp) ** 2 / (2.0_wp * sigma ** 2))
         q_prim_vf(momxb + 1)%sf(i, j, 0) = &
                 0.1_wp * sin(4.0_wp * pi * x_cc(i)) * (gauss1 + gauss2)
 
@@ -145,13 +145,13 @@
         eps = 1.0e-6_wp
         ei = 5.0_wp
         ! Smoothening function to smooth out sharp discontinuity in the interface
-        if(x_cc(i) <= 0.7 * lam) then
+        if(x_cc(i) <= 0.7_wp * lam) then
                 d = x_cc(i) - lam * (0.4_wp - 0.1_wp * sin(2.0_wp * pi * (y_cc(j) / lam + 0.25_wp)))
                 fsm = 0.5_wp * (1.0_wp + erf(d / (ei * sqrt(dx * dy))))
                 alpha_air = eps + (1.0_wp - 2.0_wp * eps) * fsm
-                alpha_sf6 = 1.0 - alpha_air
-                q_prim_vf(contxb)%sf(i, j, 0) = alpha_sf6 * 5.04
-                q_prim_vf(contxe)%sf(i, j, 0) = alpha_air * 1.0
+                alpha_sf6 = 1.0_wp - alpha_air
+                q_prim_vf(contxb)%sf(i, j, 0) = alpha_sf6 * 5.04_wp
+                q_prim_vf(contxe)%sf(i, j, 0) = alpha_air * 1.0_wp
                 q_prim_vf(advxb)%sf(i, j, 0) = alpha_sf6
                 q_prim_vf(advxe)%sf(i, j, 0) = alpha_air 
         end if