From 4d4ee6824fb2c469c086eb5fd48ded1344d93c0f Mon Sep 17 00:00:00 2001 From: Manuel Franco Date: Mon, 28 Apr 2025 16:08:44 +0200 Subject: [PATCH] Skip rotation of geo_coord_nod2D for unrotated meshes --- src/oce_mesh.F90 | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/src/oce_mesh.F90 b/src/oce_mesh.F90 index 97063e7bc..6e6a83b90 100755 --- a/src/oce_mesh.F90 +++ b/src/oce_mesh.F90 @@ -2435,6 +2435,7 @@ SUBROUTINE mesh_auxiliary_arrays(partit, mesh) USE MOD_PARSUP USE o_PARAM USE o_ARRAYS +USE g_CONFIG, only: rotated_grid, force_rotation USE g_ROTATE_grid use g_comm_auto use elem_center_interface @@ -2477,7 +2478,12 @@ SUBROUTINE mesh_auxiliary_arrays(partit, mesh) END DO DO n=1,myDim_nod2D+eDim_nod2D - call r2g(lon, lat, mesh%coord_nod2D(1,n), mesh%coord_nod2D(2,n)) + if ((.not. rotated_grid) .and. (.not. force_rotation)) then + lon = mesh%coord_nod2D(1,n) + lat = mesh%coord_nod2D(2,n) + else + call r2g(lon, lat, mesh%coord_nod2D(1,n), mesh%coord_nod2D(2,n)) + end if ! in case of numerical noise at the boundaries if (lon > 2._WP*pi) lon=lon-2._WP*pi if (lon <-2._WP*pi) lon=lon+2._WP*pi