diff --git a/inputs/namelist_remote.bdy b/inputs/namelist_remote.bdy
index 77680055fe538e986b74f6560a525ac8fd431435..48a1d8cbd2757067abd7bbe13c3ac209eadb1959 100644
--- a/inputs/namelist_remote.bdy
+++ b/inputs/namelist_remote.bdy
@@ -68,7 +68,7 @@
 !  unstructured open boundaries tidal parameters                        
 !------------------------------------------------------------------------------
     ln_tide        = .true.              !  =T : produce bdy tidal conditions
-    sn_tide_model  = 'fes'                !  Name of tidal model (fes|tpxo)
+    sn_tide_model  = 'tpxo'                !  Name of tidal model (fes|tpxo)
     clname(1)      = 'M2'                 !  constituent name
     clname(2)      = 'S2'
     ln_trans       = .true.               !  interpolate transport rather than
diff --git a/pynemo/tide/fes_extract_HC.py b/pynemo/tide/fes_extract_HC.py
index 6aaebf618de63a017d70c0453b5625b59015adb1..7f87b040d23b83092ec4d66389f70fc5984da891 100644
--- a/pynemo/tide/fes_extract_HC.py
+++ b/pynemo/tide/fes_extract_HC.py
@@ -229,7 +229,6 @@ class HcExtract(object):
         maskedpoints = interpolate.interpn((x_values, y_values), mask, lonlat)
 
         data_temp = np.zeros((data.shape[0], lon.shape[0], 2, ))
-        #check at same point in TPXO extract HC sc
         for cons_index in range(data.shape[0]):
             #interpolate real values
             data_temp[cons_index, :, 0] = interpolate_data(x_values, y_values,