def main(): # Read hcready dumpfiles after 002_anisopy and change derivatives type # Same particle placement but different diffusion ddwt = "2nw_sd" hc12 = Context() hc12.ReadDumpFile("output/dm.hcready", "output/fd.hcready") hc12.setup.ddw = hc12.esd[ddwt] hc12.Apply() hc12.ContinueRun() hc12.BackupOutput("output-002-1") # print("Done: | ", ddwt, " | ", outrest , "|", # hc12.PrintCell(hc12.setup.resultfile, -1, 3)) print("Error: ", hc12.PrintCell(hc12.setup.resultfile, -1, 3))
def main(): place = Context() resfile = place.GetDateTimeString() place.setup = defaultSetup() place.CleanAllLogs() place.SetThreadsOMP(4) place.SimpleMake() place.CleanRun() place.ReadFinalDump() place.setup.eqs = place.eeq['hydro'] place.setup.tfinish = 100.0 place.setup.dtprint = place.setup.tfinish / 100. place.setup.resultfile = resfile + ".info" place.setup.process = place.epc['fullyperiodic'] place.setup.time = 0.0 place.setup.gamma = 5. / 3. place.setup.disotropic = 1.0 place.setup.dcondconst = 1.0 # place.PrintState() def c(i): return lambda x: i rho = 1.0 mass = place.CalcUniformMass(rho) h = place.setup.hfac * place.setup.spacing prs = 1.0 u = prs / (place.setup.gamma - 1) / rho place.ModifyParticlesProperties( properties=['kappa', 'c', 'm', 'den', 'h', 'p', 'u'], value=[c(1.0), c(1.235e+8), c(mass), c(rho), c(h), c(prs), c(u)]) place.Apply() place.ContinueRun()