Code to compute hydrostatic solutions of coronal loops