Add external forces. Remove duplicated rot differentiation. Make dt a global var.
Use upsampled open/close, force, force_ext instead of resampled ones