Opened 10 years ago
Closed 7 years ago
#1212 closed Bug (wontfix)
Bug in UBS (dynamic) & partial cells
Reported by: | rbourdal | Owned by: | gm |
---|---|---|---|
Priority: | normal | Milestone: | |
Component: | OCE | Version: | trunk |
Severity: | Keywords: | ||
Cc: |
Description
In dynadv_ubs, the scale factor are not taking into account in the laplacian computation. The laplacian is OK in full step, sigma coord. but no OK in partial cells.
lines:
DO jk = 1, jpkm1 ! Laplacian of the velocity !
! ! =========================== !
! ! horizontal volume fluxes
zfu(:,:,jk) = e2u(:,:) * fse3u(:,:,jk) * un(:,:,jk)
zfv(:,:,jk) = e1v(:,:) * fse3v(:,:,jk) * vn(:,:,jk)
!
DO jj = 2, jpjm1 ! laplacian
DO ji = fs_2, fs_jpim1 ! vector opt.
zlu_uu(ji,jj,jk,1) = ( ub (ji+1,jj,jk)-2.*ub (ji,jj,jk)+ub (ji-1,jj,jk) ) * umask(ji,jj,jk)
zlv_vv(ji,jj,jk,1) = ( vb (ji,jj+1,jk)-2.*vb (ji,jj,jk)+vb (ji,jj-1,jk) ) * vmask(ji,jj,jk)
zlu_uv(ji,jj,jk,1) = ( ub (ji,jj+1,jk)-2.*ub (ji,jj,jk)+ub (ji,jj-1,jk) ) * umask(ji,jj,jk)
zlv_vu(ji,jj,jk,1) = ( vb (ji+1,jj,jk)-2.*vb (ji,jj,jk)+vb (ji-1,jj,jk) ) * vmask(ji,jj,jk)
!
zlu_uu(ji,jj,jk,2) = ( zfu(ji+1,jj,jk)-2.*zfu(ji,jj,jk)+zfu(ji-1,jj,jk) ) * umask(ji,jj,jk)
zlv_vv(ji,jj,jk,2) = ( zfv(ji,jj+1,jk)-2.*zfv(ji,jj,jk)+zfv(ji,jj-1,jk) ) * vmask(ji,jj,jk)
zlu_uv(ji,jj,jk,2) = ( zfu(ji,jj+1,jk)-2.*zfu(ji,jj,jk)+zfu(ji,jj-1,jk) ) * umask(ji,jj,jk)
zlv_vu(ji,jj,jk,2) = ( zfv(ji+1,jj,jk)-2.*zfv(ji,jj,jk)+zfv(ji-1,jj,jk) ) * vmask(ji,jj,jk)
END DO
END DO
END DO
Romain
Commit History (0)
(No commits)
Change History (5)
comment:1 Changed 10 years ago by rbourdal
- Priority changed from minor to major
comment:2 Changed 10 years ago by rbourdal
- Version changed from v3.4 to v3.6
comment:3 Changed 10 years ago by clevy
- Owner changed from NEMO team to gm
comment:4 Changed 7 years ago by clevy
- Version changed from v3.6 to trunk
comment:5 Changed 7 years ago by clevy
- Resolution set to wontfix
- Status changed from new to closed
Suppressing tickets opened on the trunk since more than 2 years... cleaning started!