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

110:         DO j=x_min,x_max+2
111: 
112:           IF(vol_flux_x(j,k).GT.0.0)THEN
113:             upwind   =j-2
114:             donor    =j-1
115:             downwind =j
116:             dif      =donor
117:           ELSE
118:             upwind   =MIN(j+1,x_max+2)
[...]
124:           sigmat=ABS(vol_flux_x(j,k))/pre_vol(donor,k)
125:           sigma3=(1.0_8+sigmat)*(vertexdx(j)/vertexdx(dif))
126:           sigma4=2.0_8-sigmat
127: 
128:           sigma=sigmat
129:           sigmav=sigmat
130: 
131:           diffuw=density1(donor,k)-density1(upwind,k)
132:           diffdw=density1(downwind,k)-density1(donor,k)
133:           wind=1.0_8
134:           IF(diffdw.LE.0.0) wind=-1.0_8
135:           IF(diffuw*diffdw.GT.0.0)THEN
136:             limiter=(1.0_8-sigmav)*wind*MIN(ABS(diffuw),ABS(diffdw)&
137:               ,one_by_six*(sigma3*ABS(diffuw)+sigma4*ABS(diffdw)))
138:           ELSE
139:             limiter=0.0
140:           ENDIF
141:           mass_flux_x(j,k)=vol_flux_x(j,k)*(density1(donor,k)+limiter)
142: 
143:           sigmam=ABS(mass_flux_x(j,k))/(density1(donor,k)*pre_vol(donor,k))
144:           diffuw=energy1(donor,k)-energy1(upwind,k)
145:           diffdw=energy1(downwind,k)-energy1(donor,k)
146:           wind=1.0_8
147:           IF(diffdw.LE.0.0) wind=-1.0_8
148:           IF(diffuw*diffdw.GT.0.0)THEN
149:             limiter=(1.0_8-sigmam)*wind*MIN(ABS(diffuw),ABS(diffdw)&
150:               ,one_by_six*(sigma3*ABS(diffuw)+sigma4*ABS(diffdw)))
151:           ELSE
152:             limiter=0.0
153:           ENDIF
154: 
155:           ener_flux(j,k)=mass_flux_x(j,k)*(energy1(donor,k)+limiter)
