[Nek5000-users] array reshape
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Thu Jul 12 10:06:18 CDT 2018
yes, this is the beauty (or ugliness) of Fortran 77. The same holds for
arrays passed in the argument list, for instance.
In fact, you could simply just increase the first index beyond lx1 even
having your first definition.
Philipp
On 2018-07-12 16:57, nek5000-users at lists.mcs.anl.gov wrote:
> Hi Philipp,
>
> If I have
>
> common /avg/ uvag(lx1,ly1,lz1,lelt)
>
>
> in one subroutine, then can I call all data stored in that array as
>
>
> common /avg/ uvag(lx1*ly1*lz1,lelt)
>
>
> in another subroutine? Will this subroutine be able to read the data stored in this
>
> (common /avg/ uvag(lx1*ly1*lz1,lelt) ) array?
>
>
>
>> On 12 Jul 2018, at 15:23, nek5000-users-request at lists.mcs.anl.gov wrote:
>>
>> 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. array reshape (nek5000-users at lists.mcs.anl.gov)
>> 2. Re: array reshape (nek5000-users at lists.mcs.anl.gov)
>> 3. Re: Number of Elements per Characteristic Length for DNS
>> (nek5000-users at lists.mcs.anl.gov)
>> 4. Re: Number of Elements per Characteristic Length for DNS
>> (nek5000-users at lists.mcs.anl.gov)
>> 5. Re: Number of Elements per Characteristic Length for DNS
>> (nek5000-users at lists.mcs.anl.gov)
>>
>>
>> ----------------------------------------------------------------------
>>
>> Message: 1
>> Date: Thu, 12 Jul 2018 12:07:42 +0000
>> From: nek5000-users at lists.mcs.anl.gov
>> To: "nek5000-users at lists.mcs.anl.gov"
>> <nek5000-users at lists.mcs.anl.gov>
>> Subject: [Nek5000-users] array reshape
>> Message-ID:
>> <mailman.6007.1531397285.86639.nek5000-users at lists.mcs.anl.gov>
>> Content-Type: text/plain; charset="iso-8859-1"
>>
>> Hi all,
>>
>>
>> How can I reshape this array uvag(lx1,ly1,lz1,lelt) into abc(lx1*ly1*lz1,lelt)?
>>
>>
>>
>> The following array reshape which works for fortran didn't work:
>>
>>
>> reshape(uavg, ( \ lx1*ly1*lz1 , lelt \ ))
>>
>>
>>
>> Any suggestions?
>>
>>
>> Thanks.
>> -------------- next part --------------
>> An HTML attachment was scrubbed...
>> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180712/85abe829/attachment-0001.html>
>>
>> ------------------------------
>>
>> Message: 2
>> Date: Thu, 12 Jul 2018 14:13:18 +0200
>> From: nek5000-users at lists.mcs.anl.gov
>> To: nek5000-users at lists.mcs.anl.gov
>> Subject: Re: [Nek5000-users] array reshape
>> Message-ID:
>> <mailman.6009.1531397608.86639.nek5000-users at lists.mcs.anl.gov>
>> Content-Type: text/plain; charset=windows-1252; format=flowed
>>
>> I dont think you need to do anything as the first index is the faster
>> varying in Fortran.
>> Philipp
>>
>> On 2018-07-12 14:07, nek5000-users at lists.mcs.anl.gov wrote:
>>> Hi all,
>>>
>>>
>>> How can I reshape this array?uvag(lx1,ly1,lz1,lelt) into
>>> abc(lx1*ly1*lz1,lelt)?
>>>
>>>
>>>
>>> The following array reshape which works for fortran didn't work:
>>>
>>>
>>> reshape(uavg, (? \ lx1*ly1*lz1?,?lelt?\?))
>>>
>>>
>>>
>>> Any suggestions?
>>>
>>>
>>> Thanks.
>>>
>>>
>>>
>>> _______________________________________________
>>> Nek5000-users mailing list
>>> Nek5000-users at lists.mcs.anl.gov
>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>
>>
>>
>> ------------------------------
>>
>> Message: 3
>> Date: Thu, 12 Jul 2018 08:56:18 -0400
>> From: nek5000-users at lists.mcs.anl.gov
>> To: nek5000-users at lists.mcs.anl.gov
>> Subject: Re: [Nek5000-users] Number of Elements per Characteristic
>> Length for DNS
>> Message-ID:
>> <mailman.6015.1531400181.86639.nek5000-users at lists.mcs.anl.gov>
>> Content-Type: text/plain; charset="utf-8"
>>
>> Hi Philipp,
>>
>> Thank you for getting back to me. Do you have any references that you could
>> provide so that I can look into this further?
>>
>> Thanks,
>>
>> Tommy
>>
>> On Thu, Jul 12, 2018 at 4:25 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
>>
>>> Hi,
>>> typically, in wall bounded turbulence you would design the mesh based on
>>> inner scaled quantities (i.e. based on the viscous length scale l*). This
>>> is what we have done in the paper by El Khoury, in accordance to studies in
>>> channel flow. TOwards the centre, one can then also use Kolmogorov units
>>> (not sure whether Taylor units are so useful). I am quite convinced that
>>> the El Khoury simulations are well-resolved DNS, which can also be seen
>>> when looking at spectra.
>>>
>>> Philipp
>>>
>>>
>>> On 2018-07-11 20:58, nek5000-users at lists.mcs.anl.gov wrote:
>>>
>>>> Hello All,
>>>>
>>>> I am currently investigating the turbulence in a pipe with periodic
>>>> boundary conditions. I am working on building my mesh and was planning on
>>>> using Re^(3/4) as the number of elements per characteristic length, but I
>>>> have found differing methods in literature detailing the number of elements
>>>> used for DNS. One paper in particular (G.K. El Khoury et al, Direct
>>>> Numerical Simulation of Turbulent Pipe Flow at Moderately High Reynolds
>>>> Numbers, Flow Turbulence Combust, 2013) uses Nek5000 to simulate turbulent
>>>> flow in a pipe at different Reynolds numbers. The number of elements
>>>> calculated in this paper seem to be in between the LES and DNS suggested
>>>> number of elements (element size between the Taylor microscale and the
>>>> Kolmogorov scale). Turbulent Flows by Pope however suggests a number of
>>>> elements that is slightly larger than Re^(3/4), and is dependent on the
>>>> Taylor microscale Reynolds number.
>>>>
>>>> Any suggestions on the correct way to determine the number of spectral
>>>> elements in a characteristic length would be appreciated.
>>>>
>>>> Regards,
>>>>
>>>> Tommy
>>>>
>>>>
>>>>
>>>> _______________________________________________
>>>> 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: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180712/3ca13d06/attachment-0001.html>
>>
>> ------------------------------
>>
>> Message: 4
>> Date: Thu, 12 Jul 2018 15:07:49 +0200
>> From: nek5000-users at lists.mcs.anl.gov
>> To: nek5000-users at lists.mcs.anl.gov
>> Subject: Re: [Nek5000-users] Number of Elements per Characteristic
>> Length for DNS
>> Message-ID:
>> <mailman.6017.1531400891.86639.nek5000-users at lists.mcs.anl.gov>
>> Content-Type: text/plain; charset=utf-8; format=flowed
>>
>> I guess the El Khoury paper might be a good starting point; we really
>> tried to make sure that the resolution is sufficient for good DNS.
>> Otherwise, there is a resolution discussion for spectral methods in the
>> paper by Coleman, Johnstone and Spalart:
>>
>> https://aip.scitation.org/doi/abs/10.1063/1.3247176
>>
>> Philipp
>>
>> On 2018-07-12 14:56, nek5000-users at lists.mcs.anl.gov wrote:
>>> Hi Philipp,
>>>
>>> Thank you for getting back to me. Do you have any references that you
>>> could provide so that I can look into this further?
>>>
>>> Thanks,
>>>
>>> Tommy
>>>
>>> On Thu, Jul 12, 2018 at 4:25 AM, <nek5000-users at lists.mcs.anl.gov
>>> <mailto:nek5000-users at lists.mcs.anl.gov>> wrote:
>>>
>>> Hi,
>>> typically, in wall bounded turbulence you would design the mesh
>>> based on inner scaled quantities (i.e. based on the viscous length
>>> scale l*). This is what we have done in the paper by El Khoury, in
>>> accordance to studies in channel flow. TOwards the centre, one can
>>> then also use Kolmogorov units (not sure whether Taylor units are so
>>> useful). I am quite convinced that the El Khoury simulations are
>>> well-resolved DNS, which can also be seen when looking at spectra.
>>>
>>> Philipp
>>>
>>>
>>> On 2018-07-11 20:58, nek5000-users at lists.mcs.anl.gov
>>> <mailto:nek5000-users at lists.mcs.anl.gov> wrote:
>>>
>>> Hello All,
>>>
>>> I am currently investigating the turbulence in a pipe with
>>> periodic boundary conditions. I am working on building my mesh
>>> and was planning on using Re^(3/4) as the number of elements per
>>> characteristic length, but I have found differing methods in
>>> literature detailing the number of elements used for DNS. One
>>> paper in particular (G.K. El Khoury et al, Direct Numerical
>>> Simulation of Turbulent Pipe Flow at Moderately High Reynolds
>>> Numbers,? Flow Turbulence Combust, 2013) uses Nek5000 to
>>> simulate turbulent flow in a pipe at different Reynolds numbers.
>>> The number of elements calculated in this paper seem to be in
>>> between the LES and DNS suggested number of elements (element
>>> size between the Taylor microscale and the Kolmogorov scale).
>>> Turbulent Flows by Pope however suggests a number of elements
>>> that is slightly larger than Re^(3/4), and is dependent on the
>>> Taylor microscale Reynolds number.
>>>
>>> Any suggestions on the correct way to determine the number of
>>> spectral elements in a characteristic length would be appreciated.
>>>
>>> Regards,
>>>
>>> Tommy
>>>
>>>
>>>
>>> _______________________________________________
>>> Nek5000-users mailing list
>>> Nek5000-users at lists.mcs.anl.gov
>>> <mailto:Nek5000-users at lists.mcs.anl.gov>
>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>> <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
>>>
>>> _______________________________________________
>>> Nek5000-users mailing list
>>> Nek5000-users at lists.mcs.anl.gov <mailto:Nek5000-users at lists.mcs.anl.gov>
>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>> <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
>>>
>>
>>
>> ------------------------------
>>
>> Message: 5
>> Date: Thu, 12 Jul 2018 10:22:15 -0400
>> From: nek5000-users at lists.mcs.anl.gov
>> To: nek5000-users at lists.mcs.anl.gov
>> Subject: Re: [Nek5000-users] Number of Elements per Characteristic
>> Length for DNS
>> Message-ID:
>> <mailman.6023.1531405402.86639.nek5000-users at lists.mcs.anl.gov>
>> Content-Type: text/plain; charset="utf-8"
>>
>> Hello Philipp,
>>
>> Thanks again for getting back to me and for sending these references.
>>
>> Regards,
>>
>> Tommy
>>
>> On Thu, Jul 12, 2018 at 9:07 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
>>
>>> I guess the El Khoury paper might be a good starting point; we really
>>> tried to make sure that the resolution is sufficient for good DNS.
>>> Otherwise, there is a resolution discussion for spectral methods in the
>>> paper by Coleman, Johnstone and Spalart:
>>>
>>> https://aip.scitation.org/doi/abs/10.1063/1.3247176
>>>
>>> Philipp
>>>
>>> On 2018-07-12 14:56, nek5000-users at lists.mcs.anl.gov wrote:
>>>
>>>> Hi Philipp,
>>>>
>>>> Thank you for getting back to me. Do you have any references that you
>>>> could provide so that I can look into this further?
>>>>
>>>> Thanks,
>>>>
>>>> Tommy
>>>>
>>>> On Thu, Jul 12, 2018 at 4:25 AM, <nek5000-users at lists.mcs.anl.gov
>>>> <mailto:nek5000-users at lists.mcs.anl.gov>> wrote:
>>>>
>>>> Hi,
>>>> typically, in wall bounded turbulence you would design the mesh
>>>> based on inner scaled quantities (i.e. based on the viscous length
>>>> scale l*). This is what we have done in the paper by El Khoury, in
>>>> accordance to studies in channel flow. TOwards the centre, one can
>>>> then also use Kolmogorov units (not sure whether Taylor units are so
>>>> useful). I am quite convinced that the El Khoury simulations are
>>>> well-resolved DNS, which can also be seen when looking at spectra.
>>>>
>>>> Philipp
>>>>
>>>>
>>>> On 2018-07-11 20:58, nek5000-users at lists.mcs.anl.gov
>>>> <mailto:nek5000-users at lists.mcs.anl.gov> wrote:
>>>>
>>>> Hello All,
>>>>
>>>> I am currently investigating the turbulence in a pipe with
>>>> periodic boundary conditions. I am working on building my mesh
>>>> and was planning on using Re^(3/4) as the number of elements per
>>>> characteristic length, but I have found differing methods in
>>>> literature detailing the number of elements used for DNS. One
>>>> paper in particular (G.K. El Khoury et al, Direct Numerical
>>>> Simulation of Turbulent Pipe Flow at Moderately High Reynolds
>>>> Numbers, Flow Turbulence Combust, 2013) uses Nek5000 to
>>>> simulate turbulent flow in a pipe at different Reynolds numbers.
>>>> The number of elements calculated in this paper seem to be in
>>>> between the LES and DNS suggested number of elements (element
>>>> size between the Taylor microscale and the Kolmogorov scale).
>>>> Turbulent Flows by Pope however suggests a number of elements
>>>> that is slightly larger than Re^(3/4), and is dependent on the
>>>> Taylor microscale Reynolds number.
>>>>
>>>> Any suggestions on the correct way to determine the number of
>>>> spectral elements in a characteristic length would be appreciated.
>>>>
>>>> Regards,
>>>>
>>>> Tommy
>>>>
>>>>
>>>>
>>>> _______________________________________________
>>>> Nek5000-users mailing list
>>>> Nek5000-users at lists.mcs.anl.gov
>>>> <mailto:Nek5000-users at lists.mcs.anl.gov>
>>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>> <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
>>>>
>>>> _______________________________________________
>>>> Nek5000-users mailing list
>>>> Nek5000-users at lists.mcs.anl.gov <mailto:Nek5000-users at lists.mc
>>>> s.anl.gov>
>>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>> <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: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180712/b5066fc5/attachment.html>
>>
>> ------------------------------
>>
>> Subject: Digest Footer
>>
>> _______________________________________________
>> 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 113, Issue 11
>> **********************************************
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
More information about the Nek5000-users
mailing list