New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
#2178 (Wrong index used in shear calculation?) – NEMO

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

     
    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 5 years ago by mathiot

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

In 10536:

Error: Failed to load processor CommitTicketReference
No macro or processor named 'CommitTicketReference' found
Note: See TracTickets for help on using tickets.