Changeset 13493


Ignore:
Timestamp:
2020-09-21T09:00:06+02:00 (5 months ago)
Author:
vancop
Message:

Some progress on VP rheology but still does not compile

Location:
NEMO/branches/2020/SI3-03_VP_rheology/src/ICE
Files:
1 added
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/SI3-03_VP_rheology/src/ICE/icedyn_rhg.F90

    r12741 r13493  
    3434   INTEGER, PARAMETER ::   np_rhgEVP = 1   ! EVP rheology 
    3535!! INTEGER, PARAMETER ::   np_rhgEAP = 2   ! EAP rheology 
     36   INTEGER, PARAMETER ::   np_rhgVP  = 3   ! VP rheology 
    3637 
    3738   ! ** namelist (namrhg) ** 
    3839   LOGICAL ::   ln_rhg_EVP       ! EVP rheology 
     40   LOGICAL ::   ln_rhg_VP        ! EVP rheology 
    3941   ! 
    4042   !! * Substitutions 
     
    7678      !--------------!    
    7779      SELECT CASE( nice_rhg ) 
    78       !                                !------------------------! 
    79       CASE( np_rhgEVP )                ! Elasto-Viscous-Plastic ! 
    80          !                             !------------------------! 
     80      !                                !---------------------------------! 
     81      CASE( np_rhgEVP )                ! Elasto-Viscous-Plastic rheology ! 
     82         !                             !---------------------------------! 
    8183         CALL ice_dyn_rhg_evp( kt, stress1_i, stress2_i, stress12_i, shear_i, divu_i, delta_i ) 
     84         !          
     85      END SELECT 
     86 
     87      !                                !---------------------------------------------! 
     88      CASE( np_rhgVP  )                ! Viscous-Plastic rheology                    ! 
     89         !                             !---------------------------------------------! 
     90         CALL ice_dyn_rhg_vp ( kt, stress1_i, stress2_i, stress12_i, shear_i, divu_i, delta_i ) 
    8291         !          
    8392      END SELECT 
     
    8594      IF( lrst_ice ) THEN                       !* write EVP fields in the restart file 
    8695         IF( ln_rhg_EVP )   CALL rhg_evp_rst( 'WRITE', kt ) 
     96         IF( ln_rhg_VP  )   CALL rhg_lsr_rst( 'WRITE', kt ) 
    8797      ENDIF 
    8898      ! 
     
    110120      INTEGER ::   ios, ioptio   ! Local integer output status for namelist read 
    111121      !! 
    112       NAMELIST/namdyn_rhg/  ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast, ln_rhg_chkcvg 
     122      NAMELIST/namdyn_rhg/  ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast, ln_rhg_chkcvg, & 
     123    &                       ln_rhg_VP  
    113124      !!------------------------------------------------------------------- 
    114125      ! 
     
    133144         WRITE(numout,*) '         ratio of elastic timescale over ice time step     rn_relast     = ', rn_relast 
    134145         WRITE(numout,*) '      check convergence of rheology                        ln_rhg_chkcvg = ', ln_rhg_chkcvg 
     146         WRITE(numout,*) ''     
     147         WRITE(numout,*) '      rheology VP     (icedyn_rhg_lsr)                     ln_rhg_VP     = ', ln_rhg_VP 
    135148      ENDIF 
    136149      ! 
     
    139152      IF( ln_rhg_EVP ) THEN   ;   ioptio = ioptio + 1   ;   nice_rhg = np_rhgEVP    ;   ENDIF 
    140153!!    IF( ln_rhg_EAP ) THEN   ;   ioptio = ioptio + 1   ;   nice_rhg = np_rhgEAP    ;   ENDIF 
     154      IF( ln_rhg_VP  ) THEN   ;   ioptio = ioptio + 1   ;   nice_rhg = np_rhgVP     ;   ENDIF 
    141155      IF( ioptio /= 1 )   CALL ctl_stop( 'ice_dyn_rhg_init: choose one and only one ice rheology' ) 
    142156      ! 
    143157      IF( ln_rhg_EVP  )   CALL rhg_evp_rst( 'READ' )  !* read or initialize all required files 
     158      IF( ln_rhg_VP   )   CALL rhg_lsr_rst( 'READ' )  !* read or initialize all required files 
    144159      ! 
    145160   END SUBROUTINE ice_dyn_rhg_init 
Note: See TracChangeset for help on using the changeset viewer.