Changeset 15493

Timestamp:

2021-11-10T17:02:17+01:00 (3 years ago)

Author:

clem

Message:

4.0-HEAD: correct output for rheology

File:

• NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rhg_evp.F90 (modified) (2 diffs)

NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rhg_evp.F90

@@ -838 +838 @@
                
                ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008)
-               zsig_I (ji,jj)   =   zsig1 * 0.5_wp                                           ! 1st stress invariant, aka average normal stress, aka negative pressure
-               zsig_II(ji,jj)   =   SQRT ( MAX( 0._wp, zsig2 * zsig2 * 0.25_wp + zsig12 ) )  ! 2nd  ''       '', aka maximum shear stress
+               zsig_I (ji,jj)   =   zsig1 * 0.5_wp                                      ! 1st stress invariant, aka average normal stress, aka negative pressure
+               zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 * zsig12 )  ! 2nd  ''       ''    , aka maximum shear stress
                
             END DO
@@ -872 +872 @@
 !!$               
 !!$               ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008), T-point
-!!$               zsig_I(ji,jj)    =   zsig1 * 0.5_wp                                           ! 1st stress invariant, aka average normal stress, aka negative pressure
-!!$               zsig_II(ji,jj)   =   SQRT ( MAX( 0._wp, zsig2 * zsig2 * 0.25_wp + zsig12 ) )  ! 2nd  ''       '', aka maximum shear stress
+!!$               zsig_I(ji,jj)    =   zsig1 * 0.5_wp                                      ! 1st stress invariant, aka average normal stress, aka negative pressure
+!!$               zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 * zsig12 )  ! 2nd  ''       ''    , aka maximum shear stress
 !!$      
 !!$               ! Normalized  principal stresses (used to display the ellipse)