Scottyujia
Scottyujia
Thanks very much for your help, Georgia! Using 'rh' for the initial condition has definitely solved the issue for the case I posted above. I just tried it for a...
Oh, I see. You mean the init_g.f90 file.
Option 2 does not seem to work for the new case, either: phi2_vs_t
In case you are thinking about up-windings. I have tried setting vpa and zed up-windings to zero and only use hyper dissipation in vpa and zed. It does not get...
The input you attached is for CBC. The plots I showed above are for the new one which is different: [input.in.zip](https://github.com/user-attachments/files/23123959/input.in.zip)
Latest update: The problem was solved by using initial condition `rh' and implicit_drifts = .true.
I can confirm that it also works. Let me summarise the findings so far for the benefit of others reading this issue: **Explicit drifts:** No numerical instability for the zonal...
Hi, Toby. This is related to the beginning of your post where you tried to do benchmark with ky = 0.0707 in stella. I tried the following inputs and I...
Entire directory as requested [entire_dir.zip](https://github.com/user-attachments/files/21674159/entire_dir.zip)
I think if you just turn off zed_upwinding, it will agree very well. I will up load some data and plots soon.