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)
Changeset | Author | Time | ChangeLog |
---|---|---|---|
12055 | jchanut | 2019-12-04T17:20:14+01:00 | |
12054 | jchanut | 2019-12-04T17:16:11+01:00 | |
11875 | jchanut | 2019-11-08T12:00:59+01:00 | Prevent from 0 vertical index, #2327 |
Change History (4)
comment:1 Changed 4 years ago by jchanut
comment:2 Changed 4 years ago by jchanut
In 12054:
comment:3 Changed 4 years ago by jchanut
In 12055:
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.
In 11875: