anchor.diff:  
       IF(rank == 0) CALL logging_solution(-1)
   
  -
       ! Determine workload and allocate communication information.
       CALL determine_workload(num_emax_slaves, states_number_period)
   
  -
       ! Calculate the systematic payoffs
       IF(rank == 0) CALL logging_solution(2)
   
  
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: 
                      ! Constructing indicator for simulation points
                      is_simulated = get_simulated_indicator(num_points_interp, num_states, period)
         
                      ! Constructing the dependent variable for all states, including the ones where simulation will take place. All information will be used in either the construction of the prediction model or the prediction step.
                      CALL get_exogenous_variables(exogenous, maxe, period, num_states, periods_payoffs_systematic, shifts, mapping_state_idx, periods_emax, states_all)
  
                      ! Initialize missing values
                      endogenous = MISSING_FLOAT
  