- Timestamp:
- 2012-11-26T11:46:39+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_NOC_MERCATOR_2012/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn3d.F90
r3294 r3651 14 14 !!---------------------------------------------------------------------- 15 15 USE timing ! Timing 16 USE wrk_nemo ! Memory Allocation 16 17 USE oce ! ocean dynamics and tracers 17 18 USE dom_oce ! ocean space and time domain … … 19 20 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 20 21 USE in_out_manager ! 22 Use phycst 21 23 22 24 IMPLICIT NONE … … 24 26 25 27 PUBLIC bdy_dyn3d ! routine called by bdy_dyn 26 28 PUBLIC bdy_dyn3d_dmp ! routine called by step 29 30 !! * Substitutions 31 # include "domzgr_substitute.h90" 27 32 !!---------------------------------------------------------------------- 28 33 !! NEMO/OPA 3.3 , NEMO Consortium (2010) … … 55 60 CASE(jp_frs) 56 61 CALL bdy_dyn3d_frs( idx_bdy(ib_bdy), dta_bdy(ib_bdy), kt ) 62 CASE(2) 63 CALL bdy_dyn3d_spe( idx_bdy(ib_bdy), dta_bdy(ib_bdy), kt ) 64 CASE(3) 65 CALL bdy_dyn3d_zro( idx_bdy(ib_bdy), dta_bdy(ib_bdy), kt ) 57 66 CASE DEFAULT 58 67 CALL ctl_stop( 'bdy_dyn3d : unrecognised option for open boundaries for baroclinic velocities' ) … … 61 70 62 71 END SUBROUTINE bdy_dyn3d 72 73 SUBROUTINE bdy_dyn3d_spe( idx, dta, kt ) 74 !!---------------------------------------------------------------------- 75 !! *** SUBROUTINE bdy_dyn3d_spe *** 76 !! 77 !! ** Purpose : - Apply a specified value for baroclinic velocities 78 !! at open boundaries. 79 !! 80 !!---------------------------------------------------------------------- 81 INTEGER :: kt 82 TYPE(OBC_INDEX), INTENT(in) :: idx ! OBC indices 83 TYPE(OBC_DATA), INTENT(in) :: dta ! OBC external data 84 !! 85 INTEGER :: jb, jk ! dummy loop indices 86 INTEGER :: ii, ij, igrd ! local integers 87 REAL(wp) :: zwgt ! boundary weight 88 !!---------------------------------------------------------------------- 89 ! 90 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_spe') 91 ! 92 igrd = 2 ! Relaxation of zonal velocity 93 DO jb = 1, idx%nblenrim(igrd) 94 DO jk = 1, jpkm1 95 ii = idx%nbi(jb,igrd) 96 ij = idx%nbj(jb,igrd) 97 ua(ii,ij,jk) = dta%u3d(jb,jk) * umask(ii,ij,jk) 98 END DO 99 END DO 100 ! 101 igrd = 3 ! Relaxation of meridional velocity 102 DO jb = 1, idx%nblenrim(igrd) 103 DO jk = 1, jpkm1 104 ii = idx%nbi(jb,igrd) 105 ij = idx%nbj(jb,igrd) 106 va(ii,ij,jk) = dta%v3d(jb,jk) * vmask(ii,ij,jk) 107 END DO 108 END DO 109 CALL lbc_lnk( ua, 'U', -1. ) ; CALL lbc_lnk( va, 'V', -1. ) ! Boundary points should be updated 110 ! 111 IF( kt .eq. nit000 ) CLOSE( unit = 102 ) 112 113 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_spe') 114 115 END SUBROUTINE bdy_dyn3d_spe 116 117 SUBROUTINE bdy_dyn3d_zro( idx, dta, kt ) 118 !!---------------------------------------------------------------------- 119 !! *** SUBROUTINE bdy_dyn3d_zro *** 120 !! 121 !! ** Purpose : - baroclinic velocities = 0. at open boundaries. 122 !! 123 !!---------------------------------------------------------------------- 124 INTEGER :: kt 125 TYPE(OBC_INDEX), INTENT(in) :: idx ! OBC indices 126 TYPE(OBC_DATA), INTENT(in) :: dta ! OBC external data 127 !! 128 INTEGER :: ib, ik ! dummy loop indices 129 INTEGER :: ii, ij, igrd, zcoef ! local integers 130 REAL(wp) :: zwgt ! boundary weight 131 !!---------------------------------------------------------------------- 132 ! 133 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_zro') 134 ! 135 igrd = 2 ! Everything is at T-points here 136 DO ib = 1, idx%nblenrim(igrd) 137 ii = idx%nbi(ib,igrd) 138 ij = idx%nbj(ib,igrd) 139 DO ik = 1, jpkm1 140 ua(ii,ij,ik) = 0._wp 141 END DO 142 END DO 143 144 igrd = 3 ! Everything is at T-points here 145 DO ib = 1, idx%nblenrim(igrd) 146 ii = idx%nbi(ib,igrd) 147 ij = idx%nbj(ib,igrd) 148 DO ik = 1, jpkm1 149 va(ii,ij,ik) = 0._wp 150 END DO 151 END DO 152 ! 153 CALL lbc_lnk( ua, 'U', -1. ) ; CALL lbc_lnk( va, 'V', -1. ) ! Boundary points should be updated 154 ! 155 IF( kt .eq. nit000 ) CLOSE( unit = 102 ) 156 157 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_zro') 158 159 END SUBROUTINE bdy_dyn3d_zro 63 160 64 161 SUBROUTINE bdy_dyn3d_frs( idx, dta, kt ) … … 111 208 END SUBROUTINE bdy_dyn3d_frs 112 209 210 SUBROUTINE bdy_dyn3d_dmp( kt ) 211 !!---------------------------------------------------------------------- 212 !! *** SUBROUTINE bdy_dyn3d_dmp *** 213 !! 214 !! ** Purpose : Apply damping for baroclinic velocities at open boundaries. 215 !! 216 !!---------------------------------------------------------------------- 217 INTEGER :: kt 218 !! 219 INTEGER :: jb, jk ! dummy loop indices 220 INTEGER :: ii, ij, igrd ! local integers 221 REAL(wp) :: zwgt ! boundary weight 222 INTEGER :: ib_bdy ! loop index 223 !!---------------------------------------------------------------------- 224 ! 225 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_dmp') 226 ! 227 !------------------------------------------------------- 228 ! Remove barotropic part from before velocity 229 !------------------------------------------------------- 230 CALL wrk_alloc(jpi,jpj,pu2d,pv2d) 231 232 pu2d(:,:) = 0.e0 233 pv2d(:,:) = 0.e0 234 235 DO jk = 1, jpkm1 236 #if defined key_vvl 237 pu2d(:,:) = pu2d(:,:) + fse3u_b(:,:,jk)* ub(:,:,jk) *umask(:,:,jk) 238 pv2d(:,:) = pv2d(:,:) + fse3v_b(:,:,jk)* vb(:,:,jk) *vmask(:,:,jk) 239 #else 240 pu2d(:,:) = pu2d(:,:) + fse3u_0(:,:,jk) * ub(:,:,jk) * umask(:,:,jk) 241 pv2d(:,:) = pv2d(:,:) + fse3v_0(:,:,jk) * vb(:,:,jk) * vmask(:,:,jk) 242 #endif 243 END DO 244 245 IF( lk_vvl ) THEN 246 pu2d(:,:) = pu2d(:,:) * umask(:,:,1) / ( hu_0(:,:) + sshu_b(:,:) + 1._wp - umask(:,:,1) ) 247 pv2d(:,:) = pv2d(:,:) * vmask(:,:,1) / ( hv_0(:,:) + sshv_b(:,:) + 1._wp - vmask(:,:,1) ) 248 ELSE 249 pu2d(:,:) = pv2d(:,:) * hur(:,:) 250 pv2d(:,:) = pu2d(:,:) * hvr(:,:) 251 ENDIF 252 253 DO ib_bdy=1, nb_bdy 254 IF ( ln_dyn3d_dmp(ib_bdy).and.nn_dyn3d(ib_bdy).gt.0 ) THEN 255 igrd = 2 ! Relaxation of zonal velocity 256 DO jb = 1, idx_bdy(ib_bdy)%nblen(igrd) 257 ii = idx_bdy(ib_bdy)%nbi(jb,igrd) 258 ij = idx_bdy(ib_bdy)%nbj(jb,igrd) 259 zwgt = idx_bdy(ib_bdy)%nbd(jb,igrd) 260 DO jk = 1, jpkm1 261 ua(ii,ij,jk) = ( ua(ii,ij,jk) + zwgt * ( dta_bdy(ib_bdy)%u3d(jb,jk) - ub(ii,ij,jk) + pu2d(ii,ij)) ) * umask(ii,ij,jk) 262 END DO 263 END DO 264 ! 265 igrd = 3 ! Relaxation of meridional velocity 266 DO jb = 1, idx_bdy(ib_bdy)%nblen(igrd) 267 ii = idx_bdy(ib_bdy)%nbi(jb,igrd) 268 ij = idx_bdy(ib_bdy)%nbj(jb,igrd) 269 zwgt = idx_bdy(ib_bdy)%nbd(jb,igrd) 270 DO jk = 1, jpkm1 271 va(ii,ij,jk) = ( va(ii,ij,jk) + zwgt * ( dta_bdy(ib_bdy)%v3d(jb,jk) - vb(ii,ij,jk) + pv2d(ii,ij)) ) * vmask(ii,ij,jk) 272 END DO 273 END DO 274 ENDIF 275 ENDDO 276 ! 277 CALL wrk_dealloc(jpi,jpj,pu2d,pv2d) 278 ! 279 CALL lbc_lnk( ua, 'U', -1. ) ; CALL lbc_lnk( va, 'V', -1. ) ! Boundary points should be updated 280 ! 281 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_dmp') 282 283 END SUBROUTINE bdy_dyn3d_dmp 113 284 114 285 #else … … 118 289 CONTAINS 119 290 SUBROUTINE bdy_dyn3d( kt ) ! Empty routine 120 WRITE(*,*) 'bdy_dyn _frs: You should not have seen this print! error?', kt291 WRITE(*,*) 'bdy_dyn3d: You should not have seen this print! error?', kt 121 292 END SUBROUTINE bdy_dyn3d 293 294 SUBROUTINE bdy_dyn3d_dmp( kt ) ! Empty routine 295 WRITE(*,*) 'bdy_dyn3d_dmp: You should not have seen this print! error?', kt 296 END SUBROUTINE bdy_dyn3d_dmp 297 122 298 #endif 123 299
Note: See TracChangeset
for help on using the changeset viewer.