New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
#1212 (Bug in UBS (dynamic) & partial cells) – NEMO

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!

Note: See TracTickets for help on using tickets.