Tidal loading computation can be performed using the program load97.exe
You have to enter the options in load97.ini which is prepared for a test run.
The output is found in load97.prn
-------------------------------------------------------------------------------
load97.ini
-------------------------------------------------------------------------------
'csr3.k1'                                      file for ocean tide model
'mascor.csr'                                   file for mass corr. coeff.
' K1'   1                                      nom onde,espece(0=lp,1=d,2=sd)
1,0                                            cor.masse/ph.greenwich (0=non)
2                                              composante (depl rad,grav,pot)
11                                             nombre de stations
PD   Potsdam                52.3810   13.0680    81.000
BR   Bruxelles              50.7986    4.3580   101.000
BH   Bad Homburg            50.2290    8.6110   190.000
ST   Strasbourg             48.6220    7.6840   184.000 
BF   BFO Schiltach          48.3306    8.3300   589.000
CA   Cantley                45.5850  284.1929   269.000
TM   Table Mountain         40.1309  254.7673  1673.680  
KY   Kyoto                  35.0280  135.7860    57.650
PF   Pinon Flat             33.5900  243.5400  1280.000
RI   Richmond               25.6130  279.6150     2.000
WU   Wuhan                  30.5770  114.3600    34.000
-------------------------------------------------------------------------------
DESCRIPTION OF INPUT PARAMETERS
enter the correct option at the beginning of each line

1.  filename of a given wave for a given ocean tide model. 
You have eight cotidal maps available
SCW80 for Q1, O1, P1, K1, N2, M2, S2, K2, Mf
CSR3  for Q1, O1, P1, K1, N2, M2, S2, K2
FES95 for SIGMA1, Q1, O1, PI1, P1, K1, PHI1, J1, 2N2, N2, M2, T2, S2, K2
ORI96 for Q1, O1, P1, K1, N2, M2, S2, K2
CSR4  for Q1, O1, P1, K1, N2, M2, S2, K2
FES02 for Q1, O1, P1, K1, N2, M2, S2, K2, 2N2
GOT00 for Q1, O1, P1, K1, N2, M2, S2, K2
NAO99 for Q1, O1, P1, K1, N2, M2, S2, K2, OO1, M1, J1, L2, 2N2, MU2, NU2, T2, SA, SSA, MF, MM, MSF, MSM, MTM
you can only compute 1 wave in one turn by entering its filename

2.filename for mass correction coefficients. Looking in the source program 
  you will notice that an existing file name is always necessary even if 
  "no mass correction" option is used.
  There are files available for SCW(mascor.scw), CSR3(mascor.csr), CSR4 (mascor.csr4), FES95(mascor.fes)
  , FES02(mascor.fe02), got00(mascor.got0), nao99(mascor.na99)and ORI96(mascor.ori). 

3. name of wave(alphanumeric) and tidal species(integer)

4. mass correction(integer) 0=no, 1=yes
   reference meridian(integer) Greenwich phase(1) or local phase(0)

5. tidal component(integer) 1= radial displacement, 2=gravity, 3=potential

6. number of stations(integer) should be less than 500

7. one line for each station according to model with latitude(N+), 
   longitude(E+) and altitude(m)
-------------------------------------------------------------------------------
load97.prn
-------------------------------------------------------------------------------
   First call of routine GEOEXT, version 1996.08.05.
 Ocean-tide loading for gravity
=========================================
PREM MODEL
PHASE LOCALE
PAS DE CORRECTION DE MASSE
RAYON DECOUPAGE ZONE PROCHE :     5.00000000000000     
 Number of stations :           11
 PD   Potsdam                 37.619    13.068    81.000
 BR   Bruxelles               39.201     4.358   101.000
 BH   Bad Homburg             39.771     8.611   190.000
 ST   Strasbourg              41.378     7.684   184.000
 BF   BFO Schiltach           41.669     8.330   589.000
 CA   Cantley                 44.415   284.193   269.000
 TM   Table Mountain          49.869   254.767  1673.680
 KY   Kyoto                   54.972   135.786    57.650
 PF   Pinon Flat              56.410   243.540  1280.000
 RI   Richmond                64.387   279.615     2.000
 WU   Wuhan                   59.423   114.360    34.000
   0.500000000000000        0.500000000000000     
 Results from tidal model:fes952.k1    

 PD   Potsdam       52.3810    13.0680    81.0000     1.8625    55.1393
 BR   Bruxelles     50.7986     4.3580   101.0000     2.9080    63.7478
 BH   Bad Hombur    50.2290     8.6110   190.0000     2.3627    60.2153
 ST   Strasbourg    48.6220     7.6840   184.0000     2.5258    62.5177
 BF   BFO Schilt    48.3306     8.3300   589.0000     2.4469    63.0595
 CA   Cantley       45.5850   284.1929   269.0000     5.8765    61.0791
 TM   Table Moun    40.1309   254.7673  1673.6800    12.4142    57.3600
 KY   Kyoto         35.0280   135.7860    57.6500    20.8100   351.7870
 PF   Pinon Flat    33.5900   243.5400  1280.0000    30.8740    76.0778
 RI   Richmond      25.6130   279.6150     2.0000     4.0508    43.7117
 WU   Wuhan         30.5770   114.3600    34.0000     7.2295   319.1849

      Routine GEOEXT. Execution time=    21.970 sec
-------------------------------------------------------------------------------

output gives amplitude in nms-2 and phase in degrees
