anchor.diff:  
       ! 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))
   
       is_head = .False.
       IF(rank == zero_int) is_head = .True.
  
positive.code:     IF(rank == 0) CALL logging_solution(-1)
  
      ! Determine workload and allocate communication information.
      ALLOCATE(num_emax_slaves(num_periods, num_slaves), num_obs_slaves(num_slaves))
  
      DO period = 1, num_periods
          CALL determine_workload(num_emax_slaves(period, :), states_number_period(period))   
  
negative.code:                  IF (any_interpolated) THEN
  
                      ! Allocate period-specific containers
                      ALLOCATE(is_simulated(num_states), endogenous(num_states), maxe(num_states), exogenous(num_states, 9), predictions(num_states))
  
                      ! Constructing indicator for simulation points
                      is_simulated = get_simulated_indicator(num_points_interp, num_states, period)
  