Changeset 11653


Ignore:
Timestamp:
2019-10-04T14:34:02+02:00 (12 months ago)
Author:
davestorkey
Message:

Fix for bug in Smagorinsky scheme #2318.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/LDF/ldfdyn.F90

    r11536 r11653  
    315315            DO jj = 1, jpj             ! Set local gridscale values 
    316316               DO ji = 1, jpi 
    317                   esqt(ji,jj) = ( e1e2t(ji,jj) / ( e1t(ji,jj) + e2t(ji,jj) ) )**2  
    318                   esqf(ji,jj) = ( e1e2f(ji,jj) / ( e1f(ji,jj) + e2f(ji,jj) ) )**2  
     317                  esqt(ji,jj) = ( 2._wp * e1e2t(ji,jj) / ( e1t(ji,jj) + e2t(ji,jj) ) )**2  
     318                  esqf(ji,jj) = ( 2._wp * e1e2f(ji,jj) / ( e1f(ji,jj) + e2f(ji,jj) ) )**2  
    319319               END DO 
    320320            END DO 
Note: See TracChangeset for help on using the changeset viewer.