Ignore:
Timestamp:
2019-11-08T17:27:43+01:00 (12 months ago)
Author:
stefryn
Message:

add initialisation and restart of sturcture tensor

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/icedyn_rhg_eap.F90

    r11867 r11877  
    5656CONTAINS 
    5757 
    58    SUBROUTINE ice_dyn_rhg_eap( kt, pstress1_i, pstress2_i, pstress12_i, pshear_i, pdivu_i, pdelta_i ) 
     58   SUBROUTINE ice_dyn_rhg_eap( kt, pstress1_i, pstress2_i, pstress12_i, pshear_i, pdivu_i, pdelta_i, aniso_11, aniso_12 ) 
    5959      !!------------------------------------------------------------------- 
    6060      !!                 ***  SUBROUTINE ice_dyn_rhg_eap  *** 
     
    111111      REAL(wp), DIMENSION(:,:), INTENT(inout) ::   pstress1_i, pstress2_i, pstress12_i   ! 
    112112      REAL(wp), DIMENSION(:,:), INTENT(  out) ::   pshear_i  , pdivu_i   , pdelta_i      ! 
     113      REAL(wp), DIMENSION(:,:), INTENT(inout) ::   aniso_11, aniso_12   ! structure tensor components 
    113114      !! 
    114115      INTEGER ::   ji, jj       ! dummy loop indices 
     
    867868      ! 
    868869      INTEGER  ::   iter            ! local integer 
    869       INTEGER  ::   id1, id2, id3   ! local integers 
     870      INTEGER  ::   id1, id2, id3, id4, id5   ! local integers 
    870871      !!---------------------------------------------------------------------- 
    871872      ! 
     
    877878            id2 = iom_varid( numrir, 'stress2_i' , ldstop = .FALSE. ) 
    878879            id3 = iom_varid( numrir, 'stress12_i', ldstop = .FALSE. ) 
    879             ! 
    880             IF( MIN( id1, id2, id3 ) > 0 ) THEN      ! fields exist 
     880            id4 = iom_varid( numrir, 'aniso_11' , ldstop = .FALSE. ) 
     881            id5 = iom_varid( numrir, 'aniso_12' , ldstop = .FALSE. ) 
     882            ! 
     883            IF( MIN( id1, id2, id3, id4, id5 ) > 0 ) THEN      ! fields exist 
    881884               CALL iom_get( numrir, jpdom_autoglo, 'stress1_i' , stress1_i  ) 
    882885               CALL iom_get( numrir, jpdom_autoglo, 'stress2_i' , stress2_i  ) 
    883886               CALL iom_get( numrir, jpdom_autoglo, 'stress12_i', stress12_i ) 
     887               CALL iom_get( numrir, jpdom_autoglo, 'aniso_11' , aniso_11 ) 
     888               CALL iom_get( numrir, jpdom_autoglo, 'aniso_12' , aniso_12 ) 
    884889            ELSE                                     ! start rheology from rest 
    885890               IF(lwp) WRITE(numout,*) 
     
    888893               stress2_i (:,:) = 0._wp 
    889894               stress12_i(:,:) = 0._wp 
     895               aniso_11 (:,:) = 0.5_wp 
     896               aniso_12 (:,:) = 0._wp 
    890897            ENDIF 
    891898         ELSE                                   !* Start from rest 
     
    895902            stress2_i (:,:) = 0._wp 
    896903            stress12_i(:,:) = 0._wp 
     904            aniso_11 (:,:) = 0.5_wp 
     905            aniso_12 (:,:) = 0._wp 
    897906         ENDIF 
    898907         ! 
     
    905914         CALL iom_rstput( iter, nitrst, numriw, 'stress2_i' , stress2_i  ) 
    906915         CALL iom_rstput( iter, nitrst, numriw, 'stress12_i', stress12_i ) 
     916         CALL iom_rstput( iter, nitrst, numriw, 'aniso_11' , aniso_11 ) 
     917         CALL iom_rstput( iter, nitrst, numriw, 'aniso_12' , aniso_12 ) 
    907918         ! 
    908919      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.