Changeset 13105 for NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/icedyn_rhg_eap.F90

Timestamp:

2020-06-12T14:42:37+02:00 (4 years ago)

Author:

stefryn

Message:

small changes to calls

File:

• NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/icedyn_rhg_eap.F90 (modified) (7 diffs)

NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/icedyn_rhg_eap.F90

@@ -67 +67 @@
 CONTAINS
 
-   SUBROUTINE ice_dyn_rhg_eap( kt, pstress1_i, pstress2_i, pstress12_i, pshear_i, pdivu_i, pdelta_i, paniso_11, paniso_12 )
+   SUBROUTINE ice_dyn_rhg_eap( kt, pstress1_i, pstress2_i, pstress12_i, pshear_i, pdivu_i, pdelta_i, paniso_11, paniso_12, prdg_conv )
       !!-------------------------------------------------------------------
       !!                 ***  SUBROUTINE ice_dyn_rhg_eap  ***
@@ -123 +123 @@
       REAL(wp), DIMENSION(:,:), INTENT(  out) ::   pshear_i  , pdivu_i   , pdelta_i      !
       REAL(wp), DIMENSION(:,:), INTENT(inout) ::   paniso_11 , paniso_12                 ! structure tensor components
+      REAL(wp), DIMENSION(:,:), INTENT(inout) ::   prdg_conv                             ! for ridging
       !!
       INTEGER ::   ji, jj       ! dummy loop indices
@@ -415 +416 @@
                   yield22(ji,jj) = 0.5_wp * (stressptmp - stressmtmp)
                   yield12(ji,jj) = stress12tmp(ji,jj)
-                  rdg_conv(ji,jj) = -min( alphar, 0._wp)    
+                  prdg_conv(ji,jj) = -min( alphar, 0._wp)    
                ENDIF
 
@@ -443 +444 @@
          END DO
          !CALL lbc_lnk( 'icedyn_rhg_eap', zp_delt, 'T', 1. )
-         CALL lbc_lnk_multi( 'icedyn_rhg_eap', stress12tmp, 'T', 1., rdg_conv, 'T', 1. )
+         CALL lbc_lnk( 'icedyn_rhg_eap', stress12tmp, 'T', 1. )
 
          DO jj = 1, jpjm1
@@ -717 +718 @@
       !                                                   ! ==================== !
       !
+      CALL lbc_lnk( 'icedyn_rhg_eap', prdg_conv, 'T', 1. )  ! only need this in ridging module after rheology completed
+      !
       !------------------------------------------------------------------------------!
       ! 4) Recompute delta, shear and div (inputs for mechanical redistribution) 
@@ -838 +841 @@
       ! --- yieldcurve --- !
       IF( iom_use('yield11') .OR. iom_use('yield12') .OR. iom_use('yield22')) THEN
+
+         CALL lbc_lnk_multi( 'icedyn_rhg_eap', yield11, 'T', 1., yield22, 'T', 1., yield12, 'T', 1. )
+
          CALL iom_put( 'yield11', yield11 * zmsk00 )
          CALL iom_put( 'yield22', yield22 * zmsk00 )
@@ -844 +850 @@
 
       ! --- anisotropy tensor --- !
-      IF( iom_use('aniso') ) THEN
+      IF( iom_use('aniso') ) THEN       
+         CALL lbc_lnk( 'icedyn_rhg_eap', aniso_11, 'T', 1. ) 
          CALL iom_put( 'aniso' , aniso_11 * zmsk00 )
       ENDIF