nc write now working for tide HC boundary generation

inputs/namelist_remote.bdy

@@ -62,7 +62,7 @@
                                           !  baroclinic velocities
     ln_tra         = .true.               !  boundary conditions for T and S
     ln_ice         = .false.              !  ice boundary condition   
-    nn_rimwidth    = 9                    !  width of the relaxation zone
+    nn_rimwidth    = 1                    !  width of the relaxation zone
 
 !------------------------------------------------------------------------------
 !  unstructured open boundaries tidal parameters                        

pynemo/tide/nemo_bdy_tide3.py

@@ -288,7 +288,10 @@ def nemo_bdy_tide_rot(setup, DstCoord, Grid_T, Grid_U, Grid_V, comp,tide_model):
     sinv = rot_rep(sinuY, sinvY, 'v', 'en to j', dst_gcos, dst_gsin)
 
     #return the values
-    return cosz, sinz, cosu, sinu, cosv, sinv
+    cons = {}
+    cons['cos'] = {'z':cosz,'u':cosu,'v':cosv}
+    cons['sin'] = {'z':sinz,'u':sinu,'v':sinv}
+    return cons
 
 
 def constituents_index(constituents, inputcons):

pynemo/tide/nemo_bdy_tide_ncgen.py

@@ -36,9 +36,11 @@ def CreateBDYTideNetcdfFile(filename, N,I,J,h,fv,grd):
         varz1ID = ncid.createVariable('z1','f4',('yb','xb',),fill_value=fv)
         varz2ID = ncid.createVariable('z2','f4',('yb','xb',),fill_value=fv)
     elif grd == 'U':
+        varmskID = ncid.createVariable('bdy_msk', 'f4', ('y', 'x',), fill_value=fv)
         varu1ID = ncid.createVariable('u1','f4',('yb','xb',),fill_value=fv)
         varu2ID = ncid.createVariable('u2','f4',('yb','xb',),fill_value=fv)
     elif grd == 'V':
+        varmskID = ncid.createVariable('bdy_msk', 'f4', ('y', 'x',), fill_value=fv)
         varv1ID = ncid.createVariable('v1','f4',('yb','xb',),fill_value=fv)
         varv2ID = ncid.createVariable('v2','f4',('yb','xb',),fill_value=fv)
     else :