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

Timestamp:

2019-11-08T17:27:43+01:00 (4 years ago)

Author:

stefryn

Message:

add initialisation and restart of sturcture tensor

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

@@ -56 +56 @@
 CONTAINS
 
-   SUBROUTINE ice_dyn_rhg_eap( kt, pstress1_i, pstress2_i, pstress12_i, pshear_i, pdivu_i, pdelta_i )
+   SUBROUTINE ice_dyn_rhg_eap( kt, pstress1_i, pstress2_i, pstress12_i, pshear_i, pdivu_i, pdelta_i, aniso_11, aniso_12 )
       !!-------------------------------------------------------------------
       !!                 ***  SUBROUTINE ice_dyn_rhg_eap  ***
@@ -111 +111 @@
       REAL(wp), DIMENSION(:,:), INTENT(inout) ::   pstress1_i, pstress2_i, pstress12_i   !
       REAL(wp), DIMENSION(:,:), INTENT(  out) ::   pshear_i  , pdivu_i   , pdelta_i      !
+      REAL(wp), DIMENSION(:,:), INTENT(inout) ::   aniso_11, aniso_12   ! structure tensor components
       !!
       INTEGER ::   ji, jj       ! dummy loop indices
@@ -867 +868 @@
       !
       INTEGER  ::   iter            ! local integer
-      INTEGER  ::   id1, id2, id3   ! local integers
+      INTEGER  ::   id1, id2, id3, id4, id5   ! local integers
       !!----------------------------------------------------------------------
       !
@@ -877 +878 @@
             id2 = iom_varid( numrir, 'stress2_i' , ldstop = .FALSE. )
             id3 = iom_varid( numrir, 'stress12_i', ldstop = .FALSE. )
-            !
-            IF( MIN( id1, id2, id3 ) > 0 ) THEN      ! fields exist
+            id4 = iom_varid( numrir, 'aniso_11' , ldstop = .FALSE. )
+            id5 = iom_varid( numrir, 'aniso_12' , ldstop = .FALSE. )
+            !
+            IF( MIN( id1, id2, id3, id4, id5 ) > 0 ) THEN      ! fields exist
                CALL iom_get( numrir, jpdom_autoglo, 'stress1_i' , stress1_i  )
                CALL iom_get( numrir, jpdom_autoglo, 'stress2_i' , stress2_i  )
                CALL iom_get( numrir, jpdom_autoglo, 'stress12_i', stress12_i )
+               CALL iom_get( numrir, jpdom_autoglo, 'aniso_11' , aniso_11 )
+               CALL iom_get( numrir, jpdom_autoglo, 'aniso_12' , aniso_12 )
             ELSE                                     ! start rheology from rest
                IF(lwp) WRITE(numout,*)
@@ -888 +893 @@
                stress2_i (:,:) = 0._wp
                stress12_i(:,:) = 0._wp
+               aniso_11 (:,:) = 0.5_wp
+               aniso_12 (:,:) = 0._wp
             ENDIF
          ELSE                                   !* Start from rest
@@ -895 +902 @@
             stress2_i (:,:) = 0._wp
             stress12_i(:,:) = 0._wp
+            aniso_11 (:,:) = 0.5_wp
+            aniso_12 (:,:) = 0._wp
          ENDIF
          !
@@ -905 +914 @@
          CALL iom_rstput( iter, nitrst, numriw, 'stress2_i' , stress2_i  )
          CALL iom_rstput( iter, nitrst, numriw, 'stress12_i', stress12_i )
+         CALL iom_rstput( iter, nitrst, numriw, 'aniso_11' , aniso_11 )
+         CALL iom_rstput( iter, nitrst, numriw, 'aniso_12' , aniso_12 )
          !
       ENDIF