Skip to content

Commit

Permalink
Patch to hide "bad lon" messages when using maskmap
Browse files Browse the repository at this point in the history
- Using a mask_table to mask out processors necessarily means that
  some geographic coordinate data is not known in all halo-regions.
- The iceberg code scans the geographic coordinate data and reports
  points that have "suspect" values. Using a mask_table caused a
  lot of messages to be written.
- I've disable this sanity check when a mask_table is being used.
- No answer changes.
  • Loading branch information
adcroft committed Sep 25, 2014
1 parent fcb1934 commit af926d5
Showing 1 changed file with 4 additions and 2 deletions.
6 changes: 4 additions & 2 deletions ice_bergs.F90
Original file line number Diff line number Diff line change
Expand Up @@ -2764,10 +2764,12 @@ subroutine icebergs_init(bergs, &
if (grd%lat(i,j).gt.900.) grd%lat(i,j)=2.*grd%lat(i,j+1)-grd%lat(i,j+2)
enddo; enddo

do j=grd%jsd,grd%jed; do i=grd%isd,grd%ied
if (.not. present(maskmap)) then ! Using a maskmap causes tickles this sanity check
do j=grd%jsd,grd%jed; do i=grd%isd,grd%ied
if (grd%lon(i,j).gt.900.) write(stderrunit,*) 'bad lon: ',mpp_pe(),i-grd%isc+1,j-grd%jsc+1,grd%lon(i,j)
if (grd%lat(i,j).gt.900.) write(stderrunit,*) 'bad lat: ',mpp_pe(),i-grd%isc+1,j-grd%jsc+1,grd%lat(i,j)
enddo; enddo
enddo; enddo
endif

! Sanitize lon for the tile (need continuous longitudes within one tile)
j=grd%jsc; do i=grd%isc+1,grd%ied
Expand Down

0 comments on commit af926d5

Please sign in to comment.