[Nek5000-users] question on subroutine sem2lex

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Mon Jun 2 00:46:16 CDT 2014


Hi,all:
    In the parallel version of sem2lex, there is a part of code wrote as following:
if (gllnid(eg).eq.nid) then  ! element eg is on this processor
           eg = eg+1    
            le = gllel(ge)            ! le is local element number
            do iz=1,nz1
            do iy=1,ny1
            do ix=1,nx1

However, ge doesn't appear anywhere except here,(le=gllel(ge)), so is this should be "le=gllel(eg)"? Or just my misunderstanding?

Xianbei
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140602/7c9f083d/attachment.html>


More information about the Nek5000-users mailing list