/home/eoseret/qaas_runs_CPU_9468/171-152-3172/intel/CloverLeafFC/build/CloverLeafFC/CloverLeaf_ref/kernels/calc_dt_kernel.f90: 89 - 133
--------------------------------------------------------------------------------

89:     !$OMP PARALLEL
90: 
91:     !$OMP DO PRIVATE(dsx,dsy,cc,dv1,dv2,div,dtct,dtut,dtvt,dtdivt) REDUCTION(MIN : dt_min_val)
92:     DO k=y_min,y_max
93:     !$OMP SIMD
94:       DO j=x_min,x_max
95: 
96:         dsx=celldx(j)
97:         dsy=celldy(k)
98: 
99:         cc=soundspeed(j,k)*soundspeed(j,k)
100:         cc=cc+2.0_8*viscosity_a(j,k)/density0(j,k)
101:         cc=MAX(SQRT(cc),g_small)
102: 
103:         dtct=dtc_safe*MIN(dsx,dsy)/cc
104: 
105:         div=0.0
106: 
107:         dv1=(xvel0(j  ,k)+xvel0(j  ,k+1))*xarea(j  ,k)
108:         dv2=(xvel0(j+1,k)+xvel0(j+1,k+1))*xarea(j+1,k)
109: 
110:         div=div+dv2-dv1
111: 
112:         dtut=dtu_safe*2.0_8*volume(j,k)/MAX(ABS(dv1),ABS(dv2),g_small*volume(j,k))
113: 
114:         dv1=(yvel0(j,k  )+yvel0(j+1,k  ))*yarea(j,k  )
115:         dv2=(yvel0(j,k+1)+yvel0(j+1,k+1))*yarea(j,k+1)
116: 
117:         div=div+dv2-dv1
118: 
119:         dtvt=dtv_safe*2.0_8*volume(j,k)/MAX(ABS(dv1),ABS(dv2),g_small*volume(j,k))
120: 
121:         div=div/(2.0_8*volume(j,k))
122: 
123:         IF(div.LT.-g_small)THEN
[...]
129:         dt_min_val=MIN(dt_min_val,dtct,dtut,dtvt,dtdivt)
130: 
131:       ENDDO
132:     ENDDO
133:     !$OMP END DO
