230
230
! Compute current maximum deviation ...........
233
!$OMP&collapse(2), private(is,ind), reduction(max:dmax)
236
234
dmax = max(dmax, abs(dmnew(ind,is) - dmold(ind,is)))
242
239
! Ensure that dmax is the same on all nodes for convergence test/output