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 are copied to .
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 NONEINPUT PARAMETERS:
integer, intent(in) :: loopREVISION HISTORY:
Original author(s): Hans Burchard & Karsten BoldingLOCAL VARIABLES:
integer :: i,j #ifdef USE_BREAKS integer :: n,break_flag,break_flags(nprocs) #endif