Skip to content

put in check for mask differences #23

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

Merged
dabail10 merged 1 commit into from
Oct 7, 2022
Merged

put in check for mask differences #23

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
56 changes: 49 additions & 7 deletions cicecore/drivers/nuopc/cmeps/ice_mesh_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -202,7 +202,7 @@ subroutine ice_mesh_set_distgrid(localpet, npes, distgrid, rc)
deallocate(gindex)

end subroutine ice_mesh_set_distgrid

!=======================================================================
subroutine ice_mesh_setmask_from_maskfile(ice_maskfile, ice_mesh, rc)

Expand Down Expand Up @@ -429,7 +429,7 @@ subroutine ice_mesh_create_scolumn(scol_lon, scol_lat, ice_mesh, rc)
! Allocate module variable ocn_gridcell_frac
allocate(ocn_gridcell_frac(2,2,1))
ocn_gridcell_frac(:,:,:) = scol_frac

end subroutine ice_mesh_create_scolumn

!===============================================================================
Expand Down Expand Up @@ -560,7 +560,7 @@ subroutine ice_mesh_check(gcomp, ice_mesh, rc)
! Check CICE mesh

use ice_constants, only : c1,c0,c360
use ice_grid , only : tlon, tlat
use ice_grid , only : tlon, tlat, hm

! input/output parameters
type(ESMF_GridComp) , intent(inout) :: gcomp
Expand All @@ -569,7 +569,8 @@ subroutine ice_mesh_check(gcomp, ice_mesh, rc)

! local variables
type(ESMF_DistGrid) :: distGrid
integer :: n,c,g,i,j,m ! indices
type(ESMF_Array) :: elemMaskArray
integer :: n,i,j ! indices
integer :: iblk, jblk ! indices
integer :: ilo, ihi, jlo, jhi ! beginning and end of physical domain
type(block) :: this_block ! block information for current block
Expand All @@ -578,11 +579,15 @@ subroutine ice_mesh_check(gcomp, ice_mesh, rc)
real(dbl_kind), pointer :: ownedElemCoords(:)
real(dbl_kind), pointer :: lat(:), latMesh(:)
real(dbl_kind), pointer :: lon(:), lonMesh(:)
integer , pointer :: model_mask(:)
real(dbl_kind) :: diff_lon
real(dbl_kind) :: diff_lat
real(dbl_kind) :: rad_to_deg
real(dbl_kind) :: tmplon, eps_imesh
logical :: isPresent, isSet
logical :: mask_error
integer :: mask_internal
integer :: mask_file
character(len=char_len_long) :: cvalue
character(len=char_len_long) :: logmsg
character(len=*), parameter :: subname = ' ice_mesh_check: '
Expand All @@ -606,7 +611,7 @@ subroutine ice_mesh_check(gcomp, ice_mesh, rc)
allocate(ownedElemCoords(spatialDim*numownedelements))
allocate(lonmesh(numOwnedElements))
allocate(latmesh(numOwnedElements))
call ESMF_MeshGet(ice_mesh, ownedElemCoords=ownedElemCoords)
call ESMF_MeshGet(ice_mesh, ownedElemCoords=ownedElemCoords, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return
do n = 1,numOwnedElements
lonMesh(n) = ownedElemCoords(2*n-1)
Expand Down Expand Up @@ -650,8 +655,45 @@ subroutine ice_mesh_check(gcomp, ice_mesh, rc)
enddo
enddo

100 format('ERROR: CICE n, lonmesh, lon, diff_lon = ',i6,2(f21.13,3x),d21.5)
101 format('ERROR: CICE n, latmesh, lat, diff_lat = ',i6,2(f21.13,3x),d21.5)
! obtain internally generated ice mask for error checks
allocate(model_mask(numOwnedElements))
call ESMF_MeshGet(ice_mesh, elementdistGrid=distGrid, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return
elemMaskArray = ESMF_ArrayCreate(distGrid, model_mask, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return
call ESMF_MeshGet(ice_mesh, elemMaskArray=elemMaskArray, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

mask_error = .false.
n=0
do iblk = 1, nblocks
this_block = get_block(blocks_ice(iblk),iblk)
do j = jlo, jhi
jlo = this_block%jlo
jhi = this_block%jhi
do i = ilo, ihi
ilo = this_block%ilo
ihi = this_block%ihi
n = n+1
mask_internal = nint(hm(i,j,iblk),kind=dbl_kind)
mask_file = model_mask(n)
if (mask_internal /= mask_file) then
write(6,102) n,mask_internal,mask_file
mask_error = .true.
end if
enddo !i
enddo !j
enddo !iblk
if (mask_error) then
call abort_ice(error_message=subname, file=__FILE__, line=__LINE__)
end if

call ESMF_ArrayDestroy(elemMaskArray, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

100 format('ERROR: CICE n, mesh_lon , lon, diff_lon = ',i8,2(f21.13,3x),d21.5)
101 format('ERROR: CICE n, mesh_lat , lat, diff_lat = ',i8,2(f21.13,3x),d21.5)
102 format('ERROR: CICE n, mesh_internal, mask_file = ',i8,2(i2,2x))

! deallocate memory
deallocate(ownedElemCoords)
Expand Down