anchor.diff:    ! -------------------------------------------------------------
     function initialize_vector(seed) result(twister)
       integer(kind=im), dimension(0:), intent(in) :: seed
  -    type(randomNumberSequence)                      :: twister 
  -    
  +    type(randomNumberSequence)                      :: twister
  +
       integer(kind=im) :: i, j, k, nFirstLoop, nWraps
  -    
  +
       nWraps  = 0
       twister = initialize_scalar(19650218_im)
  -    
  +
       nFirstLoop = max(blockSize, size(seed))
       do k = 1, nFirstLoop
          i = mod(k + nWraps, blockSize)
  
positive.code:             fac111 = fk1*fac11(lay)
              fac211 = fk2*fac11(lay)
           else if (specparm1 .gt. 0.875_rb) then
              p = -fs1 
              p4 = p**4
              fk0 = p4
              fk1 = 1 - p - 2.0_rb*p4
  
negative.code: !
  !  name     type     purpose
  !  ----   : ----   : ---------------------------------------------
  !fracrefao: real
  !fracrefbo: real
  ! kao     : real
  ! kbo     : real
  ! selfrefo: real
  ! forrefo : real
  !-----------------------------------------------------------------
  
        integer(kind=im), parameter :: no10 = 16
  