Skip to content

Lfric2lfric: interpolate W2 fields #238

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
ukmo-juan-castillo wants to merge 9 commits into
base: main
Choose a base branch
from
Open

Lfric2lfric: interpolate W2 fields #238

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
e7ca8f2
Issue 192: update partitioning namelist in lfric2lfric
ukmo-juan-castillo Jan 30, 2026
5d310ab
Revert previous change, intended for another branch
ukmo-juan-castillo Jan 30, 2026
1127c9b
Issue 192: update partitioning namelist in lfric2lfric
ukmo-juan-castillo Jan 30, 2026
7c61c2d
Revert previous change, intended for another branch
ukmo-juan-castillo Jan 30, 2026
aba975a
Merge branch 'main' into lfric2lfric_W2
ukmo-juan-castillo Feb 7, 2026
150e94e
Issue 60: lfric2lfric interpolation of W2 fields
ukmo-juan-castillo Feb 10, 2026
7e1ef9c
Issue 60: update missing KGOs
ukmo-juan-castillo Feb 11, 2026
5b22236
Issue 60: update KGOs
ukmo-juan-castillo Feb 11, 2026
0e7f71a
Issue 60: update KGO
ukmo-juan-castillo Feb 12, 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
117 changes: 111 additions & 6 deletions applications/lfric2lfric/source/driver/lfric2lfric_regrid_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -14,10 +14,17 @@ module lfric2lfric_regrid_mod
use field_collection_mod, only: field_collection_type
use field_collection_iterator_mod, only: &
field_collection_iterator_type
use fs_continuity_mod, only: W2, W3, Wtheta
use function_space_collection_mod, only: function_space_collection
use function_space_mod, only: function_space_type
use interpolation_alg_mod, only: interp_w2_to_w3wth_alg, &
interp_w3wth_to_w2_alg
use log_mod, only: log_event, &
log_level_info, &
log_scratch_space
use mesh_collection_mod, only: mesh_collection
use model_clock_mod, only: model_clock_type
use namelist_mod, only: namelist_type

!------------------------------------
! lfric2lfric modules
Expand Down Expand Up @@ -66,12 +73,57 @@ subroutine lfric2lfric_regrid( modeldb, oasis_clock, &

type(field_collection_iterator_type) :: iter

class(field_parent_type), pointer :: field => null()
type(field_type), pointer :: field_src => null()
type(field_type), pointer :: field_dst => null()
class(field_parent_type), pointer :: field => null()
type(field_type), pointer :: field_src => null()
type(field_type), pointer :: field_dst => null()

type(field_type), target :: u_in_w3_src, u_in_w3_dst
type(field_type), target :: v_in_w3_src, v_in_w3_dst
type(field_type), target :: w_in_wth_src, w_in_wth_dst
type(field_type), pointer :: field_src_ptr, field_dst_ptr

integer(kind=i_def) :: fs
type(function_space_type), pointer :: fs_w3_src, fs_w3_dst
type(function_space_type), pointer :: fs_wth_src, fs_wth_dst

type(namelist_type), pointer :: lfric2lfric_nml
type(namelist_type), pointer :: finite_element_nml

character(len=str_def) :: mesh_names(2)
integer(kind=i_def) :: element_order_h
integer(kind=i_def) :: element_order_v

character(len=str_def) :: field_name

integer(kind=i_def), parameter :: dst = 1
integer(kind=i_def), parameter :: src = 2


! Obtain namelist parameters
lfric2lfric_nml => modeldb%configuration%get_namelist('lfric2lfric')
finite_element_nml => modeldb%configuration%get_namelist('finite_element')

call lfric2lfric_nml%get_value( 'destination_mesh_name', &
mesh_names(dst) )
call lfric2lfric_nml%get_value( 'source_mesh_name', &
mesh_names(src) )
call finite_element_nml%get_value( 'element_order_h', element_order_h)
call finite_element_nml%get_value( 'element_order_v', element_order_v)

! Function spaces for creating temporary fields used in W2 interpolation
fs_w3_src => function_space_collection%get_fs( &
mesh_collection%get_mesh(mesh_names(src)), &
element_order_h, element_order_v, W3)
fs_wth_src => function_space_collection%get_fs( &
mesh_collection%get_mesh(mesh_names(src)), &
element_order_h, element_order_v, Wtheta)
fs_w3_dst => function_space_collection%get_fs( &
mesh_collection%get_mesh(mesh_names(dst)), &
element_order_h, element_order_v, W3)
fs_wth_dst => function_space_collection%get_fs( &
mesh_collection%get_mesh(mesh_names(dst)), &
element_order_h, element_order_v, Wtheta)

! Main loop over fields to be processed
call iter%initialise(source_fields)
do
Expand All @@ -87,10 +139,42 @@ subroutine lfric2lfric_regrid( modeldb, oasis_clock, &
trim(field_name)
call log_event(log_scratch_space, log_level_info)

! Convert W2 fields to a set of W3 and Wtheta fields
fs = field_src%which_function_space()
if (fs == W2) then
call u_in_w3_src%initialise(vector_space=fs_w3_src, &
name="u_in_w3_src")
call v_in_w3_src%initialise(vector_space=fs_w3_src, &
name="v_in_w3_src")
call w_in_wth_src%initialise(vector_space=fs_wth_src, &
name="w_in_wth_src")
call u_in_w3_dst%initialise(vector_space=fs_w3_dst, &
name="u_in_w3_dst")
call v_in_w3_dst%initialise(vector_space=fs_w3_dst, &
name="v_in_w3_dst")
call w_in_wth_dst%initialise(vector_space=fs_wth_dst, &
name="w_in_wth_dst")

call interp_w2_to_w3wth_alg(field_src, u_in_w3_src, &
v_in_w3_src, w_in_wth_src)
end if

! Regrid source field depending on regrid method
select case (regrid_method)
case (regrid_method_map)
call lfric2lfric_map_regrid(field_dst, field_src)
if (fs == W2) then
field_src_ptr => u_in_w3_src
field_dst_ptr => u_in_w3_dst
call lfric2lfric_map_regrid(field_dst_ptr, field_src_ptr)
field_src_ptr => v_in_w3_src
field_dst_ptr => v_in_w3_dst
call lfric2lfric_map_regrid(field_dst_ptr, field_src_ptr)
field_src_ptr => w_in_wth_src
field_dst_ptr => w_in_wth_dst
call lfric2lfric_map_regrid(field_dst_ptr, field_src_ptr)
else
call lfric2lfric_map_regrid(field_dst, field_src)
end if

case (regrid_method_lfric2lfric)
write(log_scratch_space, '(A)') &
Expand All @@ -101,10 +185,31 @@ subroutine lfric2lfric_regrid( modeldb, oasis_clock, &

case (regrid_method_oasis)
#ifdef MCT
call lfric2lfric_oasis_regrid(modeldb, oasis_clock, &
field_dst, field_src)
if (fs == W2) then
field_src_ptr => u_in_w3_src
field_dst_ptr => u_in_w3_dst
call lfric2lfric_oasis_regrid(modeldb, oasis_clock, &
field_dst_ptr, field_src_ptr)
field_src_ptr => v_in_w3_src
field_dst_ptr => v_in_w3_dst
call lfric2lfric_oasis_regrid(modeldb, oasis_clock, &
field_dst_ptr, field_src_ptr)
field_src_ptr => w_in_wth_src
field_dst_ptr => w_in_wth_dst
call lfric2lfric_oasis_regrid(modeldb, oasis_clock, &
field_dst_ptr, field_src_ptr)
else
call lfric2lfric_oasis_regrid(modeldb, oasis_clock, &
field_dst, field_src)
end if
#endif
end select

! Rebuild the W2 fields from a set of W3 and Wtheta fields
if (fs == W2) then
call interp_w3wth_to_w2_alg(field_dst, u_in_w3_dst, &
v_in_w3_dst, w_in_wth_dst)
end if
end do

end subroutine lfric2lfric_regrid
Expand Down
2 changes: 1 addition & 1 deletion ...c/azspice/checksum_lfric2lfric_oasis_ral_seuk-C32_lam_MG_azspice_gnu_fast-debug-64bit.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@ Inner product checksum soil_moist_sat = 7FF0000000000000
Inner product checksum lw_up_toa_rtsi = 416FAC72807CC071
Inner product checksum acc_sol_bc = 0
Inner product checksum lw_up_tile_rts = 41C28987D17C61C3
Inner product checksum u = 0
Inner product checksum u = 447D4D3FAF1B3E00
Inner product checksum lbc_u = 0
Inner product checksum ccw = 0
Inner product checksum tile_canopy_water = 443A39FEDE328EF1
Expand Down
2 changes: 1 addition & 1 deletion ...c2lfric/azspice/checksum_lfric2lfric_ral_seuk-C32_lam_MG_azspice_gnu_fast-debug-64bit.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@ Inner product checksum soil_moist_sat = 7FF0000000000000
Inner product checksum lw_up_toa_rtsi = 416FAC7209B088EF
Inner product checksum acc_sol_bc = 0
Inner product checksum lw_up_tile_rts = 41C28987BDD836EC
Inner product checksum u = 0
Inner product checksum u = 42083253875B88C4
Inner product checksum lbc_u = 0
Inner product checksum ccw = 0
Inner product checksum tile_canopy_water = 443A39FEDBCA7733
Expand Down
2 changes: 1 addition & 1 deletion ...c2lfric/ex1a/checksum_lfric2lfric_oasis_ral_seuk-C32_lam_MG_ex1a_cce_fast-debug-64bit.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@ Inner product checksum soil_moist_sat = 7FF0000000000000
Inner product checksum lw_up_toa_rtsi = 416FAC72807CC06E
Inner product checksum acc_sol_bc = 0
Inner product checksum lw_up_tile_rts = 41C28987D17C61B8
Inner product checksum u = 0
Inner product checksum u = 425E601913543247
Inner product checksum lbc_u = 0
Inner product checksum ccw = 0
Inner product checksum tile_canopy_water = 443A39FEDE328EC9
Expand Down
2 changes: 1 addition & 1 deletion ...s/lfric2lfric/ex1a/checksum_lfric2lfric_ral_seuk-C32_lam_MG_ex1a_cce_fast-debug-64bit.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@ Inner product checksum soil_moist_sat = 7FF0000000000000
Inner product checksum lw_up_toa_rtsi = 416FAC7209B05CCB
Inner product checksum acc_sol_bc = 0
Inner product checksum lw_up_tile_rts = 41C28987BDD7C02E
Inner product checksum u = 0
Inner product checksum u = 42083253875A2B6D
Inner product checksum lbc_u = 0
Inner product checksum ccw = 0
Inner product checksum tile_canopy_water = 443A39FEDBCBC772
Expand Down