#2178 closed Bug (fixed)

Wrong index used in shear calculation?

Reported by: jamesharle Owned by: systeam
Priority: low Milestone:
Component: SBC Version: trunk
Severity: minor Keywords: ISF
Cc:

Description

Context

Wrong vertical index of e3w_n used calculating zdku and zdkv?

Analysis

e3w at ik+1 should be used for calculation not at ik?

Fix

  • sbcisf.F90

     
    678678                  !! compute Rc number (as done in zdfric.F90) 
    679679!!gm better to do it like in the new zdfric.F90   i.e. avm weighted Ri computation 
    680680!!gm moreover, use Max(rn2,0) to take care of static instabilities.... 
    681                   zcoef = 0.5_wp / e3w_n(ji,jj,ikt) 
     681                  zcoef = 0.5_wp / e3w_n(ji,jj,ikt+1) 
    682682                  !                                            ! shear of horizontal velocity 
    683683                  zdku = zcoef * (  un(ji-1,jj  ,ikt  ) + un(ji,jj,ikt  )  & 
    684684                     &             -un(ji-1,jj  ,ikt+1) - un(ji,jj,ikt+1)  ) 

Commit History (1)

ChangesetAuthorTimeChangeLog
10536mathiot2019-01-16T20:21:09+01:00

Fix #2178

Change History (1)

comment:1 Changed 18 months ago by mathiot

  • Resolution set to fixed
  • Status changed from new to closed

In 10536:

Fix #2178

Note: See TracTickets for help on using tickets.