Merge pull request #15 from mathrack/hack_doconcurrent

Initialization of all coefficients had been removed by error.

Author: pbartholomew08

src/parameters.f90

@@ -116,7 +116,11 @@ subroutine parameter()
   dy2 = dy * dy
   dz2 = dz * dz
 
-  xnu=one/re
+  if (abs(re) > epsilon(re)) then
+     xnu = one / re
+  else
+     xnu = zero
+  endif
 
   anglex = sin(pi*angle/180._mytype)
   angley = cos(pi*angle/180._mytype)