From a571e82e02dd7ec6b7d27ab2f947f2ab05cdc13b Mon Sep 17 00:00:00 2001 From: Jean-Francois Lemieux Date: Tue, 12 Sep 2023 12:46:54 +0000 Subject: [PATCH] Modification of bugcheck condition for l_fixed_area=T --- cicecore/cicedyn/dynamics/ice_transport_remap.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/cicecore/cicedyn/dynamics/ice_transport_remap.F90 b/cicecore/cicedyn/dynamics/ice_transport_remap.F90 index 6a1294a6e..f91f96a33 100644 --- a/cicecore/cicedyn/dynamics/ice_transport_remap.F90 +++ b/cicecore/cicedyn/dynamics/ice_transport_remap.F90 @@ -63,7 +63,7 @@ module ice_transport_remap ! if false, area flux is determined internally ! and is passed out - logical (kind=log_kind), parameter :: bugcheck = .true. + logical (kind=log_kind), parameter :: bugcheck = .false. !======================================================================= ! Here is some information about how the incremental remapping scheme @@ -3008,7 +3008,7 @@ subroutine locate_triangles (nx_block, ny_block, & do ij = 1, icellsd i = indxid(ij) j = indxjd(ij) - if (abs(areasum(i,j) - edgearea(i,j)) > eps13*areafac_c(i,j)) then + if ( abs(areasum(i,j) - edgearea(i,j)) > eps13*areafac_c(i,j) .and. abs(edgearea(i,j)) > c0 ) then write(nu_diag,*) '' write(nu_diag,*) 'Areas do not add up: m, i, j, edge =', & my_task, i, j, trim(edge)