Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

hocn calculation and flushing velocity issue #504

Open
wants to merge 7 commits into
base: main
Choose a base branch
from
+34 −19
Open

hocn calculation and flushing velocity issue #504

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
33aa692
fixed hocn and flushing velocity issue, passes basic restart test
davidclemenssewall Nov 8, 2024
0b09d92
Merge branch 'main' into hocn_mushy_thermo_fix
davidclemenssewall Nov 11, 2024
5f53385
fixed indentation and relocated alvl in calling arguments
davidclemenssewall Nov 12, 2024
29f1d2d
Merge branch 'main' into hocn_mushy_thermo_fix
davidclemenssewall Dec 6, 2024
f0e1cc2
moved apond = apnd*alvl calculation to icepack_step_therm1
davidclemenssewall Dec 6, 2024
2392676
changed apnd for apond for neutral_drag_coeffs
davidclemenssewall Dec 6, 2024
f35a56d
removed extraneous tr_pond_lvl
davidclemenssewall Dec 6, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
9 changes: 3 additions & 6 deletions columnphysics/icepack_atmo.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -507,8 +507,7 @@ end subroutine atmo_boundary_const
!
! changes: Andrew Roberts, NPS (RASM/CESM coupling and documentation)

subroutine neutral_drag_coeffs (apnd, hpnd, &
ipnd, &
subroutine neutral_drag_coeffs (apondn, &
alvl, vlvl, &
aice, vice, &
vsno, aicen, &
Expand All @@ -526,9 +525,7 @@ subroutine neutral_drag_coeffs (apnd, hpnd, &
use icepack_tracers, only: tr_pond

real (kind=dbl_kind), dimension (:), intent(in) :: &
apnd ,& ! melt pond fraction of sea ice
hpnd ,& ! mean melt pond depth over sea ice
ipnd ,& ! mean ice pond depth over sea ice in cat n
eclare108213 marked this conversation as resolved.
Show resolved Hide resolved
apondn ,& ! melt pond fraction of sea ice category
alvl ,& ! level ice area fraction (of grid cell ?)
vlvl ! level ice mean thickness

Expand Down Expand Up @@ -666,7 +663,7 @@ subroutine neutral_drag_coeffs (apnd, hpnd, &
if (tr_pond) then
do n = 1,ncat
! area of pond per unit area of grid cell
apond = apond+apnd(n)*aicen(n)
apond = apond+apondn(n)*aicen(n)
enddo
endif

Expand Down
16 changes: 8 additions & 8 deletions columnphysics/icepack_therm_mushy.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@ module icepack_therm_mushy
use icepack_kinds
use icepack_parameters, only: c0, c1, c2, c8, c10
use icepack_parameters, only: p01, p05, p1, p2, p5, pi, bignum, puny
use icepack_parameters, only: viscosity_dyn, rhow, rhoi, rhos, cp_ocn, cp_ice, Lfresh, gravit
use icepack_parameters, only: viscosity_dyn, rhow, rhoi, rhos, cp_ocn, cp_ice, Lfresh, gravit, rhofresh
use icepack_parameters, only: hs_min, snwgrain
use icepack_parameters, only: a_rapid_mode, Rac_rapid_mode, tscale_pnd_drain
use icepack_parameters, only: aspect_rapid_mode, dSdt_slow_mode, phi_c_slow_mode
Expand Down Expand Up @@ -72,15 +72,15 @@ subroutine temperature_changes_salinity(dt, &
lhcoef , & ! transfer coefficient for latent heat
Tbot , & ! ice bottom surfce temperature (deg C)
sss ! sea surface salinity (PSU)

real (kind=dbl_kind), intent(inout) :: &
fswsfc , & ! SW absorbed at ice/snow surface (W m-2)
fswint ! SW absorbed in ice interior below surface (W m-2)

real (kind=dbl_kind), intent(inout) :: &
hilyr , & ! ice layer thickness (m)
hslyr , & ! snow layer thickness (m)
apond , & ! melt pond area fraction
apond , & ! melt pond area fraction of category
hpond ! melt pond depth (m)

real (kind=dbl_kind), dimension (:), intent(inout) :: &
Expand Down Expand Up @@ -3077,7 +3077,7 @@ subroutine flushing_velocity(zTin, phi, &
real(kind=dbl_kind), intent(in) :: &
hilyr , & ! ice layer thickness (m)
hpond , & ! melt pond thickness (m)
apond , & ! melt pond area (-)
apond , & ! melt pond area fraction of category (-)
hsn , & ! snow thickness (m)
hin , & ! ice thickness (m)
dt ! time step (s)
Expand Down Expand Up @@ -3138,8 +3138,8 @@ subroutine flushing_velocity(zTin, phi, &
perm_harm = real(nilyr,dbl_kind) / perm_harm

! calculate ocean surface height above bottom of ice
hocn = (ice_mass + hpond * apond * rhow + hsn * rhos) / rhow

hocn = (ice_mass + hpond * apond * rhofresh + hsn * rhos) / rhow
! calculate brine height above bottom of ice
hbrine = hin + hpond

Expand All @@ -3151,7 +3151,7 @@ subroutine flushing_velocity(zTin, phi, &

! maximum down flow to drain pond
w_down_max = (hpond * apond) / dt

! limit flow
w = min(w,w_down_max)

Expand All @@ -3178,7 +3178,7 @@ subroutine flush_pond(w, hpond, apond, dt)

real(kind=dbl_kind), intent(in) :: &
w , & ! vertical flushing Darcy flow rate (m s-1)
apond , & ! melt pond area (-)
apond , & ! melt pond area fraction of category (-)
dt ! time step (s)

real(kind=dbl_kind), intent(inout) :: &
Expand Down
28 changes: 23 additions & 5 deletions columnphysics/icepack_therm_vertical.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -122,7 +122,7 @@ subroutine thermo_vertical (dt, aicen, &
! tracers
real (kind=dbl_kind), intent(inout) :: &
Tsf , & ! ice/snow top surface temp, same as Tsfcn (deg C)
apond , & ! melt pond area fraction
apond , & ! melt pond area fraction of category
hpond ! melt pond depth (m)
! iage ! ice age (s)

Expand Down Expand Up @@ -2336,7 +2336,7 @@ subroutine icepack_step_therm1(dt, &
Tsfc , & ! ice/snow surface temperature, Tsfcn
alvl , & ! level ice area fraction
vlvl , & ! level ice volume fraction
apnd , & ! melt pond area fraction
apnd , & ! melt pond area fraction tracer
hpnd , & ! melt pond depth (m)
ipnd , & ! melt pond refrozen lid thickness (m)
iage , & ! volume-weighted ice age
Expand Down Expand Up @@ -2437,6 +2437,7 @@ subroutine icepack_step_therm1(dt, &
smliq ! tracer for mass of liquid in snow (kg/m^3)

real (kind=dbl_kind), dimension(ncat) :: &
apond , & ! melt pond area fraction of category
l_meltsliqn ! mass of snow melt local (kg/m^2)

real (kind=dbl_kind) :: &
Expand Down Expand Up @@ -2526,6 +2527,17 @@ subroutine icepack_step_therm1(dt, &
massicen(:,:) = c0
massliqn(:,:) = c0

!-----------------------------------------------------------------
! Initialize pond area fractions
!-----------------------------------------------------------------
do n= 1, ncat
if (tr_pond_lvl) then
apond(n) = apnd(n) * alvl(n)
else
apond(n) = apnd(n)
endif
enddo

!-----------------------------------------------------------------
! Initialize rate of snow loss to leads
!-----------------------------------------------------------------
Expand Down Expand Up @@ -2556,8 +2568,7 @@ subroutine icepack_step_therm1(dt, &
!-----------------------------------------------------------------

if (formdrag) then
call neutral_drag_coeffs (apnd , &
hpnd , ipnd , &
call neutral_drag_coeffs (apond , &
alvl , vlvl , &
aice , vice, &
vsno , aicen , &
Expand Down Expand Up @@ -2704,7 +2715,7 @@ subroutine icepack_step_therm1(dt, &
vicen=vicen (n), vsnon=vsnon (n), &
Tsf=Tsfc (n), zSin=zSin (:,n), &
zqin=zqin (:,n), zqsn=zqsn (:,n), &
apond=apnd (n), hpond=hpnd (n), &
apond=apond (n), hpond=hpnd (n), &
flw=flw, potT=potT, &
Qa=Qa, rhoa=rhoa, &
fsnow=fsnow, fpond=fpond, &
Expand Down Expand Up @@ -2736,6 +2747,13 @@ subroutine icepack_step_therm1(dt, &
return
endif

! Translate changes in apond into apnd tracer
if (tr_pond_lvl) then
apnd(n) = apond(n) / alvl(n)
else
apnd(n) = apond(n)
endif

if (snwgrain) then
rsnwn (:,n) = rsnw (:)
smicen(:,n) = smice(:)
Expand Down