How to apply a forcing on the boundary with only M2 and K1 internal tides?

Hello everyone,

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.

I would sincerely appreciate any help and advice.

Thanks,
Sihang Xie

Hi,

I think you only need to modify two things:

  1. The crocotools_param.m file, from

Ntides=10;
%
% Chose order from the rank in the TPXO file :
% “M2 S2 N2 K2 K1 O1 P1 Q1 Mf Mm”
% " 1 2 3 4 5 6 7 8 9 10"
%
tidalrank=[1 2 3 4 5 6 7 8 9 10];
%

to

Ntides=2;
%
% Chose order from the rank in the TPXO file :
% “M2 S2 N2 K2 K1 O1 P1 Q1 Mf Mm”
% " 1 2 3 4 5 6 7 8 9 10"
%
tidalrank=[1 5];
%

  1. In param.h change

else

  parameter (Ntides=8)

endif

to

else

  parameter (Ntides=2)

endif

Please try this.

2 Likes

Hi addres,
Thank you for your reply and the detailed explanation, now the case can be successfully run.
My Best Regards
Sihang Xie