Computing derivative information of sequentially coupled subsurface models