Hi Yang,
First of all, I think you should go step by step and add complexity progressively to your case.
Did you run a 2D model first with only tide? then a 3D with only tide? after that, you can try adding the regional circulation.
You didn't provide the execution log of your model, but according to what I see in your steering, there are a lot of things that explain your model can't correctly run, even in 2D with only tide:
- your mesh is in lat/lon real degrees, and you set GEOGRAPHIC SYSTEM = 5 which is not correct, no need to provide ZONE NUMBER IN GEOGRAPHIC SYSTEM in this case.
- COEFFICIENT TO CALIBRATE TIDAL VELOCITIES = 250 ! is it really what you want?
- Nikuradse coefficient for bottom friction law = 5m, is it really what you want?
- PRESCRIBED ELEVATIONS and VELOCITIES : 2 values provided, only one boundary..
- use_fortran: you didn't provide your modifications..
- 5 vertical levels, it's not much..
- TRACERS VERTICAL PROFILES = 0 -> you didn't provide your implementation
In a general way, your steering file is really a mess. Try to regroup the processes and related options to make it clearer.
So please, come back to simple case of tidal 2D simulation (maybe with a larger domain, yours is really small for tidal and regional circulation).
After this stage, you can try to extend it in 3D.
I hope it helps,
Regards,
Laurent