Skip to content

Commit

Permalink
GREPOS: Add check for filtered sideset with variables
Browse files Browse the repository at this point in the history
  • Loading branch information
gdsjaar committed Dec 20, 2023
1 parent 8399fb6 commit d93913d
Show file tree
Hide file tree
Showing 2 changed files with 9 additions and 6 deletions.
4 changes: 2 additions & 2 deletions packages/seacas/applications/grepos/gp_qainfo.blk
Original file line number Diff line number Diff line change
Expand Up @@ -6,5 +6,5 @@ C See packages/seacas/LICENSE for details

C -*- Mode: fortran -*-
QAINFO(1) = 'Grepos '
QAINFO(2) = '2022/10/06 '
QAINFO(3) = ' 1.92 '
QAINFO(2) = '2023/12/20 '
QAINFO(3) = ' 1.93 '
11 changes: 7 additions & 4 deletions packages/seacas/applications/grepos/grepos.f
Original file line number Diff line number Diff line change
Expand Up @@ -863,7 +863,7 @@ PROGRAM GREPOS
do i=0,numnps0-1
if (ia(kinpss+i) .eq.0) then
if (IA(KNNNS0+i) .ne. IA(KNNNS+i1)) then
write (*,900) 'Nodeset', ia(kidns0+i)
write (*,900) 'Nodeset', ia(kidns0+i), 'nodeset'
end if
i1 = i1 + 1
end if
Expand Down Expand Up @@ -946,15 +946,18 @@ PROGRAM GREPOS

call muntt(numess0, numess, nvarss,
$ ia(kssvok0), ia(kssvok), ia(kiesss))
endif

C ... check that the sidesets that are retained contain the same number
C of faces that the original sidesets contain. At the current time,
C can only map sideset variables if the sidesets are the same...
if (delel .and. nvarss .gt. 0) then
i1 = 0
do i=0,numess0-1
if (ia(kiesss+i) .eq.0) then
if (IA(KNESS0+i) .ne. IA(KNESS+i1)) then
write (*,900) 'Sideset', ia(kidss0+i)
write (*,900) 'Sideset', ia(kidss0+i), 'sideset'
stop 'Cannot Map Sideset Variables'
end if
i1 = i1 + 1
end if
Expand Down Expand Up @@ -1283,8 +1286,8 @@ PROGRAM GREPOS
& 14X,'GG GG RR RR EE PP OO OO SS'/
& 14X,' GGGGG RR RR EEEEEEE PP OOOOO SSSSSS ')
900 FORMAT(/,'WARNING: ',A,i5,' is a different size in the output',
$ /,9x,'database than in the input database. If there are',
$ /,9x,'variables on this sideset, they will be transferred',
$ /,9x,'database than in the input database. The ',
$ /,9x,'variables on this ',A,' will be transferred',
$ /,9x,'incorrectly. Contact [email protected]',
$ /,9x,'if you need this capability.')
END
Expand Down

0 comments on commit d93913d

Please sign in to comment.