Codes and computing

Here's Github-link to python/cython codes written for some of my published papers during PhD.


During my PhD work, I developed two relatively small-scale cython codes

(i) a pseudospectral code that captures all the linear instabilities and oscillations in a global atmosphere (with gradients of density, temperature, etc) using a Chebyshev basis for the multiple coupled differential equations (coming from the conservation of fluid properties).  Recently, I have added ideal MHD to this method (not yet updated on GitHub) to explore new global linear instabilities.

(ii) a 1D Lagrangian hydrodynamic code that can evolve gaseous haloes across redshift in an analytically growing dark matter potential well. This method was later extended and explored in multiD in PLUTO code by me and also independently by an undergrad student who worked in my PhD group. 

If I modify or extend these codes I will update on GitHub from time to time. 


Primarily, I have expertise in using PLUTO (hydrodynamic and MHD in C) code. I have several interesting ongoing problem set-ups on PLUTO currently. In the last one year, I have used P3D Particle-in-cell (PIC) code  in the full-particle version to study electron and whistler wave scattering process and consequences. 

In the past (PhD), I have used the ZEUS code (in Fortran) extensively.  I am fully familiar with Athena++ and one of my Part III (Masters) project students in Cambridge used it for simulating turbulent medium.

Analysis & Visualization tool 

I carry out all my analysis and visualization for the published papers using exclusively scripts developed by me in Python/Cython from scratch. Will upload some of the scripts on GitHub page soon. Stay tuned! 

Occasionally, I use Mathematica (not always a flexible workflow since it requires licence).  To create intuitive cartoons for physical phenomenon I have  often used Apple KeyNote (which I also use to prepare talk slides).

Systems I have used for computing

Sahasrat (IISc), CSD3 (Cambridge), Cori (NERSC), ARC (Oxford) - - I have used the first two systems extensively, third for running some simulations, and a new user of the fourth now (i.e. Nov 2023). Briefly used Dirac Tesseract (seedcorn) for testing a large number of 2D MHD simulations I designed.