@@ -503,15 +503,15 @@ contains
503503 ! Schlieren field is dropped, leaving exp (0 ) = 1 everywhere. Compute that term below.
504504 alpha_last = 1._wp
505505 do i = 1 , eqn_idx%adv%end - eqn_idx%E
506- q_sf(j, k, l) = q_sf(j, k, l) - schlieren_alpha(i)* q_cons_vf(i + eqn_idx%E)%sf(j, k, l) * gm_rho_sf(j, &
507- & k, l)/ gm_rho_max(1 )
506+ q_sf(j, k, l) = q_sf(j, k, l) - schlieren_alpha(i)* q_cons_vf(i + eqn_idx%E)%sf(j, k, &
507+ & l) * gm_rho_sf(j, k, l)/ gm_rho_max(1 )
508508 alpha_last = alpha_last - q_cons_vf(i + eqn_idx%E)%sf(j, k, l)
509509 end do
510510 q_sf(j, k, l) = q_sf(j, k, l) - schlieren_alpha(num_fluids)* alpha_last* gm_rho_sf(j, k, l)/ gm_rho_max(1 )
511- else
511+ else
512512 do i = 1 , eqn_idx%adv%end - eqn_idx%E
513- q_sf(j, k, l) = q_sf(j, k, l) - schlieren_alpha(i)* q_cons_vf(i + eqn_idx%E)%sf(j, k, l) * gm_rho_sf(j, &
514- & k, l)/ gm_rho_max(1 )
513+ q_sf(j, k, l) = q_sf(j, k, l) - schlieren_alpha(i)* q_cons_vf(i + eqn_idx%E)%sf(j, k, &
514+ & l) * gm_rho_sf(j, k, l)/ gm_rho_max(1 )
515515 end do
516516 end if
517517 end do
0 commit comments