Opened 3 years ago

Closed 3 years ago

#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: Branch review:
MP ready?: Task progress:



Wrong vertical index of e3w_n used calculating zdku and zdkv?


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


  • 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)


Fix #2178

Change History (1)

comment:1 Changed 3 years 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.