* Parallelize some subroutines with OpenMP.
  * The parallized subroutines contain the do loops for vertical axis or
    they are associated with ChebyShev polynominal spectral method.