Hello,
For TPXO database, you don't need to use your own harmonic file. You can directly use TIDAL DATA BASE = 2 and provide the BINARY DATABASEs 1-2 FOR TIDE. But you'll have to code something to only consider the 4 constituents you mention.
If you want to use your own harmonic file, the format of the file is not correct. In the header, you have to provide as follow:
- 1st line: NUMBER OF THE TIDAL BOUNDARIES (3 or 4 in your case, in don't really know)
- following lines (as long as there are tidal boundaries): THE FIRST AND LAST INDICES OF THE OPEN LIQUID BOUNDARY WITH TIDE TO PRESCRIBE for each tidal boundary
- and then:
"Total number of tidal boundary points" "number of constituents" "M"
Then you have to write for each of your 4 constituents, the amplitude and phase for the whole tidal boundary points.
I attach an example with 61 tidal constituents and 2 tidal boundaries. The first boundary goes from boundary points 1 to 55, and the second one from 4721 to 4778, for a total number of tidal boundary points of 113 = (55-1+1) + (4778-4721+1).
I hope it helps.
Regards,
Laurent