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

67:     !$OMP PARALLEL
68: 
69:     IF(predict)THEN
70: 
71:       !$OMP DO PRIVATE(right_flux,left_flux,top_flux,bottom_flux,total_flux,min_cell_volume, &
72:       !$OMP            energy_change,recip_volume,volume_change_s)
73:       DO k=y_min,y_max
74:       !$OMP SIMD
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
100: 
101:         ENDDO
102:       ENDDO
103:     !$OMP END DO
104: 
105:     ELSE
106: 
107:       !$OMP DO PRIVATE(right_flux,left_flux,top_flux,bottom_flux,total_flux,min_cell_volume, &
108:       !$OMP            energy_change,recip_volume,volume_change_s)
109:       DO k=y_min,y_max
110:       !$OMP SIMD
111:         DO j=x_min,x_max
112: 
113:           left_flux=  (xarea(j  ,k  )*(xvel0(j  ,k  )+xvel0(j  ,k+1)                     &
114:             +xvel1(j  ,k  )+xvel1(j  ,k+1)))*0.25_8*dt
115:           right_flux= (xarea(j+1,k  )*(xvel0(j+1,k  )+xvel0(j+1,k+1)                     &
116:             +xvel1(j+1,k  )+xvel1(j+1,k+1)))*0.25_8*dt
117:           bottom_flux=(yarea(j  ,k  )*(yvel0(j  ,k  )+yvel0(j+1,k  )                     &
118:             +yvel1(j  ,k  )+yvel1(j+1,k  )))*0.25_8*dt
119:           top_flux=   (yarea(j  ,k+1)*(yvel0(j  ,k+1)+yvel0(j+1,k+1)                     &
120:             +yvel1(j  ,k+1)+yvel1(j+1,k+1)))*0.25_8*dt
121:           total_flux=right_flux-left_flux+top_flux-bottom_flux
122: 
123:           volume_change_s=volume(j,k)/(volume(j,k)+total_flux)
[...]
131:           energy_change=(pressure(j,k)/density0(j,k)+viscosity(j,k)/density0(j,k))*total_flux*recip_volume
132: 
133:           energy1(j,k)=energy0(j,k)-energy_change
134: 
135:           density1(j,k)=density0(j,k)*volume_change_s
136: 
137:         ENDDO
138:       ENDDO
139:     !$OMP END DO
