<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 :: &
- bottomDepth, fVertex, dvEdge, dcEdge, areaCell, areaTriangle, ssh
+ bottomDepth, fVertex, dvEdge, dcEdge, areaCell, areaTriangle, ssh, seaLevelPressure
real (kind=RKIND), dimension(:,:), pointer :: &
weightsOnEdge, kiteAreasOnVertex, h_edge, h, u, v, pressure,&
circulation, vorticity, ke, ke_edge, MontPot, wTop, zMid, &
@@ -158,6 +158,8 @@
maxLevelVertexBot => grid % maxLevelVertexBot % array
kiteIndexOnCell => grid % kiteIndexOnCell % array
verticesOnCell => grid % verticesOnCell % array
+
+ seaLevelPressure => 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 &
+ pressure(1,iCell) = seaLevelPressure(iCell) + rho(1,iCell)*gravity &
* 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. &
config_vert_coord_movement.ne.'fixed'.and. &
</font>
</pre>