/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/aocc_2/CMakeFiles/clover_leaf.dir/CloverLeaf_ref/kernels/viscosity_kernel.f90-pp.f90: 56 - 92
--------------------------------------------------------------------------------

56:       DO j=x_min,x_max
57:         ugrad=(xvel0(j+1,k  )+xvel0(j+1,k+1))-(xvel0(j  ,k  )+xvel0(j  ,k+1))
58: # 58 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/viscosity_kernel.f90"
59:         vgrad=(yvel0(j  ,k+1)+yvel0(j+1,k+1))-(yvel0(j  ,k  )+yvel0(j+1,k  ))
60: # 60 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/viscosity_kernel.f90"
61:         div = (celldx(j)*(ugrad)+  celldy(k)*(vgrad))
62: # 62 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/viscosity_kernel.f90"
63:         strain2 = 0.5_8*(xvel0(j,  k+1) + xvel0(j+1,k+1)-xvel0(j  ,k  )-xvel0(j+1,k  ))/celldy(k) &
64:           + 0.5_8*(yvel0(j+1,k  ) + yvel0(j+1,k+1)-yvel0(j  ,k  )-yvel0(j  ,k+1))/celldx(j)
65: # 65 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/viscosity_kernel.f90"
66:         pgradx=(pressure(j+1,k)-pressure(j-1,k))/(celldx(j)+celldx(j+1))
67:         pgrady=(pressure(j,k+1)-pressure(j,k-1))/(celldy(k)+celldy(k+1))
68: # 68 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/viscosity_kernel.f90"
69:         pgradx2 = pgradx*pgradx
70:         pgrady2 = pgrady*pgrady
71: # 71 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/viscosity_kernel.f90"
72:         limiter = ((0.5_8*(ugrad)/celldx(j))*pgradx2+(0.5_8*(vgrad)/celldy(k))*pgrady2+strain2*pgradx*pgrady)  &
73:           /MAX(pgradx2+pgrady2,1.0e-16_8)
74: # 74 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/viscosity_kernel.f90"
75:         IF ((limiter.GT.0.0).OR.(div.GE.0.0))THEN
76:           viscosity(j,k) = 0.0
77:         ELSE
78:           dirx=1.0_8
79:           IF(pgradx.LT.0.0) dirx=-1.0_8
80:           pgradx = dirx*MAX(1.0e-16_8,ABS(pgradx))
81:           diry=1.0_8
82:           IF(pgradx.LT.0.0) diry=-1.0_8
83:           pgrady = diry*MAX(1.0e-16_8,ABS(pgrady))
84:           pgrad = SQRT(pgradx**2+pgrady**2)
85:           xgrad = ABS(celldx(j)*pgrad/pgradx)
86:           ygrad = ABS(celldy(k)*pgrad/pgrady)
87:           grad  = MIN(xgrad,ygrad)
88:           grad2 = grad*grad
89: # 89 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/viscosity_kernel.f90"
90:           viscosity(j,k)=2.0_8*density0(j,k)*grad2*limiter*limiter
91:         ENDIF
92: # 92 "/beegfs/hackathon/users/eoseret/qaas_runs_test/178-661-4072/intel/CloverLeaf1.3-FC/build/CloverLeaf1.3-FC/CloverLeaf_ref/kernels/viscosity_kernel.f90"
