anchor.diff:  
          allocate(rperms%perm(nOp,nD*n),taperms%perm(nOp, size(arrow_basis,2)),STAT=status)
          if (status/=0) stop "Allocation failed in get_rotation_perm_lists: rperms%perm"
  +       write(19,'("<< d-g Table >>")')
          do iOp = 1, nOp ! For each rotation, find the permutation
             dgp = 0 ! Initialize the (d,g) table
             dap = 0
  
positive.code: 
      ! Beginning of main routine for enumeration
      open(23,file="VERSION.enum")
      write(23,'(A)') "v2.0.4-40-gfc3d-dirty"
      close(23)
  
      ![TODO] Get rid of all the junk that crept in (writing files, making inactives table, etc. These should all be in routines so that this main routine is still readable)
  
negative.code:                 end do !jm
               end do !im
            enddo ! loop over d-vectors (each row in the table)
            if (any(dgp==0) .or. (any(dap==0) .and. arrows)) then
              do im = 1,nD; write(19,'(20(1x,i2))') dgp(im,:); enddo
              stop "(d,g)-->(d',g') mapping failed in get_rotation_perm_lists"; endif
  
            ! Now we have the (d',g') table for this rotation. Now
            ! record the permutation
  