In the past, I imposed eight major tidal forces (M2,S2,N2,K2,K1,O1,P1,Q1) on the boundary, and the case can be successfully run.
But now I just want to apply the forcing of M2 and K1 tidal components, so I set the My_correc_amp of all the tidal components except for M2 and K1 to 0 through make_tides, and then it blows up during runtime.
So I want to directly inquire about how to apply only M2 and K1 tidal forcing on the boundary.