FDCHQHP is a Fortran package to calculate the transverse momentum (p(t)) distribution of yield and polarization for heavy quarkonium hadroproduction at next-to-leading-order (NLO) based on non-relativistic QCD(NRQCD) framework. It contains the complete color-singlet and color-octet intermediate states in present theoretical level, and is available to calculate different polarization parameters in different frames. As the LHC running now and in the future, it supplies a very useful tool to obtain theoretical prediction on the heavy quarkonium hadroproduction. Program summary Program title: FDCHQHP v1.0 Catalogue identifier: AETT_v1_0 Program summary URL.: http://cpc.cs.qub.ac.uk/summaries/AETT_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 12020165 No. of bytes in distributed program, including test data, etc.: 103178384 Distribution format: tar.gz Programming language: Fortran 77. Computer: Any computer with Linux operating system, Intel Fortran Compiler and MPI library. Operating system: Linux. Has the code been vectorized or parallelized?: Parallelized with MPI. Classification: 11.1. External routines: MPI Library Nature of problem: This package is for the calculation of the heavy quarkonium hadroproduction at NRQCD NLO. Solution method: The Fortran codes of this package are generated by the FDC system  automatically. Additional comments: It is better to run the package on supercomputers or multi-core computers. !!!!! The distribution file for this program is over 100 MB and therefore is not delivered directly when download or Email is requested. Instead a html file giving details of how the program can be obtained is sent. !!!!! Running time: For an independent sub-process, it may take several seconds to several hours depending on the number of sample points if one CPU core is used. For a complete prompt production of heavy quarkonium hadroproduction (psi (Ns) or Y (Ns)) at one pt point, it may take an hour to thousands of hours depending on the number of sample points if one CPU core is used. In our test with less sample points, it takes 16 min to compile and 25 min to run the whole program with 4 processes with the CPU Intel (R) Core. 17-3770k @ 3.5 GHz x 8. References:  G.T. Bodwin, E. Braaten, G.P. Lepage, Phys. Rev. 051, 1125 (1995). (C) 2014 Elsevier B.V. All rights reserved.