From nek5000-users at lists.mcs.anl.gov Mon Aug 1 06:38:10 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 1 Aug 2016 11:38:10 +0000 Subject: [Nek5000-users] maketools prenek Message-ID: Hello Paul, I had already tried switching BIGMEM to false but the same error appears. I am using Cygwin platform on windows to provide the linux environment and compiling on gcc 5.4.0 and gfortran compilers. Shivang -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Aug 1 07:04:36 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 1 Aug 2016 12:04:36 +0000 Subject: [Nek5000-users] maketools prenek In-Reply-To: References: Message-ID: Hi Shivang, Do you have x-windows installed (properly)? Can you run any other X-based apps (e.g., gnu plot) ? Best, Paul ________________________________ 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, August 01, 2016 6:38 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] maketools prenek Hello Paul, I had already tried switching BIGMEM to false but the same error appears. I am using Cygwin platform on windows to provide the linux environment and compiling on gcc 5.4.0 and gfortran compilers. Shivang -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Aug 2 10:06:13 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 2 Aug 2016 11:06:13 -0400 Subject: [Nek5000-users] Running NekExamples Message-ID: Hello, I am trying to run a standalone Nek5000 Example. I have Nek5000 installed and I am following the Examples instructions from https://github.com/Nek5000/NekExamples When I try to run any particular example, using: python -m 'unittest' NekTests.ConjHt I face two issues: 1) It looks like some Nek5000 files keep getting rebuilt each time I execute the above command, even if they have been built before. 2) I always get an error: Successfully compiled nek5000! Running nek5000... Using command "~/github/NekExamples/conj_ht/nek5000" Using working directory "~/github/NekExamples/conj_ht" Using .rea file "conj_ht" Could not remove ~/github/NekExamples/conj_ht/conj_ht.sch: [Errno 2] No such file or directory: '~/github/NekExamples/conj_ht/conj_ht.sch' Could not successfully run nek5000! Caught error: [Errno 2] No such file or directory I tried looking through the Makefiles and the Python script but cannot see where the sch file is specified. Also, I tried creating a dummy sch file so that the build can find and remove it, but that didn't work either. Please advise on how I can resolve the errors. Thanks, K.N.Ramachandran Ph: 814-441-4279 -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Aug 2 10:53:58 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 2 Aug 2016 11:53:58 -0400 Subject: [Nek5000-users] Spectral method with Nek Message-ID: All, I was wondering if the following statement is true or false: "If a single element is used for the whole domain then one can claim the Nek simulation is spectral instead of spectral element". Based on this, I carried out DNS on cavity for relatively large Reynolds numbers and I was not able to capture the instability mechanisms with one element (even when I increased the power of polynomial to 12 or 18). While I'm not surprised, I just wanted to see if anyone else has a similar experience. As a side question, why Nek utilizes Legendre polynomials instead of Chebyshev for each element? Is this because we're also interested in modeling non-uniformly spaced nodes? Best, Saleh From nek5000-users at lists.mcs.anl.gov Tue Aug 2 10:59:59 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 2 Aug 2016 10:59:59 -0500 Subject: [Nek5000-users] Running NekExamples In-Reply-To: References: Message-ID: Good morning! 1) It looks like some Nek5000 files keep getting rebuilt each time I > execute the above command, even if they have been built before. Yes, this is intentional. NekTests.py is used primarily for software testing. We want to ensure the tests are reproducible, so Nek5000 is built from scratch every single time. Regarding (2), could you set VERBOSE_TESTS="true" in your environment and let us know what errors you see? Sorry that VERBOSE_TESTS is currently undocumented; I will document it in the README as soon as possible. Best wishes, Ron Rahaman Argonne National Laboratory, Math and Computer Science Division University of Chicago, Department of Computer Science On Tue, Aug 2, 2016 at 10:06 AM, wrote: > Hello, > > I am trying to run a standalone Nek5000 Example. I have Nek5000 installed > and I am following the Examples instructions from > https://github.com/Nek5000/NekExamples > > When I try to run any particular example, using: > > python -m 'unittest' NekTests.ConjHt > > I face two issues: > > 1) It looks like some Nek5000 files keep getting rebuilt each time I > execute the above command, even if they have been built before. > > 2) I always get an error: > > Successfully compiled nek5000! > Running nek5000... > Using command "~/github/NekExamples/conj_ht/nek5000" > Using working directory "~/github/NekExamples/conj_ht" > Using .rea file "conj_ht" > Could not remove ~/github/NekExamples/conj_ht/conj_ht.sch: [Errno 2] > No such file or directory: '~/github/NekExamples/conj_ht/conj_ht.sch' > Could not successfully run nek5000! Caught error: [Errno 2] No such file > or directory > > I tried looking through the Makefiles and the Python script but cannot see > where the sch file is specified. Also, I tried creating a dummy sch file so > that the build can find and remove it, but that didn't work either. > > Please advise on how I can resolve the errors. > > > Thanks, > K.N.Ramachandran > Ph: 814-441-4279 > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Aug 2 11:11:30 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 2 Aug 2016 16:11:30 +0000 Subject: [Nek5000-users] Spectral method with Nek In-Reply-To: References: Message-ID: This is absolutely true - one element corresponds to a spectral calculation. You do need, however, enough points to resolve the scales of motion, as is always the case, which means you need large values of N. (The first spectral results were 32^3, and then the breakthrough 100^3 by Orszag in the early 70s - the famed centicubed simulations.) Paul ________________________________________ 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: Tuesday, August 02, 2016 10:53 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] Spectral method with Nek All, I was wondering if the following statement is true or false: "If a single element is used for the whole domain then one can claim the Nek simulation is spectral instead of spectral element". Based on this, I carried out DNS on cavity for relatively large Reynolds numbers and I was not able to capture the instability mechanisms with one element (even when I increased the power of polynomial to 12 or 18). While I'm not surprised, I just wanted to see if anyone else has a similar experience. As a side question, why Nek utilizes Legendre polynomials instead of Chebyshev for each element? Is this because we're also interested in modeling non-uniformly spaced nodes? Best, Saleh _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Tue Aug 2 13:58:34 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 2 Aug 2016 14:58:34 -0400 Subject: [Nek5000-users] Running NekExamples In-Reply-To: References: Message-ID: Hello, Thank you for the reply. I had tried specifying the VERBOSE_TESTS as mentioned earlier but I don't see any new command line output. I am able to see that the VERBOSE_TESTS setting is picked up, i.e. Getting setup options... Using F77="gfortran" Using CC="gcc" Using IFMPI="false" Using VERBOSE_TESTS="true" Using SOURCE_ROOT at ~/github/Nek5000 Using EXAMPLES_ROOT at ~/github/NekExamples Using TOOLS_ROOT at ~/github/Nek5000/tools Using TOOLS_BIN at ~/github/Nek5000/bin But I still get the same error message and output as before, i.e. that the .sch file is not found and so the running process just stops. Any thoughts? On Tue, Aug 2, 2016 at 11:59 AM, wrote: > Good morning! > > 1) It looks like some Nek5000 files keep getting rebuilt each time I >> execute the above command, even if they have been built before. > > > Yes, this is intentional. NekTests.py is used primarily for software > testing. We want to ensure the tests are reproducible, so Nek5000 is built > from scratch every single time. > > Regarding (2), could you set VERBOSE_TESTS="true" in your environment and > let us know what errors you see? > > Sorry that VERBOSE_TESTS is currently undocumented; I will document it in > the README as soon as possible. > > Best wishes, > Ron Rahaman > > Argonne National Laboratory, Math and Computer Science Division > University of Chicago, Department of Computer Science > > On Tue, Aug 2, 2016 at 10:06 AM, wrote: > >> Hello, >> >> I am trying to run a standalone Nek5000 Example. I have Nek5000 installed >> and I am following the Examples instructions from >> https://github.com/Nek5000/NekExamples >> >> When I try to run any particular example, using: >> >> python -m 'unittest' NekTests.ConjHt >> >> I face two issues: >> >> 1) It looks like some Nek5000 files keep getting rebuilt each time I >> execute the above command, even if they have been built before. >> >> 2) I always get an error: >> >> Successfully compiled nek5000! >> Running nek5000... >> Using command "~/github/NekExamples/conj_ht/nek5000" >> Using working directory "~/github/NekExamples/conj_ht" >> Using .rea file "conj_ht" >> Could not remove ~/github/NekExamples/conj_ht/conj_ht.sch: [Errno 2] >> No such file or directory: '~/github/NekExamples/conj_ht/conj_ht.sch' >> Could not successfully run nek5000! Caught error: [Errno 2] No such file >> or directory >> >> I tried looking through the Makefiles and the Python script but cannot >> see where the sch file is specified. Also, I tried creating a dummy sch >> file so that the build can find and remove it, but that didn't work either. >> >> Please advise on how I can resolve the errors. >> >> >> Thanks, >> K.N.Ramachandran >> Ph: 814-441-4279 >> >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users at lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >> >> > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -- K.N.Ramachandran Ph: 814-441-4279 -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Aug 2 11:53:51 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 2 Aug 2016 12:53:51 -0400 Subject: [Nek5000-users] Spectral method with Nek In-Reply-To: References: Message-ID: Thanks Paul for the reply. In fact, in one case I increased the order of polynomial significantly and I was able to resolve the small scales and, thereby, the route to chaos. But my experience was using one element with larger order of polynomial as compared to various elements but lower order of polynomials was slower (less computationally efficient, specially when I performed parallel using openmpi). Is there any explanation for this? Thanks again. On 8/2/2016 12:11 PM, nek5000-users at lists.mcs.anl.gov wrote: > This is absolutely true - one element corresponds to a spectral calculation. > > You do need, however, enough points to resolve the scales of motion, as is always > the case, which means you need large values of N. (The first spectral results were > 32^3, and then the breakthrough 100^3 by Orszag in the early 70s - the famed centicubed > simulations.) > > Paul > > ________________________________________ > 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: Tuesday, August 02, 2016 10:53 AM > To: nek5000-users at lists.mcs.anl.gov > Subject: [Nek5000-users] Spectral method with Nek > > All, > > I was wondering if the following statement is true or false: > > "If a single element is used for the whole domain then one can claim the > Nek simulation is spectral instead of spectral element". > > Based on this, I carried out DNS on cavity for relatively large Reynolds > numbers and I was not able to capture the instability mechanisms with > one element (even when I increased the power of polynomial to 12 or 18). > While I'm not surprised, I just wanted to see if anyone else has a > similar experience. > > As a side question, why Nek utilizes Legendre polynomials instead of > Chebyshev for each element? Is this because we're also interested in > modeling non-uniformly spaced nodes? > > Best, > Saleh > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Tue Aug 2 14:01:59 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 2 Aug 2016 14:01:59 -0500 Subject: [Nek5000-users] Running NekExamples In-Reply-To: References: Message-ID: Are there any errors prior to the message "Successfully compiled Nek5000"? Equivalently, are there any errors in compiler.out? Best, Ron On Tue, Aug 2, 2016 at 1:58 PM, wrote: > Hello, > > Thank you for the reply. I had tried specifying the VERBOSE_TESTS as > mentioned earlier but I don't see any new command line output. > > I am able to see that the VERBOSE_TESTS setting is picked up, i.e. > > Getting setup options... > Using F77="gfortran" > Using CC="gcc" > Using IFMPI="false" > Using VERBOSE_TESTS="true" > Using SOURCE_ROOT at ~/github/Nek5000 > Using EXAMPLES_ROOT at ~/github/NekExamples > Using TOOLS_ROOT at ~/github/Nek5000/tools > Using TOOLS_BIN at ~/github/Nek5000/bin > > But I still get the same error message and output as before, i.e. that the > .sch file is not found and so the running process just stops. Any thoughts? > > > > > On Tue, Aug 2, 2016 at 11:59 AM, wrote: > >> Good morning! >> >> 1) It looks like some Nek5000 files keep getting rebuilt each time I >>> execute the above command, even if they have been built before. >> >> >> Yes, this is intentional. NekTests.py is used primarily for software >> testing. We want to ensure the tests are reproducible, so Nek5000 is built >> from scratch every single time. >> >> Regarding (2), could you set VERBOSE_TESTS="true" in your environment and >> let us know what errors you see? >> >> Sorry that VERBOSE_TESTS is currently undocumented; I will document it in >> the README as soon as possible. >> >> Best wishes, >> Ron Rahaman >> >> Argonne National Laboratory, Math and Computer Science Division >> University of Chicago, Department of Computer Science >> >> On Tue, Aug 2, 2016 at 10:06 AM, wrote: >> >>> Hello, >>> >>> I am trying to run a standalone Nek5000 Example. I have Nek5000 >>> installed and I am following the Examples instructions from >>> https://github.com/Nek5000/NekExamples >>> >>> When I try to run any particular example, using: >>> >>> python -m 'unittest' NekTests.ConjHt >>> >>> I face two issues: >>> >>> 1) It looks like some Nek5000 files keep getting rebuilt each time I >>> execute the above command, even if they have been built before. >>> >>> 2) I always get an error: >>> >>> Successfully compiled nek5000! >>> Running nek5000... >>> Using command "~/github/NekExamples/conj_ht/nek5000" >>> Using working directory "~/github/NekExamples/conj_ht" >>> Using .rea file "conj_ht" >>> Could not remove ~/github/NekExamples/conj_ht/conj_ht.sch: [Errno 2] >>> No such file or directory: '~/github/NekExamples/conj_ht/conj_ht.sch' >>> Could not successfully run nek5000! Caught error: [Errno 2] No such file >>> or directory >>> >>> I tried looking through the Makefiles and the Python script but cannot >>> see where the sch file is specified. Also, I tried creating a dummy sch >>> file so that the build can find and remove it, but that didn't work either. >>> >>> Please advise on how I can resolve the errors. >>> >>> >>> Thanks, >>> K.N.Ramachandran >>> Ph: 814-441-4279 >>> >>> _______________________________________________ >>> Nek5000-users mailing list >>> Nek5000-users at lists.mcs.anl.gov >>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >>> >>> >> >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users at lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >> >> > > > -- > K.N.Ramachandran > Ph: 814-441-4279 > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Aug 2 14:25:50 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 2 Aug 2016 19:25:50 +0000 Subject: [Nek5000-users] Spectral method with Nek In-Reply-To: References: , Message-ID: The number of MPI ranks, P, is <= E, the number elements. In the E=1 case you can use only one mpi rank. Also, no guarantees that, as a spectral code, E=1 would be faster than an equivalent number of spectral element points given that the code is optimized for the multi element case. The bottom line, however, is that E=1 corresponds to a spectral method. Paul ________________________________________ 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: Tuesday, August 02, 2016 11:53 AM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] Spectral method with Nek Thanks Paul for the reply. In fact, in one case I increased the order of polynomial significantly and I was able to resolve the small scales and, thereby, the route to chaos. But my experience was using one element with larger order of polynomial as compared to various elements but lower order of polynomials was slower (less computationally efficient, specially when I performed parallel using openmpi). Is there any explanation for this? Thanks again. On 8/2/2016 12:11 PM, nek5000-users at lists.mcs.anl.gov wrote: > This is absolutely true - one element corresponds to a spectral calculation. > > You do need, however, enough points to resolve the scales of motion, as is always > the case, which means you need large values of N. (The first spectral results were > 32^3, and then the breakthrough 100^3 by Orszag in the early 70s - the famed centicubed > simulations.) > > Paul > > ________________________________________ > 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: Tuesday, August 02, 2016 10:53 AM > To: nek5000-users at lists.mcs.anl.gov > Subject: [Nek5000-users] Spectral method with Nek > > All, > > I was wondering if the following statement is true or false: > > "If a single element is used for the whole domain then one can claim the > Nek simulation is spectral instead of spectral element". > > Based on this, I carried out DNS on cavity for relatively large Reynolds > numbers and I was not able to capture the instability mechanisms with > one element (even when I increased the power of polynomial to 12 or 18). > While I'm not surprised, I just wanted to see if anyone else has a > similar experience. > > As a side question, why Nek utilizes Legendre polynomials instead of > Chebyshev for each element? Is this because we're also interested in > modeling non-uniformly spaced nodes? > > Best, > Saleh > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Tue Aug 2 15:19:21 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 2 Aug 2016 16:19:21 -0400 Subject: [Nek5000-users] Running NekExamples In-Reply-To: References: Message-ID: Hi Ron, Yes, that helped! I had a closer look and noted errors in the compilation and it was related to the MPI compilers. Once I specified the right compilers, I am able to run the examples. Thanks for resolving this. On Tue, Aug 2, 2016 at 3:01 PM, wrote: > Are there any errors prior to the message "Successfully compiled > Nek5000"? Equivalently, are there any errors in compiler.out? > > Best, > Ron > > On Tue, Aug 2, 2016 at 1:58 PM, wrote: > >> Hello, >> >> Thank you for the reply. I had tried specifying the VERBOSE_TESTS as >> mentioned earlier but I don't see any new command line output. >> >> I am able to see that the VERBOSE_TESTS setting is picked up, i.e. >> >> Getting setup options... >> Using F77="gfortran" >> Using CC="gcc" >> Using IFMPI="false" >> Using VERBOSE_TESTS="true" >> Using SOURCE_ROOT at ~/github/Nek5000 >> Using EXAMPLES_ROOT at ~/github/NekExamples >> Using TOOLS_ROOT at ~/github/Nek5000/tools >> Using TOOLS_BIN at ~/github/Nek5000/bin >> >> But I still get the same error message and output as before, i.e. that >> the .sch file is not found and so the running process just stops. Any >> thoughts? >> >> >> >> >> On Tue, Aug 2, 2016 at 11:59 AM, wrote: >> >>> Good morning! >>> >>> 1) It looks like some Nek5000 files keep getting rebuilt each time I >>>> execute the above command, even if they have been built before. >>> >>> >>> Yes, this is intentional. NekTests.py is used primarily for software >>> testing. We want to ensure the tests are reproducible, so Nek5000 is built >>> from scratch every single time. >>> >>> Regarding (2), could you set VERBOSE_TESTS="true" in your environment >>> and let us know what errors you see? >>> >>> Sorry that VERBOSE_TESTS is currently undocumented; I will document it >>> in the README as soon as possible. >>> >>> Best wishes, >>> Ron Rahaman >>> >>> Argonne National Laboratory, Math and Computer Science Division >>> University of Chicago, Department of Computer Science >>> >>> On Tue, Aug 2, 2016 at 10:06 AM, >>> wrote: >>> >>>> Hello, >>>> >>>> I am trying to run a standalone Nek5000 Example. I have Nek5000 >>>> installed and I am following the Examples instructions from >>>> https://github.com/Nek5000/NekExamples >>>> >>>> When I try to run any particular example, using: >>>> >>>> python -m 'unittest' NekTests.ConjHt >>>> >>>> I face two issues: >>>> >>>> 1) It looks like some Nek5000 files keep getting rebuilt each time I >>>> execute the above command, even if they have been built before. >>>> >>>> 2) I always get an error: >>>> >>>> Successfully compiled nek5000! >>>> Running nek5000... >>>> Using command "~/github/NekExamples/conj_ht/nek5000" >>>> Using working directory "~/github/NekExamples/conj_ht" >>>> Using .rea file "conj_ht" >>>> Could not remove ~/github/NekExamples/conj_ht/conj_ht.sch: [Errno >>>> 2] No such file or directory: '~/github/NekExamples/conj_ht/conj_ht.sch' >>>> Could not successfully run nek5000! Caught error: [Errno 2] No such >>>> file or directory >>>> >>>> I tried looking through the Makefiles and the Python script but cannot >>>> see where the sch file is specified. Also, I tried creating a dummy sch >>>> file so that the build can find and remove it, but that didn't work either. >>>> >>>> Please advise on how I can resolve the errors. >>>> >>>> >>>> Thanks, >>>> K.N.Ramachandran >>>> Ph: 814-441-4279 >>>> >>>> _______________________________________________ >>>> Nek5000-users mailing list >>>> Nek5000-users at lists.mcs.anl.gov >>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >>>> >>>> >>> >>> _______________________________________________ >>> Nek5000-users mailing list >>> Nek5000-users at lists.mcs.anl.gov >>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >>> >>> >> >> >> -- >> K.N.Ramachandran >> Ph: 814-441-4279 >> >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users at lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >> >> > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -- K.N.Ramachandran Ph: 814-441-4279 -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Aug 5 14:36:17 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 5 Aug 2016 19:36:17 +0000 Subject: [Nek5000-users] Nek5000-users Digest, Vol 89, Issue 10 In-Reply-To: References: Message-ID: Patel,Mansi A has shared a?OneDrive for Business?file with you. To view it, click the link below. [https://r1.res.office365.com/owa/prem/images/dc-png_20.png] Screenshot-5.png Hello users, Thanks Ron Rahaman that was much helpful. I am having redhat linux workstation and I have downloaded openmpi in it to run it parallel I typed ./nekbmpi eddy_uv 4 I getting different error it says p= 4 has passed but not compiled with DMPI. I have no clue what does that mean. I have attached screen shot of the error(screen shot 5) Thanks, Mansi Patel ________________________________ From: nek5000-users-bounces at lists.mcs.anl.gov on behalf of nek5000-users-request at lists.mcs.anl.gov Sent: Thursday, July 14, 2016 1:00:10 PM To: nek5000-users at lists.mcs.anl.gov Subject: Nek5000-users Digest, Vol 89, Issue 10 Send Nek5000-users mailing list submissions to nek5000-users at lists.mcs.anl.gov To subscribe or unsubscribe via the World Wide Web, visit https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users or, via email, send a message with subject or body 'help' to nek5000-users-request at lists.mcs.anl.gov You can reach the person managing the list at nek5000-users-owner at lists.mcs.anl.gov When replying, please edit your Subject line so it is more specific than "Re: Contents of Nek5000-users digest..." Today's Topics: 1. problems in compiling the code for parallel run (nek5000-users at lists.mcs.anl.gov) 2. Re: problems in compiling the code for parallel run (nek5000-users at lists.mcs.anl.gov) ---------------------------------------------------------------------- Message: 1 Date: Wed, 13 Jul 2016 20:27:38 +0000 From: nek5000-users at lists.mcs.anl.gov To: "nek5000-users at lists.mcs.anl.gov" , "tanvikaushik13 at gmail.com" Subject: [Nek5000-users] problems in compiling the code for parallel run Message-ID: Content-Type: text/plain; charset="iso-8859-1" Hello Users, I have Redhat linux workstation, I have already compiled code with nekb by making changes in makenek file, but now I want to run it parallel. At first when I compiled it by making F77 = mpif77 and cc= mpicc and commenting IFMPI ="false" I got errors that mpi not found than I installed openmpi Now when I compile using makenek it gets compiled successfully but when I try to run it using nekbmpi i get errors like mv: cannot stat `eddy_uv.sch': No such file or directory and sometimes if i do ./nekbmpi eddy_uv 4 it and put enter i get no errors on first time but if i enter again for second time I get same error that is mv: cannot stat `eddy_uv.sch': No such file or directory I want to understand what mistake I am doing and what should I do to get rid of errors and also which file can help me to see whats going on in simulation Thanks, Mansi Patel -------------- next part -------------- An HTML attachment was scrubbed... URL: ------------------------------ Message: 2 Date: Thu, 14 Jul 2016 07:20:10 -0500 From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] problems in compiling the code for parallel run Message-ID: Content-Type: text/plain; charset="utf-8" Hi Mansi, I suspect that you can safely ignore the error. The error probably comes from line 7 of the "nekbmpi" script: mv $1.sch $1.sch1 This command is intended to move the .sch file from a previous run, in case there's one in your directory. If this is your first run, there will be no .sch file in your directory, and the "mv" command will give the error you saw. So it probably doesn't indicate a problem. Best, Ron Rahaman Argonne National Laboratory, Math and Computer Science Division University of Chicago, Department of Computer Science On Wed, Jul 13, 2016 at 3:27 PM, wrote: > Hello Users, > > I have Redhat linux workstation, I have already compiled code with nekb by > making changes in makenek file, but now I want to run it parallel. At first > when I compiled it by making F77 = mpif77 and cc= mpicc and commenting > IFMPI ="false" I got errors that mpi not found than I installed openmpi > > Now when I compile using makenek it gets compiled successfully but when I > try to run it using nekbmpi i get errors like > mv: cannot stat `eddy_uv.sch': No such file or directory > and sometimes if i do ./nekbmpi eddy_uv 4 it and put enter i get no errors > on first time but if i enter again for second time I get same error that is mv: > cannot stat `eddy_uv.sch': No such file or directory > I want to understand what mistake I am doing and what should I do to get > rid of errors and also which file can help me to see whats going on in > simulation > > > Thanks, > Mansi Patel > > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: ------------------------------ _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users End of Nek5000-users Digest, Vol 89, Issue 10 ********************************************* -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Aug 5 14:36:21 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 5 Aug 2016 19:36:21 +0000 Subject: [Nek5000-users] Patel, Mansi A wants to share the file Screenshot-5.png with you Message-ID: To view Screenshot-5.png, sign in or create an account. -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Aug 8 14:39:44 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 8 Aug 2016 21:39:44 +0200 Subject: [Nek5000-users] Visualization of simulation In-Reply-To: References: Message-ID: Hello Users I am facing some difficulties in visualising the result of simulation. On running " postx" a blank screen appears without any side menu, as mentioned in manual. Is this is a result of any error during building postnek ? Because i didn't encountered any ! Thanks Shivang -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Sun Aug 14 12:10:28 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sun, 14 Aug 2016 10:10:28 -0700 Subject: [Nek5000-users] vtk into prenek Message-ID: Neks, I am trying to import some mesh information into prenek via a vtk file. Can someone please post an example of one that works into the git repo? Ideally it would be nice to have a 2d and 3d case for vtk files that work in prenek. I have tried using the vtkUnstructuredGridWriter class from the vtk library and prenek does not accept the file. At first I thought it was because the class required 3d points and I am importing a 2d mesh, but after I modified the file to remove the extra coordinate and make it match the style guide ( http://www.vtk.org/wp-content/uploads/2015/04/file-formats.pdf) prenek reads the mesh without crashing, but no elements appear. The output to the terminal from prenek is: IMPORT VTK MESH Current number of elements: 0 input name of vtk file 4 0 1 2 3 inside gencen 0 F 4 done gencen 0 The contents of the file I am trying to import are below, so if you see a bug please let me know. # vtk DataFile Version 4.0 vtk output ASCII DATASET UNSTRUCTURED_GRID POINTS 7 float 0 0 2.1 0 1.48492 1.48492 0 2.1 3 0 2.12132 2.12132 0 3 CELLS 3 15 4 0 1 2 3 4 1 4 5 2 4 3 2 5 6 CELL_TYPES 3 9 9 9 Thanks, -- Phil S. -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Aug 18 06:02:43 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 18 Aug 2016 12:02:43 +0100 Subject: [Nek5000-users] MPI_TAG_UB Message-ID: Hi Neks, I was wondering if someone could tell me the implications of commenting this section of code out in Nekbone: if (nval.lt.(10000+max(lp,lelg))) then if(nid.eq.0) then write(6,*) 'ABORT: MPI_TAG_UB too small!' endif call exitt endif We are testing Nekbone with an experimental MPI library. Regards, Ken -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 490 bytes Desc: OpenPGP digital signature URL: -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: not available URL: From nek5000-users at lists.mcs.anl.gov Thu Aug 18 13:33:48 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 18 Aug 2016 14:33:48 -0400 Subject: [Nek5000-users] MPI_TAG_UB In-Reply-To: References: Message-ID: To answer your question I need to know why you want to comment this check? Am 8/18/16, 7:02 AM schrieb "nek5000-users-bounces at lists.mcs.anl.gov on behalf of nek5000-users at lists.mcs.anl.gov" unter : >Hi Neks, > >I was wondering if someone could tell me the implications of commenting >this section of code out in Nekbone: > > if (nval.lt.(10000+max(lp,lelg))) then > if(nid.eq.0) then > write(6,*) 'ABORT: MPI_TAG_UB too small!' > endif > call exitt > endif > >We are testing Nekbone with an experimental MPI library. > >Regards, >Ken > >The University of Edinburgh is a charitable body, registered in >Scotland, with registration number SC005336. >_______________________________________________ >Nek5000-users mailing list >Nek5000-users at lists.mcs.anl.gov >https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Thu Aug 18 07:55:58 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 18 Aug 2016 13:55:58 +0100 Subject: [Nek5000-users] MPI_TAG_UB In-Reply-To: References: Message-ID: Basically, the MPI library that we are testing sets MPI_TAG_UB=65534, which is quite low to be used with higher numbers of lp. On 18/08/16 19:33, nek5000-users at lists.mcs.anl.gov wrote: > To answer your question I need to know why you want to comment this check? > > > Am 8/18/16, 7:02 AM schrieb "nek5000-users-bounces at lists.mcs.anl.gov on > behalf of nek5000-users at lists.mcs.anl.gov" unter > nek5000-users at lists.mcs.anl.gov>: > >> Hi Neks, >> >> I was wondering if someone could tell me the implications of commenting >> this section of code out in Nekbone: >> >> if (nval.lt.(10000+max(lp,lelg))) then >> if(nid.eq.0) then >> write(6,*) 'ABORT: MPI_TAG_UB too small!' >> endif >> call exitt >> endif >> >> We are testing Nekbone with an experimental MPI library. >> >> Regards, >> Ken >> >> The University of Edinburgh is a charitable body, registered in >> Scotland, with registration number SC005336. >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users at lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 490 bytes Desc: OpenPGP digital signature URL: -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: not available URL: From nek5000-users at lists.mcs.anl.gov Thu Aug 18 14:07:01 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 18 Aug 2016 15:07:01 -0400 Subject: [Nek5000-users] MPI_TAG_UB In-Reply-To: References: Message-ID: I checked the source (nekbone 3.1) and it looks like the fortran part does not use any tagged messages but the C part does. In fact, there are C parts using the MPI rank to tag messages. So I am afraid commenting the check won?t help. Cheers, Stefan Am 8/18/16, 8:55 AM schrieb "nek5000-users-bounces at lists.mcs.anl.gov on behalf of nek5000-users at lists.mcs.anl.gov" unter : >Basically, the MPI library that we are testing sets MPI_TAG_UB=65534, >which is quite low to be used with higher numbers of lp. > > > >On 18/08/16 19:33, nek5000-users at lists.mcs.anl.gov wrote: >> To answer your question I need to know why you want to comment this >>check? >> >> >> Am 8/18/16, 7:02 AM schrieb "nek5000-users-bounces at lists.mcs.anl.gov on >> behalf of nek5000-users at lists.mcs.anl.gov" unter >> > nek5000-users at lists.mcs.anl.gov>: >> >>> Hi Neks, >>> >>> I was wondering if someone could tell me the implications of commenting >>> this section of code out in Nekbone: >>> >>> if (nval.lt.(10000+max(lp,lelg))) then >>> if(nid.eq.0) then >>> write(6,*) 'ABORT: MPI_TAG_UB too small!' >>> endif >>> call exitt >>> endif >>> >>> We are testing Nekbone with an experimental MPI library. >>> >>> Regards, >>> Ken >>> >>> The University of Edinburgh is a charitable body, registered in >>> Scotland, with registration number SC005336. >>> _______________________________________________ >>> Nek5000-users mailing list >>> Nek5000-users at lists.mcs.anl.gov >>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >> >> >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users at lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >> > >The University of Edinburgh is a charitable body, registered in >Scotland, with registration number SC005336. >_______________________________________________ >Nek5000-users mailing list >Nek5000-users at lists.mcs.anl.gov >https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Mon Aug 22 08:46:15 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Aug 2016 13:46:15 +0000 Subject: [Nek5000-users] General functionality of the .usr file Message-ID: Hi! I am a new user of Nek5000. I have read through the User Documentation, and ran some simulations. Now it is time for me to start editing a file that a colleague wrote fo me. I therefore need to know how the .usr file is used. For example: In the User Documentation it is written that uservp() is called "once per processor, and once per discrete point therein". Does this mean that this is called once before the simulation is executed or is it called once per time step or in some other fashion? I would guess that uservp() is called once before the simulation is run, is that correct? Does thi apply to all routines in the .usr file, or are some of them called repeatedly during the simulation? userbc() for instance? I have been looking for a definition of the variables that are accessible in the separate routines. For instance, the input arguments to useric are (ix,iy,iz,ieg). I have gotten to understand that ix, iy and iz relate to a spatial coordinate in some way, and that ieg is some kind of index for the fluid element. Is there a precise definition of these names available? In addition to the input arguments it seems as you have access to more variables than the input arguments in each routine. In useric() you may for instance change ux, uy and uz. Is there a full list of all parameters that are available available? /Johan -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Aug 22 09:38:31 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Aug 2016 14:38:31 +0000 Subject: [Nek5000-users] General functionality of the .usr file In-Reply-To: References: Message-ID: Johan, Essentially every variable is accessible through the include files, NEKUSE and TOTAL. NEKUSE variables are mapped and derferenced by Nek, point by point and step by step, except for useric(), which is called prior to start of time stepping. One critically important consideration is to not invoke any communication (or any expensive operation) in uservp, userbc, useric, userf, userq. If you try to communicate (e.g., alpha = glmax(t,n) invokes global communication) there will be trouble because different processors can have different numbers of points and not all of them will call the user routines the same number of times. The run would hang in that case, not to mention be very expensive. Should you need to do something like, the right approach is to do the global operation in userchk, store the result in a common block, and pass this to the user routine in question. I think if you look through the examples you will come to understand the basic approach. Paul ________________________________ 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, August 22, 2016 8:46 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] General functionality of the .usr file Hi! I am a new user of Nek5000. I have read through the User Documentation, and ran some simulations. Now it is time for me to start editing a file that a colleague wrote fo me. I therefore need to know how the .usr file is used. For example: In the User Documentation it is written that uservp() is called "once per processor, and once per discrete point therein". Does this mean that this is called once before the simulation is executed or is it called once per time step or in some other fashion? I would guess that uservp() is called once before the simulation is run, is that correct? Does thi apply to all routines in the .usr file, or are some of them called repeatedly during the simulation? userbc() for instance? I have been looking for a definition of the variables that are accessible in the separate routines. For instance, the input arguments to useric are (ix,iy,iz,ieg). I have gotten to understand that ix, iy and iz relate to a spatial coordinate in some way, and that ieg is some kind of index for the fluid element. Is there a precise definition of these names available? In addition to the input arguments it seems as you have access to more variables than the input arguments in each routine. In useric() you may for instance change ux, uy and uz. Is there a full list of all parameters that are available available? /Johan -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Aug 22 10:53:32 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Aug 2016 15:53:32 +0000 Subject: [Nek5000-users] General functionality of the .usr file In-Reply-To: References: , Message-ID: Thank you! I found many of the variables I recognize in the NEKUSE-file. The TOTAL-file there are lines of the type 'include' DXYZ Is DXYZ a separate file (I could not find it in such case)? Or is it a module (or similar)? Where can I find the description of what is included in TOTAL? I will be aware of not invoking communication in the .usr routines. /Johan ________________________________ From: nek5000-users-bounces at lists.mcs.anl.gov on behalf of nek5000-users at lists.mcs.anl.gov Sent: Monday, August 22, 2016 4:38 PM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] General functionality of the .usr file Johan, Essentially every variable is accessible through the include files, NEKUSE and TOTAL. NEKUSE variables are mapped and derferenced by Nek, point by point and step by step, except for useric(), which is called prior to start of time stepping. One critically important consideration is to not invoke any communication (or any expensive operation) in uservp, userbc, useric, userf, userq. If you try to communicate (e.g., alpha = glmax(t,n) invokes global communication) there will be trouble because different processors can have different numbers of points and not all of them will call the user routines the same number of times. The run would hang in that case, not to mention be very expensive. Should you need to do something like, the right approach is to do the global operation in userchk, store the result in a common block, and pass this to the user routine in question. I think if you look through the examples you will come to understand the basic approach. Paul ________________________________ 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, August 22, 2016 8:46 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] General functionality of the .usr file Hi! I am a new user of Nek5000. I have read through the User Documentation, and ran some simulations. Now it is time for me to start editing a file that a colleague wrote fo me. I therefore need to know how the .usr file is used. For example: In the User Documentation it is written that uservp() is called "once per processor, and once per discrete point therein". Does this mean that this is called once before the simulation is executed or is it called once per time step or in some other fashion? I would guess that uservp() is called once before the simulation is run, is that correct? Does thi apply to all routines in the .usr file, or are some of them called repeatedly during the simulation? userbc() for instance? I have been looking for a definition of the variables that are accessible in the separate routines. For instance, the input arguments to useric are (ix,iy,iz,ieg). I have gotten to understand that ix, iy and iz relate to a spatial coordinate in some way, and that ieg is some kind of index for the fluid element. Is there a precise definition of these names available? In addition to the input arguments it seems as you have access to more variables than the input arguments in each routine. In useric() you may for instance change ux, uy and uz. Is there a full list of all parameters that are available available? /Johan -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Aug 22 11:45:31 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Aug 2016 16:45:31 +0000 Subject: [Nek5000-users] General functionality of the .usr file In-Reply-To: References: , , Message-ID: You need to look in each of the include files mentioned in TOTAL. SOLN GEOM are two of the most important. Best, Paul ________________________________ 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, August 22, 2016 10:53 AM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] General functionality of the .usr file Thank you! I found many of the variables I recognize in the NEKUSE-file. The TOTAL-file there are lines of the type 'include' DXYZ Is DXYZ a separate file (I could not find it in such case)? Or is it a module (or similar)? Where can I find the description of what is included in TOTAL? I will be aware of not invoking communication in the .usr routines. /Johan ________________________________ From: nek5000-users-bounces at lists.mcs.anl.gov on behalf of nek5000-users at lists.mcs.anl.gov Sent: Monday, August 22, 2016 4:38 PM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] General functionality of the .usr file Johan, Essentially every variable is accessible through the include files, NEKUSE and TOTAL. NEKUSE variables are mapped and derferenced by Nek, point by point and step by step, except for useric(), which is called prior to start of time stepping. One critically important consideration is to not invoke any communication (or any expensive operation) in uservp, userbc, useric, userf, userq. If you try to communicate (e.g., alpha = glmax(t,n) invokes global communication) there will be trouble because different processors can have different numbers of points and not all of them will call the user routines the same number of times. The run would hang in that case, not to mention be very expensive. Should you need to do something like, the right approach is to do the global operation in userchk, store the result in a common block, and pass this to the user routine in question. I think if you look through the examples you will come to understand the basic approach. Paul ________________________________ 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, August 22, 2016 8:46 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] General functionality of the .usr file Hi! I am a new user of Nek5000. I have read through the User Documentation, and ran some simulations. Now it is time for me to start editing a file that a colleague wrote fo me. I therefore need to know how the .usr file is used. For example: In the User Documentation it is written that uservp() is called "once per processor, and once per discrete point therein". Does this mean that this is called once before the simulation is executed or is it called once per time step or in some other fashion? I would guess that uservp() is called once before the simulation is run, is that correct? Does thi apply to all routines in the .usr file, or are some of them called repeatedly during the simulation? userbc() for instance? I have been looking for a definition of the variables that are accessible in the separate routines. For instance, the input arguments to useric are (ix,iy,iz,ieg). I have gotten to understand that ix, iy and iz relate to a spatial coordinate in some way, and that ieg is some kind of index for the fluid element. Is there a precise definition of these names available? In addition to the input arguments it seems as you have access to more variables than the input arguments in each routine. In useric() you may for instance change ux, uy and uz. Is there a full list of all parameters that are available available? /Johan -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Aug 22 12:28:41 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Aug 2016 12:28:41 -0500 Subject: [Nek5000-users] Multiple Restart Files Message-ID: Hi, Sometimes I run into issues when Nek can't load the restart file that is specified in the .rea file. I'm not exactly sure why this is the case, but it's possible that the outputs are corrupted. Is it possible to specify multiple restart files in .rea i.e.: 2 PRESOLVE/RESTART OPTIONS ***** restart0.f00001 restart0.f00002 or do something equivalent so I can specify a backup restart file in case the first one is not readable? It's a non-issue if I'm running it on a system I own, but on queued systems, it's painful to resubmit the job and wait because of this issue. Please let me know if this is possible or there are other ways to read multiple restart files. Sincerely, Kento -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Aug 22 13:19:47 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Aug 2016 18:19:47 +0000 Subject: [Nek5000-users] Heat/passive scalar solver - Helmholtz Message-ID: Dear Neks, I am a little confused by the setup in the rea file to run a simulation with a constant cfl. My understanding (from the Nek5000 archives) is that by setting p012 to a positive value in the rea file, it should be sufficient to ensure a constant cfl condition during the run. When I do so, my time step is constant and thus my cfl varies. I tried to set p012 to a positive value in different rea file from the examples directory but none of them yield a constant cfl during the run. If someone has an answer to this problem I would really appreciate the help. Thanks, Marco From nek5000-users at lists.mcs.anl.gov Mon Aug 22 14:21:23 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Aug 2016 19:21:23 +0000 Subject: [Nek5000-users] Heat/passive scalar solver - Helmholtz In-Reply-To: <429D0B1D-0FF6-4F57-A387-5A92392D4884@ornl.gov> References: <429D0B1D-0FF6-4F57-A387-5A92392D4884@ornl.gov> Message-ID: Never mind. I got it. Marco > On Aug 22, 2016, at 2:19 PM, Delchini, Marc-Olivier G. wrote: > > Dear Neks, > > I am a little confused by the setup in the rea file to run a simulation with a constant cfl. My understanding (from the Nek5000 archives) is that by setting p012 to a positive value in the rea file, it should be sufficient to ensure a constant cfl condition during the run. When I do so, my time step is constant and thus my cfl varies. > > I tried to set p012 to a positive value in different rea file from the examples directory but none of them yield a constant cfl during the run. > > If someone has an answer to this problem I would really appreciate the help. > > Thanks, > Marco From nek5000-users at lists.mcs.anl.gov Mon Aug 22 14:28:06 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Aug 2016 19:28:06 +0000 Subject: [Nek5000-users] Heat/passive scalar solver - Helmholtz In-Reply-To: References: Message-ID: Hi Marco, If param(12) is positive, then it means that dt will not exceed param(12). and that dt will be chose such that CFL ~ CFL_target (specified in the .rea file). Thus, if you have: .001 p12 DT : : 0.5000000 p26 COURANT/NTAU then dt will be the min of (.001, dt_cfl ), where dt_cfl is such that CFL=0.5 (p26). Note that there is some tolerance in resetting dt ... i.e., I think you can get to CFL=0.55 before it would reset dt to target CFL of 0.5 Similarly, it would not increase dt until CFL < 0.45. ONE ADDITIONAL POINT: If you are using characteristics, then COURANT (p26) takes on a different meaning. In that case, it means how many substeps the characteristics scheme takes for each Navier-Stokes timestep. (1, 2 or 3 are most common values, but it can vary from case to case...) See: www.mcs.anl.gov/~fischer/new.pdf , last section on timestepping, for a description of characteristics as implemented in Nek. Paul ________________________________________ 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, August 22, 2016 1:19 PM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] Heat/passive scalar solver - Helmholtz Dear Neks, I am a little confused by the setup in the rea file to run a simulation with a constant cfl. My understanding (from the Nek5000 archives) is that by setting p012 to a positive value in the rea file, it should be sufficient to ensure a constant cfl condition during the run. When I do so, my time step is constant and thus my cfl varies. I tried to set p012 to a positive value in different rea file from the examples directory but none of them yield a constant cfl during the run. If someone has an answer to this problem I would really appreciate the help. Thanks, Marco _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Mon Aug 22 14:29:36 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Aug 2016 19:29:36 +0000 Subject: [Nek5000-users] Multiple Restart Files In-Reply-To: References: Message-ID: Hi Kento, I would suggest looking at the files in VisIt if you're worried that they may be corrupted? Paul ________________________________ 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, August 22, 2016 12:28 PM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] Multiple Restart Files Hi, Sometimes I run into issues when Nek can't load the restart file that is specified in the .rea file. I'm not exactly sure why this is the case, but it's possible that the outputs are corrupted. Is it possible to specify multiple restart files in .rea i.e.: 2 PRESOLVE/RESTART OPTIONS ***** restart0.f00001 restart0.f00002 or do something equivalent so I can specify a backup restart file in case the first one is not readable? It's a non-issue if I'm running it on a system I own, but on queued systems, it's painful to resubmit the job and wait because of this issue. Please let me know if this is possible or there are other ways to read multiple restart files. Sincerely, Kento -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Aug 22 14:42:32 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Aug 2016 19:42:32 +0000 Subject: [Nek5000-users] Heat/passive scalar solver - Helmholtz In-Reply-To: References: Message-ID: Hi Paul, thank you for your answer. What I misunderstood was the part 'dt will not exceed param(12)?. I was setting param(12) to a small value. Thanks, Marco > On Aug 22, 2016, at 3:28 PM, nek5000-users at lists.mcs.anl.gov wrote: > > > Hi Marco, > > If param(12) is positive, then it means that > > dt will not exceed param(12). > > and that > > dt will be chose such that CFL ~ CFL_target (specified in the .rea file). > > Thus, if you have: > > > .001 p12 DT > : : > 0.5000000 p26 COURANT/NTAU > > then dt will be the min of (.001, dt_cfl ), where dt_cfl is such that CFL=0.5 (p26). > > Note that there is some tolerance in resetting dt ... i.e., I think you can get to CFL=0.55 > before it would reset dt to target CFL of 0.5 > > Similarly, it would not increase dt until CFL < 0.45. > > ONE ADDITIONAL POINT: > > If you are using characteristics, then COURANT (p26) takes on a different meaning. > > In that case, it means how many substeps the characteristics scheme takes for each > Navier-Stokes timestep. (1, 2 or 3 are most common values, but it can vary from > case to case...) > > See: www.mcs.anl.gov/~fischer/new.pdf , last section on timestepping, for a > description of characteristics as implemented in Nek. > > Paul > > ________________________________________ > 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, August 22, 2016 1:19 PM > To: nek5000-users at lists.mcs.anl.gov > Subject: Re: [Nek5000-users] Heat/passive scalar solver - Helmholtz > > Dear Neks, > > I am a little confused by the setup in the rea file to run a simulation with a constant cfl. My understanding (from the Nek5000 archives) is that by setting p012 to a positive value in the rea file, it should be sufficient to ensure a constant cfl condition during the run. When I do so, my time step is constant and thus my cfl varies. > > I tried to set p012 to a positive value in different rea file from the examples directory but none of them yield a constant cfl during the run. > > If someone has an answer to this problem I would really appreciate the help. > > Thanks, > Marco > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > From nek5000-users at lists.mcs.anl.gov Tue Aug 23 09:11:45 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 23 Aug 2016 09:11:45 -0500 Subject: [Nek5000-users] Multiple Restart Files In-Reply-To: References: Message-ID: Ok, I will do that next time this issue shows up. The issues might also be caused by Cetus so I'll investigate further. Thank you for the feedback. Sincerely, Kento On Mon, Aug 22, 2016 at 2:29 PM, wrote: > > Hi Kento, > > I would suggest looking at the files in VisIt if you're worried that they > may > be corrupted? > > Paul > > ------------------------------ > *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, August 22, 2016 12:28 PM > *To:* nek5000-users at lists.mcs.anl.gov > *Subject:* [Nek5000-users] Multiple Restart Files > > Hi, > > Sometimes I run into issues when Nek can't load the restart file that is > specified in the .rea file. I'm not exactly sure why this is the case, but > it's possible that the outputs are corrupted. > > Is it possible to specify multiple restart files in .rea i.e.: > > 2 PRESOLVE/RESTART OPTIONS ***** > restart0.f00001 > restart0.f00002 > > or do something equivalent so I can specify a backup restart file in case > the first one is not readable? It's a non-issue if I'm running it on a > system I own, but on queued systems, it's painful to resubmit the job and > wait because of this issue. > > Please let me know if this is possible or there are other ways to read > multiple restart files. > > Sincerely, > Kento > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Aug 26 11:18:33 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 26 Aug 2016 12:18:33 -0400 Subject: [Nek5000-users] Error reading initial condition/drive force data aborting in routine rdicdf. EXIT: rdicdf error Message-ID: Hello all, I have been using the same .rea and .usr files in a different RESTART calculation where the input file is ASCII. But when the file is BINARY, it is showing this error. I have changed the Parameters 66 and 67 but to no avail. Any help is appreciated. Thanks, Saikat. *Saikat Mukherjee,* *Graduate** Student,* *Engineering Science and Mechanics,* *BEAM,* *Virginia Tech.* -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Aug 26 12:00:26 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 26 Aug 2016 19:00:26 +0200 Subject: [Nek5000-users] Error reading initial condition/drive force data aborting in routine rdicdf. EXIT: rdicdf error In-Reply-To: References: Message-ID: Please post your .rea file. Von: on behalf of Antworten an: Datum: Friday, August 26, 2016 at 6:18 PM An: Betreff: [Nek5000-users] Error reading initial condition/drive force data aborting in routine rdicdf. EXIT: rdicdf error Hello all, I have been using the same .rea and .usr files in a different RESTART calculation where the input file is ASCII. But when the file is BINARY, it is showing this error. I have changed the Parameters 66 and 67 but to no avail. Any help is appreciated. Thanks, Saikat. Saikat Mukherjee, Graduate Student, Engineering Science and Mechanics, BEAM, Virginia Tech. _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Aug 26 12:22:39 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 26 Aug 2016 13:22:39 -0400 Subject: [Nek5000-users] Error reading initial condition/drive force data aborting in routine rdicdf. EXIT: rdicdf error In-Reply-To: References: Message-ID: Here is the .rea file. Thanks, Saikat. *Saikat Mukherjee,* *Graduate** Student,* *Engineering Science and Mechanics,* *BEAM,* *Virginia Tech.* On Fri, Aug 26, 2016 at 1:00 PM, wrote: > Please post your .rea file. > > Von: on behalf of < > nek5000-users at lists.mcs.anl.gov> > Antworten an: > Datum: Friday, August 26, 2016 at 6:18 PM > An: > Betreff: [Nek5000-users] Error reading initial condition/drive force data > aborting in routine rdicdf. EXIT: rdicdf error > > Hello all, > > I have been using the same .rea and .usr files in a different RESTART > calculation where the input file is ASCII. But when the file is BINARY, it > is showing this error. I have changed the Parameters 66 and 67 but to no > avail. Any help is appreciated. > > Thanks, > > Saikat. > *Saikat Mukherjee,* > *Graduate** Student,* > *Engineering Science and Mechanics,* > *BEAM,* > *Virginia Tech.* > _______________________________________________ Nek5000-users mailing > list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/ > mailman/listinfo/nek5000-users > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: g10cond.rea Type: application/octet-stream Size: 558809 bytes Desc: not available URL: From nek5000-users at lists.mcs.anl.gov Fri Aug 26 12:30:43 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 26 Aug 2016 19:30:43 +0200 Subject: [Nek5000-users] Error reading initial condition/drive force data aborting in routine rdicdf. EXIT: rdicdf error In-Reply-To: References: Message-ID: I think the error is due to the blank lines (and hence incorrect line count) in the .rea file Von: on behalf of Antworten an: Datum: Friday, August 26, 2016 at 7:22 PM An: Betreff: Re: [Nek5000-users] Error reading initial condition/drive force data aborting in routine rdicdf. EXIT: rdicdf error Here is the .rea file. Thanks, Saikat. Saikat Mukherjee, Graduate Student, Engineering Science and Mechanics, BEAM, Virginia Tech. On Fri, Aug 26, 2016 at 1:00 PM, wrote: > Please post your .rea file. > > Von: on behalf of > > Antworten an: > Datum: Friday, August 26, 2016 at 6:18 PM > An: > Betreff: [Nek5000-users] Error reading initial condition/drive force data > aborting in routine rdicdf. EXIT: rdicdf error > > Hello all, > > I have been using the same .rea and .usr files in a different RESTART > calculation where the input file is ASCII. But when the file is BINARY, it is > showing this error. I have changed the Parameters 66 and 67 but to no avail. > Any help is appreciated. > > Thanks, > > Saikat. > Saikat Mukherjee, > Graduate Student, > Engineering Science and Mechanics, > BEAM, > Virginia Tech. > _______________________________________________ Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.govhttps://lists.mcs.anl.gov/mailman/listinfo/nek5 > 000-users > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Aug 26 14:20:46 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 26 Aug 2016 15:20:46 -0400 Subject: [Nek5000-users] set passive scalar initial conditions after some time t Message-ID: Hello all, I would like to apply some initial conditions for a passive scalar after some time t in the simulation as opposed to at the start of the simulation. In other words, start a simulation, let the flow field evolve, then after some time initialize the passive scalar with some prescribed initial condition, and then continue simulating the flow field and the passive scalar. Is there a simple way to accomplish this in the .usr file? What subroutine would be the best for this situation? If userchk ? how do we set the values of the passive scalars in this routine? Side note: we have found that for this case we can not first run the flow field calculation and then restart the calculation with the passive scalar values set as an initial condition. The restart affects the flow field too strongly since we are interested in a very specific state of the flow field. Thanks, Mark Mark Paul Professor of Mechanical Engineering Virginia Tech www.me.vt.edu/mpaul -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Aug 26 14:08:55 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 26 Aug 2016 15:08:55 -0400 Subject: [Nek5000-users] Error reading initial condition/drive force data aborting in routine rdicdf. EXIT: rdicdf error In-Reply-To: References: Message-ID: Thanks it seems to work now. Saikat. *Saikat Mukherjee,* *Graduate** Student,* *Engineering Science and Mechanics,* *BEAM,* *Virginia Tech.* On Fri, Aug 26, 2016 at 1:30 PM, wrote: > I think the error is due to the blank lines (and hence incorrect line > count) in the .rea file > > Von: on behalf of < > nek5000-users at lists.mcs.anl.gov> > Antworten an: > Datum: Friday, August 26, 2016 at 7:22 PM > An: > Betreff: Re: [Nek5000-users] Error reading initial condition/drive force > data aborting in routine rdicdf. EXIT: rdicdf error > > Here is the .rea file. > > Thanks, > Saikat. > > *Saikat Mukherjee,* > *Graduate** Student,* > *Engineering Science and Mechanics,* > *BEAM,* > *Virginia Tech.* > > On Fri, Aug 26, 2016 at 1:00 PM, wrote: > >> Please post your .rea file. >> >> Von: on behalf of < >> nek5000-users at lists.mcs.anl.gov> >> Antworten an: >> Datum: Friday, August 26, 2016 at 6:18 PM >> An: >> Betreff: [Nek5000-users] Error reading initial condition/drive force >> data aborting in routine rdicdf. EXIT: rdicdf error >> >> Hello all, >> >> I have been using the same .rea and .usr files in a different RESTART >> calculation where the input file is ASCII. But when the file is BINARY, it >> is showing this error. I have changed the Parameters 66 and 67 but to no >> avail. Any help is appreciated. >> >> Thanks, >> >> Saikat. >> *Saikat Mukherjee,* >> *Graduate** Student,* >> *Engineering Science and Mechanics,* >> *BEAM,* >> *Virginia Tech.* >> _______________________________________________ Nek5000-users mailing >> list Nek5000-users at lists.mcs.anl.govhttps://lists.mcs.anl.gov/mai >> lman/listinfo/nek5000-users >> >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users at lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >> >> > _______________________________________________ Nek5000-users mailing > list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/ > mailman/listinfo/nek5000-users > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Aug 30 10:50:14 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 30 Aug 2016 15:50:14 +0000 Subject: [Nek5000-users] postx not working in Ubuntu 16.04 Message-ID: Hi All, I built Nek5000 in my Ubuntu 16.04 machine. I ran the eddy_uv case and it finished running successfully. However, when I tried to visualize the results with postx, I got nothing but a black window. Have anyone experienced this before? Thanks in advance !! -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Wed Aug 31 21:33:11 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 1 Sep 2016 10:33:11 +0800 (GMT+08:00) Subject: [Nek5000-users] Failed in Helmholtz error Message-ID: Hi! nek, I met the 'failed in Helmholtz' error when I was running the example case vortex breakdown in the file 'examples/vortex'. However I still got the .fld results of which the whole velocity field was static. Obviously the results were uncorrect. Do anyone know how to deal with it? Any help is appreciated. -- Xu Bolun, University of Science and Techonology of China, Hefei,Anhui,China -------------- next part -------------- An HTML attachment was scrubbed... URL: