Opened 5 years ago
Closed 5 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: |
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
678 678 !! compute Rc number (as done in zdfric.F90) 679 679 !!gm better to do it like in the new zdfric.F90 i.e. avm weighted Ri computation 680 680 !!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) 682 682 ! ! shear of horizontal velocity 683 683 zdku = zcoef * ( un(ji-1,jj ,ikt ) + un(ji,jj,ikt ) & 684 684 & -un(ji-1,jj ,ikt+1) - un(ji,jj,ikt+1) )
Change History (1)
comment:1 Changed 5 years ago by mathiot
- Resolution set to fixed
- Status changed from new to closed
Note: See
TracTickets for help on using
tickets.
In 10536: