anchor.diff:                          endogenous_slaves(count) = emax_simulated - maxe(k + 1)
                           count = count + 1 
   
  -
  -
                       END DO
                       
  -
                       ! Distribute exogenous information
  -                    ! TODO: POLYMORPHISM
                       CALL distribute_information(num_emax_slaves, period, endogenous_slaves, endogenous)
                       
                       ! Create prediction model based on the random subset of points where the EMAX is actually simulated and thus endogenous and exogenous variables are available. For the interpolation  points, the actual values are used.
  
positive.code:     ! This part creates (or reads from disk) the draws for the Monte Carlo integration of the EMAX. For is_debugging purposes, these might  also be read in from disk or set to zero/one.   
      CALL create_draws(periods_draws_emax, num_draws_emax, seed_emax)
  
  
      ALLOCATE(draws_emax(num_draws_emax, 4))
  
      IF(rank == zero_int) THEN
          is_head = .True.
      ELSE
          is_head = .False.
      END IF
  
      DO WHILE (STAY_AVAILABLE)  
          
  
negative.code:     ! Read in model specification.
      CALL read_specification(coeffs_a, coeffs_b, coeffs_edu, coeffs_home, shocks_cholesky)
  
      ! Allocate arrays
      IF(rank == 0) CALL logging_solution(1)
  
  