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.
Changeset 15679 – NEMO

Changeset 15679


Ignore:
Timestamp:
2022-01-31T11:22:34+01:00 (2 years ago)
Author:
edblockley
Message:

Porting rheology invariant fix across to EAP

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/UKMO/NEMO_4.0.4_EAP_rheology_fix/src/ICE/icedyn_rhg_eap.F90

    r15571 r15679  
    905905               ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008) 
    906906               zsig_I (ji,jj)   =   zsig1 * 0.5_wp                                           ! 1st stress invariant, aka average normal stress, aka negative pressure 
    907                zsig_II(ji,jj)   =   SQRT ( MAX( 0._wp, zsig2 * zsig2 * 0.25_wp + zsig12 ) )  ! 2nd  ''       '', aka maximum shear stress 
     907               zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 * zsig12 )       ! 2nd  ''       '', aka maximum shear stress 
    908908                
    909909            END DO 
Note: See TracChangeset for help on using the changeset viewer.