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.
#2327 (Out of bound vertical index in 4th order vertical advection) – NEMO

Opened 4 years ago

Closed 4 years ago

#2327 closed Defect (fixed)

Out of bound vertical index in 4th order vertical advection

Reported by: jchanut Owned by: systeam
Priority: low Milestone:
Component: TRA Version: trunk
Severity: minor Keywords:
Cc:

Description

Context

Compact vertical advection

Analysis

Bottom 2nd order condition gives rise to null vertical index if mbkt=1 (over land for example).
i.e. problem is around line 661 of traadv_fct.F90:

            ikb = mbkt(ji,jj)         !     -   above the last wet point
            zwrm(ji,jj,ikb) = 0.5_wp * ( pt_in(ji,jj,ikb-1) + pt_in(ji,jj,ikb) )

Recommendation

A trivial fix is to set:

ikb = MAX(mbkt(ji,jj), 2)

Commit History (3)

ChangesetAuthorTimeChangeLog
12055jchanut2019-12-04T17:20:14+01:00

#2327

12054jchanut2019-12-04T17:16:11+01:00

#2327

11875jchanut2019-11-08T12:00:59+01:00

Prevent from 0 vertical index, #2327

Change History (4)

comment:1 Changed 4 years ago by jchanut

In 11875:

Error: Failed to load processor CommitTicketReference
No macro or processor named 'CommitTicketReference' found

comment:2 Changed 4 years ago by jchanut

In 12054:

Error: Failed to load processor CommitTicketReference
No macro or processor named 'CommitTicketReference' found

comment:3 Changed 4 years ago by jchanut

In 12055:

Error: Failed to load processor CommitTicketReference
No macro or processor named 'CommitTicketReference' found

comment:4 Changed 4 years ago by jchanut

  • Resolution set to fixed
  • Status changed from new to closed
Note: See TracTickets for help on using tickets.