Simulate matrix with 'GWAS signal'

We are evaluating methods to run GWAS and post-GWAS analyses. It would be great to be able to simulate a multi-ancestry hail matrix that causes a GWAS analysis to produce “buildings” in a Manhattan plot.
Is there a parametrization of the Balding Nichols model or the phenotype simulations in the ldscsim module to generate such a matrix?
thanks!!

Juan Felipe Ortiz