sealevel - using the cont. eq. to get the sealevel. (Source File: sealevel.F90)

INTERFACE:

    subroutine sealevel(loop)
DESCRIPTION:

Here, the sea surface elevation is iterated according to the vertically integrated continuity equation given in (57) on page [*].

When working with the option SLICE_MODEL, the elevations at $j=2$ are copied to $j=3$.

Now with consideration of fresh water fluxes (precipitation and evaporation). Positive for flux into the water. USES:

    use domain, only: imin,imax,jmin,jmax,az,H
 #if defined(SPHERICAL) || defined(CURVILINEAR)
    use domain, only : arcd1,dxv,dyu
 #else
    use domain, only : dx,dy,ard1
 #endif
    use m2d, only: dtm,sealevel_check
    use variables_2d, only: z,zo,U,V,fwf
    use getm_timers, only: tic, toc, TIM_SEALEVEL, TIM_SEALEVELH
    use halo_zones, only : update_2d_halo,wait_halo,z_TAG
 #ifdef USE_BREAKS
    use halo_zones, only : nprocs,set_flag,u_TAG,v_TAG
    use variables_2d, only: break_mask,break_stat
    use domain, only : min_depth,au,av
 #endif
  $ use omp_lib
    IMPLICIT NONE
INPUT PARAMETERS:
    integer, intent(in)                 :: loop
REVISION HISTORY:
    Original author(s): Hans Burchard & Karsten Bolding
LOCAL VARIABLES:
    integer                   :: i,j
 #ifdef USE_BREAKS
    integer                   :: n,break_flag,break_flags(nprocs)
 #endif