diff --git a/CONTRIBUTORS.md b/CONTRIBUTORS.md index 22a2e3be..9d443117 100644 --- a/CONTRIBUTORS.md +++ b/CONTRIBUTORS.md @@ -6,4 +6,5 @@ | doucla | Douglas Clark | UKCEH | 2026-07-03 | | james-bruten-mo | James Bruten | Met Office | 2026-07-09 | | maggiehendry | Maggie Hendry | Met Office | 2026-08-21 | +| DanCopsey | Dan Copsey | Met Office | 2026-08-24 | | Pierre-siddall | Pierre Siddall | Met Office | 2026-08-26 | diff --git a/src/initialisation/standalone/ancillaries/init_rivers_process_data_mod.F90 b/src/initialisation/standalone/ancillaries/init_rivers_process_data_mod.F90 index d2b50cfe..b1bc89fc 100644 --- a/src/initialisation/standalone/ancillaries/init_rivers_process_data_mod.F90 +++ b/src/initialisation/standalone/ancillaries/init_rivers_process_data_mod.F90 @@ -286,7 +286,7 @@ SUBROUTINE process_rivers_data( l_use_direction, rivers, rivers_data ) !---------------------------------------------------------------------------- IF ( l_shift_x .OR. l_reverse_y ) THEN CALL remap_ancil_control( nx_rivers, ny_rivers, rivers_dx, l_shift_x, & - l_reverse_y, rivers ) + l_reverse_y, rivers, l_use_land_fraction ) END IF !---------------------------------------------------------------------------- @@ -938,7 +938,7 @@ END SUBROUTINE process_rivers_coords !############################################################################## SUBROUTINE remap_ancil_control( nx_rivers, ny_rivers, rivers_dx, l_shift_x, & - l_reverse_y, rivers ) + l_reverse_y, rivers, l_use_land_fraction ) !------------------------------------------------------------------------------ ! Description: @@ -967,9 +967,11 @@ SUBROUTINE remap_ancil_control( nx_rivers, ny_rivers, rivers_dx, l_shift_x, & l_shift_x, & ! Flag indicating that a cyclic shift is applied in the x direction to ! river ancillaries. This can be TRUE only if the coordinate is longitude. - l_reverse_y + l_reverse_y, & ! Flag indicating if the order in the y direction of ancillary fields is ! to be reversed so that the coordinate is monotonically increasing. + l_use_land_fraction + ! Flag indicating that we are using land fractions !------------------------------------------------------------------------------ ! Array arguments with INTENT(IN) @@ -998,6 +1000,14 @@ SUBROUTINE remap_ancil_control( nx_rivers, ny_rivers, rivers_dx, l_shift_x, & var(ivar), rivers ) END DO +!------------------------------------------------------------------------------ +! Deal with the land fraction field +!------------------------------------------------------------------------------ +IF ( l_use_land_fraction ) THEN + CALL remap_ancil( nx_rivers, ny_rivers, rivers_dx, l_shift_x, l_reverse_y, & + 'land_fraction_2d', rivers ) +END IF + !------------------------------------------------------------------------------ ! Having used the river coordinates to remap all the required fields, we now ! change the coordinates themselves. This stage must be done after all other