Changeset 5951 for branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90
- Timestamp:
- 2015-11-30T12:48:01+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90
r5950 r5951 14 14 !! - ! 2013-04 (F. Roquet, G. Madec) use of eosbn2 instead of local hard coded alpha and beta 15 15 !!---------------------------------------------------------------------- 16 #if defined key_trabbl || defined key_esopa16 #if defined key_trabbl 17 17 !!---------------------------------------------------------------------- 18 18 !! 'key_trabbl' or bottom boundary layer … … 29 29 USE phycst ! physical constant 30 30 USE eosbn2 ! equation of state 31 USE trd_oce ! trends: ocean variables31 USE trd_oce ! trends: ocean variables 32 32 USE trdtra ! trends: active tracers 33 33 ! … … 198 198 DO jj = 1, jpj 199 199 DO ji = 1, jpi 200 ik = mbkt(ji,jj) 201 zptb(ji,jj) = ptb(ji,jj,ik,jn) ! bottom before T and S200 ik = mbkt(ji,jj) ! bottom T-level index 201 zptb(ji,jj) = ptb(ji,jj,ik,jn) ! bottom before T and S 202 202 END DO 203 203 END DO … … 205 205 DO jj = 2, jpjm1 ! Compute the trend 206 206 DO ji = 2, jpim1 207 ik = mbkt(ji,jj) 208 zbtr = r1_e12t(ji,jj) / fse3t(ji,jj,ik)209 pta(ji,jj,ik,jn) = pta(ji,jj,ik,jn)&210 & + ( ahu_bbl(ji ,jj ) * ( zptb(ji+1,jj ) - zptb(ji ,jj ) )&211 & - ahu_bbl(ji-1,jj ) * ( zptb(ji ,jj ) - zptb(ji-1,jj ) )&212 & + ahv_bbl(ji ,jj ) * ( zptb(ji ,jj+1) - zptb(ji ,jj ) )&213 & - ahv_bbl(ji ,jj-1) * ( zptb(ji ,jj ) - zptb(ji ,jj-1) ) ) * zbtr207 ik = mbkt(ji,jj) ! bottom T-level index 208 pta(ji,jj,ik,jn) = pta(ji,jj,ik,jn) & 209 & + ( ahu_bbl(ji ,jj ) * ( zptb(ji+1,jj ) - zptb(ji ,jj ) ) & 210 & - ahu_bbl(ji-1,jj ) * ( zptb(ji ,jj ) - zptb(ji-1,jj ) ) & 211 & + ahv_bbl(ji ,jj ) * ( zptb(ji ,jj+1) - zptb(ji ,jj ) ) & 212 & - ahv_bbl(ji ,jj-1) * ( zptb(ji ,jj ) - zptb(ji ,jj-1) ) ) & 213 & / ( e1e2t(ji,jj) * fse3t(ji,jj,ik) ) 214 214 END DO 215 215 END DO … … 263 263 ! 264 264 ! ! up -slope T-point (shelf bottom point) 265 zbtr = r1_e1 2t(iis,jj) / fse3t(iis,jj,ikus)265 zbtr = r1_e1e2t(iis,jj) / fse3t(iis,jj,ikus) 266 266 ztra = zu_bbl * ( ptb(iid,jj,ikus,jn) - ptb(iis,jj,ikus,jn) ) * zbtr 267 267 pta(iis,jj,ikus,jn) = pta(iis,jj,ikus,jn) + ztra 268 268 ! 269 269 DO jk = ikus, ikud-1 ! down-slope upper to down T-point (deep column) 270 zbtr = r1_e1 2t(iid,jj) / fse3t(iid,jj,jk)270 zbtr = r1_e1e2t(iid,jj) / fse3t(iid,jj,jk) 271 271 ztra = zu_bbl * ( ptb(iid,jj,jk+1,jn) - ptb(iid,jj,jk,jn) ) * zbtr 272 272 pta(iid,jj,jk,jn) = pta(iid,jj,jk,jn) + ztra 273 273 END DO 274 274 ! 275 zbtr = r1_e1 2t(iid,jj) / fse3t(iid,jj,ikud)275 zbtr = r1_e1e2t(iid,jj) / fse3t(iid,jj,ikud) 276 276 ztra = zu_bbl * ( ptb(iis,jj,ikus,jn) - ptb(iid,jj,ikud,jn) ) * zbtr 277 277 pta(iid,jj,ikud,jn) = pta(iid,jj,ikud,jn) + ztra … … 285 285 ! 286 286 ! up -slope T-point (shelf bottom point) 287 zbtr = r1_e1 2t(ji,ijs) / fse3t(ji,ijs,ikvs)287 zbtr = r1_e1e2t(ji,ijs) / fse3t(ji,ijs,ikvs) 288 288 ztra = zv_bbl * ( ptb(ji,ijd,ikvs,jn) - ptb(ji,ijs,ikvs,jn) ) * zbtr 289 289 pta(ji,ijs,ikvs,jn) = pta(ji,ijs,ikvs,jn) + ztra 290 290 ! 291 291 DO jk = ikvs, ikvd-1 ! down-slope upper to down T-point (deep column) 292 zbtr = r1_e1 2t(ji,ijd) / fse3t(ji,ijd,jk)292 zbtr = r1_e1e2t(ji,ijd) / fse3t(ji,ijd,jk) 293 293 ztra = zv_bbl * ( ptb(ji,ijd,jk+1,jn) - ptb(ji,ijd,jk,jn) ) * zbtr 294 294 pta(ji,ijd,jk,jn) = pta(ji,ijd,jk,jn) + ztra 295 295 END DO 296 296 ! ! down-slope T-point (deep bottom point) 297 zbtr = r1_e1 2t(ji,ijd) / fse3t(ji,ijd,ikvd)297 zbtr = r1_e1e2t(ji,ijd) / fse3t(ji,ijd,ikvd) 298 298 ztra = zv_bbl * ( ptb(ji,ijs,ikvs,jn) - ptb(ji,ijd,ikvd,jn) ) * zbtr 299 299 pta(ji,ijd,ikvd,jn) = pta(ji,ijd,ikvd,jn) + ztra … … 566 566 567 567 ! !* masked diffusive flux coefficients 568 ahu_bbl_0(:,:) = rn_ahtbbl * e2 u(:,:) * e3u_bbl_0(:,:) / e1u(:,:)* umask(:,:,1)569 ahv_bbl_0(:,:) = rn_ahtbbl * e1 v(:,:) * e3v_bbl_0(:,:) / e2v(:,:)* vmask(:,:,1)568 ahu_bbl_0(:,:) = rn_ahtbbl * e2_e1u(:,:) * e3u_bbl_0(:,:) * umask(:,:,1) 569 ahv_bbl_0(:,:) = rn_ahtbbl * e1_e2v(:,:) * e3v_bbl_0(:,:) * vmask(:,:,1) 570 570 571 571
Note: See TracChangeset
for help on using the changeset viewer.