Index | Rundeck: P1SoM40 | Created: Thu May 9 03:30:17 EDT 2024 |
icelake_fluxes | Module: SEAICE | File: SEAICE.f |
Author : Gavin Schmidt
Version:
Usage:
At the ice-lake interface the following equations hold: interface is at freezing point (Tb=0.) (1) -lam_i Ti/dh - rho_m shw g_T Tm = -m Lh(Tib) + m shw Tib (2) with Tib=Ti m>0, or Tib=0. m<0 (4)
Declaration:
subroutine icelake_fluxes(ti, tm, dh, dtsrc, mlsh, mflux, hflux)
- hflux :
real*8, intent(out)
- mass and heat fluxes at base of ice
- mflux :
real*8, intent(out)
- mass and heat fluxes at base of ice
- dtsrc :
real*8, intent(in)
- source time step (s)
- tm :
real*8, intent(in)
- temperatures in ice and upper lake layer (C)
- ti :
real*8, intent(in)
- temperatures in ice and upper lake layer (C)
- dh :
real*8, intent(in)
- distance from center of bottom ice layer to base of ice (m)
- mlsh :
real*8, intent(in)
- mixed layer specific heat capactity (J/m^2 C) (UNUSED)
- rsg :
real*8, parameter
- = rhow * shw * g_T turbulent energy flux (J/m^2 K s)
Initial Value= rhow*shw*1.3d-5
- rgs :
real*8, parameter
- = rhow * g_S turbulent tracer flux (kg/m^2 s)
Initial Value= rhow*3.2d-7
- left2 :
real*8
- lh :
real*8
- m :
real*8
- alamdh :
real*8