<p><b>mhoffman@lanl.gov</b> 2012-07-24 21:41:33 -0600 (Tue, 24 Jul 2012)</p><p>BRANCH COMMIT<br>
Adjusting the ice2sea elevation feedback branch to treat the flux correction separately from SMB.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice_projects/ice2sea_elev_feedback/src/core_land_ice/Registry
===================================================================
--- branches/land_ice_projects/ice2sea_elev_feedback/src/core_land_ice/Registry        2012-07-24 23:51:35 UTC (rev 2054)
+++ branches/land_ice_projects/ice2sea_elev_feedback/src/core_land_ice/Registry        2012-07-25 03:41:33 UTC (rev 2055)
@@ -208,6 +208,8 @@
 var persistent real    SMB_adjusted10 ( TEN nCells ) 0 r SMB_adjusted10 mesh - - 
 % a copy of the initial thickness field - used to calculate the thickness difference at each time step
 var persistent real    upperSurfaceIC ( nCells ) 0 r upperSurfaceIC mesh - - 
+% flux correction needed to maintain geometry in steady state in the absence of a SMB forcing.
+var persistent real    fluxCorrection ( nCells ) 0 ir fluxCorrection mesh - - 
 
 % STATE VARIABLES =====================
 

Modified: branches/land_ice_projects/ice2sea_elev_feedback/src/core_land_ice/mpas_land_ice_tendency.F
===================================================================
--- branches/land_ice_projects/ice2sea_elev_feedback/src/core_land_ice/mpas_land_ice_tendency.F        2012-07-24 23:51:35 UTC (rev 2054)
+++ branches/land_ice_projects/ice2sea_elev_feedback/src/core_land_ice/mpas_land_ice_tendency.F        2012-07-25 03:41:33 UTC (rev 2055)
@@ -157,7 +157,9 @@
      case default
         ! Add surface mass balance to tendency
         ! NOTE: Need to decide where to put this for layer by layer thickness tendency, and how to deal with negative SMB that eliminates top layer or all ice (check for negative thickness?)
-        layerThickness_tend(1,:) = layerThickness_tend(1,:) + mesh % sfcMassBal % array  ! (tendency in meters per year)
+        layerThickness_tend(1,:) = layerThickness_tend(1,:) + mesh % sfcMassBal % array &amp;
+                                   + mesh % fluxCorrection % array ! (tendency in meters per year)  
+        !ice2sea: adding in the flux correction here to keep it separate from the SMB calculations.
         ! THIS MIGHT RESULT IN  NEGATIVE LAYER THICKNESS!
 
         ! Do basal mass balance too?

</font>
</pre>