Changeset 9736


Ignore:
Timestamp:
2018-06-05T10:03:33+02:00 (2 years ago)
Author:
jchanut
Message:

Fix bugs in Griffies triads, #1928

File:
1 edited

Legend:

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

    r9598 r9736  
    458458                     zdks = 0._wp 
    459459                  ENDIF 
    460                   zdzrho_raw = ( - zalbet(ji,jj,jk) * zdkt + zbeta0*zdks ) / e3w_n(ji,jj,jk+kp) 
     460                  zdzrho_raw = ( - rab_b(ji,jj,jk+kp,jp_tem) * zdkt &  
     461                             &   + rab_b(ji,jj,jk+kp,jp_sal) * zdks & 
     462                             & ) / e3w_n(ji,jj,jk+kp)   
    461463                  zdzrho(ji,jj,jk,kp) = - MIN( - repsln , zdzrho_raw )    ! force zdzrho >= repsln 
    462464                 END DO 
     
    586588                     ENDIF 
    587589                     !                                      ! switching triad scheme  
    588                      zisw = (rn_sw_triad - 1._wp ) + rn_sw_triad    & 
     590                     zisw = (1._wp - rn_sw_triad ) + rn_sw_triad    & 
    589591                        &            * 2._wp * ABS( 0.5_wp - kp - ( 0.5_wp - ip ) * SIGN( 1._wp , zdxrho(ji+ip,jj,jk,1-ip) )  ) 
    590                      zjsw = (rn_sw_triad - 1._wp ) + rn_sw_triad    & 
     592                     zjsw = (1._wp - rn_sw_triad ) + rn_sw_triad    & 
    591593                        &            * 2._wp * ABS( 0.5_wp - kp - ( 0.5_wp - jp ) * SIGN( 1._wp , zdyrho(ji,jj+jp,jk,1-jp) )  ) 
    592594                     ! 
Note: See TracChangeset for help on using the changeset viewer.