[Nek5000-users] question on subroutine sem2lex

[nek5000-users at lists.mcs.anl.gov]

Hi Xianbei,

I would suggest checking out a stand-alone interpolation tool

nek5_svn/trunk/tools/int_tp/README

Best.
Aleks

________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Monday, June 02, 2014 8:05 AM
To: Nek5000
Subject: Re: [Nek5000-users] question on subroutine sem2lex

Hi,Paul:
    You are right. The code seems just out put the u in another ix,iy,iz,le way in spectral data, not what we expected as ix,iy,iz, where ix(max)=(nelx-1)*lx1+1. So is there any better way to output the data ready for the energy spectrum and two-point correlation already exists?

Xianbei







At 2014-06-02 08:22:55, nek5000-users at lists.mcs.anl.gov wrote:

Hi Xianbei,

It should be eg.

However, note that it's really not as simple as this....You'll need to understand
where your data resides, etc.

I think this code example is quite dated and definitely not tested as it is not
part of the standard source.

Paul

________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov<mailto:nek5000-users-bounces at lists.mcs.anl.gov> [nek5000-users-bounces at lists.mcs.anl.gov<mailto:nek5000-users-bounces at lists.mcs.anl.gov>] on behalf of nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> [nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>]
Sent: Monday, June 02, 2014 12:46 AM
To: Nek5000
Subject: [Nek5000-users] question on subroutine sem2lex

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


来自网易手机号码邮箱了解更多<http://shouji.163.com>


来自网易手机号码邮箱了解更多<http://shouji.163.com>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140602/e6d803dd/attachment-0001.html>