Here are some small scripts that I wrote for astrophysical purpose, please feel free to use them and send comment to me.

  1. Python script, Perform non-LTE modelling with MCMC (emcee+RADEX)
    Fork HERE: >> radex_emcee
    This code combines emcee and RADEX, performing MCMC sampling on the RADEX-parameters space WITHOUT predefined grids, which leads to faster a convergence time and a better sampling of the parameter space.
  2. GILDAS script, Extract the spectrum from UVFIT or a CLEANed image using GILDAS.
    Download HERE: >>,
    These two scripts ar for extracting the velocity, flux and flux error in each channel from the current *.UVFIT file after performing uvfit, or from the selected region under "go view" in GILDAS. The resulted spectra will be stored into an ASCII file.
  3. IDL script, Plotting radial integrated profile from FITS images using IDL.
    Download HERE: >>
    Sometimes, you need to know the total power in a particular radius of a FITS image, and the growth of this total power along the radius you choose for checking whether the object is extended or not, or contains some ring or disk structures. Therefore, this script could be very helpful for those kinds of quick checks.
  4. IDL script, Filling spectral plot (psym=10) with colors in IDL.
    Download HERE: >>
    It could be difficult when plotting the spectra with IDL with the area between the "flux" and continuum filled in some color. I firstly tried the POLYFILL, but failed with the PSYM, namely unable to set PSYM=10 if using POLYFILL directly. Then I tried to search if there exists any .pro file written by others for this purpose, but turned out to be nothing. So this pro could be useful when plotting some spectrum (e.g., CO, HCN, CS) and fill them with some color, yellow for instance. Updated on 2015-01-20, with a bug fixed.
  5. IDL script, Estimating IRAM/NOEMA time.
    Download HERE: >>
    The pro will calculate the on source and total time needed for requested SNR at NOEMA telescope based on the telescope status described in the 2016 winter IRAM call (this code will be updated according to each call).