/home/eoseret/qaas_runs_CPU_9468/171-152-3172/intel/CloverLeafFC/build/CloverLeafFC/CloverLeaf_ref/kernels/advec_cell_kernel.f90: 83 - 248
--------------------------------------------------------------------------------

83:     IF(dir.EQ.g_xdir) THEN
[...]
204:           IF(vol_flux_y(j,k).GT.0.0)THEN
[...]
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)
247: 
248:         ENDDO
