Adding nodal attribute (#251)

* adding geoid offset nodal attribute * Update README.md
CHLNDDEV · Aug 10, 2021 · cada81f · cada81f
1 parent e0bae23
commit cada81f
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 1 deletion.
diff --git a/README.md b/README.md
@@ -148,7 +148,8 @@ Changelog
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/)
 
 ### Unreleased (on current HEAD of the Projection branch)
-- None
+## Added
+- Geoid offset nodal attribute in `Calc_f13` subroutine. https://github.com/CHLNDDEV/OceanMesh2D/pull/251
 
 ### [5.0.0] - 2021-07-29
 ## Added

diff --git a/utilities/Calc_f13.m b/utilities/Calc_f13.m
@@ -14,6 +14,8 @@
 %              'Ad' ('advection_state')
 %              'Sb' ('subgrid_barrier')
 %              'Es' ('elemental_slope_limiter')
+%              'Ge' ('sea_surface_height_above_geoid')
+%
 %
 %            3) then either: 
 %              'inpoly' followed by...
@@ -59,6 +61,9 @@
 elseif strcmpi(attribute,'Es')
      attrname = 'elemental_slope_limiter';
      default_val = -1e-3;
+elseif strcmpi(attribute,'Ge')
+     attrname = 'sea_surface_height_above_geoid';
+     default_val = 0.0; 
 else
     error(['Attribute ' attribute ' not currently supported'])
 end