Changeset 15499 for NEMO/branches/2020/SI3-03_VP_rheology/src/ICE/icedyn_rhg_vp.F90

Timestamp:

2021-11-12T10:23:40+01:00 (3 years ago)

Author:

vancop

Message:

northfold and stress diagnostic issues in VP rheology

File:

• NEMO/branches/2020/SI3-03_VP_rheology/src/ICE/icedyn_rhg_vp.F90 (modified) (9 diffs)

NEMO/branches/2020/SI3-03_VP_rheology/src/ICE/icedyn_rhg_vp.F90

@@ -883 +883 @@
                      END DO
                      
-                     CALL lbc_lnk( 'icedyn_rhg_vp', zFU, 'U',  1. )
+                     CALL lbc_lnk( 'icedyn_rhg_vp', zFU, 'U',  -1. ) ! Nov 5 bug test MV zFU is a vector
                      
                      !---------------
@@ -905 +905 @@
                      END DO
 
-                     CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime, 'U',  1., zBU_prime, 'U', 1. )
+                     CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime, 'U',  -1., zBU_prime, 'U',  1. ) ! Nov 5 bug test zFU and zBU
+                                                                                                      ! are vectors
  
                      !-----------------------------
@@ -914 +915 @@
                     
                         ! --- Backward sweep 
+
                         ! last row 
                         zfac = SIGN( 1._wp , zBU_prime(jpi-1,jj) ) * MAX( 0._wp , SIGN( 1._wp , ABS( zBU_prime(jpi-1,jj) ) - epsi20 ) )
                         u_ice(jpi-1,jj)    = zfac * zFU_prime(jpi-1,jj) / MAX( ABS ( zBU_prime(jpi-1,jj) ) , epsi20 ) & 
                                            &            * umask(jpi-1,jj,1)
-                        DO ji = jpi-2 , 2, -1 ! all other rows    !  ---> original backward loop
+
+                        DO ji = jpi - 2 , 2, -1 ! all other rows    !  ---> original backward loop
                            zfac = SIGN( 1._wp , zBU_prime(ji,jj) ) * MAX( 0._wp , SIGN( 1._wp , ABS( zBU_prime(ji,jj) ) - epsi20 ) )
                            u_ice(ji,jj)    = zfac * ( zFU_prime(ji,jj) - zCU(ji,jj) * u_ice(ji+1,jj) ) * umask(ji,jj,1)   & 
@@ -985 +988 @@
                      END DO
 
-                     CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV, 'V',  1.)
+                     CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV, 'V',  -1.) ! MV Nov 5 bug test to change boundary condition to -1
+                                                                          ! since zFV is a vector!!!
                      
                      !---------------
@@ -1007 +1011 @@
                      END DO
 
-                     CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV_prime, 'V',  1., zBV_prime, 'V', 1. )
+                     CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV_prime, 'V',  -1., zBV_prime, 'V',  1. ) ! MV Nov 5 bug test
                      
                      !-----------------------------
@@ -1139 +1143 @@
 
       IF ( lwp ) WRITE(numout,*) ' We are out of outer loop '
-
-      CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU  , 'U',  1., zFV  , 'V',  1. )
-      CALL lbc_lnk_multi( 'icedyn_rhg_vp', zBU_prime  , 'U',  1., zBV_prime  , 'V',  1. )
-      CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime  , 'U',  1., zFV_prime  , 'V',  1. )
-      CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCU_prime  , 'U',  1., zCV_prime  , 'V',  1. )
 
       CALL iom_put( 'zFU'           , zFU            ) ! MV DEBUG
@@ -1368 +1367 @@
             DO ji = 2, jpi - 1
                ! Stress invariants
-               zsig_I(ji,jj)    =   zs1(ji,jj) * 0.5_wp                                        ! 1st invariant, aka average normal stress aka negative pressure
-               zsig_II(ji,jj)   =   SQRT ( zs2(ji,jj) * zs2(ji,jj) * 0.25_wp + zs12(ji,jj) )  ! 2nd invariant, aka maximum shear stress               
+               zsig_I(ji,jj)    =   zs1(ji,jj) * 0.5_wp                                                   ! 1st invariant, aka average normal stress aka negative pressure
+               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               
             END DO
          END DO
@@ -1402 +1401 @@
                zfac             =   zp_deltastar_t(ji,jj)
                zsig1            =   zfac * ( pdivu_i(ji,jj) - zdeltastar_t(ji,jj) )
-               zsig1            = 0._wp !!! FUCKING DEBUG TEST !!!
                zsig2            =   zfac * z1_ecc2 * zten_i(ji,jj)
                zsig12           =   zfac * z1_ecc2 * pshear_i(ji,jj)
@@ -1408 +1406 @@
                ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008), T-point
                zsig_I(ji,jj)    =   zsig1 * 0.5_wp                              ! 1st invariant
-               zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 )   ! 2nd invariant
+               zsig_II(ji,jj)   =   SQRT ( zsig2 * zsig2 * 0.25_wp + zsig12 * zsig12 )   ! 2nd invariant
 
                ! Normalized  principal stresses (used to display the ellipse)