The scripts and workflow are supplementary material to "3D Modelling of Vertical Gravity Gradients and the delimitation of tectonic boundaries: The Caribbean oceanic domain as a case study" (Gómez-García et al., 2019). The codes include the calculation of the VGG response of a 3D lithospheric model, in spherical coordinates, using the software Tesseroids (Uieda, 2016). The "Readme_Workflow_2019_002.pdf" file provide very detail information about the structure of this repository, as well as the step-by-step for the scripts execution, and the list of the requiered software for the correct workflow performance. All the information provided here will allow the user to reproduce the results and figures of the main paper. Detailed information are also given in the associated README.