&density_init_pars ldensity_nolog=T /
&density_run_pars diffrho_hyper3=2e-9, idiff='hyper3' /
Rmesh3(e8.1)
2e-14 512^3 2e-11 128^3 6e-9 64^3 2e-8 32^3