Opened 3 months ago

Closed 8 weeks 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 3 months ago by jchanut

In 11875:

Prevent from 0 vertical index, #2327

comment:2 Changed 8 weeks ago by jchanut

In 12054:

comment:3 Changed 8 weeks ago by jchanut

In 12055:

comment:4 Changed 8 weeks ago by jchanut

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