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 15672 – NEMO

Changeset 15672


Ignore:
Timestamp:
2022-01-26T17:03:17+01:00 (2 years ago)
Author:
edblockley
Message:

Bug fix to EVP invariant stress diagnostics

File:
1 edited

Legend:

Unmodified
Added
Removed

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

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