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. BitzINTERFACE:
module ice_vthermoUSES:
use ice_model_size use ice_domain use ice_constants use ice_state
INTERFACE:
subroutine init_thermoDESCRIPTION:
Setup the salinity profile and the melting temperature for each layer
REVISION HISTORY:
author: C.M.BitzUSES:
use ice_itd use shr_sys_mod, only : shr_sys_abortINPUT/OUTPUT PARAMETERS:
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.BitzUSES:
use ice_flux use ice_albedo use ice_atmo use ice_history use ice_calendar use ice_itd use ice_dh use ice_tstmINPUT/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)