/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/aocc_2/CMakeFiles/clover_leaf.dir/CloverLeaf_ref/kernels/advec_cell_kernel.f90-pp.f90: 202 - 246
--------------------------------------------------------------------------------

202:       DO k=y_min,y_max+2
203:         DO j=x_min,x_max
204: # 204 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/advec_cell_kernel.f90"
[...]
210:           ELSE
[...]
216: # 216 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/advec_cell_kernel.f90"
217:           sigmat=ABS(vol_flux_y(j,k))/pre_vol(j,donor)
218:           sigma3=(1.0_8+sigmat)*(vertexdy(k)/vertexdy(dif))
219:           sigma4=2.0_8-sigmat
220: # 220 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/advec_cell_kernel.f90"
221:           sigma=sigmat
222:           sigmav=sigmat
223: # 223 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/advec_cell_kernel.f90"
224:           diffuw=density1(j,donor)-density1(j,upwind)
225:           diffdw=density1(j,downwind)-density1(j,donor)
226:           wind=1.0_8
227:           IF(diffdw.LE.0.0) wind=-1.0_8
228:           IF(diffuw*diffdw.GT.0.0)THEN
229:             limiter=(1.0_8-sigmav)*wind*MIN(ABS(diffuw),ABS(diffdw)&
230:               ,one_by_six*(sigma3*ABS(diffuw)+sigma4*ABS(diffdw)))
231:           ELSE
232:             limiter=0.0
233:           ENDIF
234:           mass_flux_y(j,k)=vol_flux_y(j,k)*(density1(j,donor)+limiter)
235: # 235 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/advec_cell_kernel.f90"
236:           sigmam=ABS(mass_flux_y(j,k))/(density1(j,donor)*pre_vol(j,donor))
237:           diffuw=energy1(j,donor)-energy1(j,upwind)
238:           diffdw=energy1(j,downwind)-energy1(j,donor)
239:           wind=1.0_8
240:           IF(diffdw.LE.0.0) wind=-1.0_8
241:           IF(diffuw*diffdw.GT.0.0)THEN
242:             limiter=(1.0_8-sigmam)*wind*MIN(ABS(diffuw),ABS(diffdw)&
243:               ,one_by_six*(sigma3*ABS(diffuw)+sigma4*ABS(diffdw)))
244:           ELSE
245:             limiter=0.0
246:           ENDIF
