How to define maximum of Intermediate and another value in Python Gekko, when using sequential solver?

Question

In a modeling framework for solving systems of differential equations with GEKKO that I am writing I want to calculate the maximum of a parameter (an external forcing) and a