From cafc39b183363b64b4c88fc3cbbf77b738ca555b Mon Sep 17 00:00:00 2001
From: Denny	Tjahjanto <dennyt@kth.se>
Date: Thu, 13 Aug 2009 13:21:22 +0000
Subject: [PATCH] replaces Nast with Nreps in hybridIA

formerly, if the number of orientations requested was (far) less than the number of nonzero bins in the texture representation, you only got a subset of low phi1 values..! (that is the top corner in Euler space -- Hi Bing ;-) )
---
 code/IO.f90 | 11 ++++-------
 1 file changed, 4 insertions(+), 7 deletions(-)

diff --git a/code/IO.f90 b/code/IO.f90
index 6dae90389..0da6d9cb9 100644
--- a/code/IO.f90
+++ b/code/IO.f90
@@ -201,7 +201,7 @@
      enddo
    enddo
  enddo  
- 
+
  dV_V = dV_V/sum_dV_V  ! normalize to 1
  
 !--- now fix bounds ---
@@ -229,7 +229,7 @@
      exit
    endif
  enddo
- 
+
  allocate(binSet(Nreps))
  bin = 0 ! bin counter
  i = 1 ! set counter
@@ -247,7 +247,7 @@
  do i=1,Nast
    if (i < Nast) then
      call random_number(rnd)
-     j = nint(rnd*(Nast-i)+i+0.5_pReal,pInt)
+     j = nint(rnd*(Nreps-i)+i+0.5_pReal,pInt)
    else
      j = i
    endif
@@ -858,7 +858,7 @@ endfunction
  case (265)
    msg = 'Limit for crystallite loop too small'
  case (266)
-   msg = 'Limit for crystallite state loop too small'
+   msg = 'Limit for state loop too small'
  case (267)
    msg = 'Limit for stress loop too small'
  case (268)
@@ -884,9 +884,6 @@ endfunction
  case (277)
    msg = 'Non-positive relevant mismatch in RGC'
 
- case (279)
-   msg = 'Limit for materialpoint state loop too small'
-
  case (300)
    msg = 'This material can only be used with elements with three direct stress components'
  case (500)