4.10 Fortran: Module Interface ice_dyn_evp - elastic-viscous-plastic sea ice dynamics model (Source File: ice_dyn_evp.F)

Elastic-viscous-plastic sea ice dynamics model
Computes ice velocity

See E. C. Hunke and J. K. Dukowicz. An elastic-viscous-plastic model
for sea ice dynamics. 1997, J. Phys. Oceanogr., Vol 27, 1849-1867.

REVISION HISTORY:

   author: Elizabeth C. Hunke
           Fluid Dynamics Group, Los Alamos National Laboratory

INTERFACE:

       module ice_dyn_evp

USES:

       use ice_kinds_mod
       use ice_domain
       use ice_grid
       use ice_constants

4.10.1 evp - elastic-viscous-plastic dynamics driver

INTERFACE:

       subroutine evp(kstrngth)

DESCRIPTION:

Elastic-viscous-plastic dynamics driver

REVISION HISTORY:

   author: Elizabeth C. Hunke
           Fluid Dynamics Group, Los Alamos National Laboratory

USES:

       use ice_timers

INPUT/OUTPUT PARAMETERS:

       integer (kind=int_kind), intent(in) ::
      &   kstrngth

4.10.2 init_evp - initialize parameters needed for evp dynamics

INTERFACE:

       subroutine init_evp

DESCRIPTION:

Initialize parameters and variables needed for the evp dynamics

REVISION HISTORY:

   author: Elizabeth C. Hunke
           Fluid Dynamics Group, Los Alamos National Laboratory

USES:

       use ice_calendar

INPUT/OUTPUT PARAMETERS:

4.10.3 evp_prep - compute quantities needed for stress tensor and mom eqns

INTERFACE:

       subroutine evp_prep(kstrngth)

DESCRIPTION:

Computes quantities needed in the stress tensor (sigma)
and momentum (u) equations, but which do not change during
the thermodynamics/transport time step:
-wind stress shift to U grid,
-ice mass and ice extent masks,
-pressure (strength), and part of the forcing stresses
initializes ice velocity for new points to ocean sfc current

REVISION HISTORY:

   author: Elizabeth C. Hunke
           Fluid Dynamics Group, Los Alamos National Laboratory

USES:

       use ice_flux
       use ice_calendar
       use ice_state

INPUT/OUTPUT PARAMETERS:

       integer (kind=int_kind), intent(in) ::
      &   kstrngth

4.10.4 stress - computes strain rates and internal stress components

INTERFACE:

       subroutine stress(ksub)

DESCRIPTION:

Computes the rates of strain and internal stress components for
each of the four corners on each T-grid cell

REVISION HISTORY:

   author: Elizabeth C. Hunke
           Fluid Dynamics Group, Los Alamos National Laboratory

INPUT/OUTPUT PARAMETERS:

       integer (kind=int_kind), intent(in) :: ksub  ! subcycling step

4.10.5 stepu - surface stresses and integrates mom eqn for u,v

INTERFACE:

       subroutine stepu

DESCRIPTION:

Calculation of the surface stresses
Integration of the momentum equation to find velocity (u,v)

REVISION HISTORY:

   author: Elizabeth C. Hunke
           Fluid Dynamics Group, Los Alamos National Laboratory

USES:

       use ice_flux

INPUT/OUTPUT PARAMETERS:

4.10.6 evp_finish - calculates ice-ocean stress

INTERFACE:

       subroutine evp_finish

DESCRIPTION:

Calculation of the ice-ocean stress
...the sign will be reversed later...

REVISION HISTORY:

   author: Elizabeth C. Hunke
           Fluid Dynamics Group, Los Alamos National Laboratory

USES:

       use ice_flux

INPUT/OUTPUT PARAMETERS:

4.10.7 principal_stress - computes principal stress for yield curve

INTERFACE:

       subroutine principal_stress

DESCRIPTION:

Computes principal stresses for comparison with the theoretical yield curve; northeast values

REVISION HISTORY:

   author: Elizabeth C. Hunke
           Fluid Dynamics Group, Los Alamos National Laboratory

USES:

INPUT/OUTPUT PARAMETERS: