@@ -18,12 +18,6 @@ module m_chemistry
1818
1919 implicit none
2020
21- #:if USING_AMD
22- real (wp) :: molecular_weights_nonparameter(10 ) = (/ 2.016 , 1.008 , 15.999 , 31.998 , 17.007 , 18.015 , 33.006 , 34.014 , 39.95 , &
23- & 28.014 / )
24- $:GPU_DECLARE(create= ' [molecular_weights_nonparameter]' )
25- #:endif
26-
2721 type(int_bounds_info) :: isc1, isc2, isc3
2822 $:GPU_DECLARE(create= ' [isc1, isc2, isc3]' )
2923 integer , dimension (3 ) :: offsets
@@ -143,12 +137,7 @@ contains
143137
144138 $:GPU_LOOP(parallelism= ' [seq]' )
145139 do eqn = eqn_idx%species%beg, eqn_idx%species%end
146- #:if USING_AMD
147- omega_m = molecular_weights_nonparameter(eqn - eqn_idx%species%beg + 1 )* omega(eqn &
148- & - eqn_idx%species%beg + 1 )
149- #:else
150- omega_m = molecular_weights(eqn - eqn_idx%species%beg + 1 )* omega(eqn - eqn_idx%species%beg + 1 )
151- #:endif
140+ omega_m = molecular_weights(eqn - eqn_idx%species%beg + 1 )* omega(eqn - eqn_idx%species%beg + 1 )
152141 rhs_vf(eqn)%sf(x, y, z) = rhs_vf(eqn)%sf(x, y, z) + omega_m
153142 end do
154143 end do
@@ -259,17 +248,10 @@ contains
259248 ! Calculate species properties and gradients
260249 $:GPU_LOOP(parallelism= ' [seq]' )
261250 do i = eqn_idx%species%beg, eqn_idx%species%end
262- #:if USING_AMD
263- h_l(i - eqn_idx%species%beg + 1 ) = h_l(i - eqn_idx%species%beg + 1 ) &
264- & * gas_constant* T_L/ molecular_weights_nonparameter(i - eqn_idx%species%beg + 1 )
265- h_r(i - eqn_idx%species%beg + 1 ) = h_r(i - eqn_idx%species%beg + 1 ) &
266- & * gas_constant* T_R/ molecular_weights_nonparameter(i - eqn_idx%species%beg + 1 )
267- #:else
268- h_l(i - eqn_idx%species%beg + 1 ) = h_l(i - eqn_idx%species%beg + 1 ) &
269- & * gas_constant* T_L/ molecular_weights(i - eqn_idx%species%beg + 1 )
270- h_r(i - eqn_idx%species%beg + 1 ) = h_r(i - eqn_idx%species%beg + 1 ) &
271- & * gas_constant* T_R/ molecular_weights(i - eqn_idx%species%beg + 1 )
272- #:endif
251+ h_l(i - eqn_idx%species%beg + 1 ) = h_l(i - eqn_idx%species%beg + 1 ) &
252+ & * gas_constant* T_L/ molecular_weights(i - eqn_idx%species%beg + 1 )
253+ h_r(i - eqn_idx%species%beg + 1 ) = h_r(i - eqn_idx%species%beg + 1 ) &
254+ & * gas_constant* T_R/ molecular_weights(i - eqn_idx%species%beg + 1 )
273255 Xs_cell(i - eqn_idx%species%beg + 1 ) = 0.5_wp * (Xs_L(i - eqn_idx%species%beg + 1 ) + Xs_R(i &
274256 & - eqn_idx%species%beg + 1 ))
275257 h_k(i - eqn_idx%species%beg + 1 ) = 0.5_wp * (h_l(i - eqn_idx%species%beg + 1 ) + h_r(i &
@@ -294,15 +276,9 @@ contains
294276
295277 $:GPU_LOOP(parallelism= ' [seq]' )
296278 do eqn = eqn_idx%species%beg, eqn_idx%species%end
297- #:if USING_AMD
298- Mass_Diffu_Flux(eqn - eqn_idx%species%beg + 1 ) = rho_cell* mass_diffusivities_mixavg_Cell(eqn &
299- & - eqn_idx%species%beg + 1 )* molecular_weights_nonparameter(eqn &
300- & - eqn_idx%species%beg + 1 )/ MW_cell* dXk_dxi(eqn - eqn_idx%species%beg + 1 )
301- #:else
302- Mass_Diffu_Flux(eqn - eqn_idx%species%beg + 1 ) = rho_cell* mass_diffusivities_mixavg_Cell(eqn &
303- & - eqn_idx%species%beg + 1 )* molecular_weights(eqn - eqn_idx%species%beg + 1 ) &
304- & / MW_cell* dXk_dxi(eqn - eqn_idx%species%beg + 1 )
305- #:endif
279+ Mass_Diffu_Flux(eqn - eqn_idx%species%beg + 1 ) = rho_cell* mass_diffusivities_mixavg_Cell(eqn &
280+ & - eqn_idx%species%beg + 1 )* molecular_weights(eqn - eqn_idx%species%beg + 1 ) &
281+ & / MW_cell* dXk_dxi(eqn - eqn_idx%species%beg + 1 )
306282 rho_Vic = rho_Vic + Mass_Diffu_Flux(eqn - eqn_idx%species%beg + 1 )
307283 Mass_Diffu_Energy = Mass_Diffu_Energy + h_k(eqn - eqn_idx%species%beg + 1 )* Mass_Diffu_Flux(eqn &
308284 & - eqn_idx%species%beg + 1 )
0 commit comments