/home/eoseret/qaas/qaas_runs/178-231-1255/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/advec_cell_kernel.f90: 199 - 246
--------------------------------------------------------------------------------

199:       !$OMP DO PRIVATE(upwind,donor,downwind,dif,sigmat,sigma3,sigma4,sigmav,sigma,sigmam, &
200:       !$OMP            diffuw,diffdw,limiter,wind)
201:       DO k=y_min,y_max+2
202:         DO j=x_min,x_max
203: 
204:           IF(vol_flux_y(j,k).GT.0.0)THEN
[...]
210:             upwind   =MIN(k+1,y_max+2)
[...]
216:           sigmat=ABS(vol_flux_y(j,k))/pre_vol(j,donor)
217:           sigma3=(1.0_8+sigmat)*(vertexdy(k)/vertexdy(dif))
218:           sigma4=2.0_8-sigmat
219: 
220:           sigma=sigmat
221:           sigmav=sigmat
222: 
223:           diffuw=density1(j,donor)-density1(j,upwind)
224:           diffdw=density1(j,downwind)-density1(j,donor)
225:           wind=1.0_8
226:           IF(diffdw.LE.0.0) wind=-1.0_8
227:           IF(diffuw*diffdw.GT.0.0)THEN
228:             limiter=(1.0_8-sigmav)*wind*MIN(ABS(diffuw),ABS(diffdw)&
229:               ,one_by_six*(sigma3*ABS(diffuw)+sigma4*ABS(diffdw)))
230:           ELSE
231:             limiter=0.0
232:           ENDIF
233:           mass_flux_y(j,k)=vol_flux_y(j,k)*(density1(j,donor)+limiter)
234: 
235:           sigmam=ABS(mass_flux_y(j,k))/(density1(j,donor)*pre_vol(j,donor))
236:           diffuw=energy1(j,donor)-energy1(j,upwind)
237:           diffdw=energy1(j,downwind)-energy1(j,donor)
238:           wind=1.0_8
239:           IF(diffdw.LE.0.0) wind=-1.0_8
240:           IF(diffuw*diffdw.GT.0.0)THEN
241:             limiter=(1.0_8-sigmam)*wind*MIN(ABS(diffuw),ABS(diffdw)&
242:               ,one_by_six*(sigma3*ABS(diffuw)+sigma4*ABS(diffdw)))
243:           ELSE
244:             limiter=0.0
245:           ENDIF
246:           ener_flux(j,k)=mass_flux_y(j,k)*(energy1(j,donor)+limiter)
