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

Changeset 15499


Ignore:
Timestamp:
2021-11-12T10:23:40+01:00 (6 months ago)
Author:
vancop
Message:

northfold and stress diagnostic issues in VP rheology

Location:
NEMO/branches/2020/SI3-03_VP_rheology/src/ICE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/SI3-03_VP_rheology/src/ICE/icedyn_rhg_evp.F90

    r13591 r15499  
    831831                
    832832               ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008) 
    833                zsig_I(ji,jj)    =   zsig1 * 0.5_wp                                ! 1st stress invariant, aka average normal stress, aka negative pressure 
    834                zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 )   ! 2nd  ''       '', aka maximum shear stress 
     833               zsig_I(ji,jj)    =   zsig1 * 0.5_wp                                       ! 1st stress invariant, aka average normal stress, aka negative pressure 
     834               zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 * zsig12 )   ! 2nd  ''       '', aka maximum shear stress 
    835835                
    836836            END DO 
     
    869869               ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008), T-point 
    870870               zsig_I(ji,jj)    =   zsig1 * 0.5_wp                                ! 1st stress invariant, aka average normal stress, aka negative pressure 
    871                zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 )     ! 2nd  ''       '', aka maximum shear stress 
     871               zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 * zsig12 )     ! 2nd  ''       '', aka maximum shear stress 
    872872 
    873873               ! Normalized  principal stresses (used to display the ellipse) 
  • NEMO/branches/2020/SI3-03_VP_rheology/src/ICE/icedyn_rhg_vp.F90

    r15472 r15499  
    883883                     END DO 
    884884                      
    885                      CALL lbc_lnk( 'icedyn_rhg_vp', zFU, 'U',  1. ) 
     885                     CALL lbc_lnk( 'icedyn_rhg_vp', zFU, 'U',  -1. ) ! Nov 5 bug test MV zFU is a vector 
    886886                      
    887887                     !--------------- 
     
    905905                     END DO 
    906906 
    907                      CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime, 'U',  1., zBU_prime, 'U', 1. ) 
     907                     CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime, 'U',  -1., zBU_prime, 'U',  1. ) ! Nov 5 bug test zFU and zBU 
     908                                                                                                      ! are vectors 
    908909  
    909910                     !----------------------------- 
     
    914915                     
    915916                        ! --- Backward sweep  
     917 
    916918                        ! last row  
    917919                        zfac = SIGN( 1._wp , zBU_prime(jpi-1,jj) ) * MAX( 0._wp , SIGN( 1._wp , ABS( zBU_prime(jpi-1,jj) ) - epsi20 ) ) 
    918920                        u_ice(jpi-1,jj)    = zfac * zFU_prime(jpi-1,jj) / MAX( ABS ( zBU_prime(jpi-1,jj) ) , epsi20 ) &  
    919921                                           &            * umask(jpi-1,jj,1) 
    920                         DO ji = jpi-2 , 2, -1 ! all other rows    !  ---> original backward loop 
     922 
     923                        DO ji = jpi - 2 , 2, -1 ! all other rows    !  ---> original backward loop 
    921924                           zfac = SIGN( 1._wp , zBU_prime(ji,jj) ) * MAX( 0._wp , SIGN( 1._wp , ABS( zBU_prime(ji,jj) ) - epsi20 ) ) 
    922925                           u_ice(ji,jj)    = zfac * ( zFU_prime(ji,jj) - zCU(ji,jj) * u_ice(ji+1,jj) ) * umask(ji,jj,1)   &  
     
    985988                     END DO 
    986989 
    987                      CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV, 'V',  1.) 
     990                     CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV, 'V',  -1.) ! MV Nov 5 bug test to change boundary condition to -1 
     991                                                                          ! since zFV is a vector!!! 
    988992                      
    989993                     !--------------- 
     
    10071011                     END DO 
    10081012 
    1009                      CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV_prime, 'V',  1., zBV_prime, 'V', 1. ) 
     1013                     CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV_prime, 'V',  -1., zBV_prime, 'V',  1. ) ! MV Nov 5 bug test 
    10101014                      
    10111015                     !----------------------------- 
     
    11391143 
    11401144      IF ( lwp ) WRITE(numout,*) ' We are out of outer loop ' 
    1141  
    1142       CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU  , 'U',  1., zFV  , 'V',  1. ) 
    1143       CALL lbc_lnk_multi( 'icedyn_rhg_vp', zBU_prime  , 'U',  1., zBV_prime  , 'V',  1. ) 
    1144       CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime  , 'U',  1., zFV_prime  , 'V',  1. ) 
    1145       CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCU_prime  , 'U',  1., zCV_prime  , 'V',  1. ) 
    11461145 
    11471146      CALL iom_put( 'zFU'           , zFU            ) ! MV DEBUG 
     
    13681367            DO ji = 2, jpi - 1 
    13691368               ! Stress invariants 
    1370                zsig_I(ji,jj)    =   zs1(ji,jj) * 0.5_wp                                        ! 1st invariant, aka average normal stress aka negative pressure 
    1371                zsig_II(ji,jj)   =   SQRT ( zs2(ji,jj) * zs2(ji,jj) * 0.25_wp + zs12(ji,jj) )  ! 2nd invariant, aka maximum shear stress                
     1369               zsig_I(ji,jj)    =   zs1(ji,jj) * 0.5_wp                                                   ! 1st invariant, aka average normal stress aka negative pressure 
     1370               zsig_II(ji,jj)   =   SQRT ( zs2(ji,jj) * zs2(ji,jj) * 0.25_wp + zs12(ji,jj)*zs12(ji,jj) )  ! 2nd invariant, aka maximum shear stress                
    13721371            END DO 
    13731372         END DO 
     
    14021401               zfac             =   zp_deltastar_t(ji,jj) 
    14031402               zsig1            =   zfac * ( pdivu_i(ji,jj) - zdeltastar_t(ji,jj) ) 
    1404                zsig1            = 0._wp !!! FUCKING DEBUG TEST !!! 
    14051403               zsig2            =   zfac * z1_ecc2 * zten_i(ji,jj) 
    14061404               zsig12           =   zfac * z1_ecc2 * pshear_i(ji,jj) 
     
    14081406               ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008), T-point 
    14091407               zsig_I(ji,jj)    =   zsig1 * 0.5_wp                              ! 1st invariant 
    1410                zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 )   ! 2nd invariant 
     1408               zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 * zsig12 )   ! 2nd invariant 
    14111409 
    14121410               ! Normalized  principal stresses (used to display the ellipse) 
Note: See TracChangeset for help on using the changeset viewer.