Nonlinear reduced-order models are being investigated for use as digital twins for advanced aircraft and spacecraft. In concept, the digital twin would be tuned to accurately describe the behavior of the vehicle, and then it would be updated as the vehicle is modified or ages. Unfortunately, the linear model correlation and validation techniques that are commonly used in the aerospace industry are not valid for nonlinear models, so a new set of tools is needed. This work presents an algorithm to update geometrically nonlinear reduced-order models using nonlinear normal modes as a correlation metric. The nonlinear normal modes serve as a strong metric to correlate the numerical models because they can be extracted from experiments, and they describe the dynamics of the nonlinear system over a range of amplitudes and are independent of the loading applied to the system. This paper presents a novel method of computing analytical gradients of nonlinear normal mode solutions with respect to system parameters using the multiharmonic balance method. The procedure is first tested using numerical simulations and then applied to tune the reduced-order model of a curved beam based on nominal blueprints and idealized boundary conditions to match a nonlinear normal mode that has been measured experimentally from a three-dimensional printed test specimen.