@@ -266,9 +266,7 @@ contains
266266 rho_K = > rho_sf(i, j, k)
267267 gamma_K = > gamma_sf(i, j, k)
268268 pi_inf_K = > pi_inf_sf(i, j, k)
269- #endif
270-
271- #ifndef MFC_POST_PROCESS
269+ #else
272270 rho_K = > rho
273271 gamma_K = > gamma
274272 pi_inf_K = > pi_inf
@@ -323,9 +321,7 @@ contains
323321 rho_K = > rho_sf(j, k, l)
324322 gamma_K = > gamma_sf(j, k, l)
325323 pi_inf_K = > pi_inf_sf(j, k, l)
326- #endif
327-
328- #ifndef MFC_POST_PROCESS
324+ #else
329325 rho_K = > rho
330326 gamma_K = > gamma
331327 pi_inf_K = > pi_inf
@@ -420,9 +416,7 @@ contains
420416 rho_K = > rho_sf(k, l, r)
421417 gamma_K = > gamma_sf(k, l, r)
422418 pi_inf_K = > pi_inf_sf(k, l, r)
423- #endif
424-
425- #ifndef MFC_POST_PROCESS
419+ #else
426420 rho_K = > rho
427421 gamma_K = > gamma
428422 pi_inf_K = > pi_inf
@@ -606,7 +600,10 @@ contains
606600 integer :: i, j
607601!$acc update device(momxb, momxe, bubxb, bubxe, advxb, advxe, contxb, contxe, strxb, strxe)
608602
609- #ifndef MFC_PRE_PROCESS
603+ #ifdef MFC_PRE_PROCESS
604+ ixb = 0 ; iyb = 0 ; izb = 0 ;
605+ ixe = m; iye = n; ize = p;
606+ #else
610607 ixb = - buff_size
611608 ixe = m - ixb
612609
@@ -622,11 +619,6 @@ contains
622619 end if
623620#endif
624621
625- #ifdef MFC_PRE_PROCESS
626- ixb = 0 ; iyb = 0 ; izb = 0 ;
627- ixe = m; iye = n; ize = p;
628- #endif
629-
630622!$acc update device(ixb, ixe, iyb, iye, izb, ize)
631623
632624 allocate (gammas(1 :num_fluids))
@@ -813,17 +805,6 @@ contains
813805 end do
814806
815807 if (model_eqns /= 4 ) then
816- #ifndef MFC_SIMULATION
817- ! If pre- processing, use non acc mixture subroutines
818- if (hypoelasticity) then
819- call s_convert_to_mixture_variables(qK_cons_vf, j, k, l, &
820- rho_K, gamma_K, pi_inf_K, Re_K, G_K)
821- else
822- call s_convert_to_mixture_variables(qK_cons_vf, j, k, l, &
823- rho_K, gamma_K, pi_inf_K)
824- end if
825- #endif
826-
827808#ifdef MFC_SIMULATION
828809 ! If in simulation, use acc mixture subroutines
829810 if (hypoelasticity) then
@@ -836,6 +817,15 @@ contains
836817 call s_convert_species_to_mixture_variables_acc(rho_K, gamma_K, pi_inf_K, &
837818 alpha_K, alpha_rho_K, Re_K, j, k, l)
838819 end if
820+ #else
821+ ! If pre- processing, use non acc mixture subroutines
822+ if (hypoelasticity) then
823+ call s_convert_to_mixture_variables(qK_cons_vf, j, k, l, &
824+ rho_K, gamma_K, pi_inf_K, Re_K, G_K)
825+ else
826+ call s_convert_to_mixture_variables(qK_cons_vf, j, k, l, &
827+ rho_K, gamma_K, pi_inf_K)
828+ end if
839829#endif
840830 end if
841831
@@ -1032,9 +1022,8 @@ contains
10321022 end do
10331023 end do
10341024 end do
1035- #endif
10361025
1037- #ifdef MFC_SIMULATION
1026+ #else
10381027 if (proc_rank == 0 ) then
10391028 print ' (A)' , ' Conversion from primitive to ' // &
10401029 ' conservative variables not ' // &
0 commit comments