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

Separating sea ice melt/formation from liquid precipitation field #214

Open
wants to merge 5 commits into
base: dev/gfdl
Choose a base branch
from
Open

Separating sea ice melt/formation from liquid precipitation field #214

Show file tree
Hide file tree
Changes from 4 commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
1c8dc7e
Attempt at separating sea ice melt/formation from liquid precipitation
Aug 14, 2024
ec559ba
Added CMOR diagnostic for sea ice melt
Aug 15, 2024
5ecae71
Reset IOF%seaice_melt to zero every timestep just like old net_melt
Dec 11, 2024
49a2829
Simpler way of resetting seaice_melt to zero
Dec 17, 2024
9b49d6a
Add parentheses to preserve order of operations
Feb 4, 2025
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
33 changes: 13 additions & 20 deletions src/SIS_slow_thermo.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -573,11 +573,8 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)
tmp2d, & ! A temporary array for mass balance diagnostics [R Z yr-1 ~> kg m-2 yr-1]
qflx_lim_ice, & ! Ice limiting heat flux [Q R Z T-1 ~> W m-2]
qflx_res_ice, & ! Ice restoring heat flux [Q R Z T-1 ~> W m-2]
cool_nudge, & ! A heat flux out of the sea ice that
cool_nudge ! A heat flux out of the sea ice that
! acts to create sea-ice [Q R Z T-1 ~> W m-2].
net_melt ! The net mass flux from the ice and snow into the
! ocean due to melting and freezing integrated
! across all categories [R Z T-1 ~> kg m-2 s-1].
real, dimension(SZI_(G),SZJ_(G),1:IG%CatIce) :: &
heat_in, & ! The input heat [Q R Z ~> J m-2]
enth_prev, & ! The previous column integrated enthalpy by category [Q R Z ~> J m-2]
Expand Down Expand Up @@ -818,7 +815,7 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)

call cpu_clock_begin(iceClock5)

snow_to_ice(:,:,:) = 0.0 ; net_melt(:,:) = 0.0
snow_to_ice(:,:,:) = 0.0 ;
bsnk(:,:) = 0.0
salt_change(:,:) = 0.0
h2o_change(:,:) = 0.0
Expand Down Expand Up @@ -857,6 +854,9 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)
IOF%lprec_ocn_top(i,j) = IOF%lprec_ocn_top(i,j) + IST%part_size(i,j,k) * FIA%lprec_top(i,j,k)
enddo ; enddo ; enddo

! Reset sea ice melt rate to zero
IOF%seaice_melt(:,:) = 0.

! Add fluxes of snow and other properties to the ocean due to recent ridging or drifting events.
if (allocated(IST%snow_to_ocn)) then
!$OMP do
Expand All @@ -876,7 +876,7 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)

!$OMP parallel do default(none) shared(isc,iec,jsc,jec,ncat,G,US,IST,S_col0,NkIce,S_col,dt_slow, &
!$OMP snow_to_ice,heat_in,I_NK,enth_prev,enth_mass_in_col,bsnk, &
!$OMP Idt_slow,salt_change,net_melt,LatHtFus,LatHtVap,IG,CS,OSS, &
!$OMP Idt_slow,salt_change,LatHtFus,LatHtVap,IG,CS,OSS, &
!$OMP FIA,IOF,npassive,nb) &
!$OMP private(mass_prev,enthalpy,enthalpy_ocean,Salin, &
!$OMP heat_to_ocn,h2o_ice_to_ocn,h2o_ocn_to_ice, &
Expand Down Expand Up @@ -1010,7 +1010,7 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)
enddo
IOF%flux_sw_ocn(i,j,VIS_DIF) = IOF%flux_sw_ocn(i,j,VIS_DIF) + IST%part_size(i,j,k) * &
(sw_tot * FIA%sw_abs_ocn(i,j,k))
net_melt(i,j) = net_melt(i,j) + IST%part_size(i,j,k) * &
IOF%seaice_melt(i,j) = IOF%seaice_melt(i,j) + IST%part_size(i,j,k) * &
((h2o_ice_to_ocn-h2o_ocn_to_ice)*Idt_slow)
bsnk(i,j) = bsnk(i,j) - IST%part_size(i,j,k)*bablt*Idt_slow ! bot. melt. ablation

Expand All @@ -1019,7 +1019,7 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)

!$OMP parallel do default(none) shared(isc,iec,jsc,jec,ncat,npassive,G,US,IG,IST,S_col0,NkIce, &
!$OMP S_col,dt_slow,snow_to_ice,heat_in,I_NK,enth_mass_in_col, &
!$OMP enth_prev,Idt_slow,bsnk,salt_change,net_melt,LatHtFus, &
!$OMP enth_prev,Idt_slow,bsnk,salt_change,LatHtFus, &
!$OMP FIA,CS,OSS,IOF) &
!$OMP private(mass_prev,enthalpy,enthalpy_ocean,Salin,heat_to_ocn, &
!$OMP h2o_ice_to_ocn,h2o_ocn_to_ice,evap_from_ocn,salt_to_ice, &
Expand Down Expand Up @@ -1156,7 +1156,7 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)
! IST%part_size(i,j,k) * enth_ocn_to_ice
IOF%Enth_Mass_in_ocn(i,j) = IOF%Enth_Mass_in_ocn(i,j) + &
IST%part_size(i,j,k) * (h2o_ocn_to_ice * enthalpy_ocean)
net_melt(i,j) = net_melt(i,j) - &
IOF%seaice_melt(i,j) = IOF%seaice_melt(i,j) - &
(h2o_ocn_to_ice * IST%part_size(i,j,k)) * Idt_slow

if (CS%column_check) then
Expand Down Expand Up @@ -1222,7 +1222,7 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)
IST%mH_ice(i,j,k) = frac_keep*IST%mH_ice(i,j,k)
IST%mH_snow(i,j,k) = frac_keep*IST%mH_snow(i,j,k)
enddo
net_melt(i,j) = net_melt(i,j) + h2o_ice_to_ocn * Idt_slow
IOF%seaice_melt(i,j) = IOF%seaice_melt(i,j) + h2o_ice_to_ocn * Idt_slow
qflx_lim_ice(i,j) = enth_to_melt * Idt_slow
IOF%Enth_Mass_out_ocn(i,j) = IOF%Enth_Mass_out_ocn(i,j) - enth_ice_to_ocn
if (CS%ice_rel_salin > 0.0) then
Expand Down Expand Up @@ -1322,7 +1322,7 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)

! With transmuted ice, the ice is non-conservatively changed to match the ocean properties.
IOF%flux_salt(i,j) = IOF%flux_salt(i,j) + salt_to_ocn * (0.001*Idt_slow)
net_melt(i,j) = net_melt(i,j) + water_to_ocn * Idt_slow ! This goes to IOF%lprec_ocn_top.
IOF%seaice_melt(i,j) = IOF%seaice_melt(i,j) + water_to_ocn * Idt_slow
IOF%Enth_mass_out_ocn(i,j) = IOF%Enth_mass_out_ocn(i,j) + heat_to_ocn

! With transmuted ice, the imbalances are stored to close the heat and salt budgets.
Expand Down Expand Up @@ -1361,20 +1361,14 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)
scale=US%RZ_T_to_kg_m2s*Idt_slow)
if (CS%id_qflim>0) call post_data(CS%id_qflim, qflx_lim_ice, CS%diag)
if (CS%id_qfres>0) call post_data(CS%id_qfres, qflx_res_ice, CS%diag)
if (CS%id_net_melt>0) call post_data(CS%id_net_melt, net_melt, CS%diag)
if (CS%id_CMOR_melt>0) call post_data(CS%id_CMOR_melt, net_melt, CS%diag)
if (CS%id_net_melt>0) call post_data(CS%id_net_melt, IOF%seaice_melt, CS%diag)
if (CS%id_CMOR_melt>0) call post_data(CS%id_CMOR_melt, IOF%seaice_melt, CS%diag)
Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

to reproduce old answers, I think this should be okay, since IOF%seaice_melt and net_melt are the same. Right?


if (coupler_type_initialized(IOF%tr_flux_ocn_top)) &
call coupler_type_send_data(IOF%tr_flux_ocn_top, CS%Time)

call disable_SIS_averaging(CS%diag)

! Combine the liquid precipitation with the net melt of ice and snow for
Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think a bug fix flag here will work to get the old diagnostics, but it would be ideal (in my opinion) if there was a way to save sea ice melt and sea ice melt + lprec from the same simulation.

! passing to the ocean. These may later be kept separate.
do j=jsc,jec ; do i=isc,iec
IOF%lprec_ocn_top(i,j) = IOF%lprec_ocn_top(i,j) + net_melt(i,j)
enddo ; enddo

! Make sure TrLay is no longer allocated
if(allocated(TrLay)) deallocate(TrLay)
end subroutine SIS2_thermodynamics
Expand Down Expand Up @@ -1562,7 +1556,6 @@ subroutine SIS_slow_thermo_init(Time, G, US, IG, param_file, diag, CS, tracer_fl
CS%id_CMOR_melt = register_diag_field('ice_model','fsitherm' ,diag%axesT1, Time, &
'water_flux_into_sea_water_due_to_sea_ice_thermodynamics', &
'kg m-2 s-1', conversion=US%RZ_T_to_kg_m2s, missing_value=missing)

if (CS%do_ice_restore) then
CS%id_qfres = register_diag_field('ice_model', 'QFLX_RESTORE_ICE', diag%axesT1, Time, &
'Ice Restoring heat flux', 'W/m^2', conversion=US%QRZ_T_to_W_m2, missing_value=missing)
Expand Down
1 change: 1 addition & 0 deletions src/SIS_sum_output.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -759,6 +759,7 @@ subroutine accumulate_bottom_input(IST, OSS, FIA, IOF, dt, G, US, IG, CS)
do j=jsc,jec ; do i=isc,iec
CS%water_in_col(i,j) = CS%water_in_col(i,j) - dt * &
( ((FIA%runoff(i,j) + FIA%calving(i,j)) + &
IOF%seaice_melt(i,j) + &
Copy link
Member

@Hallberg-NOAA Hallberg-NOAA Feb 1, 2025 •
edited
Loading

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Please add parentheses to specify the order with which IOF%seaice_melt is added to the other fields so that this diagnostic does not capriciously change with the level of optimization.
My preference would be for IOF%seaice_melt to be added first to IOF%lprec_ocn_top so that the answers more closely resemble what was done previously, but other options could be justified as long as the order of arithmetic is specified unambiguously.

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thank you @Hallberg-NOAA; I've started to make the necessary changes wherever fluxes%lprec shows up.

I am not sure what to do here:
https://github.com/NOAA-GFDL/MOM6/blob/37d1385b2979a453c5b2c24d91bbea7c74c6a82e/src/core/MOM_forcing_type.F90#L2275

Is there any way to separately increment these averages for fluxes%lprec and fluxes%seaice_melt while preserving the overall order of operations?

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Similarly, in any if associated(fluxes%lprec) blocks, do I need to separately treat the two cases of if associated(fluxes%seaice_melt)?

(IOF%lprec_ocn_top(i,j) + IOF%fprec_ocn_top(i,j))) - IOF%evap_ocn_top(i,j) )

Flux_SW = 0.0
Expand Down
5 changes: 4 additions & 1 deletion src/SIS_types.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -358,6 +358,7 @@ module SIS_types
flux_lh_ocn_top, & !< The upward flux of latent heat at the ocean surface [Q R Z T-1 ~> W m-2].
lprec_ocn_top, & !< The downward flux of liquid precipitation at the ocean surface [R Z T-1 ~> kg m-2 s-1].
fprec_ocn_top, & !< The downward flux of frozen precipitation at the ocean surface [R Z T-1 ~> kg m-2 s-1].
seaice_melt, & !< The downward freshwater flux into the ocean due to sea ice melt [R Z T-1 ~> kg m-2 s-1].
flux_u_ocn, & !< The flux of x-momentum into the ocean at locations given by
!! flux_uv_stagger [R Z L T-2 ~> Pa].
!! Note that regardless of the staggering, flux_u_ocn is allocated as though on an A-grid.
Expand Down Expand Up @@ -925,6 +926,7 @@ subroutine alloc_ice_ocean_flux(IOF, HI, do_stress_mag, do_iceberg_fields, do_tr
allocate(IOF%flux_sw_ocn(SZI_(HI), SZJ_(HI), NBANDS), source=0.0)
allocate(IOF%lprec_ocn_top(SZI_(HI), SZJ_(HI)), source=0.0)
allocate(IOF%fprec_ocn_top(SZI_(HI), SZJ_(HI)), source=0.0)
allocate(IOF%seaice_melt(SZI_(HI), SZJ_(HI)), source=0.0)
allocate(IOF%flux_u_ocn(SZI_(HI), SZJ_(HI)), source=0.0)
allocate(IOF%flux_v_ocn(SZI_(HI), SZJ_(HI)), source=0.0)
if (alloc_stress_mag) then
Expand Down Expand Up @@ -2171,7 +2173,7 @@ subroutine dealloc_ice_ocean_flux(IOF)
deallocate(IOF%flux_sh_ocn_top, IOF%evap_ocn_top)
deallocate(IOF%flux_lw_ocn_top, IOF%flux_lh_ocn_top)
deallocate(IOF%flux_sw_ocn)
deallocate(IOF%lprec_ocn_top, IOF%fprec_ocn_top, IOF%flux_salt)
deallocate(IOF%lprec_ocn_top, IOF%fprec_ocn_top, IOF%seaice_melt, IOF%flux_salt)
deallocate(IOF%flux_u_ocn, IOF%flux_v_ocn, IOF%pres_ocn_top, IOF%mass_ice_sn_p)
if (allocated(IOF%stress_mag)) deallocate(IOF%stress_mag)
if (allocated(IOF%transmutation_salt_flux)) deallocate(IOF%transmutation_salt_flux)
Expand Down Expand Up @@ -2218,6 +2220,7 @@ subroutine IOF_chksum(mesg, IOF, G, US, mech_fluxes, thermo_fluxes)
call hchksum(IOF%evap_ocn_top, trim(mesg)//" IOF%evap_ocn_top", G%HI, scale=US%RZ_T_to_kg_m2s)
call hchksum(IOF%lprec_ocn_top, trim(mesg)//" IOF%lprec_ocn_top", G%HI, scale=US%RZ_T_to_kg_m2s)
call hchksum(IOF%fprec_ocn_top, trim(mesg)//" IOF%fprec_ocn_top", G%HI, scale=US%RZ_T_to_kg_m2s)
call hchksum(IOF%seaice_melt, trim(mesg)//" IOF%seaice_melt", G%HI, scale=US%RZ_T_to_kg_m2s)
endif
if (do_mech) then
call hchksum(IOF%flux_u_ocn, trim(mesg)//" IOF%flux_u_ocn", G%HI, scale=US%RZ_T_to_kg_m2s*US%L_T_to_m_s)
Expand Down
10 changes: 6 additions & 4 deletions src/combined_ice_ocean_driver.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -304,6 +304,7 @@ subroutine direct_flux_ice_to_IOB(Time, Ice, IOB, do_thermo)
if (ASSOCIATED(IOB%lw_flux)) IOB%lw_flux(:,:) = Ice%flux_lw(:,:)
if (ASSOCIATED(IOB%lprec)) IOB%lprec(:,:) = Ice%lprec(:,:)
if (ASSOCIATED(IOB%fprec)) IOB%fprec(:,:) = Ice%fprec(:,:)
if (ASSOCIATED(IOB%seaice_melt)) IOB%seaice_melt(:,:) = Ice%seaice_melt(:,:)
if (ASSOCIATED(IOB%runoff)) IOB%runoff(:,:) = Ice%runoff(:,:)
if (ASSOCIATED(IOB%calving)) IOB%calving(:,:) = Ice%calving
if (ASSOCIATED(IOB%runoff_hflx)) IOB%runoff_hflx(:,:) = Ice%runoff_hflx(:,:)
Expand All @@ -326,10 +327,11 @@ subroutine direct_flux_ice_to_IOB(Time, Ice, IOB, do_thermo)
call data_override('OCN', 'sw_flux_nir_dif', IOB%sw_flux_nir_dif, Time)
call data_override('OCN', 'sw_flux_vis_dir', IOB%sw_flux_vis_dir, Time)
call data_override('OCN', 'sw_flux_vis_dif', IOB%sw_flux_vis_dif, Time)
call data_override('OCN', 'lprec', IOB%lprec , Time)
call data_override('OCN', 'fprec', IOB%fprec , Time)
call data_override('OCN', 'runoff', IOB%runoff , Time)
call data_override('OCN', 'calving', IOB%calving , Time)
call data_override('OCN', 'lprec', IOB%lprec , Time)
call data_override('OCN', 'fprec', IOB%fprec , Time)
call data_override('OCN', 'seaice_melt', IOB%seaice_melt , Time)
call data_override('OCN', 'runoff', IOB%runoff , Time)
call data_override('OCN', 'calving', IOB%calving , Time)
call data_override('OCN', 'runoff_hflx', IOB%runoff_hflx , Time)
call data_override('OCN', 'calving_hflx', IOB%calving_hflx , Time)
endif
Expand Down
6 changes: 3 additions & 3 deletions src/ice_boundary_types.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -37,9 +37,9 @@ module ice_boundary_types
frazil => NULL(), & !< The frazil heat rejected by the ocean [J m-2].
sea_level => NULL(), & !< The sea level after adjustment for any surface
!! pressure that the ocean allows to be expressed [m].
calving => NULL(), & !< The mass flux per unit area of the ice shelf to convert to
!!bergs [RZ_T ~> kg m-2 s-1].
calving_hflx => NULL() !< Calving heat flux [Q R Z T-1 ~> W m-2].
calving => NULL(), & !< The mass flux per unit area of the ice shelf to convert to
!!bergs [RZ_T ~> kg m-2 s-1].
calving_hflx => NULL() !< Calving heat flux [Q R Z T-1 ~> W m-2].
real, dimension(:,:,:), pointer :: data =>NULL() !< S collective field for "named" fields above
integer :: stagger = BGRID_NE !< A flag indicating how the velocities are staggered.
integer :: xtype !< A flag indicating the exchange type, which may be set to
Expand Down
4 changes: 3 additions & 1 deletion src/ice_model.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -623,6 +623,7 @@ subroutine set_ocean_top_fluxes(Ice, IST, IOF, FIA, OSS, G, US, IG, sCS)
Ice%flux_sw_vis_dir(:,:) = 0.0 ; Ice%flux_sw_vis_dif(:,:) = 0.0
Ice%flux_lw(:,:) = 0.0 ; Ice%flux_lh(:,:) = 0.0
Ice%fprec(:,:) = 0.0 ; Ice%lprec(:,:) = 0.0
Ice%seaice_melt(:,:) = 0.0
call coupler_type_rescale_data(Ice%ocean_fluxes, 0.0)

i_off = LBOUND(Ice%flux_t,1) - G%isc ; j_off = LBOUND(Ice%flux_t,2) - G%jsc
Expand All @@ -640,6 +641,7 @@ subroutine set_ocean_top_fluxes(Ice, IST, IOF, FIA, OSS, G, US, IG, sCS)
Ice%flux_lh(i2,j2) = US%QRZ_T_to_W_m2*IOF%flux_lh_ocn_top(i,j)
Ice%fprec(i2,j2) = US%RZ_T_to_kg_m2s*IOF%fprec_ocn_top(i,j)
Ice%lprec(i2,j2) = US%RZ_T_to_kg_m2s*IOF%lprec_ocn_top(i,j)
Ice%seaice_melt(i2,j2) = US%RZ_T_to_kg_m2s*IOF%seaice_melt(i,j)
Ice%runoff(i2,j2) = US%RZ_T_to_kg_m2s*FIA%runoff(i,j)
Ice%calving(i2,j2) = US%RZ_T_to_kg_m2s*FIA%calving(i,j)
Ice%runoff_hflx(i2,j2) = US%QRZ_T_to_W_m2*FIA%runoff_hflx(i,j)
Expand All @@ -659,7 +661,7 @@ subroutine set_ocean_top_fluxes(Ice, IST, IOF, FIA, OSS, G, US, IG, sCS)
if (allocated(IOF%melt_nudge)) then
do j=jsc,jec ; do i=isc,iec
i2 = i+i_off ; j2 = j+j_off! Use these to correct for indexing differences.
Ice%lprec(i2,j2) = Ice%lprec(i2,j2) + US%RZ_T_to_kg_m2s*IOF%melt_nudge(i,j)
Ice%seaice_melt(i2,j2) = Ice%seaice_melt(i2,j2) + US%RZ_T_to_kg_m2s*IOF%melt_nudge(i,j)
hdrake marked this conversation as resolved.
Show resolved Hide resolved
enddo ; enddo
endif

Expand Down
12 changes: 9 additions & 3 deletions src/ice_type.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -97,9 +97,10 @@ module ice_type_mod
flux_sw_vis_dif => NULL(), & !< The diffuse visible shortwave heat flux into the ocean [W m-2].
flux_sw_nir_dir => NULL(), & !< The direct near-infrared heat flux into the ocean [W m-2].
flux_sw_nir_dif => NULL(), & !< The diffuse near-infrared heat flux into the ocean [W m-2].
flux_lh => NULL(), & !< The latent heat flux out of the ocean [W m-2].
lprec => NULL(), & !< The liquid precipitation flux into the ocean [kg m-2].
fprec => NULL(), & !< The frozen precipitation flux into the ocean [kg m-2].
flux_lh => NULL(), & !< The latent heat flux out of the ocean [W m-2].
lprec => NULL(), & !< The liquid precipitation flux into the ocean [kg m-2].
fprec => NULL(), & !< The frozen precipitation flux into the ocean [kg m-2].
seaice_melt => NULL(), & !< The sea ice meltwater flux into the ocean [kg m-2].
p_surf => NULL(), & !< The pressure at the ocean surface [Pa]. This may
!! or may not include atmospheric pressure.
runoff => NULL(), & !< Liquid runoff into the ocean [kg m-2].
Expand Down Expand Up @@ -198,6 +199,7 @@ subroutine ice_type_slow_reg_restarts(domain, CatIce, param_file, Ice, &
call safe_alloc_ptr(Ice%flux_lh, isc, iec, jsc, jec) !NI
call safe_alloc_ptr(Ice%lprec, isc, iec, jsc, jec)
call safe_alloc_ptr(Ice%fprec, isc, iec, jsc, jec)
call safe_alloc_ptr(Ice%seaice_melt, isc, iec, jsc, jec)
call safe_alloc_ptr(Ice%p_surf, isc, iec, jsc, jec)
call safe_alloc_ptr(Ice%runoff, isc, iec, jsc, jec)
call safe_alloc_ptr(Ice%calving, isc, iec, jsc, jec)
Expand Down Expand Up @@ -233,6 +235,7 @@ subroutine ice_type_slow_reg_restarts(domain, CatIce, param_file, Ice, &
call register_restart_field(Ice_restart, 'flux_lw', Ice%flux_lw)
call register_restart_field(Ice_restart, 'lprec', Ice%lprec)
call register_restart_field(Ice_restart, 'fprec', Ice%fprec)
call register_restart_field(Ice_restart, 'seaice_melt', Ice%seaice_melt)
call register_restart_field(Ice_restart, 'runoff', Ice%runoff)
call register_restart_field(Ice_restart, 'calving', Ice%calving)
call register_restart_field(Ice_restart, 'runoff_hflx', Ice%runoff_hflx, mandatory=.false.)
Expand Down Expand Up @@ -338,6 +341,7 @@ subroutine dealloc_Ice_arrays(Ice)
if (associated(Ice%flux_lh)) deallocate(Ice%flux_lh)
if (associated(Ice%lprec)) deallocate(Ice%lprec)
if (associated(Ice%fprec)) deallocate(Ice%fprec)
if (associated(Ice%seaice_melt)) deallocate(Ice%seaice_melt)
if (associated(Ice%p_surf)) deallocate(Ice%p_surf)
if (associated(Ice%runoff)) deallocate(Ice%runoff)
if (associated(Ice%calving)) deallocate(Ice%calving)
Expand Down Expand Up @@ -419,6 +423,7 @@ subroutine Ice_public_type_chksum(mesg, Ice, check_fast, check_slow, check_rough
call chksum(Ice%flux_lh, trim(mesg)//" Ice%flux_lh")
call chksum(Ice%lprec, trim(mesg)//" Ice%lprec")
call chksum(Ice%fprec, trim(mesg)//" Ice%fprec")
call chksum(Ice%seaice_melt, trim(mesg)//" Ice%seaice_melt")
call chksum(Ice%p_surf, trim(mesg)//" Ice%p_surf")
call chksum(Ice%calving, trim(mesg)//" Ice%calving")
call chksum(Ice%runoff, trim(mesg)//" Ice%runoff")
Expand Down Expand Up @@ -680,6 +685,7 @@ subroutine ice_data_type_chksum(mesg, timestep, Ice, init_call)
chks = SIS_chksum(Ice%flux_lh ) ; if (root) write(outunit,100) 'ice_data_type%flux_lh ', chks
chks = SIS_chksum(Ice%lprec ) ; if (root) write(outunit,100) 'ice_data_type%lprec ', chks
chks = SIS_chksum(Ice%fprec ) ; if (root) write(outunit,100) 'ice_data_type%fprec ', chks
chks = SIS_chksum(Ice%seaice_melt ) ; if (root) write(outunit,100) 'ice_data_type%seaice_melt ', chks
chks = SIS_chksum(Ice%p_surf ) ; if (root) write(outunit,100) 'ice_data_type%p_surf ', chks
chks = SIS_chksum(Ice%runoff ) ; if (root) write(outunit,100) 'ice_data_type%runoff ', chks
chks = SIS_chksum(Ice%calving ) ; if (root) write(outunit,100) 'ice_data_type%calving ', chks
Expand Down