diff --git a/Interface.py b/Interface.py
index b340478862e10f99c5786c226842784a5b260195..9ef45cb585002eb3ba888beedfabb70f706145c4 100644
--- a/Interface.py
+++ b/Interface.py
@@ -407,14 +407,10 @@ class HydroSuper :
         nbpt_glo = part.l2glandindex(nbpt_loc)
         self.add_variable(outnf, procgrid, NCFillValue, part, ('y','x'), "nbpt_glo", "Grid point Global", "-", nbpt_glo[:,0], vtyp)
         #
-        # gridarea
+        # contfrac
         contfrac = np.array(procgrid.contfrac)
         self.add_variable(outnf, procgrid, NCFillValue, part, ('y','x'), "contfrac", "Land fraction", "-", np.array(procgrid.contfrac), vtyp)
         #
-        # gridarea
-        nbptarea = np.array(procgrid.area)
-        self.add_variable(outnf, procgrid, NCFillValue, part, ('y','x'), "gridarea", "Grid area", "-", nbptarea, vtyp)
-        #
         # basin_id
         self.add_variable(outnf, procgrid, NCFillValue, part, ('htuext','y','x'), "basin_id", "ID for each HTU", "-", self.basin_id, vtyp)
         #
@@ -462,9 +458,9 @@ class HydroSuper :
         # Inflow Grid -> convert to global
         gingrid = part.l2glandindex(self.inflow_grid[:,:,:inflow_size])
         self.add_variable(outnf, procgrid, NCFillValue, part, ('in','htu','y','x'), "HTUingrid", "Inflow grid", "-", gingrid, vtyp)
-
+        #
+        # Inflow Basin
         self.add_variable(outnf, procgrid, NCFillValue, part, ('in','htu','y','x'), "HTUinbas", "Inflow basin", "-", self.inflow_basin[:,:,:inflow_size], vtyp)
-
         #
         # Save the fetch of each basin
         #