docs: def PROFILE_HT(sg0,GamD,Gam0,Gam2,Shift0,Shift2,anuVC,eta,sg):
      """
      #-------------------------------------------------
      #      "pCqSDHC": partially-Correlated quadratic-Speed-Dependent Hard-Collision
      #      Subroutine to Compute the complex normalized spectral shape of an 
      #      isolated line by the pCqSDHC model
      #
      #      References:
      #
      #      1) N.H. Ngo, D. Lisak, H. Tran, J.-M. Hartmann.
      #         An isolated line-shape model to go beyond the Voigt profile in 
      #         spectroscopic databases and radiative transfer codes.
      #         JQSRT, Volume 129, November 2013, Pages 89–100
      #         http://dx.doi.org/10.1016/j.jqsrt.2013.05.034
      #
      #      2) H. Tran, N.H. Ngo, J.-M. Hartmann.
      #         Efficient computation of some speed-dependent isolated line profiles.
      #         JQSRT, Volume 129, November 2013, Pages 199–203
      #         http://dx.doi.org/10.1016/j.jqsrt.2013.06.015
      #
      #      3) H. Tran, N.H. Ngo, J.-M. Hartmann.
      #         Erratum to “Efficient computation of some speed-dependent isolated line profiles”.
      #         JQSRT, Volume 134, February 2014, Pages 104
      #         http://dx.doi.org/10.1016/j.jqsrt.2013.10.015
      #
      #      Input/Output Parameters of Routine (Arguments or Common)
      #      ---------------------------------
      #      T       : Temperature in Kelvin (Input).
      #      amM1    : Molar mass of the absorber in g/mol(Input).
      #      sg0     : Unperturbed line position in cm-1 (Input).
      #      GamD    : Doppler HWHM in cm-1 (Input)
      #      Gam0    : Speed-averaged line-width in cm-1 (Input).       
      #      Gam2    : Speed dependence of the line-width in cm-1 (Input).
      #      anuVC   : Velocity-changing frequency in cm-1 (Input).
      #      eta     : Correlation parameter, No unit (Input).
      #      Shift0  : Speed-averaged line-shift in cm-1 (Input).
      #      Shift2  : Speed dependence of the line-shift in cm-1 (Input)       
      #      sg      : Current WaveNumber of the Computation in cm-1 (Input).
      #
      #      The function has two outputs:
      #      -----------------
      #      (1): Real part of the normalized spectral shape (cm)
      #      (2): Imaginary part of the normalized spectral shape (cm)
      #
      #      Called Routines: 'CPF'       (Complex Probability Function)
      #      ---------------  'CPF3'      (Complex Probability Function for the region 3)
      #
      #      Based on a double precision Fortran version
      #
      #-------------------------------------------------
      """
      return pcqsdhc(sg0,GamD,Gam0,Gam2,Shift0,Shift2,anuVC,eta,sg)