/home/eoseret/qaas_runs_CPU_9468/171-152-3172/intel/CloverLeafFC/build/CloverLeafFC/CloverLeaf_ref/kernels/PdV_kernel.f90: 69 - 99
--------------------------------------------------------------------------------

69:     IF(predict)THEN
[...]
75:         DO j=x_min,x_max
76: 
77:           left_flux=  (xarea(j  ,k  )*(xvel0(j  ,k  )+xvel0(j  ,k+1)                     &
78:             +xvel0(j  ,k  )+xvel0(j  ,k+1)))*0.25_8*dt*0.5
79:           right_flux= (xarea(j+1,k  )*(xvel0(j+1,k  )+xvel0(j+1,k+1)                     &
80:             +xvel0(j+1,k  )+xvel0(j+1,k+1)))*0.25_8*dt*0.5
81:           bottom_flux=(yarea(j  ,k  )*(yvel0(j  ,k  )+yvel0(j+1,k  )                     &
82:             +yvel0(j  ,k  )+yvel0(j+1,k  )))*0.25_8*dt*0.5
83:           top_flux=   (yarea(j  ,k+1)*(yvel0(j  ,k+1)+yvel0(j+1,k+1)                     &
84:             +yvel0(j  ,k+1)+yvel0(j+1,k+1)))*0.25_8*dt*0.5
85:           total_flux=right_flux-left_flux+top_flux-bottom_flux
86: 
87:           volume_change_s=volume(j,k)/(volume(j,k)+total_flux)
[...]
95:           energy_change=(pressure(j,k)/density0(j,k)+viscosity(j,k)/density0(j,k))*total_flux*recip_volume
96: 
97:           energy1(j,k)=energy0(j,k)-energy_change
98: 
99:           density1(j,k)=density0(j,k)*volume_change_s
