Skip to content

Diagnostics required by lfric2um for PS49 #658

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
iboutle wants to merge 21 commits into
base: main
Choose a base branch
from
Open

Diagnostics required by lfric2um for PS49 #658

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
21 commits
Select commit Hold shift + click to select a range
88faf33
new diagnostics needed
iboutle May 29, 2026
68ca1a7
add snow prognostics
iboutle Jun 4, 2026
946049d
forgot doxygen comments
iboutle Jun 4, 2026
0090382
sea and sea-ice temperature
iboutle Jul 21, 2026
51c312b
bugfixes
iboutle Jul 22, 2026
417d266
fix sea-ice area fraction for coupled models
iboutle Jul 23, 2026
ecb42cc
sea_ice surface temperature
iboutle Jul 23, 2026
829edc8
Merge branch 'stable' into new_diags
iboutle Jul 23, 2026
3d06e07
xios doesn't work for sea-ice fraction
iboutle Jul 24, 2026
2930e26
Add open sea sensible heat flux and evaporation diagnostics
Aug 10, 2026
b451ebb
fix names
iboutle Aug 10, 2026
2ba3ad8
add to OS file
iboutle Aug 10, 2026
0f9215d
Add sea-ice top melt heat flux diagnostic
Aug 11, 2026
ac57990
name updates
iboutle Aug 11, 2026
b484e19
fixes to code
iboutle Aug 11, 2026
3c340f2
Add land-mean surface heat and latent heat flux diagnostics
Aug 11, 2026
db25b61
change order
iboutle Aug 11, 2026
8ace74f
Add sea and sea-ice mean surface heat and latent heat flux diagnostics
Aug 11, 2026
14cd336
Derive the number of sea and sea-ice tiles in jules_control_init
Aug 11, 2026
6211ca4
Add the open sea roughness length diagnostic
Aug 11, 2026
15088db
Add the mean over ice sea-ice thickness and temperature diagnostics
Aug 11, 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
34 changes: 29 additions & 5 deletions applications/lfric_atm/metadata/field_def_diags.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -500,18 +500,19 @@
<!-- Jules surface exchange diagnostics -->
<field id="surface__tile_fraction" name="tile_fraction" long_name="surface_tile_fraction" unit="1" domain_ref="face" axis_ref="surface_tiles" />
<!-- Extracts tile types for sea and sea-ice, removes pseudo-level coordinate, then calculate sea-ice fraction -->
<field id="zoomed__sea_ice_tile" field_ref="surface__tile_fraction" domain_ref="face" axis_ref="surface_tiles_sea_ice_zoom" />
<field id="surface__sea_ice_tile" name="sea_ice_total_fraction" long_name="sea_ice_total_fraction_of_grid_box" unit="1" domain_ref="face"/>
<field id="zoomed__sea_tile" field_ref="surface__tile_fraction" domain_ref="face" axis_ref="surface_tiles_sea_zoom" />
<field id="surface__sea_ice_tile" domain_ref="face"> zoomed__sea_ice_tile </field>
<field id="surface__sea_tile" domain_ref="face"> zoomed__sea_tile </field>
<field id="surface__sea_tile" name="sea_tile_fraction" long_name="sea_tile_fraction_of_grid_box" unit="1" domain_ref="face"> zoomed__sea_tile </field>
<field id="surface__sea_ice_fraction" name="sea_ice_fraction" long_name="sea_ice_area_fraction" standard_name="sea_ice_area_fraction" unit="1" domain_ref="face"> surface__sea_ice_tile / (surface__sea_ice_tile + surface__sea_tile + 1.0e-16) </field>
<field id="zoomed__urban_tile" field_ref="surface__tile_fraction" domain_ref="face" axis_ref="surface_tiles_urban_zoom"/>
<field id="surface__urban_tile" name="urban_tile_fraction" long_name="urban_tile_fraction" domain_ref="face"> zoomed__urban_tile </field>
<field id="surface__urban_tile" name="urban_tile_fraction" long_name="urban_tile_fraction" unit="1" domain_ref="face"> zoomed__urban_tile </field>
<field id="zoomed__land_ice_tile" field_ref="surface__tile_fraction" domain_ref="face" axis_ref="surface_tiles_land_ice_zoom"/>
<field id="surface__land_ice_tile" name="land_ice_tile_fraction" long_name="land_ice_tile_fraction" domain_ref="face"> zoomed__land_ice_tile </field>
<field id="surface__land_ice_tile" name="land_ice_tile_fraction" long_name="land_ice_tile_fraction" unit="1" domain_ref="face"> zoomed__land_ice_tile </field>
<field id="surface__land_fraction" name="land_fraction" long_name="land_area_fraction" standard_name="land_area_fraction" unit="1" domain_ref="face"> 1.0 - (surface__sea_ice_tile + surface__sea_tile) </field>
<field id="surface__z0m" name="z0m" long_name="surface_roughness_length" unit="m" domain_ref="face" />
<field id="surface__z0m_tile" name="z0m_tile" long_name="surface_roughness_length_on_tiles" unit="m" domain_ref="face" axis_ref="surface_tiles" />
<field id="zoomed__z0m_sea" field_ref="surface__z0m_tile" domain_ref="face" axis_ref="surface_tiles_sea_zoom" />
<field id="surface__z0m_sea" name="z0m_sea" long_name="sea_surface_roughness_length" unit="m" domain_ref="face"> zoomed__z0m_sea </field>
<field id="surface__z0m_eff" name="z0m_eff" long_name="effective_surface_roughness_length_for_momentum_in_air" unit="m" domain_ref="face" />
<field id="surface__z0h_eff" name="z0h_eff" long_name="effective_surface_roughness_length_for_heat_in_air" unit="m" domain_ref="face" />
<field id="surface__canopy_height" name="canopy_height" long_name="canopy_height_of_vegetated_tiles" unit="m" domain_ref="face" axis_ref="plant_func_types" />
Expand All @@ -531,8 +532,11 @@
<field id="surface__pseudotaux" name="pseudotaux" long_name="surface_eastward_wind_pseudo_stress" unit="N m-2" domain_ref="face" />
<field id="surface__pseudotauy" name="pseudotauy" long_name="surface_northward_wind_pseudo_stress" unit="N m-2" domain_ref="face" />
<field id="surface__tile_temperature" name="tile_temperature" long_name="surface_tile_temperature" unit="K" domain_ref="face" axis_ref="surface_tiles" />
<field id="zoomed__sea_temperature" field_ref="surface__tile_temperature" domain_ref="face" axis_ref="surface_tiles_sea_zoom" />
<field id="surface__sea_temperature" name="sea_surface_temperature" long_name="sea_surface_temperature" unit="K" domain_ref="face"> zoomed__sea_temperature </field>
<field id="surface__grid_surface_temperature" name="grid_surface_temperature" long_name="grid_surface_temperature" unit="K" domain_ref="face" />
<field id="surface__land_surface_temperature" name="land_surface_temperature" long_name="land_surface_temperature" unit="K" domain_ref="face" />
<field id="surface__sea_ice_surface_temperature" name="sea_ice_surface_temperature" long_name="sea_ice_surface_temperature" unit="K" domain_ref="face" />
<field id="surface__t1p5m_surft" name="t1p5m_surft" long_name="land_tile_screen_level_temperature" unit="K" domain_ref="face" axis_ref="surface_tiles" />
<field id="surface__t1p5m" name="t1p5m" long_name="temperature_at_screen_level" unit="K" domain_ref="face" />
<field id="surface__t1p5m_ssi" name="t1p5m_ssi" long_name="temperature_at_screen_level_over_sea_and_sea_ice" unit="K" domain_ref="face" />
Expand All @@ -558,9 +562,22 @@
<field id="surface__visibility_no_precip" name="visibility_no_precip" long_name="visibility_excluding precipitation_at_screen_level" standard_name="visibility_in_air" unit="m" domain_ref="face" />
<field id="surface__visibility_with_precip" name="visibility_with_precip" long_name="visibility_including_precipitation_at_screen_level" standard_name="visibility_in_air" unit="m" domain_ref="face" />
<field id="surface__tile_heat_flux" name="tile_heat_flux" long_name="surface_sensible_heat_flux" unit="W m-2" domain_ref="face" axis_ref="surface_tiles" />
<!-- First line extracts the open-sea tile, second line removes pseudo-level coordinate, cannot do as one -->
<!-- The tiled flux is per unit open water, so the third line rescales it to a mean over the whole sea portion of the grid-box -->
<field id="zoomed__sea_heat_flux" field_ref="surface__tile_heat_flux" domain_ref="face" axis_ref="surface_tiles_sea_zoom" />
<field id="surface__sea_tile_heat_flux" name="sea_tile_heat_flux" long_name="surface_sensible_heat_flux_of_sea_tile" unit="W m-2" domain_ref="face"> zoomed__sea_heat_flux </field>
<field id="surface__sea_mean_heat_flux" name="sea_mean_heat_flux" long_name="surface_sensible_heat_flux_sea_area_mean" unit="W m-2" domain_ref="face"> surface__sea_tile_heat_flux * (1.0 - surface__sea_ice_fraction) </field>
<field id="surface__land_mean_heat_flux" name="land_mean_heat_flux" long_name="surface_sensible_heat_flux_land_area_mean" unit="W m-2" domain_ref="face" />
<field id="surface__ssi_mean_heat_flux" name="ssi_mean_heat_flux" long_name="surface_sensible_heat_flux_mean_over_sea_seaice" unit="W m-2" domain_ref="face" />
<field id="surface__tile_moisture_flux" name="tile_moisture_flux" long_name="surface_moisture_flux" unit="kg m-2 s-1" domain_ref="face" axis_ref="surface_tiles" />
<!-- As above, extract the open-sea tile, remove the pseudo-level coordinate, then rescale to a sea mean -->
<field id="zoomed__sea_moisture_flux" field_ref="surface__tile_moisture_flux" domain_ref="face" axis_ref="surface_tiles_sea_zoom" />
<field id="surface__sea_tile_moisture_flux" name="sea_tile_moisture_flux" long_name="surface_moisture_flux_of_sea_tile" unit="kg m-2 s-1" domain_ref="face"> zoomed__sea_moisture_flux </field>
<field id="surface__sea_mean_moisture_flux" name="sea_mean_moisture_flux" long_name="surface_moisture_flux_sea_area_mean" unit="kg m-2 s-1" domain_ref="face"> surface__sea_tile_moisture_flux * (1.0 - surface__sea_ice_fraction) </field>
<field id="surface__latent_heat" name="latent_heat" long_name="surface_upward_latent_heat_flux" unit="W m-2" domain_ref="face" axis_ref="surface_tiles" />
<field id="surface__grid_latent_heat" name="grid_latent_heat" long_name="grid_surface_upward_latent_heat_flux" unit="W m-2" domain_ref="face" />
<field id="surface__land_mean_latent_heat" name="land_mean_latent_heat" long_name="surface_upward_latent_heat_flux_land_area_mean" unit="W m-2" domain_ref="face" />
<field id="surface__ssi_mean_latent_heat" name="ssi_mean_latent_heat" long_name="surface_upward_latent_heat_flux_mean_over_sea_seaice" unit="W m-2" domain_ref="face" />
<field id="surface__canopy_water" name="canopy_water_content" long_name="canopy_water_amount" unit="kg m-2" domain_ref="face" axis_ref="surface_tiles" />
<field id="surface__grid_canopy_water" name="grid_canopy_water_content" long_name="grid_canopy_water_amount" unit="kg m-2" domain_ref="face" />
<field id="surface__throughfall" name="throughfall_rate" long_name="canopy_throughfall_flux" unit="kg m-2 s-1" domain_ref="face" axis_ref="surface_tiles" />
Expand Down Expand Up @@ -610,6 +627,7 @@
<field id="snow__snow_depth" name="snow_depth" long_name="snow_depth" unit="m" domain_ref="face" axis_ref="surface_tiles" />
<field id="snow__tile_snow_rgrain" name="tile_snow_rgrain" long_name="tile_snow_grain_size" unit="microns" domain_ref="face" axis_ref="surface_tiles" />
<field id="snow__snowpack_density" name="snowpack_density" long_name="snowpack_density" unit="kg m-3" domain_ref="face" axis_ref="surface_tiles" />
<field id="snow__snow_under_canopy" name="snow_under_canopy" long_name="snow_under_canopy_on_tiles" unit="m" domain_ref="face" axis_ref="surface_tiles" />
<field id="snow__snow_layer_temp" name="snow_layer_temp" long_name="snow_layer_temperature" unit="K" domain_ref="face" axis_ref="snow_layers_and_tiles" />
<field id="snow__snow_layer_thickness" name="snow_layer_thickness" long_name="snow_layer_thickness" unit="m" domain_ref="face" axis_ref="snow_layers_and_tiles" />
<field id="snow__snow_layer_ice_mass" name="snow_layer_ice_mass" long_name="snow_layer_ice_mass" unit="kg m-2" domain_ref="face" axis_ref="snow_layers_and_tiles" />
Expand All @@ -626,6 +644,9 @@
<field id="processed__sea_ice_pensolar" name="sea_ice_pensolar_gbm" long_name="sea_ice_pensolar_grid_box_mean" unit="W m-2" domain_ref="face" axis_ref="sea_ice_categories">seaice__sea_ice_pensolar * seaice__sea_ice_fraction</field>
<field id="seaice__melt_pond_fraction" name="melt_pond_fraction" long_name="sea ice melt pond fraction" unit="1" domain_ref="face" axis_ref="sea_ice_categories" />
<field id="seaice__melt_pond_depth" name="melt_pond_depth" long_name="sea ice melt pond depth" unit="m" domain_ref="face" axis_ref="sea_ice_categories" />
<field id="seaice__ssi_mean_topmelt_heat_flux" name="ssi_mean_topmelt_heat_flux" long_name="sea_ice_surface_melt_heat_flux_mean_over_sea_seaice" unit="W m-2" domain_ref="face" />
<field id="seaice__ice_mean_thickness" name="ice_mean_thickness" long_name="sea_ice_thickness_mean_over_ice" unit="m" domain_ref="face" />
<field id="seaice__ice_mean_temperature" name="ice_mean_temperature" long_name="sea_ice_temperature_mean_over_ice" unit="K" domain_ref="face" />
<!-- Jules soil diagnostics -->
<field id="soil__surface_runoff" name="surface_runoff" long_name="surface_runoff_flux" unit="kg m-2 s-1" domain_ref="face" />
<field id="soil__sub_surface_runoff" name="sub_surface_runoff" long_name="subsurface_runoff_flux" unit="kg m-2 s-1" domain_ref="face" />
Expand All @@ -641,6 +662,8 @@
<field id="soil__soil_moist_wilt" name="soil_moist_wilt" long_name="soil_moisture_at_wilting_point" unit="m3 m-3" domain_ref="face" />
<field id="soil__soil_moist_crit" name="soil_moist_crit" long_name="soil_moisture_at_critical_point" unit="m3 m-3" domain_ref="face" />
<field id="soil__soil_moist_sat" name="soil_moist_sat" long_name="soil_moisture_at_saturation" unit="m3 m-3" domain_ref="face" />
<field id="soil__soil_sat_frac" name="soil_sat_frac" long_name="soil_saturated_fraction" unit="1" domain_ref="face" />
<field id="soil__wetness_under_soil" name="wetness_under_soil" long_name="wetness_below_soil_column" unit="1" domain_ref="face" />
<field id="soil__inland_basin_flow" name="inland_basin_flow" long_name="Inland basin flow" unit="kg m-2 s-1" domain_ref="face" />
<!-- Boundary layer turbulence diagnostics -->
<field id="turbulence__dtheta_bl" name="dtheta_bl" long_name="potential_temperature_increment_from_bl_scheme" unit="K" grid_ref="full_level_face_grid" />
Expand Down Expand Up @@ -1023,6 +1046,7 @@
<field id="convection__lowest_cv_top" name="lowest_cv_top" long_name="model_level_number_at_lowest_conv_cloud_top" unit="1" domain_ref="face" />
<field id="convection__cv_base" name="cv_base" long_name="model_level_number_at_conv_cloud_base" unit="1" domain_ref="face" />
<field id="convection__cv_top" name="cv_top" long_name="model_level_number_at_conv_cloud_top" unit="1" domain_ref="face" />
<field id="snow__n_snow_layers" name="n_snow_layers" long_name="number_of_snow_layers" unit="1" domain_ref="face" axis_ref="surface_tiles" />
<field id="processed__trop_level" name="trop_level" long_name="tropopause_level" unit="1" grid_ref="var_face" />
</field_group>
</field_definition>
8 changes: 7 additions & 1 deletion interfaces/jules_interface/source/algorithm/jules_extra_alg_mod.x90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,6 +3,9 @@
! The file LICENCE, distributed with this code, contains details of the terms
! under which the code may be used.
!-------------------------------------------------------------------------------
! Some of the content of this file has been produced with the assistance of
! Anthropic Claude Opus 5 (Claude Code).
Comment thread
iboutle marked this conversation as resolved.
!-------------------------------------------------------------------------------
!> @brief Interface to the JULES surf_couple_extra scheme

module jules_extra_alg_mod
Expand Down Expand Up @@ -285,6 +288,7 @@ contains
call output_diags_for_jules_snow(tile_fraction, tile_snow_mass, &
snowice_melt, &
snow_depth, tile_snow_rgrain, &
snow_under_canopy, n_snow_layers, &
snowpack_density, snow_layer_temp, &
snow_layer_thickness, snow_layer_ice_mass,&
snow_layer_liq_mass, snow_layer_rgrain, &
Expand All @@ -299,10 +303,12 @@ contains
water_table, soil_moist_wilt, &
soil_moist_crit, soil_moist_sat, &
grid_canopy_water,grid_throughfall, &
soil_sat_frac, wetness_under_soil, &
inland_basin_flow)

! Output sea ice diagnostics
call output_diags_for_jules_seaice(tile_fraction, sea_ice_thickness, &
call output_diags_for_jules_seaice(tile_fraction, snowice_melt, &
sea_ice_thickness, &
sea_ice_temperature, sea_ice_conductivity, &
sea_ice_pensolar, melt_pond_fraction, &
melt_pond_depth)
Expand Down
Loading
Loading