anchor.diff:              if (op) dmqmc_calc_type = dmqmc_calc_type + dmqmc_full_r2
               call aot_get_val(op, err, lua_state, table, 'energy', default=.false.)
               if (op) dmqmc_calc_type = dmqmc_calc_type + dmqmc_energy
  +            call aot_get_val(op, err, lua_state, table, 'energy2', default=.false.)
  +            if (op) dmqmc_calc_type = dmqmc_calc_type + dmqmc_energy_squared
               call aot_get_val(op, err, lua_state, table, 'staggered_magnetisation', default=.false.)
               if (op) dmqmc_calc_type = dmqmc_calc_type + dmqmc_staggered_magnetisation
               if (aot_exists(lua_state, table, 'correlation')) then
  
positive.code:         !     beta_loops = Nb,
          !     sampling_weights = { ... },
          !     find_weights = true/false,
          ! }
          ! ipdmqmc = { -- sets propagate_to_beta to true
          !     initial_beta = b,
  
negative.code:             if (err /= 0) call stop_all('read_dmqmc_in', 'spawned_state_size not specified in rdm table.')
              if (aot_exists(lua_state, table, 'rdms')) then
                  dmqmc_in%rdm%doing_rdm = .true.
                  call aot_table_open(lua_state, table, subtable, 'rdms')
                  dmqmc_in%rdm%nrdms = aot_table_length(lua_state, subtable)
                  allocate(rdm_info(dmqmc_in%rdm%nrdms), stat=err)
                  call check_allocate('rdm_info', dmqmc_in%rdm%nrdms, err)
  