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 15493 for NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rhg_evp.F90 – NEMO

Ignore:
Timestamp:
2021-11-10T17:02:17+01:00 (8 months ago)
Author:
clem
Message:

4.0-HEAD: correct output for rheology

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rhg_evp.F90

    r13646 r15493  
    838838                
    839839               ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008) 
    840                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 
     840               zsig_I (ji,jj)   =   zsig1 * 0.5_wp                                      ! 1st stress invariant, aka average normal stress, aka negative pressure 
     841               zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 * zsig12 )  ! 2nd  ''       ''    , aka maximum shear stress 
    842842                
    843843            END DO 
     
    872872!!$                
    873873!!$               ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008), T-point 
    874 !!$               zsig_I(ji,jj)    =   zsig1 * 0.5_wp                                           ! 1st stress invariant, aka average normal stress, aka negative pressure 
    875 !!$               zsig_II(ji,jj)   =   SQRT ( MAX( 0._wp, zsig2 * zsig2 * 0.25_wp + zsig12 ) )  ! 2nd  ''       '', aka maximum shear stress 
     874!!$               zsig_I(ji,jj)    =   zsig1 * 0.5_wp                                      ! 1st stress invariant, aka average normal stress, aka negative pressure 
     875!!$               zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 * zsig12 )  ! 2nd  ''       ''    , aka maximum shear stress 
    876876!!$       
    877877!!$               ! Normalized  principal stresses (used to display the ellipse) 
Note: See TracChangeset for help on using the changeset viewer.