anchor.diff:      INTEGER(our_int)                :: k
       INTEGER(our_int)                :: j
   
  -    REAL(our_dble)                  :: draws_prob_raw(num_draws_prob, 4)
  -    REAL(our_dble)                  :: payoffs_systematic(4)
  -    REAL(our_dble)                  :: crit_val_contrib
  -    REAL(our_dble)                  :: total_payoffs(4)
  -    REAL(our_dble)                  :: draws_cond(4)
  -    REAL(our_dble)                  :: draws_stan(4)
  -    REAL(our_dble)                  :: prob_choice
  -    REAL(our_dble)                  :: prob_wage
  -    REAL(our_dble)                  :: prob_obs
  -    REAL(our_dble)                  :: draws(4)
  -    REAL(our_dble)                  :: dist
  -
       LOGICAL                         :: is_deterministic
       LOGICAL                         :: is_working
   
  
positive.code:     REAL(our_dble), INTENT(OUT)     :: rslt
  
      REAL(our_dble), INTENT(IN)      :: periods_payoffs_systematic(num_periods, max_states_period, 4)
      REAL(our_dble), INTENT(IN)      :: data_array(num_agents_est * num_periods, 8)
      REAL(our_dble), INTENT(IN)      :: periods_draws_prob(:, :, :)
      REAL(our_dble), INTENT(IN)      :: shocks_cholesky(4, 4)
      REAL(our_dble), INTENT(IN)      :: periods_emax(:, :)
  
  
negative.code:     !/* internal objects        */
  
      REAL(our_dble)                  :: crit_val(num_agents_est * num_periods)
      REAL(our_dble)                  :: draws_prob_raw(num_draws_prob, 4)
      REAL(our_dble)                  :: payoffs_systematic(4)
      REAL(our_dble)                  :: crit_val_contrib
      REAL(our_dble)                  :: shocks_cov(4, 4)
      REAL(our_dble)                  :: total_payoffs(4)
      REAL(our_dble)                  :: draws_cond(4)
      REAL(our_dble)                  :: draws_stan(4)
      REAL(our_dble)                  :: prob_choice
      REAL(our_dble)                  :: prob_wage
      REAL(our_dble)                  :: prob_obs
      REAL(our_dble)                  :: draws(4)
      REAL(our_dble)                  :: dist
  
      INTEGER(our_int)                :: edu_lagged
      INTEGER(our_int)                :: counts(4)
  