The Sky3D code has been widely used to describe nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The approach is based on Skyrme forces or related energy density functionals. The static and dynamic equations are solved on a three-dimensional grid, and pairing is been implemented in the BCS approximation. This updated version of the code aims to facilitate the calculation of nuclear strength functions in the regime of linear response theory, while retaining all existing functionality and use cases. The strength functions are benchmarked against available RPA codes, and the user has the freedom of choice when selecting the nature of external excitation (from monopole to hexadecapole and more). Some utility programs are also provided that calculate the strength function from the time-dependent output of the dynamic calculations of the Sky3D code. New version program summaryProgram Title: Sky3DCPC Library link to program files:https://doi.org/10.17632/vzbrzvyrn4.2Developer's repository link:https://github.com/manybody/sky3dLicensing provisions: GPLv3Programming language: Fortran, with one post-processing utility in Python.Journal reference of previous version:: Schuetrumpf, B., Reinhard, P.G., Stevenson, P.D., Umar, A.S., and Maruhn, J.A. (2018). The TDHF code Sky3D version 1.1. Comput. Phys. Commun. 229 (2018) 211–213.Does the new version supersede the previous version?: Yes.Reasons for the new version: The capability of reproducing the nuclear strength function for a variety of newly-coded external boosts has been added.Nature of problem: Calculating nuclear multipole strength functions is a crucial probe that can help model the nuclear system and its structure properties. A variety of models exist for this task, such as QRPA (Quasiparticle Random Phase Approximation) and its variants, but such approaches are often limited due to symmetry constraints. Time-dependent Hartree Fock (TDHF) has been used to simulate nuclear vibrations and collisions between nuclei for low energies without assuming any symmetry in the system. This code extends the TDHF to calculate the multi-pole strength functions of atomic nuclei. We showcase its reliability by comparing it with the established RPA codes for the calculation of such strength functions.Solution method: We extended previous versions of the Sky3D code [1,2] to include an external boost of multipole type where the user can provide custom input that decides the nature of the multipole (monopole, quadrupole, octupole, and so on) boost. The principal aim is to calculate the multipole strength function, which is the Fourier transform of the time-dependent expectation value of the multipole operator, which has the same form as the external boost. The proper unit conversion is done so that one can extract the exact unit of the thus calculated strength function, which is comparable to the available literature in the field. The boundary conditions are chosen such that a Woods-Saxon-like function cuts off the external field, driving it to zero at the boundary.Summary of revisions: The ability to give a more general external field for excitation has been implemented, which can be of a chosen multipolarity and isospin nature. Then, as a function of time, moments for both isospins and for multipolarities up to L=5 are followed. A new analysis program is included to calculate strength functions and energy weighted sum rules. The ability to specify a mixing parameter in the pairing force is included. This allows a mixed surface-volume pairing to be selected in a way more compatable with QRPA codes for comparison. An updated Makefile is included which makes compilation on Apple silicon computers easier.