Skip to content

Commit ec251c0

Browse files
Merge pull request #1885 from MetOffice/1759_branch_off_import_nimrod_branch
Enable transverse mercator plotting in spatial plot operator
2 parents fe82709 + 318c5fb commit ec251c0

File tree

1 file changed

+10
-0
lines changed

1 file changed

+10
-0
lines changed

src/CSET/operators/plot.py

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -407,6 +407,16 @@ def _setup_spatial_map(
407407
central_rotated_longitude=central_longitude,
408408
)
409409
crs = projection
410+
elif isinstance(coord_system, iris.coord_systems.TransverseMercator):
411+
# Define Transverse Mercator projection for TM inputs.
412+
projection = ccrs.TransverseMercator(
413+
central_longitude=coord_system.longitude_of_central_meridian,
414+
central_latitude=coord_system.latitude_of_projection_origin,
415+
false_easting=coord_system.false_easting,
416+
false_northing=coord_system.false_northing,
417+
scale_factor=coord_system.scale_factor_at_central_meridian,
418+
)
419+
crs = projection
410420
else:
411421
# Define regular map projection for non-rotated pole inputs.
412422
# Alternatives might include e.g. for global model outputs:

0 commit comments

Comments
 (0)