<p><b>dwj07@fsu.edu</b> 2013-02-20 13:32:54 -0700 (Wed, 20 Feb 2013)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Adding an optional seaLevelPressure input field.<br>
<br>
        This is added to the pressure field in the top layer to depress the surface of the ocean.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/sea_level_pressure/src/core_ocean/Registry
===================================================================
--- branches/ocean_projects/sea_level_pressure/src/core_ocean/Registry        2013-02-20 18:47:45 UTC (rev 2492)
+++ branches/ocean_projects/sea_level_pressure/src/core_ocean/Registry        2013-02-20 20:32:54 UTC (rev 2493)
@@ -371,3 +371,6 @@
 var persistent integer edgeSignOnCell ( maxEdges nCells ) 0 - edgeSignOnCell mesh - -
 var persistent integer edgeSignOnVertex ( maxEdges nVertices ) 0 - edgeSignOnVertex mesh - -
 var persistent integer kiteIndexOnCell ( maxEdges nCells ) 0 - kiteIndexOnCell mesh - -
+
+% Sea level pressure, for coupling
+var persistent real    seaLevelPressure ( nCells Time ) 0 iro seaLevelPressure mesh - -

Modified: branches/ocean_projects/sea_level_pressure/src/core_ocean/mpas_ocn_diagnostics.F
===================================================================
--- branches/ocean_projects/sea_level_pressure/src/core_ocean/mpas_ocn_diagnostics.F        2013-02-20 18:47:45 UTC (rev 2492)
+++ branches/ocean_projects/sea_level_pressure/src/core_ocean/mpas_ocn_diagnostics.F        2013-02-20 20:32:54 UTC (rev 2493)
@@ -97,7 +97,7 @@
       real (kind=RKIND), dimension(:), allocatable:: pTop, div_hu
 
       real (kind=RKIND), dimension(:), pointer :: &amp;
-        bottomDepth, fVertex, dvEdge, dcEdge, areaCell, areaTriangle, ssh
+        bottomDepth, fVertex, dvEdge, dcEdge, areaCell, areaTriangle, ssh, seaLevelPressure
       real (kind=RKIND), dimension(:,:), pointer :: &amp;
         weightsOnEdge, kiteAreasOnVertex, h_edge, h, u, v, pressure,&amp;
         circulation, vorticity, ke, ke_edge, MontPot, wTop, zMid, &amp;
@@ -158,6 +158,8 @@
       maxLevelVertexBot =&gt; grid % maxLevelVertexBot % array
       kiteIndexOnCell =&gt; grid % kiteIndexOnCell % array
       verticesOnCell =&gt; grid % verticesOnCell % array
+
+      seaLevelPressure =&gt; grid % seaLevelPressure % array
                   
       nCells      = grid % nCells
       nEdges      = grid % nEdges
@@ -424,7 +426,7 @@
         do iCell=1,nCells
            ! pressure for generalized coordinates
            ! assume atmospheric pressure at the surface is zero for now.
-           pressure(1,iCell) = rho(1,iCell)*gravity &amp;
+           pressure(1,iCell) = seaLevelPressure(iCell) + rho(1,iCell)*gravity &amp;
               * 0.5*h(1,iCell)
 
            do k=2,maxLevelCell(iCell)

Modified: branches/ocean_projects/sea_level_pressure/src/core_ocean/mpas_ocn_mpas_core.F
===================================================================
--- branches/ocean_projects/sea_level_pressure/src/core_ocean/mpas_ocn_mpas_core.F        2013-02-20 18:47:45 UTC (rev 2492)
+++ branches/ocean_projects/sea_level_pressure/src/core_ocean/mpas_ocn_mpas_core.F        2013-02-20 20:32:54 UTC (rev 2493)
@@ -126,7 +126,7 @@
 
       if (.not.config_do_restart) call ocn_init_split_timestep(domain)
 
-      write (0,'(a,a10)') ' Vertical coordinate movement is: ',config_vert_coord_movement
+      write (0,'(a,a)') ' Vertical coordinate movement is: ',trim(config_vert_coord_movement)
 
       if (config_vert_coord_movement.ne.'isopycnal'.and. &amp;
           config_vert_coord_movement.ne.'fixed'.and. &amp;

</font>
</pre>