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.
icerhg.F90 in branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3 – NEMO

source: branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icerhg.F90 @ 8517

Last change on this file since 8517 was 8517, checked in by clem, 7 years ago

changes in style - part6 - one more round
File size: 7.2 KB
RevLine 
[8407]1MODULE icerhg
2   !!======================================================================
3   !!                     ***  MODULE  icerhg  ***
4   !!   Sea-Ice dynamics : 
5   !!======================================================================
6   !! history :  1.0  ! 2002-08  (C. Ethe, G. Madec)  original VP code
7   !!            3.0  ! 2007-03  (MA Morales Maqueda, S. Bouillon, M. Vancoppenolle)  LIM3: EVP-Cgrid
8   !!            3.5  ! 2011-02  (G. Madec) dynamical allocation
9   !!----------------------------------------------------------------------
10#if defined key_lim3
11   !!----------------------------------------------------------------------
[8486]12   !!   'key_lim3'                                       LIM3 sea-ice model
[8407]13   !!----------------------------------------------------------------------
14   !!    ice_rhg      : computes ice velocities
15   !!    ice_rhg_init : initialization and namelist read
16   !!----------------------------------------------------------------------
[8486]17   USE phycst         ! physical constants
18   USE dom_oce        ! ocean space and time domain
19   USE ice            ! sea-ice: variables
20   USE icerhg_evp     ! sea-ice: EVP rheology
21   USE icectl         ! sea-ice: control prints
[8407]22   !
[8486]23   USE lbclnk         ! lateral boundary conditions - MPP exchanges
24   USE lib_mpp        ! MPP library
25   USE in_out_manager ! I/O manager
26   USE lib_fortran    ! glob_sum
27   USE timing         ! Timing
[8407]28
29   IMPLICIT NONE
30   PRIVATE
31
[8512]32   PUBLIC   ice_rhg        ! called by icestp.F90
33   PUBLIC   ice_rhg_init   ! called by icestp.F90
[8517]34
35   INTEGER ::              nice_rhg   ! choice of the type of rheology
36   !                                        ! associated indices:
37   INTEGER, PARAMETER ::   np_rhgEVP = 1   ! EVP rheology
38!! INTEGER, PARAMETER ::   np_rhgEAP = 2   ! EAP rheology
[8512]39   
[8407]40   !! * Substitutions
41#  include "vectopt_loop_substitute.h90"
42   !!----------------------------------------------------------------------
[8486]43   !! NEMO/ICE 4.0 , NEMO Consortium (2017)
[8407]44   !! $Id: icerhg.F90 8378 2017-07-26 13:55:59Z clem $
45   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
46   !!----------------------------------------------------------------------
47CONTAINS
48
49   SUBROUTINE ice_rhg( kt )
50      !!-------------------------------------------------------------------
51      !!               ***  ROUTINE ice_rhg  ***
52      !!               
53      !! ** Purpose :   compute ice velocity
54      !!
55      !! ** Action  : comupte - ice velocity (u_ice, v_ice)
56      !!                      - 3 components of the stress tensor (stress1_i, stress2_i, stress12_i)
57      !!                      - shear, divergence and delta (shear_i, divu_i, delta_i)
58      !!--------------------------------------------------------------------
[8486]59      INTEGER, INTENT(in) ::   kt     ! ice time step
[8407]60      !!--------------------------------------------------------------------
[8517]61      ! controls
62      IF( nn_timing == 1 )   CALL timing_start('icerhg')                                                             ! timing
63      IF( ln_icediachk   )   CALL ice_cons_hsm(0, 'icerhg', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation
[8486]64      !
[8426]65      IF( kt == nit000 .AND. lwp ) THEN
66         WRITE(numout,*)
[8512]67         WRITE(numout,*)'ice_rhg: sea-ice rheology'
68         WRITE(numout,*)'~~~~~~~'
[8426]69      ENDIF
70
[8517]71      ! --------
72      ! Rheology
73      ! --------   
74      SELECT CASE( nice_rhg )
75      !                                !------------------------!
76      CASE( np_rhgEVP )                ! Elasto-Viscous-Plastic !
77         !                             !------------------------!
78         CALL ice_rhg_evp( kt, stress1_i, stress2_i, stress12_i, u_ice, v_ice, shear_i, divu_i, delta_i )
79
80      END SELECT
[8407]81      !
[8517]82      ! controls
83      IF( ln_icediachk   )   CALL ice_cons_hsm(1, 'icerhg', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation
84      IF( ln_ctl         )   CALL ice_prt3D   ('icerhg')                                                             ! prints
85      IF( nn_timing == 1 )   CALL timing_stop ('icerhg')                                                             ! timing
[8407]86      !
87   END SUBROUTINE ice_rhg
88
89
90   SUBROUTINE ice_rhg_init
91      !!-------------------------------------------------------------------
92      !!                  ***  ROUTINE ice_rhg_init  ***
93      !!
94      !! ** Purpose : Physical constants and parameters linked to the ice
95      !!      dynamics
96      !!
[8512]97      !! ** Method  :  Read the namice_rhg namelist and check the ice-dynamic
[8407]98      !!       parameter values called at the first timestep (nit000)
99      !!
[8512]100      !! ** input   :   Namelist namice_rhg
[8407]101      !!-------------------------------------------------------------------
[8517]102      INTEGER ::   ios, ioptio   ! Local integer output status for namelist read
[8486]103      !!
[8516]104      NAMELIST/namice_rhg/  ln_rhg_EVP, rn_creepl, rn_ecc , nn_nevp, rn_relast
[8407]105      !!-------------------------------------------------------------------
[8486]106      !
[8512]107      REWIND( numnam_ice_ref )         ! Namelist namice_rhg in reference namelist : Ice dynamics
108      READ  ( numnam_ice_ref, namice_rhg, IOSTAT = ios, ERR = 901)
109901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_rhg in reference namelist', lwp )
[8486]110      !
[8512]111      REWIND( numnam_ice_cfg )         ! Namelist namice_rhg in configuration namelist : Ice dynamics
112      READ  ( numnam_ice_cfg, namice_rhg, IOSTAT = ios, ERR = 902 )
113902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_rhg in configuration namelist', lwp )
114      IF(lwm) WRITE ( numoni, namice_rhg )
[8486]115      !
116      IF(lwp) THEN                     ! control print
[8407]117         WRITE(numout,*)
[8514]118         WRITE(numout,*) 'ice_rhg_init: ice parameters for ice dynamics '
[8407]119         WRITE(numout,*) '~~~~~~~~~~~~'
[8512]120         WRITE(numout,*) '   Namelist namice_rhg'
[8516]121         WRITE(numout,*) '      rheology EVP (icerhg_evp)                                   ln_rhg_EVP    = ', ln_rhg_EVP
[8486]122         WRITE(numout,*) '         creep limit                                              rn_creepl     = ', rn_creepl
123         WRITE(numout,*) '         eccentricity of the elliptical yield curve               rn_ecc        = ', rn_ecc
124         WRITE(numout,*) '         number of iterations for subcycling                      nn_nevp       = ', nn_nevp
125         WRITE(numout,*) '         ratio of elastic timescale over ice time step            rn_relast     = ', rn_relast
[8407]126      ENDIF
127      !
[8517]128      !                             !== set the choice of ice advection ==!
129      ioptio = 0 
130      IF( ln_rhg_EVP ) THEN   ;   ioptio = ioptio + 1   ;   nice_rhg = np_rhgEVP    ;   ENDIF
131!!    IF( ln_rhg_EAP ) THEN   ;   ioptio = ioptio + 1   ;   nice_rhg = np_rhgEAP    ;   ENDIF
132      IF( ioptio /= 1 )   CALL ctl_stop( 'ice_rhg_init: choose one and only one ice rheology' )
133
134     !                              ! allocate tke arrays
135!!clem example      IF( zdf_tke_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'zdf_tke_init : unable to allocate arrays' )
136      !
[8407]137   END SUBROUTINE ice_rhg_init
138
[8486]139#else
140   !!----------------------------------------------------------------------
141   !!   Default option         Empty module          NO LIM-3 sea-ice model
142   !!----------------------------------------------------------------------
[8407]143#endif 
144
145   !!======================================================================
146END MODULE icerhg
Note: See TracBrowser for help on using the repository browser.