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()