convection in annulus geometry

[yanick ricard]

Hello and happy new year!

I have been running dedalus for some time in cartesian and spherical geometry and

now I try to simulate convection (with infinite Prandlt number) in annulus geometry.

I just get:

****

annulus = d3.AnnulusBasis(coords, shape=(Nphi, Nr), radii=(Ri, Ro), dealias=dealias, dtype=dtype)

  File "/Users/ricard/anaconda3/envs/dedalus3/lib/python3.11/site-packages/dedalus/core/basis.py", line 1976, in __init__

    super().__init__(coordsystem, shape, dtype, k=k, dealias=tuple(dealias), azimuth_library=azimuth_library)

TypeError: 'float' object is not iterable

****

Even when I tried to get other 2D codes in annulus geometry found in the dedalus-Users mails (that should work), I get a similar answer...

Can somebody help me?

By the way Dedalus is a great tool!

Yan Ricard (ENS de Lyon)

[Keaton Burns]

Hi Yan,

This may be hitting some recent changes based on how we handle dealiasing for the annulus. Please try with the latest released version (v3.0.0, which is now on conda and pip).  Otherwise I think the tau fields need to be modified to use either one of the edge bases.

Best,

-Keaton

--
You received this message because you are subscribed to the Google Groups "Dedalus Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dedalus-user...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/dedalus-users/09dd4f2c-196a-4e3b-8331-b9690fce6620n%40googlegroups.com.

[Daniel Lecoanet]

Hi Yanick,

The following script runs for me on the latest version of Dedalus.

Daniel

To view this discussion on the web visit https://groups.google.com/d/msgid/dedalus-users/CADZXxBgAEwyZfva-YHnnY8rruZGujMir-yTfXjgmK0sDZa%2B-og%40mail.gmail.com.

[yanick ricard]

Thanks Keaton and Daniel for your prompt answer !

I must not have the latest released version as the code (the version sent back by Daniel) still, does not run. I will update my version tomorrow.

Best