Skip to content

Configure LFRic to receive meltpond fields as grid box mean #431

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

Open
Tim Graham (timgraham-Met) wants to merge 6 commits into
base: main
Choose a base branch
from
Open

Configure LFRic to receive meltpond fields as grid box mean #431

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
6e40a0a
Copied changes from old branch version
timgraham-Met Apr 8, 2026
99f31b2
Update NEMO to make use of new functionality
timgraham-Met Apr 14, 2026
ad09584
Added Alex to contributors list
timgraham-Met Apr 15, 2026
a29b200
Updated KGOs
timgraham-Met Apr 15, 2026
c2502ab
Merge branch 'main' into meltponds_as_gbm
timgraham-Met Apr 15, 2026
c6a4618
Update CONTRIBUTORS.md
timgraham-Met Apr 15, 2026
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
11 changes: 11 additions & 0 deletions interfaces/coupled_interface/source/kernel/process_o2a_kernel_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -134,6 +134,7 @@ subroutine process_o2a_code(nlayers, sea_surf_temp, sea_surf_temp_pert, &
integer(kind=i_def) :: i
real(kind=r_def) :: total_sea_ice_fraction
real(kind=r_def) :: recip_sea_ice_fraction
real(kind=r_def) :: recip_gbm_melt_pond_fraction

if (l_couple_ocean) then
if (sea_surf_temp(map_ocn(1))>1.0_r_def) then
Expand Down Expand Up @@ -202,6 +203,7 @@ subroutine process_o2a_code(nlayers, sea_surf_temp, sea_surf_temp_pert, &
(therm_cond_sice/therm_cond_sice_snow)

! Undo sea ice fraction scaling done on ocean side of coupler
! Melt pond thickness needs to be divided by grid box mean melt pond fraction
if ( sea_ice_fraction(map_ice(1) + i) > 0.0_r_def ) then
recip_sea_ice_fraction = 1.0_r_def / sea_ice_fraction(map_ice(1) + i)
sea_ice_thickness(map_ice(1) + i) = sea_ice_thickness(map_ice(1) + i) &
Expand All @@ -212,6 +214,15 @@ subroutine process_o2a_code(nlayers, sea_surf_temp, sea_surf_temp_pert, &
* recip_sea_ice_fraction
sea_ice_conductivity(map_ice(1)+i) = sea_ice_conductivity(map_ice(1)+i)&
* recip_sea_ice_fraction

if (melt_pond_fraction(map_ice(1) + i) > 0.0_r_def ) then
recip_gbm_melt_pond_fraction = 1.0_r_def / &
melt_pond_fraction(map_ice(1) + i)
melt_pond_depth(map_ice(1)+i) = melt_pond_depth(map_ice(1)+i) &
* recip_gbm_melt_pond_fraction
end if
melt_pond_fraction(map_ice(1)+i) = melt_pond_fraction(map_ice(1)+i) &
* recip_sea_ice_fraction
end if

! Apply bounds to sea ice layer temperatures
Expand Down
2 changes: 1 addition & 1 deletion rose-stem/app/fcm_make_ocean/rose-app.conf
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -28,7 +28,7 @@ nemo_base=branches/UKMO/NEMO_4.0.4_mirror
nemo_path_incl=src/OCE cfgs/SHARED src/ICE cfgs/METO_GO
nemo_project_location=fcm:nemo.xm/NEMO
nemo_rev=14075
nemo_sources=branches/UKMO/NEMO_4.0.4_GOSI9_package@16897
nemo_sources=branches/UKMO/NEMO_4.0.4_GOSI9_package@17142
nemo_version_file=nemo4.0-version.cfg
netcdf_inc_path=.
netcdf_lib_path=.
Expand Down
2 changes: 1 addition & 1 deletion rose-stem/app/lfric_coupled_nemo/rose-app.conf
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -629,7 +629,7 @@ sn_snd_cond='weighted ice','yes','','',''
sn_snd_crt='mixed oce-ice','no','spherical','eastward-northward','U,V'
sn_snd_crtw='none','no','','','U,V'
sn_snd_ifrac='none','no','','',''
sn_snd_mpnd='ice only','yes','','',''
sn_snd_mpnd='weighted ice','yes','','',''
sn_snd_sstfrz='none','no','','',''
sn_snd_temp='oce only','yes','','',''
sn_snd_thick='weighted ice and snow','yes','','',''
Expand Down
14 changes: 7 additions & 7 deletions ...kgos/lfric_coupled/ex1a/checksum_lfric_coupled_nwp_gal9-C48_ex1a_cce_fast-debug-64bit.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,9 +1,9 @@
Inner product checksum rho = 411B540A101D14D3
Inner product checksum theta = 42735B179026DDD2
Inner product checksum u = 456F4424668465D2
Inner product checksum mr1 = 403684A084256805
Inner product checksum mr2 = 3F3FD0C0751D3CC1
Inner product checksum mr3 = 3F012DD8506F6E14
Inner product checksum mr4 = 3F3693B42B44FFA1
Inner product checksum rho = 411B540D01244EA2
Inner product checksum theta = 42735B1FDDE2D7A5
Inner product checksum u = 456F41816571CCD1
Inner product checksum mr1 = 403684FBC201727A
Inner product checksum mr2 = 3F3F886F3FCF3DBA
Inner product checksum mr3 = 3F006A5ED4A55E92
Inner product checksum mr4 = 3F36AE3943C84DC6
Inner product checksum mr5 = 0
Inner product checksum mr6 = 0
Loading