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

Changeset 11867


Ignore:
Timestamp:
2019-11-06T13:08:50+01:00 (5 years ago)
Author:
stefryn
Message:

add EAP option to namelist and changes in icedyn_rhg.F90, icedyn_rhg_eap.F90 is just a copy of icedyn_rhg_evp.F90

Location:
NEMO/branches/2019/dev_r11842_SI3-10_EAP
Files:
1 added
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11842_SI3-10_EAP/cfgs/SHARED/namelist_ice_ref

    r11731 r11867  
    9797      rn_relast     =   0.333         !     ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
    9898                                      !        advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 
     99   ln_rhg_EAP       = .false.          !  EVP rheology 
    99100/ 
    100101!------------------------------------------------------------------------------ 
  • NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/icedyn_rhg.F90

    r11536 r11867  
    1717   USE ice            ! sea-ice: variables 
    1818   USE icedyn_rhg_evp ! sea-ice: EVP rheology 
     19   USE icedyn_rhg_eap ! sea-ice: EAP rheology 
    1920   USE icectl         ! sea-ice: control prints 
    2021   ! 
     
    3334   !                                        ! associated indices: 
    3435   INTEGER, PARAMETER ::   np_rhgEVP = 1   ! EVP rheology 
    35 !! INTEGER, PARAMETER ::   np_rhgEAP = 2   ! EAP rheology 
     36  INTEGER, PARAMETER ::   np_rhgEAP = 2   ! EAP rheology 
    3637 
    3738   ! ** namelist (namrhg) ** 
    3839   LOGICAL ::   ln_rhg_EVP       ! EVP rheology 
     40   LOGICAL ::   ln_rhg_EAP       ! EVP rheology 
    3941   ! 
    4042   !! * Substitutions 
     
    8183         CALL ice_dyn_rhg_evp( kt, stress1_i, stress2_i, stress12_i, shear_i, divu_i, delta_i ) 
    8284         !          
     85         !                             !----------------------------! 
     86      CASE( np_rhgEAP )                ! Elasto-Anisotropic-Plastic ! 
     87         !                             !----------------------------! 
     88         CALL ice_dyn_rhg_eap( kt, stress1_i, stress2_i, stress12_i, shear_i, divu_i, delta_i ) 
    8389      END SELECT 
    8490      ! 
    85       IF( lrst_ice ) THEN                       !* write EVP fields in the restart file 
    86          IF( ln_rhg_EVP )   CALL rhg_evp_rst( 'WRITE', kt ) 
     91      IF( lrst_ice ) THEN                       
     92         IF( ln_rhg_EVP )   CALL rhg_evp_rst( 'WRITE', kt ) !* write EVP fields in the restart file 
     93         IF( ln_rhg_EAP )   CALL rhg_eap_rst( 'WRITE', kt ) !* write EAP fields in the restart file 
    8794      ENDIF 
    8895      ! 
     
    110117      INTEGER ::   ios, ioptio   ! Local integer output status for namelist read 
    111118      !! 
    112       NAMELIST/namdyn_rhg/  ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast 
     119      NAMELIST/namdyn_rhg/  ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast, ln_rhg_EAP 
    113120      !!------------------------------------------------------------------- 
    114121      ! 
     
    132139         WRITE(numout,*) '         number of iterations for subcycling               nn_nevp    = ', nn_nevp 
    133140         WRITE(numout,*) '         ratio of elastic timescale over ice time step     rn_relast  = ', rn_relast 
     141         WRITE(numout,*) '      rheology EAP (icedyn_rhg_eap)                        ln_rhg_EAP = ', ln_rhg_EAP 
    134142      ENDIF 
    135143      ! 
     
    137145      ioptio = 0  
    138146      IF( ln_rhg_EVP ) THEN   ;   ioptio = ioptio + 1   ;   nice_rhg = np_rhgEVP    ;   ENDIF 
    139 !!    IF( ln_rhg_EAP ) THEN   ;   ioptio = ioptio + 1   ;   nice_rhg = np_rhgEAP    ;   ENDIF 
     147      IF( ln_rhg_EAP ) THEN   ;   ioptio = ioptio + 1   ;   nice_rhg = np_rhgEAP    ;   ENDIF 
    140148      IF( ioptio /= 1 )   CALL ctl_stop( 'ice_dyn_rhg_init: choose one and only one ice rheology' ) 
    141149      ! 
    142150      IF( ln_rhg_EVP  )   CALL rhg_evp_rst( 'READ' )  !* read or initialize all required files 
     151      IF( ln_rhg_EAP  )   CALL rhg_eap_rst( 'READ' )  !* read or initialize all required files 
    143152      ! 
    144153   END SUBROUTINE ice_dyn_rhg_init 
Note: See TracChangeset for help on using the changeset viewer.