Changeset 10481
- Timestamp:
- 2019-01-09T17:31:56+01:00 (6 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 4 edited
-
doc/latex/NEMO/subfiles/chap_LBC.tex (modified) (1 diff)
-
src/OCE/BDY/bdyini.F90 (modified) (4 diffs)
-
src/OCE/BDY/bdyvol.F90 (modified) (7 diffs)
-
src/OCE/DYN/dynspg_ts.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/doc/latex/NEMO/subfiles/chap_LBC.tex
r10442 r10481 613 613 \label{subsec:BDY_vol_corr} 614 614 615 There is an option to force the total volume in the regional model to be constant ,616 similar to the option in the OBC module.617 This is controlled by the \np{nn\_volctl} parameter in the namelist.618 A value of \np{nn\_volctl}\forcode{ = 0} indicates that this option is not used. 619 If \np{nn\_volctl}\forcode{ = 1} then a correction is applied to the normalvelocities around the boundary at615 There is an option to force the total volume in the regional model to be constant. 616 This is controlled by the \np{ln\_vol} parameter in the namelist. 617 A value of \np{ln\_vol}\forcode{ = .false.} indicates that this option is not used. 618 Two options to control the volume are available (\np{nn\_volctl}). 619 If \np{nn\_volctl}\forcode{ = 0} then a correction is applied to the normal barotropic velocities around the boundary at 620 620 each timestep to ensure that the integrated volume flow through the boundary is zero. 621 If \np{nn\_volctl}\forcode{ = 2} then the calculation of the volume change on621 If \np{nn\_volctl}\forcode{ = 1} then the calculation of the volume change on 622 622 the timestep includes the change due to the freshwater flux across the surface and 623 623 the correction velocity corrects for this as well. -
NEMO/trunk/src/OCE/BDY/bdyini.F90
r10425 r10481 133 133 INTEGER :: i_offset, j_offset ! - - 134 134 INTEGER , POINTER :: nbi, nbj, nbr ! short cuts 135 REAL(wp), POINTER :: flagu, flagv ! - -136 135 REAL(wp), POINTER, DIMENSION(:,:) :: pmask ! pointer to 2D mask fields 137 136 REAL(wp) :: zefl, zwfl, znfl, zsfl ! local scalars … … 385 384 END SELECT 386 385 IF(lwp) WRITE(numout,*) 386 ! 387 ! sanity check if used with tides 388 IF( ln_tide ) THEN 389 IF(lwp) WRITE(numout,*) ' The total volume correction is not working with tides. ' 390 IF(lwp) WRITE(numout,*) ' Set ln_vol to .FALSE. ' 391 IF(lwp) WRITE(numout,*) ' or ' 392 IF(lwp) WRITE(numout,*) ' equilibriate your bdy input files ' 393 CALL ctl_stop( 'The total volume correction is not working with tides.' ) 394 END IF 387 395 ELSE 388 396 IF(lwp) WRITE(numout,*) 'No volume correction applied at open boundaries' … … 1244 1252 ! 1245 1253 END DO 1246 1247 ! Compute total lateral surface for volume correction:1248 ! ----------------------------------------------------1249 ! JC: this must be done at each time step with non-linear free surface1250 bdysurftot = 0._wp1251 IF( ln_vol ) THEN1252 igrd = 2 ! Lateral surface at U-points1253 DO ib_bdy = 1, nb_bdy1254 DO ib = 1, idx_bdy(ib_bdy)%nblenrim(igrd)1255 nbi => idx_bdy(ib_bdy)%nbi(ib,igrd)1256 nbj => idx_bdy(ib_bdy)%nbj(ib,igrd)1257 flagu => idx_bdy(ib_bdy)%flagu(ib,igrd)1258 bdysurftot = bdysurftot + hu_n (nbi , nbj) &1259 & * e2u (nbi , nbj) * ABS( flagu ) &1260 & * tmask_i(nbi , nbj) &1261 & * tmask_i(nbi+1, nbj)1262 END DO1263 END DO1264 1265 igrd=3 ! Add lateral surface at V-points1266 DO ib_bdy = 1, nb_bdy1267 DO ib = 1, idx_bdy(ib_bdy)%nblenrim(igrd)1268 nbi => idx_bdy(ib_bdy)%nbi(ib,igrd)1269 nbj => idx_bdy(ib_bdy)%nbj(ib,igrd)1270 flagv => idx_bdy(ib_bdy)%flagv(ib,igrd)1271 bdysurftot = bdysurftot + hv_n (nbi, nbj ) &1272 & * e1v (nbi, nbj ) * ABS( flagv ) &1273 & * tmask_i(nbi, nbj ) &1274 & * tmask_i(nbi, nbj+1)1275 END DO1276 END DO1277 !1278 CALL mpp_sum( 'bdyini', bdysurftot ) ! sum over the global domain1279 END IF1280 1254 ! 1281 1255 ! Tidy up … … 1284 1258 ! 1285 1259 END SUBROUTINE bdy_segs 1286 1287 1260 1288 1261 SUBROUTINE bdy_ctl_seg -
NEMO/trunk/src/OCE/BDY/bdyvol.F90
r10425 r10481 9 9 !! 3.0 ! 2008-04 (NEMO team) add in the reference version 10 10 !! 3.4 ! 2011 (D. Storkey) rewrite in preparation for OBC-BDY merge 11 !! 4.0 ! 2019-01 (P. Mathiot) adapted to time splitting 11 12 !!---------------------------------------------------------------------- 12 13 USE oce ! ocean dynamics and tracers … … 24 25 PRIVATE 25 26 26 PUBLIC bdy_vol ! called by ???27 PUBLIC bdy_vol2d ! called by dynspg_ts 27 28 28 29 !!---------------------------------------------------------------------- … … 33 34 CONTAINS 34 35 35 SUBROUTINE bdy_vol ( kt)36 SUBROUTINE bdy_vol2d( kt, kc, pua2d, pva2d, phu, phv ) 36 37 !!---------------------------------------------------------------------- 37 38 !! *** ROUTINE bdyvol *** … … 40 41 !! A correction velocity is calculated to correct the total transport 41 42 !! through the unstructured OBC. 42 !! The total depth used is constant (H0) to be consistent with the43 !! linear free surface coded in OPA 8.2 <<<=== !!gm ???? true ????44 43 !! 45 44 !! ** Method : The correction velocity (zubtpecor here) is defined calculating … … 65 64 !! (set nn_volctl to 1 in tne namelist for this option) 66 65 !!---------------------------------------------------------------------- 67 INTEGER, INTENT(in) :: kt ! ocean time-step index66 INTEGER, INTENT(in) :: kt, kc ! ocean time-step index, cycle time-step 68 67 ! 69 68 INTEGER :: ji, jj, jk, jb, jgrd 70 69 INTEGER :: ib_bdy, ii, ij 71 REAL(wp) :: zubtpecor, z_cflxemp, ztranst 70 REAL(wp) :: zubtpecor, ztranst 71 REAL(wp), SAVE :: z_cflxemp ! cumulated emp flux 72 REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: pua2d, pva2d ! Barotropic velocities 73 REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: phu, phv ! Ocean depth at U- and V-points 72 74 TYPE(OBC_INDEX), POINTER :: idx 73 75 !!----------------------------------------------------------------------------- 74 76 ! 75 IF( ln_vol ) THEN76 !77 IF( kt == nit000 ) THEN78 IF(lwp) WRITE(numout,*)79 IF(lwp) WRITE(numout,*)'bdy_vol : Correction of velocities along unstructured OBC'80 IF(lwp) WRITE(numout,*)'~~~~~~~'81 END IF82 83 77 ! Calculate the cumulate surface Flux z_cflxemp (m3/s) over all the domain 84 78 ! ----------------------------------------------------------------------- 85 !!gm replace these lines : 86 z_cflxemp = SUM ( ( emp(:,:) - rnf(:,:) + fwfisf(:,:) ) * bdytmask(:,:) * e1e2t(:,:) ) / rau0 87 CALL mpp_sum( 'bdyvol', z_cflxemp ) ! sum over the global domain 88 !!gm by : 89 !!gm z_cflxemp = glob_sum( 'bdyvol', ( emp(:,:)-rnf(:,:)+fwfisf(:,:) ) * bdytmask(:,:) * e1e2t(:,:) ) / rau0 90 !!gm 79 IF ( kc == 1 ) z_cflxemp = glob_sum( 'bdyvol', ( emp(:,:) - rnf(:,:) + fwfisf(:,:) ) * bdytmask(:,:) * e1e2t(:,:) ) / rau0 80 81 ! Compute bdy surface each cycle if non linear free surface 82 ! --------------------------------------------------------- 83 IF ( .NOT. ln_linssh ) THEN 84 ! compute area each time step 85 bdysurftot = bdy_segs_surf( phu, phv ) 86 ELSE 87 ! compute area only the first time 88 IF ( ( kt == nit000 ) .AND. ( kc == 1 ) ) bdysurftot = bdy_segs_surf( phu, phv ) 89 END IF 91 90 92 91 ! Transport through the unstructured open boundary … … 98 97 jgrd = 2 ! cumulate u component contribution first 99 98 DO jb = 1, idx%nblenrim(jgrd) 100 DO jk = 1, jpkm1 99 ii = idx%nbi(jb,jgrd) 100 ij = idx%nbj(jb,jgrd) 101 zubtpecor = zubtpecor + idx%flagu(jb,jgrd) * pua2d(ii,ij) * e2u(ii,ij) * phu(ii,ij) * tmask_i(ii,ij) * tmask_i(ii+1,ij) 102 END DO 103 jgrd = 3 ! then add v component contribution 104 DO jb = 1, idx%nblenrim(jgrd) 105 ii = idx%nbi(jb,jgrd) 106 ij = idx%nbj(jb,jgrd) 107 zubtpecor = zubtpecor + idx%flagv(jb,jgrd) * pva2d(ii,ij) * e1v(ii,ij) * phv(ii,ij) * tmask_i(ii,ij) * tmask_i(ii,ij+1) 108 END DO 109 ! 110 END DO 111 IF( lk_mpp ) CALL mpp_sum( 'bdyvol', zubtpecor ) ! sum over the global domain 112 113 ! The normal velocity correction 114 ! ------------------------------ 115 IF( nn_volctl==1 ) THEN ; zubtpecor = ( zubtpecor - z_cflxemp ) / bdysurftot ! maybe should be apply only once at the end 116 ELSE ; zubtpecor = zubtpecor / bdysurftot 117 END IF 118 119 ! Correction of the total velocity on the unstructured boundary to respect the mass flux conservation 120 ! ------------------------------------------------------------- 121 DO ib_bdy = 1, nb_bdy 122 idx => idx_bdy(ib_bdy) 123 ! 124 jgrd = 2 ! correct u component 125 DO jb = 1, idx%nblenrim(jgrd) 101 126 ii = idx%nbi(jb,jgrd) 102 127 ij = idx%nbj(jb,jgrd) 103 zubtpecor = zubtpecor + idx%flagu(jb,jgrd) * ua(ii,ij, jk) * e2u(ii,ij) * e3u_n(ii,ij,jk) 104 END DO 105 END DO 106 jgrd = 3 ! then add v component contribution 107 DO jb = 1, idx%nblenrim(jgrd) 108 DO jk = 1, jpkm1 128 pua2d(ii,ij) = pua2d(ii,ij) - idx%flagu(jb,jgrd) * zubtpecor * tmask_i(ii,ij) * tmask_i(ii+1,ij) 129 END DO 130 jgrd = 3 ! correct v component 131 DO jb = 1, idx%nblenrim(jgrd) 109 132 ii = idx%nbi(jb,jgrd) 110 133 ij = idx%nbj(jb,jgrd) 111 zubtpecor = zubtpecor + idx%flagv(jb,jgrd) * va(ii,ij, jk) * e1v(ii,ij) * e3v_n(ii,ij,jk) 112 END DO 113 END DO 114 ! 115 END DO 116 CALL mpp_sum( 'bdyvol', zubtpecor ) ! sum over the global domain 117 118 ! The normal velocity correction 119 ! ------------------------------ 120 IF( nn_volctl==1 ) THEN ; zubtpecor = ( zubtpecor - z_cflxemp ) / bdysurftot 121 ELSE ; zubtpecor = zubtpecor / bdysurftot 122 END IF 123 124 ! Correction of the total velocity on the unstructured boundary to respect the mass flux conservation 125 ! ------------------------------------------------------------- 126 ztranst = 0._wp 127 DO ib_bdy = 1, nb_bdy 128 idx => idx_bdy(ib_bdy) 129 ! 130 jgrd = 2 ! correct u component 131 DO jb = 1, idx%nblenrim(jgrd) 132 DO jk = 1, jpkm1 133 ii = idx%nbi(jb,jgrd) 134 ij = idx%nbj(jb,jgrd) 135 ua(ii,ij,jk) = ua(ii,ij,jk) - idx%flagu(jb,jgrd) * zubtpecor * umask(ii,ij,jk) 136 ztranst = ztranst + idx%flagu(jb,jgrd) * ua(ii,ij,jk) * e2u(ii,ij) * e3u_n(ii,ij,jk) 137 END DO 138 END DO 139 jgrd = 3 ! correct v component 140 DO jb = 1, idx%nblenrim(jgrd) 141 DO jk = 1, jpkm1 142 ii = idx%nbi(jb,jgrd) 143 ij = idx%nbj(jb,jgrd) 144 va(ii,ij,jk) = va(ii,ij,jk) -idx%flagv(jb,jgrd) * zubtpecor * vmask(ii,ij,jk) 145 ztranst = ztranst + idx%flagv(jb,jgrd) * va(ii,ij,jk) * e1v(ii,ij) * e3v_n(ii,ij,jk) 146 END DO 147 END DO 148 ! 149 END DO 150 CALL mpp_sum( 'bdyvol', ztranst ) ! sum over the global domain 151 134 pva2d(ii,ij) = pva2d(ii,ij) - idx%flagv(jb,jgrd) * zubtpecor * tmask_i(ii,ij) * tmask_i(ii,ij+1) 135 END DO 136 ! 137 END DO 138 ! 152 139 ! Check the cumulated transport through unstructured OBC once barotropic velocities corrected 153 140 ! ------------------------------------------------------ 154 IF( lwp .AND. MOD( kt, nwrite ) == 0 ) THEN 141 IF( MOD( kt, nwrite ) == 0 .AND. ( kc == 1 ) ) THEN 142 ! 143 ! compute residual transport across boundary 144 ztranst = 0._wp 145 DO ib_bdy = 1, nb_bdy 146 idx => idx_bdy(ib_bdy) 147 ! 148 jgrd = 2 ! correct u component 149 DO jb = 1, idx%nblenrim(jgrd) 150 ii = idx%nbi(jb,jgrd) 151 ij = idx%nbj(jb,jgrd) 152 ztranst = ztranst + idx%flagu(jb,jgrd) * pua2d(ii,ij) * e2u(ii,ij) * phu(ii,ij) * tmask_i(ii,ij) * tmask_i(ii+1,ij) 153 END DO 154 jgrd = 3 ! correct v component 155 DO jb = 1, idx%nblenrim(jgrd) 156 ii = idx%nbi(jb,jgrd) 157 ij = idx%nbj(jb,jgrd) 158 ztranst = ztranst + idx%flagv(jb,jgrd) * pva2d(ii,ij) * e1v(ii,ij) * phv(ii,ij) * tmask_i(ii,ij) * tmask_i(ii,ij+1) 159 END DO 160 ! 161 END DO 162 IF( lk_mpp ) CALL mpp_sum('bdyvol', ztranst ) ! sum over the global domain 163 164 155 165 IF(lwp) WRITE(numout,*) 156 166 IF(lwp) WRITE(numout,*)'bdy_vol : time step :', kt … … 162 172 END IF 163 173 ! 164 END IF ! ln_vol 165 ! 166 END SUBROUTINE bdy_vol 167 174 END SUBROUTINE bdy_vol2d 175 ! 176 REAL(wp) FUNCTION bdy_segs_surf(phu, phv) 177 !!---------------------------------------------------------------------- 178 !! *** ROUTINE bdy_ctl_seg *** 179 !! 180 !! ** Purpose : Compute total lateral surface for volume correction 181 !! 182 !!---------------------------------------------------------------------- 183 184 REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: phu, phv ! water column thickness at U and V points 185 INTEGER :: igrd, ib_bdy, ib ! loop indexes 186 INTEGER , POINTER :: nbi, nbj ! short cuts 187 REAL(wp), POINTER :: zflagu, zflagv ! - - 188 189 ! Compute total lateral surface for volume correction: 190 ! ---------------------------------------------------- 191 bdy_segs_surf = 0._wp 192 igrd = 2 ! Lateral surface at U-points 193 DO ib_bdy = 1, nb_bdy 194 DO ib = 1, idx_bdy(ib_bdy)%nblenrim(igrd) 195 nbi => idx_bdy(ib_bdy)%nbi(ib,igrd) 196 nbj => idx_bdy(ib_bdy)%nbj(ib,igrd) 197 zflagu => idx_bdy(ib_bdy)%flagu(ib,igrd) 198 bdy_segs_surf = bdy_segs_surf + phu(nbi, nbj) & 199 & * e2u(nbi, nbj) * ABS( zflagu ) & 200 & * tmask_i(nbi, nbj) * tmask_i(nbi+1, nbj) 201 END DO 202 END DO 203 204 igrd=3 ! Add lateral surface at V-points 205 DO ib_bdy = 1, nb_bdy 206 DO ib = 1, idx_bdy(ib_bdy)%nblenrim(igrd) 207 nbi => idx_bdy(ib_bdy)%nbi(ib,igrd) 208 nbj => idx_bdy(ib_bdy)%nbj(ib,igrd) 209 zflagv => idx_bdy(ib_bdy)%flagv(ib,igrd) 210 bdy_segs_surf = bdy_segs_surf + phv(nbi, nbj) & 211 & * e1v(nbi, nbj) * ABS( zflagv ) & 212 & * tmask_i(nbi, nbj) * tmask_i(nbi, nbj+1) 213 END DO 214 END DO 215 ! 216 ! redirect the time to bdyvol as this variable is only used by bdyvol 217 IF( lk_mpp ) CALL mpp_sum( 'bdyvol', bdy_segs_surf ) ! sum over the global domain 218 ! 219 END FUNCTION bdy_segs_surf 168 220 !!====================================================================== 169 221 END MODULE bdyvol -
NEMO/trunk/src/OCE/DYN/dynspg_ts.F90
r10425 r10481 41 41 USE wet_dry ! wetting/drying flux limter 42 42 USE bdy_oce ! open boundary 43 USE bdyvol ! open boundary volume conservation 43 44 USE bdytides ! open boundary condition data 44 45 USE bdydyn2d ! open boundary conditions on barotropic variables … … 792 793 ! !* after ssh 793 794 ! ! ----------- 794 ! One should enforce volume conservation at open boundaries here 795 ! considering fluxes below: 795 ! 796 ! Enforce volume conservation at open boundaries: 797 IF( ln_bdy .AND. ln_vol ) CALL bdy_vol2d( kt, jn, ua_e, va_e, zhup2_e, zhvp2_e ) 796 798 ! 797 799 zwx(:,:) = e2u(:,:) * ua_e(:,:) * zhup2_e(:,:) ! fluxes at jn+0.5
Note: See TracChangeset
for help on using the changeset viewer.
