anchor.diff:  
       !/* internal objects        */
   
  -    REAL(our_dble), ALLOCATABLE    :: crit_val(:)
  +    REAL(our_dble), ALLOCATABLE     :: crit_val(:)
   
       REAL(our_dble)                  :: draws_prob_raw(num_draws_prob, 4)
       REAL(our_dble)                  :: payoffs_systematic(4)
  
positive.code: 
      USE shared_constants
  
      USE shared_containers
  
  	!/*	setup	*/
  
      IMPLICIT NONE
  
negative.code:             IF (is_working) THEN
  
                  ! Calculate the disturbance, which follows a normal distribution.
                  dist = clip_value(LOG(data_evaluate(j, 4)), -HUGE_FLOAT, HUGE_FLOAT) - clip_value(LOG(payoffs_systematic(idx)), -HUGE_FLOAT, HUGE_FLOAT) 
  
                  ! If there is no random variation in payoffs, then the observed wages need to be identical their systematic components. The discrepancy between the observed wages and their systematic components might be small due to the reading in of the dataset.
                  IF (is_deterministic) THEN
  