Skip to content

Allow outputting of vertical height in W0 and W2H function spaces #165

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

Allow outputting of vertical height in W0 and W2H function spaces #165

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
03640d9
implement ability to write out initial height fields
tommbendall Jan 20, 2026
078ccec
Fix fringe cases, and use a loop to simplify code
tommbendall Jan 20, 2026
cfa0246
force update of various iodef files
tommbendall Jan 27, 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
1 change: 1 addition & 0 deletions CONTRIBUTORS.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,3 +5,4 @@
| james-bruten-mo | James Bruten | Met Office | 2025-12-09 |
| jennyhickson | Jenny Hickson | Met Office | 2025-12-10 |
| mo-marqh | mark Hedley | Met Office | 2025-12-11 |
| tommbendall | Thomas Bendall | Met Office | 2025-01-20 |
2 changes: 2 additions & 0 deletions applications/gungho_model/example/file_def_initial.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13,6 +13,8 @@
<field field_ref="init_exner"/>
<field field_ref="init_height_w3"/>
<field field_ref="init_height_wth"/>
<field field_ref="init_height_w2h"/>
<field field_ref="init_height_w0"/>
</field_group>
</file>

Expand Down
4 changes: 4 additions & 0 deletions applications/jedi_lfric_tests/example/iodef.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -140,6 +140,8 @@
<field id="exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="half_level_edge_grid" />
<field id="height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="node_grid" />
<!-- Processed diagnostics -->
<field id="processed__tot_col_dry_air_mass" name="tot_col_dry_air_mass" long_name="atmosphere_mass_of_air_per_unit_area" unit="kg m-2" domain_ref="face" />
<field id="processed__tot_col_m_v" name="tot_col_m_v" long_name="atmosphere_mass_content_of_water_vapor" unit="kg m-2" domain_ref="face" />
Expand Down Expand Up @@ -186,6 +188,8 @@
<field id="init_exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="init_height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="init_height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="init_height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" unit="m" grid_ref="node_grid" />
<field id="init_height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" unit="m" grid_ref="half_level_edge_grid" />
<field id="init_area_at_msl" name="cell_area" long_name="grid_cell_area" standard_name="grid_cell_area" unit="m2" domain_ref="face"/>

<field id="init_ls_rho" name="ls_rho" long_name="LS values for air_density" standard_name="air_density" unit="kg m-3" grid_ref="half_level_face_grid" />
Expand Down
2 changes: 2 additions & 0 deletions applications/jedi_lfric_tests/example_forecast/iodef.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -66,6 +66,8 @@
<field id="init_exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="init_height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="init_height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="init_height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" unit="m" grid_ref="node_grid" />
<field id="init_height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" unit="m" grid_ref="half_level_edge_grid" />
<field id="init_ageofair" name="ageofair" long_name="ageofair" unit="s" grid_ref="half_level_face_grid" />
<field id="init_area_at_msl" name="cell_area" long_name="grid_cell_area" standard_name="grid_cell_area" unit="m2" domain_ref="face"/>
</field_group>
Expand Down
4 changes: 4 additions & 0 deletions applications/jedi_lfric_tests/example_id_tlm_tests/iodef.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -234,6 +234,8 @@
<field id="exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="half_level_edge_grid" />
<field id="height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="node_grid" />
<!-- Processed diagnostics -->
<field id="processed__tot_col_dry_air_mass" name="tot_col_dry_air_mass" long_name="atmosphere_mass_of_air_per_unit_area" unit="kg m-2" domain_ref="face" />
<field id="processed__tot_col_m_v" name="tot_col_m_v" long_name="atmosphere_mass_content_of_water_vapor" unit="kg m-2" domain_ref="face" />
Expand Down Expand Up @@ -280,6 +282,8 @@
<field id="init_exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="init_height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="init_height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="init_height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" unit="m" grid_ref="node_grid" />
<field id="init_height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" unit="m" grid_ref="half_level_edge_grid" />
<field id="init_area_at_msl" name="cell_area" long_name="grid_cell_area" standard_name="grid_cell_area" unit="m2" domain_ref="face"/>

<field id="init_ls_rho" name="ls_rho" long_name="LS values for air_density" standard_name="air_density" unit="kg m-3" grid_ref="half_level_face_grid" />
Expand Down
4 changes: 4 additions & 0 deletions applications/jedi_lfric_tests/example_tlm_forecast_tl/iodef.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -234,6 +234,8 @@
<field id="exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="half_level_edge_grid" />
<field id="height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="node_grid" />
<!-- Processed diagnostics -->
<field id="processed__tot_col_dry_air_mass" name="tot_col_dry_air_mass" long_name="atmosphere_mass_of_air_per_unit_area" unit="kg m-2" domain_ref="face" />
<field id="processed__tot_col_m_v" name="tot_col_m_v" long_name="atmosphere_mass_content_of_water_vapor" unit="kg m-2" domain_ref="face" />
Expand Down Expand Up @@ -280,6 +282,8 @@
<field id="init_exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="init_height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="init_height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="init_height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" unit="m" grid_ref="node_grid" />
<field id="init_height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" unit="m" grid_ref="half_level_edge_grid" />
<field id="init_area_at_msl" name="cell_area" long_name="grid_cell_area" standard_name="grid_cell_area" unit="m2" domain_ref="face"/>

<field id="init_ls_rho" name="ls_rho" long_name="LS values for air_density" standard_name="air_density" unit="kg m-3" grid_ref="half_level_face_grid" />
Expand Down
4 changes: 4 additions & 0 deletions applications/jedi_lfric_tests/example_tlm_tests/iodef.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -234,6 +234,8 @@
<field id="exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="half_level_edge_grid" />
<field id="height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="node_grid" />
<!-- Processed diagnostics -->
<field id="processed__tot_col_dry_air_mass" name="tot_col_dry_air_mass" long_name="atmosphere_mass_of_air_per_unit_area" unit="kg m-2" domain_ref="face" />
<field id="processed__tot_col_m_v" name="tot_col_m_v" long_name="atmosphere_mass_content_of_water_vapor" unit="kg m-2" domain_ref="face" />
Expand Down Expand Up @@ -280,6 +282,8 @@
<field id="init_exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="init_height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="init_height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="init_height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" unit="m" grid_ref="node_grid" />
<field id="init_height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" unit="m" grid_ref="half_level_edge_grid" />
<field id="init_area_at_msl" name="cell_area" long_name="grid_cell_area" standard_name="grid_cell_area" unit="m2" domain_ref="face"/>

<field id="init_ls_rho" name="ls_rho" long_name="LS values for air_density" standard_name="air_density" unit="kg m-3" grid_ref="half_level_face_grid" />
Expand Down
2 changes: 2 additions & 0 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 @@ -43,6 +43,8 @@
<field id="w_in_w3_star" name="w_in_w3_star" long_name="physics_predictor_upward_air_velocity_at_cell_centres" unit="ms-1" grid_ref="half_level_face_grid" />
<field id="height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="half_level_face_grid" />
<field id="height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="full_level_face_grid" />
<field id="height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="half_level_edge_grid" />
<field id="height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="node_grid" />
<!-- First line extracts level 0, second line removes vertical coordinate, cannot do as one -->
<field id="zoomed__orography" field_ref="height_wth" grid_ref="full_level_0_face_zoom" />
<field id="orography" name="orography" long_name="surface_altitude" standard_name="surface_altitude" unit="m" domain_ref="face" > zoomed__orography </field>
Expand Down
2 changes: 2 additions & 0 deletions applications/lfric_atm/metadata/field_def_initial_diags.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -21,6 +21,8 @@
<field id="init_exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="init_height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="init_height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="init_height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" unit="m" grid_ref="node_grid" />
<field id="init_height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" unit="m" grid_ref="half_level_edge_grid" />
<field id="init_ageofair" name="ageofair" long_name="ageofair" unit="s" grid_ref="half_level_face_grid" />
<field id="init_area_at_msl" name="cell_area" long_name="grid_cell_area" standard_name="grid_cell_area" unit="m2" domain_ref="face"/>
</field_group>
Expand Down
10 changes: 7 additions & 3 deletions applications/linear_model/example/iodef.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -80,9 +80,9 @@
<field id="exner_in_wth" name="exner_in_wth" long_name="Exner pressure on wtheta points" standard_name="air_pressure" unit="1" grid_ref="full_level_face_grid" />
<field id="exner_wth_n" name="exner_wth_n" long_name="start_of_timestep_exner_pressure_at_cell_interfaces" unit="1" grid_ref="full_level_face_grid" />
<field id="theta_star" name="theta_star" long_name="air potential temperature: physics predictor" standard_name="air_potential_temperature" unit="K" grid_ref="full_level_face_grid" />
<field id="u_star1" name="u_star1" long_name="eastward_wind: physics predictor" standard_name="eastward_wind" unit="m s-1" grid_ref="half_level_edge_grid" />
<field id="u_star2" name="u_star2" long_name="northward_wind: physics predictor" standard_name="northward_wind" unit="m s-1" grid_ref="half_level_edge_grid" />
<field id="u_star3" name="u_star3" long_name="upward_air_velocity: physics predictor" standard_name="upward_air_velocity" unit="m s-1" grid_ref="full_level_face_grid" />
<field id="u_star1" name="u_star1" long_name="eastward_wind: physics predictor" standard_name="eastward_wind" unit="m s-1" grid_ref="half_level_edge_grid" />
<field id="u_star2" name="u_star2" long_name="northward_wind: physics predictor" standard_name="northward_wind" unit="m s-1" grid_ref="half_level_edge_grid" />
<field id="u_star3" name="u_star3" long_name="upward_air_velocity: physics predictor" standard_name="upward_air_velocity" unit="m s-1" grid_ref="full_level_face_grid" />
<field id="shear" name="shear" long_name="3D wind shear" unit="s-1" grid_ref="full_level_face_grid" />
<field id="u_in_w3" name="u_in_w3" long_name="physics U wind on W3 points" standard_name="eastward_wind" unit="ms-1" grid_ref="half_level_face_grid" />
<field id="v_in_w3" name="v_in_w3" long_name="physics V wind on W3 points" standard_name="northward_wind" unit="ms-1" grid_ref="half_level_face_grid" />
Expand Down Expand Up @@ -134,6 +134,8 @@
<field id="exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="half_level_edge_grid" />
<field id="height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="node_grid" />
<!-- Processed diagnostics -->
<field id="processed__tot_col_dry_air_mass" name="tot_col_dry_air_mass" long_name="atmosphere_mass_of_air_per_unit_area" unit="kg m-2" domain_ref="face" />
<field id="processed__tot_col_m_v" name="tot_col_m_v" long_name="atmosphere_mass_content_of_water_vapor" unit="kg m-2" domain_ref="face" />
Expand Down Expand Up @@ -180,6 +182,8 @@
<field id="init_exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="init_height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="init_height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="init_height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" unit="m" grid_ref="node_grid" />
<field id="init_height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" unit="m" grid_ref="half_level_edge_grid" />
<field id="init_area_at_msl" name="cell_area" long_name="grid_cell_area" standard_name="grid_cell_area" unit="m2" domain_ref="face"/>

<field id="init_ls_rho" name="ls_rho" long_name="LS values for air_density" standard_name="air_density" unit="kg m-3" grid_ref="half_level_face_grid" />
Expand Down
4 changes: 2 additions & 2 deletions dependencies.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -30,8 +30,8 @@ lfric_apps:
ref:

lfric_core:
source: git@github.com:MetOffice/lfric_core.git
ref: 2025.12.1
source: git@github.com:tommbendall/lfric_core.git
ref: TBendall/HeightW0

moci:
source: git@github.com:MetOffice/moci.git
Expand Down
2 changes: 2 additions & 0 deletions rose-stem/app/gungho_model/file/file_def_initial.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13,6 +13,8 @@
<field field_ref="init_exner"/>
<field field_ref="init_height_w3"/>
<field field_ref="init_height_wth"/>
<field field_ref="init_height_w2h"/>
<field field_ref="init_height_w0"/>
</field_group>
</file>

Expand Down
2 changes: 2 additions & 0 deletions rose-stem/app/jedi_forecast/file/iodef.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -66,6 +66,8 @@
<field id="init_exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="init_height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="init_height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="init_height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" unit="m" grid_ref="node_grid" />
<field id="init_height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" unit="m" grid_ref="half_level_edge_grid" />
<field id="init_ageofair" name="ageofair" long_name="ageofair" unit="s" grid_ref="half_level_face_grid" />
<field id="init_area_at_msl" name="cell_area" long_name="grid_cell_area" standard_name="grid_cell_area" unit="m2" domain_ref="face"/>
</field_group>
Expand Down
4 changes: 4 additions & 0 deletions rose-stem/app/jedi_id_tlm_tests/file/iodef.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -234,6 +234,8 @@
<field id="exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="half_level_edge_grid" />
<field id="height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="node_grid" />
<!-- Processed diagnostics -->
<field id="processed__tot_col_dry_air_mass" name="tot_col_dry_air_mass" long_name="atmosphere_mass_of_air_per_unit_area" unit="kg m-2" domain_ref="face" />
<field id="processed__tot_col_m_v" name="tot_col_m_v" long_name="atmosphere_mass_content_of_water_vapor" unit="kg m-2" domain_ref="face" />
Expand Down Expand Up @@ -280,6 +282,8 @@
<field id="init_exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
<field id="init_height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
<field id="init_height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
<field id="init_height_w0" name="height_w0" long_name="height_of_cell_vertices_above_earth_radius" unit="m" grid_ref="node_grid" />
<field id="init_height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" unit="m" grid_ref="half_level_edge_grid" />
<field id="init_area_at_msl" name="cell_area" long_name="grid_cell_area" standard_name="grid_cell_area" unit="m2" domain_ref="face"/>

<field id="init_ls_rho" name="ls_rho" long_name="LS values for air_density" standard_name="air_density" unit="kg m-3" grid_ref="half_level_face_grid" />
Expand Down
Loading
Loading