4.31 Fortran: Module Interface ice_vthermo - main routines for energy-conserving vertical thermodynamics (Source File: ice_vthermo.F)

Energy-conserving sea ice model
Main routines for vertical thermodynamics

See Bitz, C.M., and W.H. Lipscomb, 1999: An energy-conserving thermodynamic model of sea ice, J. Geophys. Res., 104, 15,669-15,677.

See Bitz, C.M., M.M. Holland, A.J. Weaver, M. Eby, 2001: Simulating the ice-thickness distribution in a climate model, J. Geophys. Res., 106, 2441-2464.

REVISION HISTORY:

   author C. M. Bitz

INTERFACE:

       module ice_vthermo

USES:

       use ice_model_size
       use ice_domain 
       use ice_constants
       use ice_state

4.31.1 init_thermo - setup salinity profile and tmelt for each layer

INTERFACE:

       subroutine init_thermo

DESCRIPTION:

Setup the salinity profile and the melting temperature for each layer

REVISION HISTORY:

   author: C.M.Bitz

USES:

       use ice_itd
       use shr_sys_mod, only : shr_sys_abort

INPUT/OUTPUT PARAMETERS:

4.31.2 thermo_vertical - heat budget over open water and ice

INTERFACE:

       subroutine thermo_vertical(i,j,hin,hsn,tiz,Fbot,
      &         dhi1,dhin,dhsn,hfatm,hfocn,qin)

DESCRIPTION:

Heat budget over open water and ice
NOTE the wind stress is also computed here for later use

REVISION HISTORY:

   author:  C.M.Bitz

USES:

       use ice_flux
       use ice_albedo
       use ice_atmo
       use ice_history
       use ice_calendar
       use ice_itd
       use ice_dh
       use ice_tstm

INPUT/OUTPUT PARAMETERS:

       integer (kind=int_kind), intent(in) :: i, j
       real (kind=dbl_kind), intent(in) ::
      &   Fbot
 
       real (kind=dbl_kind), intent(inout) ::
      &   hin(ncat)          ! ice thickness for each cat               (m)
      &,  hsn(ncat)          ! snow thickness for each cat              (m)
      &,  tiz(0:nmax,ncat)   ! temp of each layer for each cat          (C)
 
       real (kind=dbl_kind), intent(out) ::
      &   dhi1               ! melt at bottom and top plus sublimation  (m)
      &,  dhin(ncat)         ! ice thickness change                     (m)
      &,  dhsn(ncat)         ! snow thickness change                    (m)
      &,  hfatm              ! heat flx  to ice-snow from atm           (J/m**2)
      &,  hfocn              ! heat flx  to ice-snow from ocn           (J/m**2)
      &,  qin(nmax,ncat)     ! enthalpy per unit volume                 (J/m**3)