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

81: # 81 "/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"
82:     !$OMP PARALLEL
83: # 83 "/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"
84:     IF(dir.EQ.g_xdir) THEN
85: # 85 "/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"
86:       IF(sweep_number.EQ.1)THEN
87:         !$OMP DO
88:         DO k=y_min-2,y_max+2
89:         !$OMP SIMD
90:           DO j=x_min-2,x_max+2
91:             pre_vol(j,k)=volume(j,k)+(vol_flux_x(j+1,k  )-vol_flux_x(j,k)+vol_flux_y(j  ,k+1)-vol_flux_y(j,k))
92:             post_vol(j,k)=pre_vol(j,k)-(vol_flux_x(j+1,k  )-vol_flux_x(j,k))
93:           ENDDO
94:         ENDDO
95:       !$OMP END DO
96:       ELSE
97:         !$OMP DO
98:         DO k=y_min-2,y_max+2
99:         !$OMP SIMD
100:           DO j=x_min-2,x_max+2
101:             pre_vol(j,k)=volume(j,k)+vol_flux_x(j+1,k)-vol_flux_x(j,k)
102:             post_vol(j,k)=volume(j,k)
103:           ENDDO
[...]
109:       !$OMP            diffuw,diffdw,limiter,wind)
110:       DO k=y_min,y_max
111:         DO j=x_min,x_max+2
112: # 112 "/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"
113:           IF(vol_flux_x(j,k).GT.0.0)THEN
114:             upwind   =j-2
115:             donor    =j-1
116:             downwind =j
117:             dif      =donor
118:           ELSE
119:             upwind   =MIN(j+1,x_max+2)
120:             donor    =j
121:             downwind =j-1
122:             dif      =upwind
123:           ENDIF
124: # 124 "/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"
125:           sigmat=ABS(vol_flux_x(j,k))/pre_vol(donor,k)
126:           sigma3=(1.0_8+sigmat)*(vertexdx(j)/vertexdx(dif))
127:           sigma4=2.0_8-sigmat
128: # 128 "/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"
129:           sigma=sigmat
130:           sigmav=sigmat
131: # 131 "/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"
132:           diffuw=density1(donor,k)-density1(upwind,k)
133:           diffdw=density1(downwind,k)-density1(donor,k)
134:           wind=1.0_8
135:           IF(diffdw.LE.0.0) wind=-1.0_8
136:           IF(diffuw*diffdw.GT.0.0)THEN
137:             limiter=(1.0_8-sigmav)*wind*MIN(ABS(diffuw),ABS(diffdw)&
138:               ,one_by_six*(sigma3*ABS(diffuw)+sigma4*ABS(diffdw)))
139:           ELSE
140:             limiter=0.0
141:           ENDIF
142:           mass_flux_x(j,k)=vol_flux_x(j,k)*(density1(donor,k)+limiter)
143: # 143 "/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"
144:           sigmam=ABS(mass_flux_x(j,k))/(density1(donor,k)*pre_vol(donor,k))
145:           diffuw=energy1(donor,k)-energy1(upwind,k)
146:           diffdw=energy1(downwind,k)-energy1(donor,k)
147:           wind=1.0_8
148:           IF(diffdw.LE.0.0) wind=-1.0_8
149:           IF(diffuw*diffdw.GT.0.0)THEN
[...]
155: # 155 "/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"
156:           ener_flux(j,k)=mass_flux_x(j,k)*(energy1(donor,k)+limiter)
157: # 157 "/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"
158:         ENDDO
159:       ENDDO
160:       !$OMP END DO
161: # 161 "/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"
162:       !$OMP DO PRIVATE(pre_mass_s,post_mass_s,post_ener_s,advec_vol_s)
163:       DO k=y_min,y_max
164:       !$OMP SIMD
165:         DO j=x_min,x_max
166:           pre_mass_s=density1(j,k)*pre_vol(j,k)
167:           post_mass_s=pre_mass_s+mass_flux_x(j,k)-mass_flux_x(j+1,k)
168:           post_ener_s=(energy1(j,k)*pre_mass_s+ener_flux(j,k)-ener_flux(j+1,k))/post_mass_s
169:           advec_vol_s=pre_vol(j,k)+vol_flux_x(j,k)-vol_flux_x(j+1,k)
170:           density1(j,k)=post_mass_s/advec_vol_s
171:           energy1(j,k)=post_ener_s
172:         ENDDO
173:       ENDDO
174:     !$OMP END DO
175: # 175 "/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"
176:     ELSEIF(dir.EQ.g_ydir) THEN
177: # 177 "/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"
178:       IF(sweep_number.EQ.1)THEN
179:         !$OMP DO
180:         DO k=y_min-2,y_max+2
181:         !$OMP SIMD
182:           DO j=x_min-2,x_max+2
183:             pre_vol(j,k)=volume(j,k)+(vol_flux_y(j  ,k+1)-vol_flux_y(j,k)+vol_flux_x(j+1,k  )-vol_flux_x(j,k))
184:             post_vol(j,k)=pre_vol(j,k)-(vol_flux_y(j  ,k+1)-vol_flux_y(j,k))
185:           ENDDO
186:         ENDDO
187:       !$OMP END DO
188:       ELSE
189:         !$OMP DO
190:         DO k=y_min-2,y_max+2
191:         !$OMP SIMD
192:           DO j=x_min-2,x_max+2
193:             pre_vol(j,k)=volume(j,k)+vol_flux_y(j  ,k+1)-vol_flux_y(j,k)
194:             post_vol(j,k)=volume(j,k)
195:           ENDDO
[...]
201:       !$OMP            diffuw,diffdw,limiter,wind)
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
247:           ener_flux(j,k)=mass_flux_y(j,k)*(energy1(j,donor)+limiter)
248: # 248 "/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"
249:         ENDDO
250:       ENDDO
251:       !$OMP END DO
252: # 252 "/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"
253:       !$OMP DO PRIVATE(pre_mass_s,post_mass_s,post_ener_s,advec_vol_s)
254:       DO k=y_min,y_max
255:       !$OMP SIMD
256:         DO j=x_min,x_max
257:           pre_mass_s=density1(j,k)*pre_vol(j,k)
258:           post_mass_s=pre_mass_s+mass_flux_y(j,k)-mass_flux_y(j,k+1)
259:           post_ener_s=(energy1(j,k)*pre_mass_s+ener_flux(j,k)-ener_flux(j,k+1))/post_mass_s
260:           advec_vol_s=pre_vol(j,k)+vol_flux_y(j,k)-vol_flux_y(j,k+1)
261:           density1(j,k)=post_mass_s/advec_vol_s
262:           energy1(j,k)=post_ener_s
263:         ENDDO
264:       ENDDO
265:     !$OMP END DO
266: # 266 "/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"
267:     ENDIF
268: # 268 "/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"
