Post-Processing

[Lorenzo Campana]

Dear Users,

I am trying to do some post-processing with Dedalus 3. I have the problem to go back to dealiased grid:

1. Read Data:

write, dt = solver.load_state('snapshots_s39.h5', -1, allow_missing=True)

2. Compute stuff with dealiasing as the main problem (vorticity)

w = dx(v1)-dy(u1)

3. Go back to original resolution?

... I have tried w = w.change_scales(1,1) doesn't work with error.

Someone could help me?

Many tanks.

Lorenzo

[Miguel Beneitez]

Hi Lorenzo,

I think that you might be missing the evaluation of the operation. See here: https://dedalus-project.readthedocs.io/en/latest/notebooks/dedalus_tutorial_2.html

I also attach a slightly modified version of the tutorial for Rayleigh-Benard with some dummy post processing at the end trying to do the same thing that you mention. This works well for me without throwing any errors.

Cheers,

Miguel