From nek5000-users at lists.mcs.anl.gov Sun Mar 2 12:42:23 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sun, 2 Mar 2014 18:42:23 +0000 Subject: [Nek5000-users] Pipe example - Strange velocity/pressure signals Message-ID: Dear Neks, We are interested to do a smoothly deformed axisymmetric geometry, much like the pipe/stenosis example in the repository. So, we added hpts probes in the pipe example, and modified nothing else, but this turns out to give strange values initially. The first probe sits 3 diameters downstream of the stenosis throat (z=3), and gives U=0.48 at the first time step, despite zero initial condition. The pressure in this position shoots up to very high values too, so maybe a spurious pressure peak is the cause? Could you please help us to understand/prevent this phenomenon? The .usr, hpts.in and essential parts of .rea are attached. We are using revision 1008. This would be important for us as we aim to use Nek5000 to study transient as well as steady state dynamics of similar flows. Thank you in advance! Outi -------------- next part -------------- A non-text attachment was scrubbed... Name: stenosis.usr Type: application/octet-stream Size: 3967 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: hpts.in Type: application/octet-stream Size: 66 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: stenosis.rea Type: application/octet-stream Size: 4649 bytes Desc: not available URL: -------------- next part -------------- From nek5000-users at lists.mcs.anl.gov Sun Mar 2 15:02:25 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sun, 2 Mar 2014 15:02:25 -0600 (CST) Subject: [Nek5000-users] Pipe example - Strange velocity/pressure signals In-Reply-To: References: Message-ID: Dear Outi, What value of lx1 and lx2 were you using? I'm not observing the phenomenon you report. I have: lx1=6 lxd=9 lx2=lx1-2 lx3=lx2 ... Once I know your SIZE parameters I can retry. Best, Paul On Sun, 2 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: > Dear Neks, > > We are interested to do a smoothly deformed axisymmetric geometry, much > like the pipe/stenosis example in the repository. > So, we added hpts probes in the pipe example, and modified nothing else, > but this turns out to give strange values initially. The first probe sits > 3 diameters downstream of the stenosis throat (z=3), and gives U=0.48 at > the first time step, despite zero initial condition. The pressure in this > position shoots up to very high values too, so maybe a spurious pressure > peak is the cause? > > Could you please help us to understand/prevent this phenomenon? The .usr, > hpts.in and essential parts of .rea are attached. We are using revision > 1008. > This would be important for us as we aim to use Nek5000 to study > transient as well as steady state dynamics of similar flows. > > Thank you in advance! > Outi > > > > > _______________________________________________ > 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 Sun Mar 2 17:29:43 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sun, 2 Mar 2014 23:29:43 +0000 Subject: [Nek5000-users] Pipe example - Strange velocity/pressure signals In-Reply-To: References: Message-ID: Dear Paul, Thank you very much for looking into this. My parameters are as in the repository example, I didn't change them: lx1=6, lxd=10, lx2=lx1-2 lx3=lx1 The output after the first time step (dt=2.5e-2) of hpts.out is attached below. Please note the third column (z-velocity). 2.5000000E-02 -4.7626817E-05 -4.4872405E-05 4.8233404E-01 1.3486090E+02 2.5000000E-02 -6.0842050E-05 -1.8870027E-02 5.9476094E-01 1.7381205E+02 2.5000000E-02 2.0295515E-07 -1.5249089E-02 4.8242873E-01 1.3487912E+02 2.5000000E-02 -1.5349154E-02 -1.1545977E-05 4.8241771E-01 1.3487899E+02 2.5000000E-02 9.8682769E-03 -1.7779070E-02 5.7201214E-01 1.7367640E+02 Best regards, Outi On Mar 2, 2014, at 9:02 PM, nek5000-users at lists.mcs.anl.gov wrote: > > Dear Outi, > > What value of lx1 and lx2 were you using? > > I'm not observing the phenomenon you report. > I have: > > lx1=6 > lxd=9 > lx2=lx1-2 > lx3=lx2 > > ... > > Once I know your SIZE parameters I can retry. > > Best, Paul > > > On Sun, 2 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: > >> Dear Neks, >> We are interested to do a smoothly deformed axisymmetric geometry, much >> like the pipe/stenosis example in the repository. >> So, we added hpts probes in the pipe example, and modified nothing else, >> but this turns out to give strange values initially. The first probe sits >> 3 diameters downstream of the stenosis throat (z=3), and gives U=0.48 at >> the first time step, despite zero initial condition. The pressure in this >> position shoots up to very high values too, so maybe a spurious pressure >> peak is the cause? >> Could you please help us to understand/prevent this phenomenon? The .usr, >> hpts.in and essential parts of .rea are attached. We are using revision >> 1008. >> This would be important for us as we aim to use Nek5000 to study >> transient as well as steady state dynamics of similar flows. >> Thank you in advance! >> Outi >> _______________________________________________ >> 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 Sun Mar 2 18:05:01 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sun, 2 Mar 2014 18:05:01 -0600 (CST) Subject: [Nek5000-users] Pipe example - Strange velocity/pressure signals In-Reply-To: References: Message-ID: Hi Outi, This looks ok to me... I see in your userbc file: xp = x/0.5 yp = y/0.5 rr = xp*xp + yp*yp ux = 0.0 uy = 0.0 uz = 1 - rr ! Hagen-Poiseiulle flow, pipe diamter = 1. which would give a nonzero velocity after Step 1 ? Best, Paul On Sun, 2 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: > Dear Paul, > > Thank you very much for looking into this. My parameters are as in the repository example, I didn't change them: > lx1=6, > lxd=10, > lx2=lx1-2 > lx3=lx1 > > The output after the first time step (dt=2.5e-2) of hpts.out is attached below. Please note the third column (z-velocity). > > 2.5000000E-02 -4.7626817E-05 -4.4872405E-05 4.8233404E-01 1.3486090E+02 > 2.5000000E-02 -6.0842050E-05 -1.8870027E-02 5.9476094E-01 1.7381205E+02 > 2.5000000E-02 2.0295515E-07 -1.5249089E-02 4.8242873E-01 1.3487912E+02 > 2.5000000E-02 -1.5349154E-02 -1.1545977E-05 4.8241771E-01 1.3487899E+02 > 2.5000000E-02 9.8682769E-03 -1.7779070E-02 5.7201214E-01 1.7367640E+02 > > Best regards, > Outi > > > > On Mar 2, 2014, at 9:02 PM, nek5000-users at lists.mcs.anl.gov wrote: > >> >> Dear Outi, >> >> What value of lx1 and lx2 were you using? >> >> I'm not observing the phenomenon you report. >> I have: >> >> lx1=6 >> lxd=9 >> lx2=lx1-2 >> lx3=lx2 >> >> ... >> >> Once I know your SIZE parameters I can retry. >> >> Best, Paul >> >> >> On Sun, 2 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: >> >>> Dear Neks, >>> We are interested to do a smoothly deformed axisymmetric geometry, much >>> like the pipe/stenosis example in the repository. >>> So, we added hpts probes in the pipe example, and modified nothing else, >>> but this turns out to give strange values initially. The first probe sits >>> 3 diameters downstream of the stenosis throat (z=3), and gives U=0.48 at >>> the first time step, despite zero initial condition. The pressure in this >>> position shoots up to very high values too, so maybe a spurious pressure >>> peak is the cause? >>> Could you please help us to understand/prevent this phenomenon? The .usr, >>> hpts.in and essential parts of .rea are attached. We are using revision >>> 1008. >>> This would be important for us as we aim to use Nek5000 to study >>> transient as well as steady state dynamics of similar flows. >>> Thank you in advance! >>> Outi >>> _______________________________________________ >>> 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 Mon Mar 3 03:29:56 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 03 Mar 2014 10:29:56 +0100 Subject: [Nek5000-users] Pipe example - Strange velocity/pressure signals In-Reply-To: References: Message-ID: Hi Paul, As far as I can see, the userbc is set only at the inlet. The first probe point is at z=3, that is, 6 diameters downstream of the inlet (after transformation in usrdat2 from -3 to 10). What concerns me is that the z-velocity goes up to 0.48e-1 at z=3 at the first time step, when the inflow hasn't reached near this point yet. When I output the whole velocity field at the first time step (as .f00001) and look at it in paraview, the z-velocity has oscillations near the throat. Maybe at the end this is some transient phenomenon which doesn't matter, but I thought it was safest to check with you. Best, Outi Quoting nek5000-users at lists.mcs.anl.gov: > Hi Outi, > > This looks ok to me... I see in your userbc file: > > xp = x/0.5 > yp = y/0.5 > rr = xp*xp + yp*yp > > ux = 0.0 > uy = 0.0 > uz = 1 - rr ! Hagen-Poiseiulle flow, pipe diamter = 1. > > which would give a nonzero velocity after Step 1 ? > > Best, Paul > > > > > > > On Sun, 2 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: > >> Dear Paul, >> >> Thank you very much for looking into this. My parameters are as in >> the repository example, I didn't change them: >> lx1=6, >> lxd=10, >> lx2=lx1-2 >> lx3=lx1 >> >> The output after the first time step (dt=2.5e-2) of hpts.out is >> attached below. Please note the third column (z-velocity). >> >> 2.5000000E-02 -4.7626817E-05 -4.4872405E-05 4.8233404E-01 1.3486090E+02 >> 2.5000000E-02 -6.0842050E-05 -1.8870027E-02 5.9476094E-01 1.7381205E+02 >> 2.5000000E-02 2.0295515E-07 -1.5249089E-02 4.8242873E-01 1.3487912E+02 >> 2.5000000E-02 -1.5349154E-02 -1.1545977E-05 4.8241771E-01 1.3487899E+02 >> 2.5000000E-02 9.8682769E-03 -1.7779070E-02 5.7201214E-01 1.7367640E+02 >> >> Best regards, >> Outi >> >> >> >> On Mar 2, 2014, at 9:02 PM, nek5000-users at lists.mcs.anl.gov wrote: >> >>> >>> Dear Outi, >>> >>> What value of lx1 and lx2 were you using? >>> >>> I'm not observing the phenomenon you report. >>> I have: >>> >>> lx1=6 >>> lxd=9 >>> lx2=lx1-2 >>> lx3=lx2 >>> >>> ... >>> >>> Once I know your SIZE parameters I can retry. >>> >>> Best, Paul >>> >>> >>> On Sun, 2 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: >>> >>>> Dear Neks, >>>> We are interested to do a smoothly deformed axisymmetric geometry, much >>>> like the pipe/stenosis example in the repository. >>>> So, we added hpts probes in the pipe example, and modified nothing else, >>>> but this turns out to give strange values initially. The first probe sits >>>> 3 diameters downstream of the stenosis throat (z=3), and gives U=0.48 at >>>> the first time step, despite zero initial condition. The pressure in this >>>> position shoots up to very high values too, so maybe a spurious pressure >>>> peak is the cause? >>>> Could you please help us to understand/prevent this phenomenon? The .usr, >>>> hpts.in and essential parts of .rea are attached. We are using revision >>>> 1008. >>>> This would be important for us as we aim to use Nek5000 to study >>>> transient as well as steady state dynamics of similar flows. >>>> Thank you in advance! >>>> Outi >>>> _______________________________________________ >>>> 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 >> > _______________________________________________ > 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 Mar 3 06:16:57 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 3 Mar 2014 12:16:57 +0000 Subject: [Nek5000-users] Pipe example - Strange velocity/pressure signals In-Reply-To: References: , Message-ID: Hi Outi, Generally, the incompressibility constraint will enforce the flow to be nonzero after the first step, i.e., when you start to push water through a pipe it comes out the other end immediately with the incompressible model. Is this not what you're expecting for this case? 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, March 03, 2014 3:29 AM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] Pipe example - Strange velocity/pressure signals Hi Paul, As far as I can see, the userbc is set only at the inlet. The first probe point is at z=3, that is, 6 diameters downstream of the inlet (after transformation in usrdat2 from -3 to 10). What concerns me is that the z-velocity goes up to 0.48e-1 at z=3 at the first time step, when the inflow hasn't reached near this point yet. When I output the whole velocity field at the first time step (as .f00001) and look at it in paraview, the z-velocity has oscillations near the throat. Maybe at the end this is some transient phenomenon which doesn't matter, but I thought it was safest to check with you. Best, Outi Quoting nek5000-users at lists.mcs.anl.gov: > Hi Outi, > > This looks ok to me... I see in your userbc file: > > xp = x/0.5 > yp = y/0.5 > rr = xp*xp + yp*yp > > ux = 0.0 > uy = 0.0 > uz = 1 - rr ! Hagen-Poiseiulle flow, pipe diamter = 1. > > which would give a nonzero velocity after Step 1 ? > > Best, Paul > > > > > > > On Sun, 2 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: > >> Dear Paul, >> >> Thank you very much for looking into this. My parameters are as in >> the repository example, I didn't change them: >> lx1=6, >> lxd=10, >> lx2=lx1-2 >> lx3=lx1 >> >> The output after the first time step (dt=2.5e-2) of hpts.out is >> attached below. Please note the third column (z-velocity). >> >> 2.5000000E-02 -4.7626817E-05 -4.4872405E-05 4.8233404E-01 1.3486090E+02 >> 2.5000000E-02 -6.0842050E-05 -1.8870027E-02 5.9476094E-01 1.7381205E+02 >> 2.5000000E-02 2.0295515E-07 -1.5249089E-02 4.8242873E-01 1.3487912E+02 >> 2.5000000E-02 -1.5349154E-02 -1.1545977E-05 4.8241771E-01 1.3487899E+02 >> 2.5000000E-02 9.8682769E-03 -1.7779070E-02 5.7201214E-01 1.7367640E+02 >> >> Best regards, >> Outi >> >> >> >> On Mar 2, 2014, at 9:02 PM, nek5000-users at lists.mcs.anl.gov wrote: >> >>> >>> Dear Outi, >>> >>> What value of lx1 and lx2 were you using? >>> >>> I'm not observing the phenomenon you report. >>> I have: >>> >>> lx1=6 >>> lxd=9 >>> lx2=lx1-2 >>> lx3=lx2 >>> >>> ... >>> >>> Once I know your SIZE parameters I can retry. >>> >>> Best, Paul >>> >>> >>> On Sun, 2 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: >>> >>>> Dear Neks, >>>> We are interested to do a smoothly deformed axisymmetric geometry, much >>>> like the pipe/stenosis example in the repository. >>>> So, we added hpts probes in the pipe example, and modified nothing else, >>>> but this turns out to give strange values initially. The first probe sits >>>> 3 diameters downstream of the stenosis throat (z=3), and gives U=0.48 at >>>> the first time step, despite zero initial condition. The pressure in this >>>> position shoots up to very high values too, so maybe a spurious pressure >>>> peak is the cause? >>>> Could you please help us to understand/prevent this phenomenon? The .usr, >>>> hpts.in and essential parts of .rea are attached. We are using revision >>>> 1008. >>>> This would be important for us as we aim to use Nek5000 to study >>>> transient as well as steady state dynamics of similar flows. >>>> Thank you in advance! >>>> Outi >>>> _______________________________________________ >>>> 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 >> > _______________________________________________ > 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 Mon Mar 3 10:27:03 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 03 Mar 2014 17:27:03 +0100 Subject: [Nek5000-users] Pipe example - Strange velocity/pressure signals In-Reply-To: <20140303102956.Horde.ny7OJQ-FGe_ftOCXPICi8g8@webmail.mech.kth.se> References: <20140303102956.Horde.ny7OJQ-FGe_ftOCXPICi8g8@webmail.mech.kth.se> Message-ID: Thanks! Indeed this is what happens, we seem to have mixed up the mean flow development and the disturbance development.... So, I am sorry for bothering you but happy that all is correct. Outi Quoting Outi Tammisola : > Hi Paul, > > As far as I can see, the userbc is set only at the inlet. The first > probe point is at z=3, that is, 6 diameters downstream of the inlet > (after transformation in usrdat2 from -3 to 10). What concerns me is > that the z-velocity goes up to 0.48e-1 at z=3 at the first time > step, when the inflow hasn't reached near this point yet. When I > output the whole velocity field at the first time step (as .f00001) > and look at it in paraview, the z-velocity has oscillations near the > throat. Maybe at the end this is some transient phenomenon which > doesn't matter, but I thought it was safest to check with you. > > Best, > Outi > > Quoting nek5000-users at lists.mcs.anl.gov: > >> Hi Outi, >> >> This looks ok to me... I see in your userbc file: >> >> xp = x/0.5 >> yp = y/0.5 >> rr = xp*xp + yp*yp >> >> ux = 0.0 >> uy = 0.0 >> uz = 1 - rr ! Hagen-Poiseiulle flow, pipe diamter = 1. >> >> which would give a nonzero velocity after Step 1 ? >> >> Best, Paul >> >> >> >> >> >> >> On Sun, 2 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: >> >>> Dear Paul, >>> >>> Thank you very much for looking into this. My parameters are as in >>> the repository example, I didn't change them: >>> lx1=6, >>> lxd=10, >>> lx2=lx1-2 >>> lx3=lx1 >>> >>> The output after the first time step (dt=2.5e-2) of hpts.out is >>> attached below. Please note the third column (z-velocity). >>> >>> 2.5000000E-02 -4.7626817E-05 -4.4872405E-05 4.8233404E-01 1.3486090E+02 >>> 2.5000000E-02 -6.0842050E-05 -1.8870027E-02 5.9476094E-01 1.7381205E+02 >>> 2.5000000E-02 2.0295515E-07 -1.5249089E-02 4.8242873E-01 1.3487912E+02 >>> 2.5000000E-02 -1.5349154E-02 -1.1545977E-05 4.8241771E-01 1.3487899E+02 >>> 2.5000000E-02 9.8682769E-03 -1.7779070E-02 5.7201214E-01 1.7367640E+02 >>> >>> Best regards, >>> Outi >>> >>> >>> >>> On Mar 2, 2014, at 9:02 PM, nek5000-users at lists.mcs.anl.gov wrote: >>> >>>> >>>> Dear Outi, >>>> >>>> What value of lx1 and lx2 were you using? >>>> >>>> I'm not observing the phenomenon you report. >>>> I have: >>>> >>>> lx1=6 >>>> lxd=9 >>>> lx2=lx1-2 >>>> lx3=lx2 >>>> >>>> ... >>>> >>>> Once I know your SIZE parameters I can retry. >>>> >>>> Best, Paul >>>> >>>> >>>> On Sun, 2 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: >>>> >>>>> Dear Neks, >>>>> We are interested to do a smoothly deformed axisymmetric geometry, much >>>>> like the pipe/stenosis example in the repository. >>>>> So, we added hpts probes in the pipe example, and modified nothing else, >>>>> but this turns out to give strange values initially. The first probe sits >>>>> 3 diameters downstream of the stenosis throat (z=3), and gives U=0.48 at >>>>> the first time step, despite zero initial condition. The pressure in this >>>>> position shoots up to very high values too, so maybe a spurious pressure >>>>> peak is the cause? >>>>> Could you please help us to understand/prevent this phenomenon? The .usr, >>>>> hpts.in and essential parts of .rea are attached. We are using revision >>>>> 1008. >>>>> This would be important for us as we aim to use Nek5000 to study >>>>> transient as well as steady state dynamics of similar flows. >>>>> Thank you in advance! >>>>> Outi >>>>> _______________________________________________ >>>>> 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 >>> >> _______________________________________________ >> 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 Mar 4 09:01:32 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 04 Mar 2014 10:01:32 -0500 Subject: [Nek5000-users] Which version of NEK to download Message-ID: Hi, If I download the trunk, shall I get all the features in the branches: adapt_mesh moab_shared_tag_memory mvmesh Thanks, Tania From nek5000-users at lists.mcs.anl.gov Tue Mar 4 14:49:25 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 4 Mar 2014 20:49:25 +0000 Subject: [Nek5000-users] Which version of NEK to download In-Reply-To: References: Message-ID: Hi Tania, No there is no need -- I usually update only trunk/nek/ 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: Tuesday, March 04, 2014 9:01 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] Which version of NEK to download Hi, If I download the trunk, shall I get all the features in the branches: adapt_mesh moab_shared_tag_memory mvmesh Thanks, Tania _______________________________________________ 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 Mar 4 08:51:36 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 4 Mar 2014 15:51:36 +0100 Subject: [Nek5000-users] 2D axisymmetric problems Message-ID: Hi Nek's, I was playing with the example vortex2 and I was wondering if, in addition to the forcing terms ffy= temp*temp/y in userf and qvol = -uy*temp/y in userq , one should also add in userq the term: -param(2)*temp/(y*y) . Do we have to add this forcing term or the laplacian term is already modified by the flag IFAZIV ? Best regards, Francesco -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Mar 4 17:01:31 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 5 Mar 2014 00:01:31 +0100 Subject: [Nek5000-users] how to dump drag and lift in a file Message-ID: Hi, nek users, Instead of letting the function torque_calc dumping the drag lift and torque on screen, how can I have the them dumped to a file? Thank you. Best regards Tubois -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Mar 4 18:53:07 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 4 Mar 2014 18:53:07 -0600 (CST) Subject: [Nek5000-users] 2D axisymmetric problems In-Reply-To: References: Message-ID: Hi Francesco, Those terms are already in axhelm. Paul On Tue, 4 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: > Hi Nek's, > > I was playing with the example vortex2 and I was wondering if, in addition to the forcing terms > ffy= temp*temp/y in userf and qvol = -uy*temp/y in userq , > one should also add in userq the term: -param(2)*temp/(y*y) . > > Do we have to add this forcing term or the laplacian term is already modified by the flag IFAZIV ? > > Best regards, > > Francesco From nek5000-users at lists.mcs.anl.gov Tue Mar 4 19:01:20 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 4 Mar 2014 19:01:20 -0600 (CST) Subject: [Nek5000-users] how to dump drag and lift in a file In-Reply-To: References: Message-ID: Hi Tubois, The normal way would be to not have your results go to the screen (which is in any case slow), but to redirect stdout to a file, as the nekb and nekbmpi scripts do. ("b" is for background). I would normally do the following: nekbmpi blah 8 ... grep dragx logfile > t.x and then plot the contents of t.x. Here, blah is my session name. If you don't like that option (but, really, you should in anycase redirect stdout), you can copy the common block ctorq from navier5.f into your userchk subroutine and then have something like: common /ctorq/ dragx(0:maxobj),dragpx(0:maxobj),dragvx(0:maxobj) call torque_calc(scale,x0,.false.,.false.) if (nid.eq.0) write(14,1) istep,time,dragx(1) 1 format(i9,1p2e17.9) Best, Paul On Tue, 4 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: > Hi, nek users, > > Instead of letting the function torque_calc dumping the drag lift and > torque on screen, how can I have the them dumped to a file? Thank you. > > Best regards > > Tubois > From nek5000-users at lists.mcs.anl.gov Wed Mar 5 12:29:29 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 5 Mar 2014 19:29:29 +0100 Subject: [Nek5000-users] boundary condition issue? Message-ID: Hi nek users, I am encountering a strange issue. Build on the example of cylinder flow, I write a matlab script to generate the Cartesian mesh and a genbox script to generate the mesh around the cylinder. When the domain is small, the code runs. When I increases the domain size and cell numbers, the calculation blows up after several steps? See the picture attached here, there is a singular point on the boundary for the pressure. But the velociy is ok, since I imposed a uniform flow for inlet and top bottom. It also sometimes shows another error, the problem size is too large, even I run them in 32 procs instead of 4 as usual. Here attached the output of nek. Does someone have an idea? Best regards Tubois -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: out Type: application/octet-stream Size: 41510 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Capture du 2014-03-05 19:28:22.png Type: image/png Size: 20559 bytes Desc: not available URL: From nek5000-users at lists.mcs.anl.gov Wed Mar 5 15:08:33 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 5 Mar 2014 22:08:33 +0100 Subject: [Nek5000-users] boundary condition issue? In-Reply-To: References: Message-ID: I check the first time step, the pressure is more or less correct. From the second step, the singular point starts to show up. I changed the Cartesian mesh of the example back, I have the same problem. So the problem looks like to produced during the generation of the mesh around the cylinder. I used the import.rea from the ext_cyl example, so that I won't introduced any other problem. After several tests, I isolated the problem to two potential issues: 1, genbox to generate the mesh around cylinder, the script is below 2, prex, after importing the geometry and boundaries conditions, there is only the boundary condition of the cylinder left, which I imposed as 'wall' by using the script below I feel that it is very probably caused by prex, but I do not see where is the bug? I attached all scripts in the mail as well. Thank you! Tubois pretex << EOF new 1 READ PREVIOUS PARAMETERS box 1 BUILD FROM FILE box 10 IMPORT MESH toto y Would you like to displace existing elements in box? 1 END ELEMENTS 1 ACCEPT MATL,QVOL 1 SET BCs 11 SET ENTIRE LEVEL 3 WALL n Is this a lattice for hex-close-packed spheres?: 1 END LEVEL 1 ACCEPT B.C.'s 1 EXIT EOF genmap << EOF new .1 EOF import.rea #======================================================== 2 spatial dimension 1 number of fields # # comments # # #======================================================== # Y cylinder 4 -40 1 nelr,nel_theta,nelz .0 .0 x0,y0 - center of cylinder cccbb descriptors: c-cyl, o-oct, b-box (1 character + space) .5 .6 .7 .8 1. r0 r1 ... r_nelr 0 1 1 theta0/2pi theta1/2pi ratio E ,E ,E ,E , bc's (3 characters + comma) On Wed, Mar 5, 2014 at 7:29 PM, wrote: > Hi nek users, > > I am encountering a strange issue. Build on the example of cylinder flow, > I write a matlab script to generate the Cartesian mesh and a genbox script > to generate the mesh around the cylinder. When the domain is small, the > code runs. When I increases the domain size and cell numbers, the > calculation blows up after several steps? > > See the picture attached here, there is a singular point on the boundary > for the pressure. But the velociy is ok, since I imposed a uniform flow for > inlet and top bottom. > > It also sometimes shows another error, the problem size is too large, even > I run them in 32 procs instead of 4 as usual. > > Here attached the output of nek. > > Does someone have an idea? > > Best regards > > Tubois > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -- Universit? Pierre et Marie Curie CNRS-LIMSI Campus d'Orsay 91403 ORSAY -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Screenshot from 2014-03-05 21:35:52.png Type: image/png Size: 16555 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Screenshot from 2014-03-05 21:36:28.png Type: image/png Size: 17698 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: test.tar.gz Type: application/x-gzip Size: 11504 bytes Desc: not available URL: From nek5000-users at lists.mcs.anl.gov Fri Mar 7 09:12:34 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 7 Mar 2014 15:12:34 +0000 Subject: [Nek5000-users] Unstructured Mesh Message-ID: Dear All, I imported a .vtk mesh (unstructured, created in gmsh) in prex with no problems but I get the following error at run-time: " done :: read .rea file 0.74162E-01 sec setup mesh topology WARNINGa: Detected non-right-handed element. Number 1 C1-4: -0.6454E-03 -0.4845E-03 -0.4873E-03 -0.3265E-03 call exitt: dying ... " Does NEK5000 support unstructured meshes? I tried various geometries to make sure elements were not too distorted but I got similar errors. Thank you, JP -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Mar 7 12:46:39 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 7 Mar 2014 12:46:39 -0600 Subject: [Nek5000-users] Unstructured Mesh In-Reply-To: References: Message-ID: Hi JP, Is your mesh made up of fully hexahedral elements (quadrilaterals in 2D)? Nek supports unstructured meshes, but only these element types. If your mesh is fully hexahedral, then my suspicion is something went wrong with converting the .vtk to the .rea in terms of ordering the nodes for each element (Nek expects a certain ordering) Josh On Fri, Mar 7, 2014 at 9:12 AM, wrote: > Dear All, > > I imported a .vtk mesh (unstructured, created in gmsh) in prex with no > problems but I get the following error at run-time: > > " > done :: read .rea file 0.74162E-01 sec > > setup mesh topology > > WARNINGa: Detected non-right-handed element. > Number 1 C1-4: -0.6454E-03 -0.4845E-03 -0.4873E-03 -0.3265E-03 > > call exitt: dying ... > " > > Does NEK5000 support unstructured meshes? > > I tried various geometries to make sure elements were not too distorted > but I got similar errors. > > Thank you, > JP > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -- Josh Camp "All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Mar 7 17:18:15 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 7 Mar 2014 17:18:15 -0600 Subject: [Nek5000-users] Latest Nek5000 update. Message-ID: Hi, The latest update to the Nek5000 svn repository (version #1010) sets the default for all outpost files to be the binary f0000? format, unless otherwise specified in the .usr file. The .rea file parameter 66 and parameter 67 will ONLY be used by postx/postex to determine file format. That is, if you want a file format that is not f0000?, it MUST be changed in the .usr file, not the .rea file. As a reminder, the format is set according to the param(66)/param(67) values: 0 - ASCII ( fld format) +/- 6 - Binary ( .f00000 format) all others - Binary ( fld format) These values are now consistent in throughout postx and Nek5000. Best, Katherine -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Mar 7 19:11:39 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sat, 8 Mar 2014 02:11:39 +0100 Subject: [Nek5000-users] BlueGene assembly code Message-ID: Hi, I have tried to use the "BG" compilation flag in the PPLIST option to use the BlueGene specific code, both with "xl" and gnu compiler, but all of them fail saying that some opcodes are not defined fxcpmadd fxcsmadd fxpmul lfpdux stfpdux I guess they are bluegene specific codes, but it looks like the compilers that I used don't recognize them. Does anyone have a hint on what to use for BlueGene specific code? Thanks in advance! Best regards, -- Ismael From nek5000-users at lists.mcs.anl.gov Sat Mar 8 12:05:22 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sat, 8 Mar 2014 19:05:22 +0100 Subject: [Nek5000-users] Unstructured Mesh Message-ID: Hi Josh, Yes, my mesh is made up of quadrilaterals. I am importing a 2D mesh. (Importing a 3D .vtk file is not possible in prex from what I can see..) Any suggestion on how I can try fixing/controlling the ordering? Thank you, JP -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Mar 10 12:34:16 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 10 Mar 2014 12:34:16 -0500 Subject: [Nek5000-users] Unstructured Mesh In-Reply-To: References: Message-ID: Hi JP, Since you are importing your .vtk file into prex, I would assume that prex would order the nodes correctly when it outputs to a .rea/.re2 file. Perhaps there is something wrong with the .vtk file itself? I've never used this process before myself so maybe another user/developer could be more helpful. Josh On Sat, Mar 8, 2014 at 12:05 PM, wrote: > Hi Josh, > > Yes, my mesh is made up of quadrilaterals. I am importing a 2D mesh. > (Importing a 3D .vtk file is not possible in prex from what I can see..) > > Any suggestion on how I can try fixing/controlling the ordering? > > Thank you, > JP > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -- Josh Camp "All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Mar 11 15:16:46 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 11 Mar 2014 21:16:46 +0100 Subject: [Nek5000-users] Mesh generation Message-ID: Hallo, alle nek5000-users I want to use Nek5000 to make a channelflow-DNS. So I need to make cuboid-mesh first. Can I just use genbox to make 3D-mesh? Or I need to make 2D-mesh first and then use n2to3 to make 3D-mesh? Is there any tutorial for this case? Thank you very much. Best regards xwjtc -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Mar 11 15:41:40 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 11 Mar 2014 15:41:40 -0500 (CDT) Subject: [Nek5000-users] Mesh generation In-Reply-To: References: Message-ID: nek5_svn/examples/turbChannel/turbChannel.box gives an example. On Tue, 11 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: > Hallo, alle nek5000-users > I want to use Nek5000 to make a channelflow-DNS. So I need to make cuboid-mesh first. Can I just use genbox to make 3D-mesh? Or I need to make 2D-mesh first and then use n2to3 to make 3D-mesh? Is there any tutorial for this case? Thank you very much. > Best regards > xwjtc > From nek5000-users at lists.mcs.anl.gov Wed Mar 12 03:17:52 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 12 Mar 2014 09:17:52 +0100 Subject: [Nek5000-users] Mesh generation In-Reply-To: References: Message-ID: On 11/03/2014 21:16, nek5000-users at lists.mcs.anl.gov wrote: > Hallo, alle nek5000-users > I want to use Nek5000 to make a channelflow-DNS. So I need to make > cuboid-mesh first. Can I just use genbox to make 3D-mesh? Or I need to > make 2D-mesh first and then use n2to3 to make 3D-mesh? Is there any > tutorial for this case? Thank you very much. > Best regards > xwjtc > Hi Xwjtc, You can take a look at TurbChannel example in Nek. The explanation for the mesh is given in turbChannel.box file. Best regards, Kamal > > _______________________________________________ > 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 Wed Mar 12 08:31:48 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 12 Mar 2014 08:31:48 -0500 (CDT) Subject: [Nek5000-users] Unstructured Mesh In-Reply-To: References: Message-ID: Hi JP, Unstructured all-quad meshes have always been supported in nek. If you send me a compressed version of your .vtk, .rea, .usr, and SIZE files off list I'll take a look. Paul On Fri, 7 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: > Dear All, > > I imported a .vtk mesh (unstructured, created in gmsh) in prex with no > problems but I get the following error at run-time: > > " > done :: read .rea file 0.74162E-01 sec > > setup mesh topology > > WARNINGa: Detected non-right-handed element. > Number 1 C1-4: -0.6454E-03 -0.4845E-03 -0.4873E-03 -0.3265E-03 > > call exitt: dying ... > " > > Does NEK5000 support unstructured meshes? > > I tried various geometries to make sure elements were not too distorted but > I got similar errors. > > Thank you, > JP > From nek5000-users at lists.mcs.anl.gov Thu Mar 13 20:43:51 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 13 Mar 2014 21:43:51 -0400 Subject: [Nek5000-users] Specifying boundary condition for Annulus Geometry in box file Message-ID: Hi All, We are trying to build a mesh to simulate external cylinder flow, similar to the one given in the example "ext_cyl". Searching though the documentation, we understand that it is possible to generate several boxes of mesh using genbox and combine using pretex to get the final mesh. We are having difficulty in interpreting the boundary condition (BC) strings in the input to files to genbox. It is clear for a rectangular box, the BC string follows the order West, East, South, and North in 2D. What is the analogous interpretation for a circular annulus geometry? such as for the one defined by the box files given in "primer.pdf" ( http://www.mcs.anl.gov/~fischer/primer.pdf), which is repeated below: x2d.rea 2 spatial dimension 1 number of fields # # comments # # #======================================================== # Y cYlinder 3 -24 1 nelr,nel_theta,nelz .5 .3 x0,y0 - center of cylinder ccbb descriptors: c-cyl, o-oct, b-box (1 character + space) .5 .55 .7 .8 r0 r1 ... r_nelr 0 1 1 theta0/2pi theta1/2pi ratio v ,W ,E ,E , bc's (3 characters + comma) Does "v" corresponds to BC at r0 and W correspond to BC at r_nelr in r-direction and the two E's correspond to BC at theta0 and theta1? Thanks in advance for your help. Regards, Ravi -- Ravi Kumar T R University of Illinois at Urbana-Champaign ph : +1 217 778 7538 http://lndvl.mechse.illinois.edu/ -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Mar 14 03:21:03 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 14 Mar 2014 09:21:03 +0100 Subject: [Nek5000-users] Specifying boundary condition for Annulus Geometry in box file In-Reply-To: References: Message-ID: Hi Ravi, I came across this issue recently as well. I found that you are correct about the first and second BC. However, the third and the fourth, I do not know what they are. I tried to find it with prex, but nothing turns out. However, when I merged it with another mesh, it played a role, which modified the whole mesh BC, so I have a bad global BC. I do not see where is the problem. I think there is some isssue with pretex. Cheers, Tubois On Fri, Mar 14, 2014 at 2:43 AM, wrote: > Hi All, > > We are trying to build a mesh to simulate external cylinder flow, similar > to the one given in the example "ext_cyl". Searching though the > documentation, we understand that it is possible to generate several boxes > of mesh using genbox and combine using pretex to get the final mesh. > > We are having difficulty in interpreting the boundary condition (BC) > strings in the input to files to genbox. It is clear for a rectangular box, > the BC string follows the order West, East, South, and North in 2D. What is > the analogous interpretation for a circular annulus geometry? such as for > the one defined by the box files given in "primer.pdf" ( > http://www.mcs.anl.gov/~fischer/primer.pdf), which is repeated below: > > x2d.rea > 2 spatial dimension > 1 number of fields > # > # comments > # > # > #======================================================== > # > Y cYlinder > 3 -24 1 nelr,nel_theta,nelz > .5 .3 x0,y0 - center of cylinder > ccbb descriptors: c-cyl, o-oct, b-box (1 character + space) > .5 .55 .7 .8 r0 r1 ... r_nelr > 0 1 1 theta0/2pi theta1/2pi ratio > v ,W ,E ,E , bc's (3 characters + comma) > > > Does "v" corresponds to BC at r0 and W correspond to BC at r_nelr in > r-direction and the two E's correspond to BC at theta0 and theta1? > > Thanks in advance for your help. > > Regards, > Ravi > > > > > > > > > > -- > Ravi Kumar T R > > University of Illinois at Urbana-Champaign > ph : +1 217 778 7538 > http://lndvl.mechse.illinois.edu/ > > _______________________________________________ > 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 Mar 14 06:32:49 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 14 Mar 2014 11:32:49 +0000 Subject: [Nek5000-users] Unstructured Mesh Message-ID: Hi Paul, I have attached the compressed .vtk, .rea, .usr and SIZE files of a simple test with simple geometry and default 2D settings that I am trying to carry out where I get the same errors. Just in case, you can also download the file from the following link: https://www.wetransfer.com/downloads/c2cf981cf76c40ad1d942e6f4aaf51c920140314112847/31b0fb8abb3b3f369d390512487c88d920140314112847/9e4167 Thank you very much, JP -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: simpletest.tar.gz Type: application/x-gzip Size: 54148 bytes Desc: not available URL: From nek5000-users at lists.mcs.anl.gov Fri Mar 14 06:53:58 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 14 Mar 2014 17:23:58 +0530 Subject: [Nek5000-users] What happens when we use p093, p094, p095 Message-ID: Hello I want to understand what these parameters do p093, p094, p095 Where can I read the theory behind these ? Thanks praveen -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Mar 14 09:54:38 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 14 Mar 2014 14:54:38 +0000 Subject: [Nek5000-users] Non conforming mesh Message-ID: Dear users and developers, are non conforming elements currently supported by Nek5000? I have built a rea file with pretex starting from a non conforming vtk mesh. Everything goes well until I try to run the solver, then I get the "Error in generalev" message. Looking at the rea file, I found that the corresponding non conforming edges are marked with the 'J' and 'SP' labels, whose meaning is not reported at the wiki page: http://nek5000.mcs.anl.gov/index.php/Boundary_Definition Thank you in advance. Best regards, Giuseppe From nek5000-users at lists.mcs.anl.gov Fri Mar 14 13:22:03 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 14 Mar 2014 14:22:03 -0400 Subject: [Nek5000-users] Conjugate heat transfer problem Message-ID: Hi, Neks I want to set up a 3-D conjugate heat transfer problem. The procedure I used is as below: 1) generate a mesh for fluid domain. 2) generate a mesh for solid domain. 3)merge the element from 1) into 2). I have one question during the steps 1 and 2. What the boundary condition I need to set in step 1 and 2 at the interface between fluid and solid? Mu Xu -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Mar 14 14:42:01 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 14 Mar 2014 14:42:01 -0500 Subject: [Nek5000-users] Specifying boundary condition for Annulus Geometry in box file In-Reply-To: References: Message-ID: Hi Ravi, You are correct, for annulus geometry, the BC's used correspond to r_0, r_nelr, theta0, theta1 You don't see theta0 and theta1 when the geometry is a full cylinder (theta0=0 to theta1=1) since the theta edges will connect to one another and essential be 'E ' BC. However, if you were working with a half or quarter cylinder, these would be relevant and visible in prex. Similarly, if you run genbox for a full cylinder and put a BC that is not 'E ' or blank (' ') for theta0 or theta1, genbox will put this BC in the .rea file. This could then result in an incorrect BC. I can look at setting up a warning message within genbox for this situation. Best, Katie On Fri, Mar 14, 2014 at 3:21 AM, wrote: > Hi Ravi, > > I came across this issue recently as well. I found that you are correct > about the first and second BC. However, the third and the fourth, I do not > know what they are. I tried to find it with prex, but nothing turns out. > > However, when I merged it with another mesh, it played a role, which > modified the whole mesh BC, so I have a bad global BC. I do not see where > is the problem. I think there is some isssue with pretex. > > Cheers, > > Tubois > > On Fri, Mar 14, 2014 at 2:43 AM, wrote: > >> Hi All, >> >> We are trying to build a mesh to simulate external cylinder flow, similar >> to the one given in the example "ext_cyl". Searching though the >> documentation, we understand that it is possible to generate several boxes >> of mesh using genbox and combine using pretex to get the final mesh. >> >> We are having difficulty in interpreting the boundary condition (BC) >> strings in the input to files to genbox. It is clear for a rectangular box, >> the BC string follows the order West, East, South, and North in 2D. What is >> the analogous interpretation for a circular annulus geometry? such as for >> the one defined by the box files given in "primer.pdf" ( >> http://www.mcs.anl.gov/~fischer/primer.pdf), which is repeated below: >> >> x2d.rea >> 2 spatial dimension >> 1 number of fields >> # >> # comments >> # >> # >> #======================================================== >> # >> Y cYlinder >> 3 -24 1 nelr,nel_theta,nelz >> .5 .3 x0,y0 - center of cylinder >> ccbb descriptors: c-cyl, o-oct, b-box (1 character + space) >> .5 .55 .7 .8 r0 r1 ... r_nelr >> 0 1 1 theta0/2pi theta1/2pi ratio >> v ,W ,E ,E , bc's (3 characters + comma) >> >> >> Does "v" corresponds to BC at r0 and W correspond to BC at r_nelr in >> r-direction and the two E's correspond to BC at theta0 and theta1? >> >> Thanks in advance for your help. >> >> Regards, >> Ravi >> >> >> >> >> >> >> >> >> >> -- >> Ravi Kumar T R >> >> University of Illinois at Urbana-Champaign >> ph : +1 217 778 7538 >> http://lndvl.mechse.illinois.edu/ >> >> _______________________________________________ >> 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 Mar 14 16:14:07 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 14 Mar 2014 16:14:07 -0500 (CDT) Subject: [Nek5000-users] What happens when we use p093, p094, p095 In-Reply-To: References: Message-ID: Hi Praveen, Fischer, P.F., "Projection techniques for iterative solution of Ax=b with successive right-hand sides." Comp. Meth. in Appl. Mech 163 (1998) which is on my webpage. Paul On Fri, 14 Mar 2014, nek5000-users at lists.mcs.anl.gov wrote: > Hello > > I want to understand what these parameters do > > p093, p094, p095 > > Where can I read the theory behind these ? > > Thanks > praveen > From nek5000-users at lists.mcs.anl.gov Fri Mar 14 17:15:05 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 14 Mar 2014 22:15:05 +0000 Subject: [Nek5000-users] Conjugate heat transfer problem In-Reply-To: References: Message-ID: Hi Xu, in Conjugate Heat Transfer problems there is no boundary condition for the energy equation at the interface between fluid and solid. The interface behaves like the one of any other internal element, which means that the boundary is just marked with 'E '. I worked a lot with Conjugate Heat Transfer. If you need help don't hesitate to contact me. Cheers, Andrea. ________________________________ 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: Friday, March 14, 2014 7:22 PM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] Conjugate heat transfer problem Hi, Neks I want to set up a 3-D conjugate heat transfer problem. The procedure I used is as below: 1) generate a mesh for fluid domain. 2) generate a mesh for solid domain. 3)merge the element from 1) into 2). I have one question during the steps 1 and 2. What the boundary condition I need to set in step 1 and 2 at the interface between fluid and solid? Mu Xu From nek5000-users at lists.mcs.anl.gov Mon Mar 17 22:24:30 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 17 Mar 2014 23:24:30 -0400 Subject: [Nek5000-users] Specifying boundary condition for Annulus Geometry in box file In-Reply-To: References: Message-ID: Thanks Katie and Tubois for your response. That clarifies my doubts. I am able to generate (although fairly coarse) mesh using boxes with annulus geometry for cylinder and merge them with rectangular mesh boxes. I am working on to get the desired mesh refinement. Regards Ravi On Thu, Mar 13, 2014 at 9:43 PM, Ravi Tumkur wrote: > Hi All, > > We are trying to build a mesh to simulate external cylinder flow, similar > to the one given in the example "ext_cyl". Searching though the > documentation, we understand that it is possible to generate several boxes > of mesh using genbox and combine using pretex to get the final mesh. > > We are having difficulty in interpreting the boundary condition (BC) > strings in the input to files to genbox. It is clear for a rectangular box, > the BC string follows the order West, East, South, and North in 2D. What is > the analogous interpretation for a circular annulus geometry? such as for > the one defined by the box files given in "primer.pdf" ( > http://www.mcs.anl.gov/~fischer/primer.pdf), which is repeated below: > > x2d.rea > 2 spatial dimension > 1 number of fields > # > # comments > # > # > #======================================================== > # > Y cYlinder > 3 -24 1 nelr,nel_theta,nelz > .5 .3 x0,y0 - center of cylinder > ccbb descriptors: c-cyl, o-oct, b-box (1 character + space) > .5 .55 .7 .8 r0 r1 ... r_nelr > 0 1 1 theta0/2pi theta1/2pi ratio > v ,W ,E ,E , bc's (3 characters + comma) > > > Does "v" corresponds to BC at r0 and W correspond to BC at r_nelr in > r-direction and the two E's correspond to BC at theta0 and theta1? > > Thanks in advance for your help. > > Regards, > Ravi > > > > > > > > > > -- > Ravi Kumar T R > > University of Illinois at Urbana-Champaign > ph : +1 217 778 7538 > http://lndvl.mechse.illinois.edu/ > -- Ravi Kumar T R University of Illinois at Urbana-Champaign ph : +1 217 778 7538 http://lndvl.mechse.illinois.edu/ -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Mar 24 10:20:33 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 24 Mar 2014 11:20:33 -0400 Subject: [Nek5000-users] Conjugate heat transfer problem In-Reply-To: References: Message-ID: Hi, Andrea Thank you for your help. With your instruction, I can build a domain for conjugate heat transfer problem now. However the problem is always crash in the middle of running. I think the problem is during the Helmholtz solver for velocity which is controlled by flag P94 in .rea file. Flag P95 controls the Helmholtz solver for pressure. If I set both P94 and P95 to zero, the case can run but the result is not right. If I set P94 to positive integer and P95 to zero which is I used to, the case is crash. If I set P95 to positive integer and P94 to zero, the case can run and the result looks fine. I have never use this setting before. I want to know is it fine to run the conjugate heat transfer problem with this setting which is P94 to zero and P95 to positive integer. Mu Xu On Fri, Mar 14, 2014 at 6:15 PM, wrote: > Hi Xu, > > in Conjugate Heat Transfer problems there is no boundary condition for the > energy equation at the interface between fluid and solid. The interface > behaves like the one of any other internal element, which means that the > boundary is just marked with 'E '. > I worked a lot with Conjugate Heat Transfer. If you need help don't > hesitate to contact me. > > Cheers, > Andrea. > ________________________________ > 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: Friday, March 14, 2014 7:22 PM > To: nek5000-users at lists.mcs.anl.gov > Subject: [Nek5000-users] Conjugate heat transfer problem > > Hi, Neks > > I want to set up a 3-D conjugate heat transfer problem. The procedure I > used is as below: > 1) generate a mesh for fluid domain. > 2) generate a mesh for solid domain. > 3)merge the element from 1) into 2). > > I have one question during the steps 1 and 2. What the boundary condition > I need to set in step 1 and 2 at the interface between fluid and solid? > > Mu Xu > _______________________________________________ > 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 Mon Mar 24 10:33:43 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 24 Mar 2014 15:33:43 +0000 Subject: [Nek5000-users] Conjugate heat transfer problem In-Reply-To: References: , Message-ID: Mu, Set p94 and 95 to zero. That will eliminate that issue. The result not being right is a separate issue. 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, March 24, 2014 10:20 AM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] Conjugate heat transfer problem Hi, Andrea Thank you for your help. With your instruction, I can build a domain for conjugate heat transfer problem now. However the problem is always crash in the middle of running. I think the problem is during the Helmholtz solver for velocity which is controlled by flag P94 in .rea file. Flag P95 controls the Helmholtz solver for pressure. If I set both P94 and P95 to zero, the case can run but the result is not right. If I set P94 to positive integer and P95 to zero which is I used to, the case is crash. If I set P95 to positive integer and P94 to zero, the case can run and the result looks fine. I have never use this setting before. I want to know is it fine to run the conjugate heat transfer problem with this setting which is P94 to zero and P95 to positive integer. Mu Xu On Fri, Mar 14, 2014 at 6:15 PM, > wrote: Hi Xu, in Conjugate Heat Transfer problems there is no boundary condition for the energy equation at the interface between fluid and solid. The interface behaves like the one of any other internal element, which means that the boundary is just marked with 'E '. I worked a lot with Conjugate Heat Transfer. If you need help don't hesitate to contact me. Cheers, Andrea. ________________________________ 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: Friday, March 14, 2014 7:22 PM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] Conjugate heat transfer problem Hi, Neks I want to set up a 3-D conjugate heat transfer problem. The procedure I used is as below: 1) generate a mesh for fluid domain. 2) generate a mesh for solid domain. 3)merge the element from 1) into 2). I have one question during the steps 1 and 2. What the boundary condition I need to set in step 1 and 2 at the interface between fluid and solid? Mu Xu _______________________________________________ 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 Mon Mar 24 10:28:36 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 24 Mar 2014 15:28:36 +0000 Subject: [Nek5000-users] Problem with makenek for a moving bc Message-ID: Hi all, We are trying to compute 2-D laminar unsteady flow past a cylinder that moves perpendicularly to the mean flow direction, with the cylinder velocity being sinusoidal in time. We referred to the moving cylinder case in the attachment in this post: https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2012-February/001702.html . We used our own mesh file and edited parameters and functions in rea and usr files according to above example. The problem is, "makenek" exited with error. Please find in attachment the file "compiler.out". The last part of the info are: ... obj/mvmesh.o: In function `u_meshv_': mvmesh.f:(.text+0x0): multiple definition of `u_meshv_' obj/mvmesh.o:mvmesh.f:(.text+0x0): first defined here obj/mvmesh.o: In function `uservp_': mvmesh.f:(.text+0x70): multiple definition of `uservp_' obj/mvmesh.o:mvmesh.f:(.text+0x70): first defined here ... ( And many other multiple definition of functions. ) drive2.f:(.text+0x56bc): undefined reference to `updcoor_' obj/drive2.o: In function `meshv_': drive2.f:(.text+0x6671): undefined reference to `mvbdry_' ... ( And a few more undefined reference.) collect2: error: ld returned 1 exit status make: *** [nek5000] Error 1 Does anybody know why this happened or how to solve it? Thank you all in advance for your help. Regards, Li -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: compiler.out Type: application/octet-stream Size: 152385 bytes Desc: compiler.out URL: From nek5000-users at lists.mcs.anl.gov Mon Mar 24 13:12:03 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 24 Mar 2014 18:12:03 +0000 Subject: [Nek5000-users] Problem with makenek for a moving bc In-Reply-To: References: Message-ID: Hi Li, My usual approach to resolving this is to retry makenek with a .usr file that is known to work -- i.e., any of the .usr files from any of the example directories. Also, it never hurts to first type makenek clean ...just in case an earlier compile messed up and there are remnants of previously generated .o files in the obj directory. 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, March 24, 2014 10:28 AM To: nek5000-users at lists.mcs.anl.gov Cc: Pearlstein, Arne J; tumkur.ravikumar at gmail.com Subject: [Nek5000-users] Problem with makenek for a moving bc Hi all, We are trying to compute 2-D laminar unsteady flow past a cylinder that moves perpendicularly to the mean flow direction, with the cylinder velocity being sinusoidal in time. We referred to the moving cylinder case in the attachment in this post: https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2012-February/001702.html . We used our own mesh file and edited parameters and functions in rea and usr files according to above example. The problem is, "makenek" exited with error. Please find in attachment the file "compiler.out". The last part of the info are: ... obj/mvmesh.o: In function `u_meshv_': mvmesh.f:(.text+0x0): multiple definition of `u_meshv_' obj/mvmesh.o:mvmesh.f:(.text+0x0): first defined here obj/mvmesh.o: In function `uservp_': mvmesh.f:(.text+0x70): multiple definition of `uservp_' obj/mvmesh.o:mvmesh.f:(.text+0x70): first defined here ... ( And many other multiple definition of functions. ) drive2.f:(.text+0x56bc): undefined reference to `updcoor_' obj/drive2.o: In function `meshv_': drive2.f:(.text+0x6671): undefined reference to `mvbdry_' ... ( And a few more undefined reference.) collect2: error: ld returned 1 exit status make: *** [nek5000] Error 1 Does anybody know why this happened or how to solve it? Thank you all in advance for your help. Regards, Li -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Mar 25 10:04:26 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 25 Mar 2014 15:04:26 +0000 Subject: [Nek5000-users] Problem with makenek for a moving bc In-Reply-To: References: , Message-ID: Hi Paul, Thank you so much for your solution. I used the same usr file from the previous example and it worked. Regards, Li ________________________________ 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, March 24, 2014 13:12 To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] Problem with makenek for a moving bc Hi Li, My usual approach to resolving this is to retry makenek with a .usr file that is known to work -- i.e., any of the .usr files from any of the example directories. Also, it never hurts to first type makenek clean ...just in case an earlier compile messed up and there are remnants of previously generated .o files in the obj directory. 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, March 24, 2014 10:28 AM To: nek5000-users at lists.mcs.anl.gov Cc: Pearlstein, Arne J; tumkur.ravikumar at gmail.com Subject: [Nek5000-users] Problem with makenek for a moving bc Hi all, We are trying to compute 2-D laminar unsteady flow past a cylinder that moves perpendicularly to the mean flow direction, with the cylinder velocity being sinusoidal in time. We referred to the moving cylinder case in the attachment in this post: https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2012-February/001702.html . We used our own mesh file and edited parameters and functions in rea and usr files according to above example. The problem is, "makenek" exited with error. Please find in attachment the file "compiler.out". The last part of the info are: ... obj/mvmesh.o: In function `u_meshv_': mvmesh.f:(.text+0x0): multiple definition of `u_meshv_' obj/mvmesh.o:mvmesh.f:(.text+0x0): first defined here obj/mvmesh.o: In function `uservp_': mvmesh.f:(.text+0x70): multiple definition of `uservp_' obj/mvmesh.o:mvmesh.f:(.text+0x70): first defined here ... ( And many other multiple definition of functions. ) drive2.f:(.text+0x56bc): undefined reference to `updcoor_' obj/drive2.o: In function `meshv_': drive2.f:(.text+0x6671): undefined reference to `mvbdry_' ... ( And a few more undefined reference.) collect2: error: ld returned 1 exit status make: *** [nek5000] Error 1 Does anybody know why this happened or how to solve it? Thank you all in advance for your help. Regards, Li -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Wed Mar 26 01:31:28 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 26 Mar 2014 07:31:28 +0100 Subject: [Nek5000-users] Unstructured Mesh Message-ID: Dear Paul/All, Any possible updates on the error? I have tried several other geometries but unfortunately I get the same error. Thank you, JP -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Mar 27 09:34:33 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 27 Mar 2014 15:34:33 +0100 Subject: [Nek5000-users] Problem with passive scalar fields - subroutine outfld Message-ID: Dear users and developers, I'm simulating a 3D conjugate heat transfer case, considering temperature as a passive scalar. When I simulate only one temperature field everything is ok. If I add more passive scalar fields (each one identical to the first one), after the subroutine prepost.f/outfld is called, then some of these fields are not good. I ran up to 8 additional fields, and results are the followings ("- -" means that the field is not good): ldimt T PS1 PS2 PS3 PS4 PS5 PS6 PS7 PS8 1 ok 2 ok -- 3 ok -- ok 4 ok -- ok ok 5 ok -- -- -- ok 6 ok -- -- -- ok ok 7 ok -- -- -- ok ok ok 8 ok -- -- -- ok ok ok ok 9 ok -- -- -- ok ok ok ok ok I'm using PN-PN. Nothing changes with the number of mesh points or the number of steps or the spatial polynomial order. I saw differences with Visit, with mean temperature profiles, and with average temperature in fluid and solid domain. Differences are noticeable also in logfile. Consider the following example, in which the first field file is dumped at step 250; at step 251 the changes are noticeable in the number of iterations per field only for PS1, PS2 and PS3. :: Step 250, t= 3.1936627E+00, DT= 1.0935650E-02, C= 0.571 7.8515E+01 2.6463E-01 Solving for heat 250 Hmholtz TEMP: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 1: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 2: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 3: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 4: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 5: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 6: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 7: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 8: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 3.1937E+00 2.0457E-01 Heat done Solving for fluid :: call outfld: ifpsco: T 250 3.1937E+00 Write checkpoint: 0 250 OPEN: cjh_tc0.f00001 250 3.1937E+00 done :: Write checkpoint file size = 5.6 MB avg data-throughput = 367.7MB/s io-nodes = 1 Step 251, t= 3.2045983E+00, DT= 1.0935650E-02, C= 0.566 7.8809E+01 2.9378E-01 Solving for heat 251 Hmholtz TEMP: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 1: 74 9.1166E-04 1.1888E+02 1.0000E-03 251 Hmholtz PS 2: 85 9.4452E-04 9.1977E+02 1.0000E-03 251 Hmholtz PS 3: 91 8.5561E-04 2.3562E+03 1.0000E-03 251 Hmholtz PS 4: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 5: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 6: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 7: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 8: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 3.2046E+00 3.6832E-01 Heat done :: Looking to the outfld subroutine I wasn't able to figure out the reason for that. Does anyone know how to solve this problem? Thanks in advance. Relevant parts of logfile, SIZE, .usr, .box and .rea files are attached at the end of the mail. Best Regards, Donato ________________________________________________________________________ ________ SIZE file parameter (ldimt= 9) ! passive scalars + T ________________________________________________________________________ ________ *.box file base.rea 3 spatial dimension ( < 0 --> generate .rea/.re2 pair) 10 number of fields #====== Box -10 14 -5 nelx,nely,nelz for Box 0. 6. 1. x0,x1,gain -1.00 -0.98 -0.94 -0.86 -0.70 -0.50 -0.25 -0.00 0.25 0.50 0.70 0.86 0.94 0.98 1.00 y0,y1,...,yn 0. 3. 1. z0,z1,gain P ,P ,W ,W ,P ,P vel BCs P ,P ,E ,E ,P ,P temp BCs P ,P ,E ,E ,P ,P PS1 BCs P ,P ,E ,E ,P ,P PS2 BCs P ,P ,E ,E ,P ,P PS3 BCs P ,P ,E ,E ,P ,P PS4 BCs P ,P ,E ,E ,P ,P PS5 BCs P ,P ,E ,E ,P ,P PS6 BCs P ,P ,E ,E ,P ,P PS7 BCs P ,P ,E ,E ,P ,P PS8 BCs Box_up -10 7 -5 nelx,nely,nelz for Box_up 0. 6. 1. x0,x1,gain 1.00 1.02 1.06 1.14 1.30 1.50 1.75 2.00 y0,y1,...,yn 0. 3. 1. z0,z1,gain , , , , , vel BCs P ,P ,I ,E ,P ,P temp BCs P ,P ,I ,E ,P ,P PS1 BCs P ,P ,I ,E ,P ,P PS2 BCs P ,P ,I ,E ,P ,P PS3 BCs P ,P ,I ,E ,P ,P PS4 BCs P ,P ,I ,E ,P ,P PS5 BCs P ,P ,I ,E ,P ,P PS6 BCs P ,P ,I ,E ,P ,P PS7 BCs P ,P ,I ,E ,P ,P PS8 BCs Box_down -10 7 -5 nelx,nely,nelz for Box_down 0. 6. 1. x0,x1,gain -2.00 -1.75 -1.50 -1.30 -1.14 -1.06 -1.02 -1.00 y0,y1,...,yn 0. 3. 1. z0,z1,gain , , , , , vel BCs P ,P ,E ,I ,P ,P temp BCs P ,P ,E ,I ,P ,P PS1 BCs P ,P ,E ,I ,P ,P PS2 BCs P ,P ,E ,I ,P ,P PS3 BCs P ,P ,E ,I ,P ,P PS4 BCs P ,P ,E ,I ,P ,P PS5 BCs P ,P ,E ,I ,P ,P PS6 BCs P ,P ,E ,I ,P ,P PS7 BCs P ,P ,E ,I ,P ,P PS8 BCs ________________________________________________________________________ ________ logfile 103 Parameters 1 1.00000 p1 DENSITY 2 -2812.94 p2 Re_b 7 1.00000 p7 RHOCP 8 1.00000 p8 CONDUCT 11 500. p11 NSTEPS 12 0.05 p12 DT 15 250. p15 IOSTEP 21 1.000000E-06 p21 DIVERGENCE 22 1.000000E-03 p22 HELMHOLTZ 23 8 p23 NPSCAL 24 0.100000E-01 p24 TOLREL 25 0.100000E-01 p25 TOLABS 26 0.50000 p26 TARGET COURANT NUMBER 27 3.00000 p27 TORDER 28 0.000000E+00 p28 30 1.000000E+00 p30 1: user-def props via USERVP subr. 52 1. p52 if >1, history point dump frequency (in steps) 54 1.000000E+00 p54 direction of flow rate -- 1-X, 2-Y, 3-Z 55 3.000E+00 p55 volumetric flow rate 65 1.00000 p65 66 0.00000 p66 write format 67 0.00000 p67 read format 68 1. p68 iostep for avg_all (if ^= 0, else = iostep) 93 20.000 p93 if>0, # res-proj (adjust MXPREF in SIZEu accordingly) 94 5.00000 p94 if>0, start proj vel and p.s. after P094 steps 95 5.00000 p95 if>0, start proj pressure after p095 steps 99 4.000000E+00 p99 dealiasing 102 1.00000 p102 dump out divergence at each time step IFTRAN = T IFFLOW = T IFHEAT = T IFSPLIT = T IFLOMACH = F IFUSERVP = T IFUSERMV = F IFSTRS = F IFCHAR = F IFCYCLIC = F IFAXIS = F IFMVBD = F IFMELT = F IFMODEL = F IFKEPS = F IFMOAB = F IFNEKNEK = F IFSYNC = T IFVCOR = T IFINTQ = F IFCWUZ = F IFSWALL = F IFGEOM = F IFSURT = F IFWCNO = F IFTMSH for field 1 = F IFADVC for field 1 = T IFNONL for field 1 = F IFTMSH for field 2 = T IFADVC for field 2 = T IFNONL for field 2 = F IFTMSH for field 3 = T IFADVC for field 3 = T IFNONL for field 3 = F IFTMSH for field 4 = T IFADVC for field 4 = T IFNONL for field 4 = F IFTMSH for field 5 = T IFADVC for field 5 = T IFNONL for field 5 = F IFTMSH for field 6 = T IFADVC for field 6 = T IFNONL for field 6 = F IFTMSH for field 7 = T IFADVC for field 7 = T IFNONL for field 7 = F IFTMSH for field 8 = T IFADVC for field 8 = T IFNONL for field 8 = F IFTMSH for field 9 = T IFADVC for field 9 = T IFNONL for field 9 = F IFTMSH for field 10 = T IFADVC for field 10 = T IFNONL for field 10 = F Dealiasing enabled, lxd= 6 ________________________________________________________________________ ________ *.rea file ***** OUTPUT FIELD SPECIFICATION ***** 5 SPECIFICATIONS FOLLOW T COORDINATES T VELOCITY T PRESSURE T TEMPERATURE F TEMPERATURE GRADIENT 8 PASSIVE SCALARS T PS 1 T PS 2 T PS 3 T PS 4 T PS 5 T PS 6 T PS 7 T PS 8 ________________________________________________________________________ ________ *.usr file subroutine uservp (ix,iy,iz,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' if (ifield.eq.1) then udiff = param(2) ! fluid kinematic viscosity utrans = param(1) ! fluid density if (ieg .gt. nelgv) then ! properties in the solid utrans = param(1) ! solid density endif else utrans = param(7) ! fluid rhocp udiff = param(8) ! fluid conductivity if (ieg .gt. nelgv) then ! properties in the solid udiff = param(8) ! solid conductivity utrans = param(7) ! solid rhocp endif endif return end c------------------------------------------ subroutine userq (ix,iy,iz,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' qvol = 1. if (ieg.gt.nelgv) qvol = -1. return end c------------------------------------------ subroutine userbc (ix,iy,iz,iside,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' ux=0.0 uy=0.0 uz=0.0 temp=0.0 return end c------------------------------------------ subroutine useric (ix,iy,iz,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' one = 1. pi = 4.*atan(one) ux = (9.*(1.-y**8))/8.-0.6283185307179586*Cos(x)*Sin(Pi*y)*Sin(z) uy = 0.1*(1. + Cos(Pi*y))*Sin(x)*Sin(z) uz = 0.3141592653589793*Cos(z)*Sin(x)*Sin(Pi*y) temp = 0.d0 return end -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Mar 27 11:49:09 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 27 Mar 2014 17:49:09 +0100 Subject: [Nek5000-users] Problem with n2to3 Message-ID: Hi Neks, When I use n2to3 command I get this error ' input gain (0=custom,1=uniform,other=geometric spacing): At line 1583 of file n2to3.f Fortran runtime error: Bad integer for item 1 in list input -Adobe-Helvetica-Medium-R-Normal--12-120-75-75-P-67-ISO8859-1 NEKTON Version 2.6 ' could some please tell me what is the reason for this error ? Thank you, Regards, Kamal From nek5000-users at lists.mcs.anl.gov Thu Mar 27 12:19:20 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 27 Mar 2014 12:19:20 -0500 Subject: [Nek5000-users] Problem with n2to3 In-Reply-To: References: Message-ID: Hi Kamal, It looks like the variable 'mhd' was not defined as a real in a recent update. I have commited a change to n2to3.f, so if you updated and re-make the tool, it should work now. Thanks Katherine On Thu, Mar 27, 2014 at 11:49 AM, wrote: > Hi Neks, > > When I use n2to3 command I get this error > > ' input gain (0=custom,1=uniform,other=geometric spacing): > At line 1583 of file n2to3.f > Fortran runtime error: Bad integer for item 1 in list input > -Adobe-Helvetica-Medium-R-Normal--12-120-75-75-P-67-ISO8859-1 > NEKTON Version 2.6 ' > > > could some please tell me what is the reason for this error ? > > Thank you, > > > Regards, > Kamal > _______________________________________________ > 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 Mar 28 04:34:59 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 28 Mar 2014 10:34:59 +0100 Subject: [Nek5000-users] Problem with passive scalar fields - subroutine outfld Message-ID: Dear Neks, Sorry for my previous mail. I hope that the layout will be good now. _________________________________________________________________ Dear users and developers, I'm simulating a 3D conjugate heat transfer case, considering temperature as a passive scalar. When I simulate only one temperature field everything is ok. If I add more passive scalar fields (each one identical to the first one), after the subroutine prepost.f/outfld is called, then some of these fields are not good. I ran up to 8 additional fields, and results are the followings ("--" means that the field is not good): ldimt T PS1 PS2 PS3 PS4 PS5 PS6 PS7 PS8 1 ok 2 ok -- 3 ok -- ok 4 ok -- ok ok 5 ok -- -- -- ok 6 ok -- -- -- ok ok 7 ok -- -- -- ok ok ok 8 ok -- -- -- ok ok ok ok 9 ok -- -- -- ok ok ok ok ok I'm using PN-PN. Nothing changes with the number of mesh points or the number of steps or the spatial polynomial order. I saw differences with Visit, with mean temperature profiles, and with average temperature in fluid and solid domain. Differences are noticeable also in logfile. Consider the following example, in which the first field file is dumped at step 250; at step 251 the changes are noticeable in the number of iterations per field only for PS1, PS2 and PS3. :: Step 250, t= 3.19E+00, DT= 1.09E-02, C= 0.571 7.85E+01 2.6463E-01 Solving for heat 250 Hmholtz TEMP: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 1: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 2: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 3: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 4: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 5: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 6: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 7: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 8: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 3.1937E+00 2.0457E-01 Heat done Solving for fluid :: call outfld: ifpsco: T 250 3.1937E+00 Write checkpoint: 0 250 OPEN: cjh_tc0.f00001 250 3.1937E+00 done :: Write checkpoint file size = 5.6 MB avg data-throughput = 367.7MB/s io-nodes = 1 Step 251, t= 3.20E+00, DT= 1.09E-02, C= 0.566 7.88E+01 2.93E-01 Solving for heat 251 Hmholtz TEMP: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 1: 74 9.1166E-04 1.1888E+02 1.0000E-03 251 Hmholtz PS 2: 85 9.4452E-04 9.1977E+02 1.0000E-03 251 Hmholtz PS 3: 91 8.5561E-04 2.3562E+03 1.0000E-03 251 Hmholtz PS 4: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 5: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 6: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 7: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 8: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 3.2046E+00 3.6832E-01 Heat done :: Looking to the outfld subroutine I wasn't able to figure out the reason for that. Does anyone know how to solve this problem? Thanks in advance. Relevant parts of logfile, SIZE, .usr, .box and .rea files are attached at the end of the mail. Best Regards, Donato ______________________________________________________________________ SIZE file parameter (ldimt= 9) ! passive scalars + T ______________________________________________________________________ *.box file base.rea 3 spatial dimension ( < 0 --> generate .rea/.re2 pair) 10 number of fields #====== Box -10 14 -5 nelx,nely,nelz for Box 0. 6. 1. x0,x1,gain -1.00 -0.98 -0.94 -0.86 -0.70 -0.50 -0.25 -0.00 0.25 0.50 0.70 0.86 0.94 0.98 1.00 y0,y1,...,yn 0. 3. 1. z0,z1,gain P ,P ,W ,W ,P ,P vel BCs P ,P ,E ,E ,P ,P temp BCs P ,P ,E ,E ,P ,P PS1 BCs P ,P ,E ,E ,P ,P PS2 BCs P ,P ,E ,E ,P ,P PS3 BCs P ,P ,E ,E ,P ,P PS4 BCs P ,P ,E ,E ,P ,P PS5 BCs P ,P ,E ,E ,P ,P PS6 BCs P ,P ,E ,E ,P ,P PS7 BCs P ,P ,E ,E ,P ,P PS8 BCs Box_up -10 7 -5 nelx,nely,nelz for Box_up 0. 6. 1. x0,x1,gain 1.00 1.02 1.06 1.14 1.30 1.50 1.75 2.00 y0,y1,...,yn 0. 3. 1. z0,z1,gain , , , , , vel BCs P ,P ,I ,E ,P ,P temp BCs P ,P ,I ,E ,P ,P PS1 BCs P ,P ,I ,E ,P ,P PS2 BCs P ,P ,I ,E ,P ,P PS3 BCs P ,P ,I ,E ,P ,P PS4 BCs P ,P ,I ,E ,P ,P PS5 BCs P ,P ,I ,E ,P ,P PS6 BCs P ,P ,I ,E ,P ,P PS7 BCs P ,P ,I ,E ,P ,P PS8 BCs Box_down -10 7 -5 nelx,nely,nelz for Box_down 0. 6. 1. x0,x1,gain -2.00 -1.75 -1.50 -1.30 -1.14 -1.06 -1.02 -1.00 y0,y1,...,yn 0. 3. 1. z0,z1,gain , , , , , vel BCs P ,P ,E ,I ,P ,P temp BCs P ,P ,E ,I ,P ,P PS1 BCs P ,P ,E ,I ,P ,P PS2 BCs P ,P ,E ,I ,P ,P PS3 BCs P ,P ,E ,I ,P ,P PS4 BCs P ,P ,E ,I ,P ,P PS5 BCs P ,P ,E ,I ,P ,P PS6 BCs P ,P ,E ,I ,P ,P PS7 BCs P ,P ,E ,I ,P ,P PS8 BCs ______________________________________________________________________ logfile 103 Parameters 1 1.00000 p1 DENSITY 2 -2812.94 p2 Re_b 7 1.00000 p7 RHOCP 8 1.00000 p8 CONDUCT 11 500. p11 NSTEPS 12 0.05 p12 DT 15 250. p15 IOSTEP 21 1.000000E-06 p21 DIVERGENCE 22 1.000000E-03 p22 HELMHOLTZ 23 8 p23 NPSCAL 24 0.100000E-01 p24 TOLREL 25 0.100000E-01 p25 TOLABS 26 0.50000 p26 TARGET COURANT NUMBER 27 3.00000 p27 TORDER 28 0.000000E+00 p28 30 1.000000E+00 p30 1: user-def props via USERVP subr. 52 1. p52 if >1, history point dump freq.(steps) 54 1.000000E+00 p54 direction of flow rate -- 1-X, 2-Y, 3-Z 55 3.000E+00 p55 volumetric flow rate 65 1.00000 p65 66 0.00000 p66 write format 67 0.00000 p67 read format 68 1. p68 iostep for avg_all 93 20.000 p93 if>0, # res-proj 94 5.00000 p94 if>0, start proj vel and p.s. 95 5.00000 p95 if>0, start proj pressure 99 4.000000E+00 p99 dealiasing 102 1.00000 p102 dump out divergence at each time step IFTRAN = T IFFLOW = T IFHEAT = T IFSPLIT = T IFLOMACH = F IFUSERVP = T IFUSERMV = F IFSTRS = F IFCHAR = F IFCYCLIC = F IFAXIS = F IFMVBD = F IFMELT = F IFMODEL = F IFKEPS = F IFMOAB = F IFNEKNEK = F IFSYNC = T IFVCOR = T IFINTQ = F IFCWUZ = F IFSWALL = F IFGEOM = F IFSURT = F IFWCNO = F IFTMSH for field 1 = F IFADVC for field 1 = T IFNONL for field 1 = F IFTMSH for field 2 = T IFADVC for field 2 = T IFNONL for field 2 = F IFTMSH for field 3 = T IFADVC for field 3 = T IFNONL for field 3 = F IFTMSH for field 4 = T IFADVC for field 4 = T IFNONL for field 4 = F IFTMSH for field 5 = T IFADVC for field 5 = T IFNONL for field 5 = F IFTMSH for field 6 = T IFADVC for field 6 = T IFNONL for field 6 = F IFTMSH for field 7 = T IFADVC for field 7 = T IFNONL for field 7 = F IFTMSH for field 8 = T IFADVC for field 8 = T IFNONL for field 8 = F IFTMSH for field 9 = T IFADVC for field 9 = T IFNONL for field 9 = F IFTMSH for field 10 = T IFADVC for field 10 = T IFNONL for field 10 = F Dealiasing enabled, lxd= 6 ______________________________________________________________________ *.rea file ***** OUTPUT FIELD SPECIFICATION ***** 5 SPECIFICATIONS FOLLOW T COORDINATES T VELOCITY T PRESSURE T TEMPERATURE F TEMPERATURE GRADIENT 8 PASSIVE SCALARS T PS 1 T PS 2 T PS 3 T PS 4 T PS 5 T PS 6 T PS 7 T PS 8 ______________________________________________________________________ *.usr file subroutine uservp (ix,iy,iz,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' if (ifield.eq.1) then udiff = param(2) ! fluid kinematic viscosity utrans = param(1) ! fluid density if (ieg .gt. nelgv) then ! properties in the solid utrans = param(1) ! solid density endif else utrans = param(7) ! fluid rhocp udiff = param(8) ! fluid conductivity if (ieg .gt. nelgv) then ! properties in the solid udiff = param(8) ! solid conductivity utrans = param(7) ! solid rhocp endif endif return end c------------------------------------------ subroutine userq (ix,iy,iz,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' qvol = 1. if (ieg.gt.nelgv) qvol = -1. return end c------------------------------------------ subroutine userbc (ix,iy,iz,iside,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' ux=0.0 uy=0.0 uz=0.0 temp=0.0 return end c------------------------------------------ subroutine useric (ix,iy,iz,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' one = 1. pi = 4.*atan(one) ux = (9.*(1.-y**8))/8.-0.6283*Cos(x)*Sin(Pi*y)*Sin(z) uy = 0.1*(1. + Cos(Pi*y))*Sin(x)*Sin(z) uz = 0.3141*Cos(z)*Sin(x)*Sin(Pi*y) temp = 0.d0 return end -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Mar 28 05:13:14 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 28 Mar 2014 11:13:14 +0100 Subject: [Nek5000-users] Emergency Exit error Message-ID: Hello Neks, I launched a simulation on cluster, but it crashes with the error 1 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 2 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** could some please help me out with fixing this. I have attached the logfile. Thank you, Kamal -------------- next part -------------- /----------------------------------------------------------\\ | _ __ ______ __ __ ______ ____ ____ ____ | | / | / // ____// //_/ / ____/ / __ \\ / __ \\ / __ \\ | | / |/ // __/ / ,< /___ \\ / / / // / / // / / / | | / /| // /___ / /| | ____/ / / /_/ // /_/ // /_/ / | | /_/ |_//_____//_/ |_|/_____/ \\____/ \\____/ \\____/ | | | |----------------------------------------------------------| | | | NEK5000: Open Source Spectral Element Solver | | COPYRIGHT (c) 2008-2010 UCHICAGO ARGONNE, LLC | | Version: 1.0rc1 / SVN r1016 | | Web: http://nek5000.mcs.anl.gov | | | \\----------------------------------------------------------/ Number of processors: 135 REAL wdsize : 8 INTEGER wdsize : 4 Beginning session: /home/2014012/jpeixi01/neksamples/expansion/expansion.rea timer accuracy: 9.5367432E-08 sec read .rea file nelgt/nelgv/lelt: 25760 25760 200 lx1 /lx2 /lx3 : 5 5 5 mapping elements to processors RANK 0 IEG 11201 11203 11217 11254 11292 11300 11337 11338 11340 11369 11373 11401 11409 11410 11411 11412 11413 11414 11415 11416 11417 11418 11419 11420 11421 11422 11423 11424 11425 11426 11427 11428 11462 11500 11508 11509 11541 11545 11546 11548 11565 11573 11577 11581 11609 11613 11617 11618 11619 11620 11621 11623 11624 11625 11626 11629 11630 11631 11632 11633 11634 11635 11636 11669 11670 11707 11708 11717 11749 11753 11755 11756 11773 11781 11785 11789 11813 11817 11821 11825 11826 11827 11828 11829 11832 11833 11834 11837 11838 11839 11840 11841 11842 11844 11877 11878 11915 11916 11925 11957 11961 11962 11964 11981 11989 11993 12021 12025 12029 12033 12034 12035 12036 12045 12046 12047 12048 12049 12052 12085 12086 12133 12165 12169 12170 12172 12181 12189 12197 12201 12205 12229 12233 12237 12241 12242 12243 12244 12253 12254 12255 12257 12293 12294 12341 12373 12377 12378 12380 12389 12397 12405 12409 12413 12437 12441 12449 12450 12451 12452 12463 12501 12502 12529 12581 12585 12586 12587 12588 12597 12605 12613 12617 12621 12657 12659 12660 12671 12709 12710 12737 12793 12794 12795 12796 12805 12813 12821 12825 12829 element load imbalance: 1 190 191 done :: mapping elements to processors 0 objects found done :: read .rea file 2.9936 sec Reset the target Courant number to .5 setup mesh topology Right-handed check complete for 25760 elements. OK. setvert3d: 5 973641 1669161 973641 973641 call usrsetvert done :: usrsetvert gs_setup: 289964 unique labels shared pairwise times (avg, min, max): 8.12765e-05 6.88791e-05 9.08136e-05 crystal router : 0.000431342 0.000419593 0.000444794 used all_to_all method: pairwise handle bytes (avg, min, max): 287033 249460 309884 buffer bytes (avg, min, max): 82103.5 47360 106528 setupds time 5.7415E-02 seconds 0 5 973641 25760 8 max multiplicity done :: setup mesh topology call usrdat done :: usrdat generate geometry data vol_t,vol_v: 1916.3691909032111 1916.3691909032111 done :: generate geometry data call usrdat2 done :: usrdat2 regenerate geometry data 1 vol_t,vol_v: 479.09229772580278 479.09229772580278 done :: regenerate geometry data 1 verify mesh topology -1.0000000000000000 1.0000000000000000 Xrange -1.0000000000000000 1.0000000000000000 Yrange -10.000000000000000 150.00000000000000 Zrange done :: verify mesh topology 118 Parameters from file:/home/2014012/jpeixi01/neksamples/expansion/expansion.rea 1 1.00000 p001 DENSITY 2 -2000.000 p002 VISCOS 8 1.00000 p008 CONDUCT 11 1000.000 p011 NSTEPS 12 -0.100000E-02 p012 DT 15 100.0000 p015 IOSTEP 21 0.100000E-05 p021 DIVERGENCE 22 0.100000E-06 p022 HELMHOLTZ 24 0.100000E-01 p024 TOLREL 25 0.100000E-01 p025 TOLABS 26 1.00000 p026 COURANT/NTAU 27 2.00000 p027 TORDER 28 0.00000 p028 TORDER: mesh velocity (0: p28=p27) 65 1.00000 p065 #iofiles (eg, 0 or 64); <0 --> sep. dirs 66 6.00000 p066 output : <0=ascii, else binary 67 6.00000 p067 restart: <0=ascii, else binary 93 20.0000 p093 Number of previous pressure solns saved 94 3.00000 p094 start projecting velocity after p94 step 95 5.00000 p095 start projecting pressure after p95 step 99 0.00000 p099 dealiasing: <0--> off/3--> old/4--> new 103 0.100000E-01 p103 weight of stabilizing filter (.01) IFTRAN = T IFFLOW = T IFHEAT = F IFSPLIT = T IFLOMACH = F IFUSERVP = F IFUSERMV = F IFSTRS = F IFCHAR = F IFCYCLIC = F IFAXIS = F IFMVBD = F IFMELT = F IFMODEL = F IFKEPS = F IFMOAB = F IFNEKNEK = F IFSYNC = T IFVCOR = F IFINTQ = F IFCWUZ = F IFSWALL = F IFGEOM = F IFSURT = F IFWCNO = F IFTMSH for field 1 = F IFADVC for field 1 = T IFNONL for field 1 = F Dealiasing enabled, lxd= 8 Estimated eigenvalues EIGAA = 1.6450625773221985 EIGGA = 6170411880.1129513 EIGAE = 3.85531421917553079E-004 EIGAS = 1.95297242402937292E-005 EIGGE = 6170411880.1129513 EIGGS = 2.0000000000000000 verify mesh topology -1.0000000000000000 1.0000000000000000 Xrange -1.0000000000000000 1.0000000000000000 Yrange -10.000000000000000 150.00000000000000 Zrange done :: verify mesh topology E-solver strategy: 0 itr mg_nx: 1 2 4 mg_ny: 1 2 4 mg_nz: 1 2 4 call usrsetvert done :: usrsetvert gs_setup: 16049 unique labels shared pairwise times (avg, min, max): 2.49049e-05 2.21014e-05 2.77042e-05 crystal router : 5.43299e-05 5.24998e-05 5.58138e-05 all reduce : 0.000609596 0.000607705 0.000614309 used all_to_all method: pairwise handle bytes (avg, min, max): 23388.4 19756 25868 buffer bytes (avg, min, max): 7389.16 3680 10112 setupds time 1.2838E-02 seconds 1 2 27147 25760 setvert3d: 4 501087 707167 501087 501087 call usrsetvert done :: usrsetvert gs_setup: 160989 unique labels shared pairwise times (avg, min, max): 5.01423e-05 4.22955e-05 5.88179e-05 crystal router : 0.000263465 0.000255799 0.000272799 used all_to_all method: pairwise handle bytes (avg, min, max): 165421 143180 179100 buffer bytes (avg, min, max): 48269.5 27328 62976 setupds time 2.9434E-02 seconds 2 4 501087 25760 setvert3d: 3 185589 211349 185589 185589 call usrsetvert done :: usrsetvert gs_setup: 69684 unique labels shared pairwise times (avg, min, max): 3.20203e-05 2.79903e-05 3.60012e-05 crystal router : 0.000139336 0.000135708 0.00014441 all reduce : 0.00236143 0.00234699 0.0023746 used all_to_all method: pairwise handle bytes (avg, min, max): 77539.4 66612 84508 buffer bytes (avg, min, max): 23364.7 12768 30912 setupds time 4.6377E-02 seconds 3 3 185589 25760 setvert3d: 5 973641 1669161 973641 973641 call usrsetvert done :: usrsetvert gs_setup: 289964 unique labels shared pairwise times (avg, min, max): 7.52006e-05 6.23941e-05 8.84056e-05 crystal router : 0.000432244 0.000421095 0.000446296 used all_to_all method: pairwise handle bytes (avg, min, max): 287033 249460 309884 buffer bytes (avg, min, max): 82103.5 47360 106528 setupds time 4.9351E-02 seconds 4 5 973641 25760 regenerate geometry data 1 vol_t,vol_v: 479.09229772580278 479.09229772580278 NOTE: All elements deformed , param(59) ^=0 done :: regenerate geometry data 1 h1_mg_nx: 1 2 4 h1_mg_ny: 1 2 4 h1_mg_nz: 1 2 4 call usrsetvert done :: usrsetvert gs_setup: 16049 unique labels shared pairwise times (avg, min, max): 2.42703e-05 2.10047e-05 2.7895e-05 crystal router : 4.95451e-05 4.82082e-05 5.20945e-05 all reduce : 0.000627152 0.00062499 0.000628805 used all_to_all method: pairwise handle bytes (avg, min, max): 23388.4 19756 25868 buffer bytes (avg, min, max): 7389.16 3680 10112 setupds time 1.2090E-02 seconds 5 2 27147 25760 setvert3d: 4 501087 707167 501087 501087 call usrsetvert done :: usrsetvert gs_setup: 160989 unique labels shared pairwise times (avg, min, max): 4.90579e-05 4.06027e-05 5.79119e-05 crystal router : 0.000273124 0.000265503 0.000281787 used all_to_all method: pairwise handle bytes (avg, min, max): 165421 143180 179100 buffer bytes (avg, min, max): 48269.5 27328 62976 setupds time 2.9640E-02 seconds 6 4 501087 25760 setvert3d: 3 185589 211349 185589 185589 call usrsetvert done :: usrsetvert gs_setup: 69684 unique labels shared pairwise times (avg, min, max): 3.21588e-05 2.67982e-05 3.56913e-05 crystal router : 0.000132228 0.000128198 0.000136113 all reduce : 0.00225109 0.0022393 0.0022593 used all_to_all method: pairwise handle bytes (avg, min, max): 77539.4 66612 84508 buffer bytes (avg, min, max): 23364.7 12768 30912 setupds time 4.4819E-02 seconds 7 3 185589 25760 setvert3d: 5 973641 1669161 973641 973641 call usrsetvert done :: usrsetvert gs_setup: 289964 unique labels shared pairwise times (avg, min, max): 7.4347e-05 6.1202e-05 8.7595e-05 crystal router : 0.000428043 0.000416994 0.000440192 used all_to_all method: pairwise handle bytes (avg, min, max): 287033 249460 309884 buffer bytes (avg, min, max): 82103.5 47360 106528 setupds time 4.9358E-02 seconds 8 5 973641 25760 setvert3d: 5 973641 1669161 973641 973641 call usrsetvert done :: usrsetvert gs_setup: 289964 unique labels shared pairwise times (avg, min, max): 7.4325e-05 5.95808e-05 8.74996e-05 crystal router : 0.00042732 0.000415993 0.000440502 used all_to_all method: pairwise handle bytes (avg, min, max): 287033 249460 309884 buffer bytes (avg, min, max): 82103.5 47360 106528 setupds time 4.9174E-02 seconds 9 5 973641 25760 setvert3d: 7 2389917 5609917 2389917 2389917 call usrsetvert done :: usrsetvert gs_setup: 660924 unique labels shared pairwise times (avg, min, max): 0.00013483 0.000116491 0.000150299 crystal router : 0.000920938 0.000892711 0.000946999 used all_to_all method: pairwise handle bytes (avg, min, max): 631448 551156 680380 buffer bytes (avg, min, max): 176559 103840 229920 setupds time 1.1020E-01 seconds 10 7 2389917 25760 setup h1 coarse grid, nx_crs= 2 call usrsetvert done :: usrsetvert gs_setup: 16049 unique labels shared pairwise times (avg, min, max): 2.6669e-05 2.19107e-05 3.05891e-05 crystal router : 5.02219e-05 4.90904e-05 5.14984e-05 all reduce : 0.000600592 0.000597286 0.000603104 used all_to_all method: pairwise handle bytes (avg, min, max): 23388.4 19756 25868 buffer bytes (avg, min, max): 7389.16 3680 10112 done :: setup h1 coarse grid 1.3920550346374512 sec call usrdat3 done :: usrdat3 set initial conditions nekuic (1) for ifld 1 call nekuic for vel xyz min -1.0000 -1.0000 -10.000 uvwpt min -0.11602E-19 -0.11602E-19 -6.0000 -0.11602E-19 -0.11602E-19 PS min 0.0000 0.0000 0.99000E+22 xyz max 1.0000 1.0000 150.00 uvwpt max 0.80000E-19 0.80000E-19 2.0000 0.80000E-19 0.80000E-19 PS max 0.0000 0.0000 -0.99000E+22 done :: set initial conditions call userchk 0 0.0000E+00 7.8540E-01 1.0000E+00 2.0000E+00 ubar 0 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1dragx 0 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1dragy 0 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1dragz 0 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1torqx 0 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1torqy 0 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1torqz done :: userchk gridpoints unique/tot: 1669161 3220000 dofs: 1631271 1664487 Initialization successfully completed 5.1737 sec Starting time loop ... DT/DTCFL/DTFS/DTINIT 0.100E-02 0.303-316 0.494-323 0.100E-02 CFL, Ctarg! 214.75997258291918 0.50000000000000000 Step 1, t= 1.0000000E-03, DT= 1.0000000E-03, C=214.760 0.0000E+00 0.0000E+00 Solving for fluid T T T 9.99999999999999955E-008 p22 1 1 New CG1-tolerance (RINIT*epsm) = 1.08207767419886776E-006 9.99999999999999955E-007 New CG1-tolerance (RINIT*epsm) = 1.08207767419886776E-006 9.99999999999999955E-007 1 1.00000E-06 3.22803E+02 1.75695E+03 1.83728E-01 1 Divergence 2 1.00000E-06 1.57519E+02 1.75695E+03 8.96545E-02 1 Divergence 3 1.00000E-06 6.28753E+01 1.75695E+03 3.57865E-02 1 Divergence 4 1.00000E-06 3.25760E+01 1.75695E+03 1.85412E-02 1 Divergence 5 1.00000E-06 1.88748E+01 1.75695E+03 1.07429E-02 1 Divergence 6 1.00000E-06 1.16254E+01 1.75695E+03 6.61678E-03 1 Divergence 7 1.00000E-06 5.93552E+00 1.75695E+03 3.37830E-03 1 Divergence 8 1.00000E-06 3.15576E+00 1.75695E+03 1.79616E-03 1 Divergence 9 1.00000E-06 1.95012E+00 1.75695E+03 1.10994E-03 1 Divergence 10 1.00000E-06 1.12041E+00 1.75695E+03 6.37698E-04 1 Divergence 11 1.00000E-06 8.50816E-01 1.75695E+03 4.84256E-04 1 Divergence 12 1.00000E-06 7.06272E-01 1.75695E+03 4.01987E-04 1 Divergence 13 1.00000E-06 6.39486E-01 1.75695E+03 3.63974E-04 1 Divergence 14 1.00000E-06 5.67954E-01 1.75695E+03 3.23261E-04 1 Divergence 15 1.00000E-06 5.16824E-01 1.75695E+03 2.94159E-04 1 Divergence 16 1.00000E-06 4.57075E-01 1.75695E+03 2.60152E-04 1 Divergence 17 1.00000E-06 4.28050E-01 1.75695E+03 2.43632E-04 1 Divergence 18 1.00000E-06 4.03813E-01 1.75695E+03 2.29837E-04 1 Divergence 19 1.00000E-06 3.84796E-01 1.75695E+03 2.19013E-04 1 Divergence 20 1.00000E-06 3.64836E-01 1.75695E+03 2.07653E-04 1 Divergence 21 1.00000E-06 3.43890E-01 1.75695E+03 1.95731E-04 1 Divergence 22 1.00000E-06 3.30000E-01 1.75695E+03 1.87825E-04 1 Divergence 23 1.00000E-06 3.17040E-01 1.75695E+03 1.80449E-04 1 Divergence 24 1.00000E-06 3.04429E-01 1.75695E+03 1.73271E-04 1 Divergence 25 1.00000E-06 2.91701E-01 1.75695E+03 1.66026E-04 1 Divergence 26 1.00000E-06 2.80348E-01 1.75695E+03 1.59565E-04 1 Divergence 27 1.00000E-06 2.71172E-01 1.75695E+03 1.54342E-04 1 Divergence 28 1.00000E-06 2.61264E-01 1.75695E+03 1.48702E-04 1 Divergence 29 1.00000E-06 2.53621E-01 1.75695E+03 1.44353E-04 1 Divergence 30 1.00000E-06 2.46152E-01 1.75695E+03 1.40102E-04 1 Divergence 31 1.00000E-06 2.42107E-01 1.75695E+03 1.37799E-04 1 Divergence 32 1.00000E-06 2.38468E-01 1.75695E+03 1.35728E-04 1 Divergence 33 1.00000E-06 2.34723E-01 1.75695E+03 1.33597E-04 1 Divergence 34 1.00000E-06 2.31380E-01 1.75695E+03 1.31694E-04 1 Divergence 35 1.00000E-06 2.27228E-01 1.75695E+03 1.29331E-04 1 Divergence 36 1.00000E-06 2.24432E-01 1.75695E+03 1.27739E-04 1 Divergence 37 1.00000E-06 2.22329E-01 1.75695E+03 1.26542E-04 1 Divergence 38 1.00000E-06 2.20257E-01 1.75695E+03 1.25363E-04 1 Divergence 39 1.00000E-06 2.17840E-01 1.75695E+03 1.23987E-04 1 Divergence 40 1.00000E-06 2.14726E-01 1.75695E+03 1.22215E-04 1 Divergence 41 1.00000E-06 2.11561E-01 1.75695E+03 1.20414E-04 1 Divergence 42 1.00000E-06 2.08644E-01 1.75695E+03 1.18753E-04 1 Divergence 43 1.00000E-06 2.05786E-01 1.75695E+03 1.17127E-04 1 Divergence 44 1.00000E-06 2.02654E-01 1.75695E+03 1.15344E-04 1 Divergence 45 1.00000E-06 1.97932E-01 1.75695E+03 1.12656E-04 1 Divergence 46 1.00000E-06 1.92774E-01 1.75695E+03 1.09720E-04 1 Divergence 47 1.00000E-06 1.87224E-01 1.75695E+03 1.06562E-04 1 Divergence 48 1.00000E-06 1.81402E-01 1.75695E+03 1.03248E-04 1 Divergence 49 1.00000E-06 1.74384E-01 1.75695E+03 9.92536E-05 1 Divergence 50 1.00000E-06 1.68176E-01 1.75695E+03 9.57201E-05 1 Divergence 51 1.00000E-06 1.61313E-01 1.75695E+03 9.18140E-05 1 Divergence 52 1.00000E-06 1.57329E-01 1.75695E+03 8.95462E-05 1 Divergence 53 1.00000E-06 1.54409E-01 1.75695E+03 8.78842E-05 1 Divergence 54 1.00000E-06 1.52259E-01 1.75695E+03 8.66606E-05 1 Divergence 55 1.00000E-06 1.50049E-01 1.75695E+03 8.54031E-05 1 Divergence 56 1.00000E-06 1.47270E-01 1.75695E+03 8.38214E-05 1 Divergence 57 1.00000E-06 1.43941E-01 1.75695E+03 8.19263E-05 1 Divergence 58 1.00000E-06 1.40918E-01 1.75695E+03 8.02056E-05 1 Divergence 59 1.00000E-06 1.38270E-01 1.75695E+03 7.86989E-05 1 Divergence 60 1.00000E-06 1.36257E-01 1.75695E+03 7.75530E-05 1 Divergence 61 1.00000E-06 1.35121E-01 1.75695E+03 7.69061E-05 1 Divergence 62 1.00000E-06 1.33722E-01 1.75695E+03 7.61100E-05 1 Divergence 63 1.00000E-06 1.32399E-01 1.75695E+03 7.53572E-05 1 Divergence 64 1.00000E-06 1.30981E-01 1.75695E+03 7.45497E-05 1 Divergence 65 1.00000E-06 1.29756E-01 1.75695E+03 7.38528E-05 1 Divergence 66 1.00000E-06 1.28722E-01 1.75695E+03 7.32643E-05 1 Divergence 67 1.00000E-06 1.27514E-01 1.75695E+03 7.25765E-05 1 Divergence 68 1.00000E-06 1.26196E-01 1.75695E+03 7.18263E-05 1 Divergence 69 1.00000E-06 1.24699E-01 1.75695E+03 7.09747E-05 1 Divergence 70 1.00000E-06 1.23081E-01 1.75695E+03 7.00538E-05 1 Divergence 71 1.00000E-06 1.21761E-01 1.75695E+03 6.93022E-05 1 Divergence 72 1.00000E-06 1.20440E-01 1.75695E+03 6.85502E-05 1 Divergence 73 1.00000E-06 1.19462E-01 1.75695E+03 6.79937E-05 1 Divergence 74 1.00000E-06 1.18755E-01 1.75695E+03 6.75911E-05 1 Divergence 75 1.00000E-06 1.18181E-01 1.75695E+03 6.72645E-05 1 Divergence 76 1.00000E-06 1.17775E-01 1.75695E+03 6.70335E-05 1 Divergence 77 1.00000E-06 1.17406E-01 1.75695E+03 6.68235E-05 1 Divergence 78 1.00000E-06 1.17017E-01 1.75695E+03 6.66022E-05 1 Divergence 79 1.00000E-06 1.16638E-01 1.75695E+03 6.63865E-05 1 Divergence 80 1.00000E-06 1.16307E-01 1.75695E+03 6.61978E-05 1 Divergence 81 1.00000E-06 1.16055E-01 1.75695E+03 6.60545E-05 1 Divergence 82 1.00000E-06 1.15873E-01 1.75695E+03 6.59511E-05 1 Divergence 83 1.00000E-06 1.15646E-01 1.75695E+03 6.58218E-05 1 Divergence 84 1.00000E-06 1.15369E-01 1.75695E+03 6.56641E-05 1 Divergence 85 1.00000E-06 1.15022E-01 1.75695E+03 6.54664E-05 1 Divergence 86 1.00000E-06 1.14685E-01 1.75695E+03 6.52747E-05 1 Divergence 87 1.00000E-06 1.14367E-01 1.75695E+03 6.50939E-05 1 Divergence 88 1.00000E-06 1.13986E-01 1.75695E+03 6.48769E-05 1 Divergence 89 1.00000E-06 1.13554E-01 1.75695E+03 6.46313E-05 1 Divergence 90 1.00000E-06 1.12810E-01 1.75695E+03 6.42074E-05 1 Divergence 91 1.00000E-06 1.12276E-01 1.75695E+03 6.39036E-05 1 Divergence 92 1.00000E-06 1.11878E-01 1.75695E+03 6.36773E-05 1 Divergence 93 1.00000E-06 1.11541E-01 1.75695E+03 6.34855E-05 1 Divergence 94 1.00000E-06 1.11214E-01 1.75695E+03 6.32993E-05 1 Divergence 95 1.00000E-06 1.10811E-01 1.75695E+03 6.30699E-05 1 Divergence 96 1.00000E-06 1.10345E-01 1.75695E+03 6.28049E-05 1 Divergence 97 1.00000E-06 1.09846E-01 1.75695E+03 6.25205E-05 1 Divergence 98 1.00000E-06 1.09432E-01 1.75695E+03 6.22851E-05 1 Divergence 99 1.00000E-06 1.08940E-01 1.75695E+03 6.20049E-05 1 Divergence 100 1.00000E-06 1.08255E-01 1.75695E+03 6.16149E-05 1 Divergence 101 1.00000E-06 1.07131E-01 1.75695E+03 6.09755E-05 1 Divergence 102 1.00000E-06 1.06032E-01 1.75695E+03 6.03500E-05 1 Divergence 103 1.00000E-06 1.05002E-01 1.75695E+03 5.97638E-05 1 Divergence 104 1.00000E-06 1.03914E-01 1.75695E+03 5.91441E-05 1 Divergence 105 1.00000E-06 1.02512E-01 1.75695E+03 5.83463E-05 1 Divergence 106 1.00000E-06 1.00626E-01 1.75695E+03 5.72729E-05 1 Divergence 107 1.00000E-06 9.82260E-02 1.75695E+03 5.59069E-05 1 Divergence 108 1.00000E-06 9.49935E-02 1.75695E+03 5.40671E-05 1 Divergence 109 1.00000E-06 9.14758E-02 1.75695E+03 5.20650E-05 1 Divergence 110 1.00000E-06 8.75772E-02 1.75695E+03 4.98460E-05 1 Divergence 111 1.00000E-06 8.33393E-02 1.75695E+03 4.74340E-05 1 Divergence 112 1.00000E-06 7.90651E-02 1.75695E+03 4.50012E-05 1 Divergence 113 1.00000E-06 7.47196E-02 1.75695E+03 4.25279E-05 1 Divergence 114 1.00000E-06 7.03656E-02 1.75695E+03 4.00498E-05 1 Divergence 115 1.00000E-06 6.70951E-02 1.75695E+03 3.81883E-05 1 Divergence 116 1.00000E-06 6.43758E-02 1.75695E+03 3.66405E-05 1 Divergence 117 1.00000E-06 6.19525E-02 1.75695E+03 3.52613E-05 1 Divergence 118 1.00000E-06 6.00907E-02 1.75695E+03 3.42016E-05 1 Divergence 119 1.00000E-06 5.86494E-02 1.75695E+03 3.33813E-05 1 Divergence 120 1.00000E-06 5.79274E-02 1.75695E+03 3.29704E-05 1 Divergence 121 1.00000E-06 5.74350E-02 1.75695E+03 3.26901E-05 1 Divergence 122 1.00000E-06 5.67243E-02 1.75695E+03 3.22856E-05 1 Divergence 123 1.00000E-06 5.59673E-02 1.75695E+03 3.18547E-05 1 Divergence 124 1.00000E-06 5.51477E-02 1.75695E+03 3.13882E-05 1 Divergence 125 1.00000E-06 5.42863E-02 1.75695E+03 3.08979E-05 1 Divergence 126 1.00000E-06 5.34865E-02 1.75695E+03 3.04427E-05 1 Divergence 127 1.00000E-06 5.26256E-02 1.75695E+03 2.99527E-05 1 Divergence 128 1.00000E-06 5.17645E-02 1.75695E+03 2.94626E-05 1 Divergence 129 1.00000E-06 5.09915E-02 1.75695E+03 2.90227E-05 1 Divergence 130 1.00000E-06 5.02108E-02 1.75695E+03 2.85783E-05 1 Divergence 131 1.00000E-06 4.94175E-02 1.75695E+03 2.81268E-05 1 Divergence 132 1.00000E-06 4.85935E-02 1.75695E+03 2.76578E-05 1 Divergence 133 1.00000E-06 4.77428E-02 1.75695E+03 2.71736E-05 1 Divergence 134 1.00000E-06 4.69242E-02 1.75695E+03 2.67077E-05 1 Divergence 135 1.00000E-06 4.59654E-02 1.75695E+03 2.61620E-05 1 Divergence 136 1.00000E-06 4.50510E-02 1.75695E+03 2.56415E-05 1 Divergence 137 1.00000E-06 4.41995E-02 1.75695E+03 2.51569E-05 1 Divergence 138 1.00000E-06 4.30613E-02 1.75695E+03 2.45090E-05 1 Divergence 139 1.00000E-06 4.19256E-02 1.75695E+03 2.38627E-05 1 Divergence 140 1.00000E-06 4.06598E-02 1.75695E+03 2.31422E-05 1 Divergence 141 1.00000E-06 3.96012E-02 1.75695E+03 2.25397E-05 1 Divergence 142 1.00000E-06 3.83529E-02 1.75695E+03 2.18292E-05 1 Divergence 143 1.00000E-06 3.73504E-02 1.75695E+03 2.12586E-05 1 Divergence 144 1.00000E-06 3.64765E-02 1.75695E+03 2.07612E-05 1 Divergence 145 1.00000E-06 3.57075E-02 1.75695E+03 2.03235E-05 1 Divergence 146 1.00000E-06 3.50544E-02 1.75695E+03 1.99518E-05 1 Divergence 147 1.00000E-06 3.44393E-02 1.75695E+03 1.96017E-05 1 Divergence 148 1.00000E-06 3.38714E-02 1.75695E+03 1.92784E-05 1 Divergence 149 1.00000E-06 3.32778E-02 1.75695E+03 1.89406E-05 1 Divergence 150 1.00000E-06 3.27883E-02 1.75695E+03 1.86620E-05 1 Divergence 151 1.00000E-06 3.24845E-02 1.75695E+03 1.84891E-05 1 Divergence 152 1.00000E-06 3.22196E-02 1.75695E+03 1.83383E-05 1 Divergence 153 1.00000E-06 3.19538E-02 1.75695E+03 1.81870E-05 1 Divergence 154 1.00000E-06 3.16804E-02 1.75695E+03 1.80314E-05 1 Divergence 155 1.00000E-06 3.14392E-02 1.75695E+03 1.78941E-05 1 Divergence 156 1.00000E-06 3.12178E-02 1.75695E+03 1.77682E-05 1 Divergence 157 1.00000E-06 3.09832E-02 1.75695E+03 1.76346E-05 1 Divergence 158 1.00000E-06 3.07586E-02 1.75695E+03 1.75068E-05 1 Divergence 159 1.00000E-06 3.05367E-02 1.75695E+03 1.73805E-05 1 Divergence 160 1.00000E-06 3.03635E-02 1.75695E+03 1.72819E-05 1 Divergence 161 1.00000E-06 3.01779E-02 1.75695E+03 1.71762E-05 1 Divergence 162 1.00000E-06 3.00244E-02 1.75695E+03 1.70889E-05 1 Divergence 163 1.00000E-06 2.98808E-02 1.75695E+03 1.70071E-05 1 Divergence 164 1.00000E-06 2.97599E-02 1.75695E+03 1.69383E-05 1 Divergence 165 1.00000E-06 2.96434E-02 1.75695E+03 1.68720E-05 1 Divergence 166 1.00000E-06 2.95418E-02 1.75695E+03 1.68142E-05 1 Divergence 167 1.00000E-06 2.94575E-02 1.75695E+03 1.67662E-05 1 Divergence 168 1.00000E-06 2.93790E-02 1.75695E+03 1.67216E-05 1 Divergence 169 1.00000E-06 2.93042E-02 1.75695E+03 1.66790E-05 1 Divergence 170 1.00000E-06 2.92303E-02 1.75695E+03 1.66369E-05 1 Divergence 171 1.00000E-06 2.91602E-02 1.75695E+03 1.65970E-05 1 Divergence 172 1.00000E-06 2.91020E-02 1.75695E+03 1.65639E-05 1 Divergence 173 1.00000E-06 2.90389E-02 1.75695E+03 1.65280E-05 1 Divergence 174 1.00000E-06 2.89804E-02 1.75695E+03 1.64947E-05 1 Divergence 175 1.00000E-06 2.89310E-02 1.75695E+03 1.64666E-05 1 Divergence 176 1.00000E-06 2.88809E-02 1.75695E+03 1.64380E-05 1 Divergence 177 1.00000E-06 2.88287E-02 1.75695E+03 1.64083E-05 1 Divergence 178 1.00000E-06 2.87651E-02 1.75695E+03 1.63721E-05 1 Divergence 179 1.00000E-06 2.86966E-02 1.75695E+03 1.63331E-05 1 Divergence 180 1.00000E-06 2.86328E-02 1.75695E+03 1.62968E-05 1 Divergence 181 1.00000E-06 2.85850E-02 1.75695E+03 1.62696E-05 1 Divergence 182 1.00000E-06 2.85365E-02 1.75695E+03 1.62420E-05 1 Divergence 183 1.00000E-06 2.84866E-02 1.75695E+03 1.62136E-05 1 Divergence 184 1.00000E-06 2.84386E-02 1.75695E+03 1.61863E-05 1 Divergence 185 1.00000E-06 2.83915E-02 1.75695E+03 1.61595E-05 1 Divergence 186 1.00000E-06 2.83403E-02 1.75695E+03 1.61303E-05 1 Divergence 187 1.00000E-06 2.82740E-02 1.75695E+03 1.60926E-05 1 Divergence 188 1.00000E-06 2.81864E-02 1.75695E+03 1.60428E-05 1 Divergence 189 1.00000E-06 2.80892E-02 1.75695E+03 1.59874E-05 1 Divergence 190 1.00000E-06 2.79528E-02 1.75695E+03 1.59098E-05 1 Divergence 191 1.00000E-06 2.77710E-02 1.75695E+03 1.58063E-05 1 Divergence 192 1.00000E-06 2.75358E-02 1.75695E+03 1.56725E-05 1 Divergence 193 1.00000E-06 2.72558E-02 1.75695E+03 1.55131E-05 1 Divergence 194 1.00000E-06 2.69550E-02 1.75695E+03 1.53419E-05 1 Divergence 195 1.00000E-06 2.65998E-02 1.75695E+03 1.51397E-05 1 Divergence 196 1.00000E-06 2.61926E-02 1.75695E+03 1.49079E-05 1 Divergence 197 1.00000E-06 2.58063E-02 1.75695E+03 1.46881E-05 1 Divergence 198 1.00000E-06 2.53231E-02 1.75695E+03 1.44130E-05 1 Divergence 199 1.00000E-06 2.48249E-02 1.75695E+03 1.41295E-05 1 Divergence 200 1.00000E-06 2.42284E-02 1.75695E+03 1.37900E-05 1 Divergence 201 1.00000E-06 2.37263E-02 1.75695E+03 1.35042E-05 1 Divergence 202 1.00000E-06 2.31411E-02 1.75695E+03 1.31711E-05 1 Divergence 203 1.00000E-06 2.26490E-02 1.75695E+03 1.28911E-05 1 Divergence 204 1.00000E-06 2.22343E-02 1.75695E+03 1.26550E-05 1 Divergence 205 1.00000E-06 2.18832E-02 1.75695E+03 1.24552E-05 1 Divergence 206 1.00000E-06 2.15959E-02 1.75695E+03 1.22916E-05 1 Divergence 207 1.00000E-06 2.13310E-02 1.75695E+03 1.21409E-05 1 Divergence 208 1.00000E-06 2.11127E-02 1.75695E+03 1.20166E-05 1 Divergence 209 1.00000E-06 2.08897E-02 1.75695E+03 1.18897E-05 1 Divergence 210 1.00000E-06 2.07127E-02 1.75695E+03 1.17890E-05 1 Divergence 211 1.00000E-06 2.05803E-02 1.75695E+03 1.17136E-05 1 Divergence 212 1.00000E-06 2.04314E-02 1.75695E+03 1.16289E-05 1 Divergence 213 1.00000E-06 2.02645E-02 1.75695E+03 1.15339E-05 1 Divergence 214 1.00000E-06 2.00689E-02 1.75695E+03 1.14225E-05 1 Divergence 215 1.00000E-06 1.98315E-02 1.75695E+03 1.12874E-05 1 Divergence 216 1.00000E-06 1.95598E-02 1.75695E+03 1.11328E-05 1 Divergence 217 1.00000E-06 1.92118E-02 1.75695E+03 1.09347E-05 1 Divergence 218 1.00000E-06 1.88482E-02 1.75695E+03 1.07278E-05 1 Divergence 219 1.00000E-06 1.84414E-02 1.75695E+03 1.04963E-05 1 Divergence 220 1.00000E-06 1.80565E-02 1.75695E+03 1.02772E-05 1 Divergence 221 1.00000E-06 1.76029E-02 1.75695E+03 1.00190E-05 1 Divergence 222 1.00000E-06 1.71749E-02 1.75695E+03 9.77535E-06 1 Divergence 223 1.00000E-06 1.67546E-02 1.75695E+03 9.53613E-06 1 Divergence 224 1.00000E-06 1.63396E-02 1.75695E+03 9.29993E-06 1 Divergence 225 1.00000E-06 1.58904E-02 1.75695E+03 9.04428E-06 1 Divergence 226 1.00000E-06 1.54637E-02 1.75695E+03 8.80141E-06 1 Divergence 227 1.00000E-06 1.50774E-02 1.75695E+03 8.58155E-06 1 Divergence 228 1.00000E-06 1.47259E-02 1.75695E+03 8.38151E-06 1 Divergence 229 1.00000E-06 1.44075E-02 1.75695E+03 8.20026E-06 1 Divergence 230 1.00000E-06 1.41102E-02 1.75695E+03 8.03108E-06 1 Divergence 231 1.00000E-06 1.38545E-02 1.75695E+03 7.88553E-06 1 Divergence 232 1.00000E-06 1.36502E-02 1.75695E+03 7.76923E-06 1 Divergence 233 1.00000E-06 1.34692E-02 1.75695E+03 7.66623E-06 1 Divergence 234 1.00000E-06 1.33166E-02 1.75695E+03 7.57937E-06 1 Divergence 235 1.00000E-06 1.32119E-02 1.75695E+03 7.51978E-06 1 Divergence 236 1.00000E-06 1.31105E-02 1.75695E+03 7.46208E-06 1 Divergence 237 1.00000E-06 1.30219E-02 1.75695E+03 7.41162E-06 1 Divergence 238 1.00000E-06 1.29362E-02 1.75695E+03 7.36283E-06 1 Divergence 239 1.00000E-06 1.28544E-02 1.75695E+03 7.31629E-06 1 Divergence 240 1.00000E-06 1.27885E-02 1.75695E+03 7.27880E-06 1 Divergence 241 1.00000E-06 1.27382E-02 1.75695E+03 7.25015E-06 1 Divergence 242 1.00000E-06 1.26839E-02 1.75695E+03 7.21923E-06 1 Divergence 243 1.00000E-06 1.26311E-02 1.75695E+03 7.18918E-06 1 Divergence 244 1.00000E-06 1.25808E-02 1.75695E+03 7.16056E-06 1 Divergence 245 1.00000E-06 1.25278E-02 1.75695E+03 7.13041E-06 1 Divergence 246 1.00000E-06 1.24829E-02 1.75695E+03 7.10483E-06 1 Divergence 247 1.00000E-06 1.24300E-02 1.75695E+03 7.07476E-06 1 Divergence 248 1.00000E-06 1.23817E-02 1.75695E+03 7.04724E-06 1 Divergence 249 1.00000E-06 1.23360E-02 1.75695E+03 7.02122E-06 1 Divergence 250 1.00000E-06 1.22864E-02 1.75695E+03 6.99299E-06 1 Divergence 251 1.00000E-06 1.22339E-02 1.75695E+03 6.96310E-06 1 Divergence 252 1.00000E-06 1.21825E-02 1.75695E+03 6.93386E-06 1 Divergence 253 1.00000E-06 1.21395E-02 1.75695E+03 6.90940E-06 1 Divergence 254 1.00000E-06 1.21014E-02 1.75695E+03 6.88769E-06 1 Divergence 255 1.00000E-06 1.20701E-02 1.75695E+03 6.86989E-06 1 Divergence 256 1.00000E-06 1.20385E-02 1.75695E+03 6.85191E-06 1 Divergence 257 1.00000E-06 1.20159E-02 1.75695E+03 6.83902E-06 1 Divergence 258 1.00000E-06 1.19942E-02 1.75695E+03 6.82672E-06 1 Divergence 259 1.00000E-06 1.19738E-02 1.75695E+03 6.81510E-06 1 Divergence 260 1.00000E-06 1.19537E-02 1.75695E+03 6.80367E-06 1 Divergence 261 1.00000E-06 1.19375E-02 1.75695E+03 6.79443E-06 1 Divergence 262 1.00000E-06 1.19217E-02 1.75695E+03 6.78545E-06 1 Divergence 263 1.00000E-06 1.19090E-02 1.75695E+03 6.77821E-06 1 Divergence 264 1.00000E-06 1.18975E-02 1.75695E+03 6.77168E-06 1 Divergence 265 1.00000E-06 1.18871E-02 1.75695E+03 6.76573E-06 1 Divergence 266 1.00000E-06 1.18742E-02 1.75695E+03 6.75841E-06 1 Divergence 267 1.00000E-06 1.18601E-02 1.75695E+03 6.75037E-06 1 Divergence 268 1.00000E-06 1.18417E-02 1.75695E+03 6.73990E-06 1 Divergence 269 1.00000E-06 1.18150E-02 1.75695E+03 6.72470E-06 1 Divergence 270 1.00000E-06 1.17826E-02 1.75695E+03 6.70628E-06 1 Divergence 271 1.00000E-06 1.17576E-02 1.75695E+03 6.69204E-06 1 Divergence 272 1.00000E-06 1.17368E-02 1.75695E+03 6.68021E-06 1 Divergence 273 1.00000E-06 1.17172E-02 1.75695E+03 6.66903E-06 1 Divergence 274 1.00000E-06 1.17008E-02 1.75695E+03 6.65968E-06 1 Divergence 275 1.00000E-06 1.16852E-02 1.75695E+03 6.65085E-06 1 Divergence 276 1.00000E-06 1.16713E-02 1.75695E+03 6.64293E-06 1 Divergence 277 1.00000E-06 1.16537E-02 1.75695E+03 6.63288E-06 1 Divergence 278 1.00000E-06 1.16316E-02 1.75695E+03 6.62030E-06 1 Divergence 279 1.00000E-06 1.16020E-02 1.75695E+03 6.60345E-06 1 Divergence 280 1.00000E-06 1.15648E-02 1.75695E+03 6.58229E-06 1 Divergence 281 1.00000E-06 1.15223E-02 1.75695E+03 6.55808E-06 1 Divergence 282 1.00000E-06 1.14718E-02 1.75695E+03 6.52937E-06 1 Divergence 283 1.00000E-06 1.14193E-02 1.75695E+03 6.49951E-06 1 Divergence 284 1.00000E-06 1.13620E-02 1.75695E+03 6.46688E-06 1 Divergence 285 1.00000E-06 1.12771E-02 1.75695E+03 6.41855E-06 1 Divergence 286 1.00000E-06 1.11857E-02 1.75695E+03 6.36654E-06 1 Divergence 287 1.00000E-06 1.10547E-02 1.75695E+03 6.29197E-06 1 Divergence 288 1.00000E-06 1.09118E-02 1.75695E+03 6.21063E-06 1 Divergence 289 1.00000E-06 1.07453E-02 1.75695E+03 6.11586E-06 1 Divergence 290 1.00000E-06 1.05493E-02 1.75695E+03 6.00433E-06 1 Divergence 291 1.00000E-06 1.03549E-02 1.75695E+03 5.89368E-06 1 Divergence 292 1.00000E-06 1.01575E-02 1.75695E+03 5.78129E-06 1 Divergence 293 1.00000E-06 9.95734E-03 1.75695E+03 5.66738E-06 1 Divergence 294 1.00000E-06 9.71237E-03 1.75695E+03 5.52796E-06 1 Divergence 295 1.00000E-06 9.52330E-03 1.75695E+03 5.42035E-06 1 Divergence 296 1.00000E-06 9.32154E-03 1.75695E+03 5.30551E-06 1 Divergence 297 1.00000E-06 9.15521E-03 1.75695E+03 5.21084E-06 1 Divergence 298 1.00000E-06 9.01234E-03 1.75695E+03 5.12952E-06 1 Divergence 299 1.00000E-06 8.89643E-03 1.75695E+03 5.06355E-06 1 Divergence 300 1.00000E-06 8.82203E-03 1.75695E+03 5.02121E-06 1 Divergence 301 1.00000E-06 8.76187E-03 1.75695E+03 4.98696E-06 1 Divergence 302 1.00000E-06 8.68829E-03 1.75695E+03 4.94508E-06 1 Divergence 303 1.00000E-06 8.59156E-03 1.75695E+03 4.89003E-06 1 Divergence 304 1.00000E-06 8.47419E-03 1.75695E+03 4.82323E-06 1 Divergence 305 1.00000E-06 8.32673E-03 1.75695E+03 4.73930E-06 1 Divergence 306 1.00000E-06 8.18628E-03 1.75695E+03 4.65936E-06 1 Divergence 307 1.00000E-06 8.01945E-03 1.75695E+03 4.56440E-06 1 Divergence 308 1.00000E-06 7.89201E-03 1.75695E+03 4.49187E-06 1 Divergence 309 1.00000E-06 7.77079E-03 1.75695E+03 4.42287E-06 1 Divergence 310 1.00000E-06 7.63619E-03 1.75695E+03 4.34627E-06 1 Divergence 311 1.00000E-06 7.48812E-03 1.75695E+03 4.26199E-06 1 Divergence 312 1.00000E-06 7.33717E-03 1.75695E+03 4.17607E-06 1 Divergence 313 1.00000E-06 7.20455E-03 1.75695E+03 4.10059E-06 1 Divergence 314 1.00000E-06 7.06833E-03 1.75695E+03 4.02306E-06 1 Divergence 315 1.00000E-06 6.94050E-03 1.75695E+03 3.95030E-06 1 Divergence 316 1.00000E-06 6.80141E-03 1.75695E+03 3.87114E-06 1 Divergence 317 1.00000E-06 6.69247E-03 1.75695E+03 3.80913E-06 1 Divergence 318 1.00000E-06 6.58354E-03 1.75695E+03 3.74713E-06 1 Divergence 319 1.00000E-06 6.46429E-03 1.75695E+03 3.67926E-06 1 Divergence 320 1.00000E-06 6.33482E-03 1.75695E+03 3.60557E-06 1 Divergence 321 1.00000E-06 6.22624E-03 1.75695E+03 3.54377E-06 1 Divergence 322 1.00000E-06 6.12717E-03 1.75695E+03 3.48738E-06 1 Divergence 323 1.00000E-06 6.05269E-03 1.75695E+03 3.44499E-06 1 Divergence 324 1.00000E-06 5.99721E-03 1.75695E+03 3.41341E-06 1 Divergence 325 1.00000E-06 5.95430E-03 1.75695E+03 3.38899E-06 1 Divergence 326 1.00000E-06 5.91114E-03 1.75695E+03 3.36442E-06 1 Divergence 327 1.00000E-06 5.86821E-03 1.75695E+03 3.33999E-06 1 Divergence 328 1.00000E-06 5.82303E-03 1.75695E+03 3.31427E-06 1 Divergence 329 1.00000E-06 5.77782E-03 1.75695E+03 3.28854E-06 1 Divergence 330 1.00000E-06 5.72939E-03 1.75695E+03 3.26098E-06 1 Divergence 331 1.00000E-06 5.70100E-03 1.75695E+03 3.24482E-06 1 Divergence 332 1.00000E-06 5.67989E-03 1.75695E+03 3.23280E-06 1 Divergence 333 1.00000E-06 5.66022E-03 1.75695E+03 3.22161E-06 1 Divergence 334 1.00000E-06 5.64448E-03 1.75695E+03 3.21265E-06 1 Divergence 335 1.00000E-06 5.63048E-03 1.75695E+03 3.20468E-06 1 Divergence 336 1.00000E-06 5.61885E-03 1.75695E+03 3.19806E-06 1 Divergence 337 1.00000E-06 5.60472E-03 1.75695E+03 3.19002E-06 1 Divergence 338 1.00000E-06 5.58556E-03 1.75695E+03 3.17912E-06 1 Divergence 339 1.00000E-06 5.56325E-03 1.75695E+03 3.16642E-06 1 Divergence 340 1.00000E-06 5.54194E-03 1.75695E+03 3.15429E-06 1 Divergence 341 1.00000E-06 5.51877E-03 1.75695E+03 3.14110E-06 1 Divergence 342 1.00000E-06 5.49758E-03 1.75695E+03 3.12904E-06 1 Divergence 343 1.00000E-06 5.47636E-03 1.75695E+03 3.11696E-06 1 Divergence 344 1.00000E-06 5.45553E-03 1.75695E+03 3.10511E-06 1 Divergence 345 1.00000E-06 5.43374E-03 1.75695E+03 3.09270E-06 1 Divergence 346 1.00000E-06 5.41500E-03 1.75695E+03 3.08204E-06 1 Divergence 347 1.00000E-06 5.39544E-03 1.75695E+03 3.07090E-06 1 Divergence 348 1.00000E-06 5.37843E-03 1.75695E+03 3.06122E-06 1 Divergence 349 1.00000E-06 5.36231E-03 1.75695E+03 3.05205E-06 1 Divergence 350 1.00000E-06 5.34706E-03 1.75695E+03 3.04337E-06 1 Divergence 351 1.00000E-06 5.33288E-03 1.75695E+03 3.03530E-06 1 Divergence 352 1.00000E-06 5.32093E-03 1.75695E+03 3.02849E-06 1 Divergence 353 1.00000E-06 5.30999E-03 1.75695E+03 3.02227E-06 1 Divergence 354 1.00000E-06 5.29709E-03 1.75695E+03 3.01493E-06 1 Divergence 355 1.00000E-06 5.28614E-03 1.75695E+03 3.00870E-06 1 Divergence 356 1.00000E-06 5.27221E-03 1.75695E+03 3.00077E-06 1 Divergence 357 1.00000E-06 5.26091E-03 1.75695E+03 2.99433E-06 1 Divergence 358 1.00000E-06 5.24874E-03 1.75695E+03 2.98741E-06 1 Divergence 359 1.00000E-06 5.23629E-03 1.75695E+03 2.98032E-06 1 Divergence 360 1.00000E-06 5.22478E-03 1.75695E+03 2.97377E-06 1 Divergence 361 1.00000E-06 5.21529E-03 1.75695E+03 2.96837E-06 1 Divergence 362 1.00000E-06 5.20636E-03 1.75695E+03 2.96328E-06 1 Divergence 363 1.00000E-06 5.19504E-03 1.75695E+03 2.95685E-06 1 Divergence 364 1.00000E-06 5.18415E-03 1.75695E+03 2.95064E-06 1 Divergence 365 1.00000E-06 5.16950E-03 1.75695E+03 2.94231E-06 1 Divergence 366 1.00000E-06 5.15805E-03 1.75695E+03 2.93579E-06 1 Divergence 367 1.00000E-06 5.14532E-03 1.75695E+03 2.92855E-06 1 Divergence 368 1.00000E-06 5.13192E-03 1.75695E+03 2.92092E-06 1 Divergence 369 1.00000E-06 5.11827E-03 1.75695E+03 2.91315E-06 1 Divergence 370 1.00000E-06 5.10116E-03 1.75695E+03 2.90341E-06 1 Divergence 371 1.00000E-06 5.07831E-03 1.75695E+03 2.89041E-06 1 Divergence 372 1.00000E-06 5.04929E-03 1.75695E+03 2.87389E-06 1 Divergence 373 1.00000E-06 5.01936E-03 1.75695E+03 2.85685E-06 1 Divergence 374 1.00000E-06 4.98139E-03 1.75695E+03 2.83524E-06 1 Divergence 375 1.00000E-06 4.93605E-03 1.75695E+03 2.80943E-06 1 Divergence 376 1.00000E-06 4.88347E-03 1.75695E+03 2.77951E-06 1 Divergence 377 1.00000E-06 4.82636E-03 1.75695E+03 2.74700E-06 1 Divergence 378 1.00000E-06 4.76060E-03 1.75695E+03 2.70958E-06 1 Divergence 379 1.00000E-06 4.68986E-03 1.75695E+03 2.66931E-06 1 Divergence 380 1.00000E-06 4.60181E-03 1.75695E+03 2.61920E-06 1 Divergence 381 1.00000E-06 4.51462E-03 1.75695E+03 2.56957E-06 1 Divergence 382 1.00000E-06 4.43191E-03 1.75695E+03 2.52249E-06 1 Divergence 383 1.00000E-06 4.36393E-03 1.75695E+03 2.48381E-06 1 Divergence 384 1.00000E-06 4.31547E-03 1.75695E+03 2.45622E-06 1 Divergence 385 1.00000E-06 4.28315E-03 1.75695E+03 2.43782E-06 1 Divergence 386 1.00000E-06 4.25815E-03 1.75695E+03 2.42360E-06 1 Divergence 387 1.00000E-06 4.23395E-03 1.75695E+03 2.40982E-06 1 Divergence 388 1.00000E-06 4.20916E-03 1.75695E+03 2.39571E-06 1 Divergence 389 1.00000E-06 4.18351E-03 1.75695E+03 2.38111E-06 1 Divergence 390 1.00000E-06 4.15404E-03 1.75695E+03 2.36434E-06 1 Divergence 391 1.00000E-06 4.13351E-03 1.75695E+03 2.35265E-06 1 Divergence 392 1.00000E-06 4.11218E-03 1.75695E+03 2.34052E-06 1 Divergence 393 1.00000E-06 4.08949E-03 1.75695E+03 2.32760E-06 1 Divergence 394 1.00000E-06 4.06497E-03 1.75695E+03 2.31365E-06 1 Divergence 395 1.00000E-06 4.04221E-03 1.75695E+03 2.30069E-06 1 Divergence 396 1.00000E-06 4.01416E-03 1.75695E+03 2.28472E-06 1 Divergence 397 1.00000E-06 3.97866E-03 1.75695E+03 2.26452E-06 1 Divergence 398 1.00000E-06 3.93274E-03 1.75695E+03 2.23838E-06 1 Divergence 399 1.00000E-06 3.88299E-03 1.75695E+03 2.21007E-06 1 Divergence 400 1.00000E-06 3.82987E-03 1.75695E+03 2.17983E-06 1 Divergence 401 1.00000E-06 3.78744E-03 1.75695E+03 2.15568E-06 1 Divergence 402 1.00000E-06 3.75404E-03 1.75695E+03 2.13667E-06 1 Divergence 403 1.00000E-06 3.71984E-03 1.75695E+03 2.11721E-06 1 Divergence 404 1.00000E-06 3.69557E-03 1.75695E+03 2.10339E-06 1 Divergence 405 1.00000E-06 3.67474E-03 1.75695E+03 2.09154E-06 1 Divergence 406 1.00000E-06 3.65320E-03 1.75695E+03 2.07928E-06 1 Divergence 407 1.00000E-06 3.63623E-03 1.75695E+03 2.06962E-06 1 Divergence 408 1.00000E-06 3.62637E-03 1.75695E+03 2.06401E-06 1 Divergence 409 1.00000E-06 3.61506E-03 1.75695E+03 2.05757E-06 1 Divergence 410 1.00000E-06 3.60602E-03 1.75695E+03 2.05242E-06 1 Divergence 411 1.00000E-06 3.59819E-03 1.75695E+03 2.04797E-06 1 Divergence 412 1.00000E-06 3.59115E-03 1.75695E+03 2.04396E-06 1 Divergence 413 1.00000E-06 3.58501E-03 1.75695E+03 2.04047E-06 1 Divergence 414 1.00000E-06 3.57953E-03 1.75695E+03 2.03735E-06 1 Divergence 415 1.00000E-06 3.57368E-03 1.75695E+03 2.03402E-06 1 Divergence 416 1.00000E-06 3.56558E-03 1.75695E+03 2.02941E-06 1 Divergence 417 1.00000E-06 3.55934E-03 1.75695E+03 2.02586E-06 1 Divergence 418 1.00000E-06 3.55178E-03 1.75695E+03 2.02155E-06 1 Divergence 419 1.00000E-06 3.54440E-03 1.75695E+03 2.01736E-06 1 Divergence 420 1.00000E-06 3.53778E-03 1.75695E+03 2.01359E-06 1 Divergence 421 1.00000E-06 3.53256E-03 1.75695E+03 2.01062E-06 1 Divergence 422 1.00000E-06 3.52759E-03 1.75695E+03 2.00779E-06 1 Divergence 423 1.00000E-06 3.52057E-03 1.75695E+03 2.00379E-06 1 Divergence 424 1.00000E-06 3.51513E-03 1.75695E+03 2.00069E-06 1 Divergence 425 1.00000E-06 3.50793E-03 1.75695E+03 1.99660E-06 1 Divergence 426 1.00000E-06 3.50219E-03 1.75695E+03 1.99333E-06 1 Divergence 427 1.00000E-06 3.49736E-03 1.75695E+03 1.99058E-06 1 Divergence 428 1.00000E-06 3.49168E-03 1.75695E+03 1.98735E-06 1 Divergence 429 1.00000E-06 3.48574E-03 1.75695E+03 1.98397E-06 1 Divergence 430 1.00000E-06 3.47931E-03 1.75695E+03 1.98031E-06 1 Divergence 431 1.00000E-06 3.47224E-03 1.75695E+03 1.97628E-06 1 Divergence 432 1.00000E-06 3.45984E-03 1.75695E+03 1.96922E-06 1 Divergence 433 1.00000E-06 3.44724E-03 1.75695E+03 1.96206E-06 1 Divergence 434 1.00000E-06 3.43180E-03 1.75695E+03 1.95327E-06 1 Divergence 435 1.00000E-06 3.40695E-03 1.75695E+03 1.93912E-06 1 Divergence 436 1.00000E-06 3.37490E-03 1.75695E+03 1.92088E-06 1 Divergence 437 1.00000E-06 3.34238E-03 1.75695E+03 1.90237E-06 1 Divergence 438 1.00000E-06 3.29822E-03 1.75695E+03 1.87724E-06 1 Divergence 439 1.00000E-06 3.23322E-03 1.75695E+03 1.84024E-06 1 Divergence 440 1.00000E-06 3.15361E-03 1.75695E+03 1.79493E-06 1 Divergence 441 1.00000E-06 3.03777E-03 1.75695E+03 1.72900E-06 1 Divergence 442 1.00000E-06 2.89501E-03 1.75695E+03 1.64774E-06 1 Divergence 443 1.00000E-06 2.77372E-03 1.75695E+03 1.57871E-06 1 Divergence 444 1.00000E-06 2.68749E-03 1.75695E+03 1.52963E-06 1 Divergence 445 1.00000E-06 2.62286E-03 1.75695E+03 1.49284E-06 1 Divergence 446 1.00000E-06 2.58394E-03 1.75695E+03 1.47069E-06 1 Divergence 447 1.00000E-06 2.55169E-03 1.75695E+03 1.45233E-06 1 Divergence 448 1.00000E-06 2.51895E-03 1.75695E+03 1.43370E-06 1 Divergence 449 1.00000E-06 2.49348E-03 1.75695E+03 1.41921E-06 1 Divergence 450 1.00000E-06 2.47441E-03 1.75695E+03 1.40835E-06 1 Divergence 451 1.00000E-06 2.46123E-03 1.75695E+03 1.40085E-06 1 Divergence 452 1.00000E-06 2.44257E-03 1.75695E+03 1.39023E-06 1 Divergence 453 1.00000E-06 2.41728E-03 1.75695E+03 1.37583E-06 1 Divergence 454 1.00000E-06 2.38686E-03 1.75695E+03 1.35852E-06 1 Divergence 455 1.00000E-06 2.35450E-03 1.75695E+03 1.34010E-06 1 Divergence 456 1.00000E-06 2.30435E-03 1.75695E+03 1.31156E-06 1 Divergence 457 1.00000E-06 2.24507E-03 1.75695E+03 1.27782E-06 1 Divergence 458 1.00000E-06 2.16961E-03 1.75695E+03 1.23487E-06 1 Divergence 459 1.00000E-06 2.09149E-03 1.75695E+03 1.19041E-06 1 Divergence 460 1.00000E-06 2.03993E-03 1.75695E+03 1.16106E-06 1 Divergence 461 1.00000E-06 1.99946E-03 1.75695E+03 1.13803E-06 1 Divergence 462 1.00000E-06 1.96684E-03 1.75695E+03 1.11946E-06 1 Divergence 463 1.00000E-06 1.94314E-03 1.75695E+03 1.10597E-06 1 Divergence 464 1.00000E-06 1.92383E-03 1.75695E+03 1.09498E-06 1 Divergence 465 1.00000E-06 1.90324E-03 1.75695E+03 1.08326E-06 1 Divergence 466 1.00000E-06 1.88995E-03 1.75695E+03 1.07570E-06 1 Divergence 467 1.00000E-06 1.87939E-03 1.75695E+03 1.06969E-06 1 Divergence 468 1.00000E-06 1.87068E-03 1.75695E+03 1.06473E-06 1 Divergence 469 1.00000E-06 1.86355E-03 1.75695E+03 1.06067E-06 1 Divergence 470 1.00000E-06 1.85819E-03 1.75695E+03 1.05762E-06 1 Divergence 471 1.00000E-06 1.85183E-03 1.75695E+03 1.05400E-06 1 Divergence 472 1.00000E-06 1.84748E-03 1.75695E+03 1.05153E-06 1 Divergence 473 1.00000E-06 1.84399E-03 1.75695E+03 1.04954E-06 1 Divergence 474 1.00000E-06 1.83976E-03 1.75695E+03 1.04713E-06 1 Divergence 475 1.00000E-06 1.83691E-03 1.75695E+03 1.04551E-06 1 Divergence 476 1.00000E-06 1.83291E-03 1.75695E+03 1.04323E-06 1 Divergence 477 1.00000E-06 1.82897E-03 1.75695E+03 1.04099E-06 1 Divergence 478 1.00000E-06 1.82506E-03 1.75695E+03 1.03877E-06 1 Divergence 479 1.00000E-06 1.82168E-03 1.75695E+03 1.03684E-06 1 Divergence 480 1.00000E-06 1.81766E-03 1.75695E+03 1.03455E-06 1 Divergence 481 1.00000E-06 1.81462E-03 1.75695E+03 1.03282E-06 1 Divergence 482 1.00000E-06 1.81236E-03 1.75695E+03 1.03153E-06 1 Divergence 483 1.00000E-06 1.80902E-03 1.75695E+03 1.02964E-06 1 Divergence 484 1.00000E-06 1.80582E-03 1.75695E+03 1.02781E-06 1 Divergence 485 1.00000E-06 1.80274E-03 1.75695E+03 1.02606E-06 1 Divergence 486 1.00000E-06 1.80053E-03 1.75695E+03 1.02480E-06 1 Divergence 487 1.00000E-06 1.79758E-03 1.75695E+03 1.02313E-06 1 Divergence 488 1.00000E-06 1.79521E-03 1.75695E+03 1.02178E-06 1 Divergence 489 1.00000E-06 1.79236E-03 1.75695E+03 1.02015E-06 1 Divergence 490 1.00000E-06 1.78822E-03 1.75695E+03 1.01780E-06 1 Divergence 491 1.00000E-06 1.78490E-03 1.75695E+03 1.01591E-06 1 Divergence 492 1.00000E-06 1.77942E-03 1.75695E+03 1.01279E-06 1 Divergence 493 1.00000E-06 1.77175E-03 1.75695E+03 1.00842E-06 1 Divergence 494 1.00000E-06 1.76417E-03 1.75695E+03 1.00410E-06 1 Divergence 495 1.00000E-06 1.75475E-03 1.75695E+03 9.98743E-07 1 Divergence 496 1.00000E-06 1.74125E-03 1.75695E+03 9.91059E-07 1 Divergence 497 1.00000E-06 1.72832E-03 1.75695E+03 9.83699E-07 1 Divergence 498 1.00000E-06 1.71485E-03 1.75695E+03 9.76033E-07 1 Divergence 499 1.00000E-06 1.69437E-03 1.75695E+03 9.64381E-07 1 Divergence 500 1.00000E-06 1.67027E-03 1.75695E+03 9.50661E-07 1 Divergence 501 1.00000E-06 1.64364E-03 1.75695E+03 9.35507E-07 1 Divergence 502 1.00000E-06 1.60946E-03 1.75695E+03 9.16050E-07 1 Divergence 503 1.00000E-06 1.58011E-03 1.75695E+03 8.99344E-07 1 Divergence 504 1.00000E-06 1.56100E-03 1.75695E+03 8.88471E-07 1 Divergence 505 1.00000E-06 1.54409E-03 1.75695E+03 8.78844E-07 1 Divergence 506 1.00000E-06 1.53153E-03 1.75695E+03 8.71698E-07 1 Divergence 507 1.00000E-06 1.52045E-03 1.75695E+03 8.65387E-07 1 Divergence 508 1.00000E-06 1.50767E-03 1.75695E+03 8.58116E-07 1 Divergence 509 1.00000E-06 1.49572E-03 1.75695E+03 8.51315E-07 1 Divergence 510 1.00000E-06 1.48375E-03 1.75695E+03 8.44499E-07 1 Divergence 1 PRES gmres: 510 1.4837E-03 1.0000E-06 1.7570E+03 4.4790E+00 6.5689E+00 F 9.99999999999999955E-008 p22 1 1 1 1 Helmholtz VELX F: 5.0372E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 2 Helmholtz VELX F: 3.7548E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 3 Helmholtz VELX F: 2.9811E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 4 Helmholtz VELX F: 3.4244E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 5 Helmholtz VELX F: 2.5126E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 6 Helmholtz VELX F: 3.1289E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 7 Helmholtz VELX F: 2.2675E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 8 Helmholtz VELX F: 2.9705E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 9 Helmholtz VELX F: 2.0600E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 10 Helmholtz VELX F: 2.5881E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 11 Helmholtz VELX F: 1.8571E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 12 Helmholtz VELX F: 2.3955E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 13 Helmholtz VELX F: 1.6764E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 14 Helmholtz VELX F: 2.1144E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 15 Helmholtz VELX F: 1.5109E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 16 Helmholtz VELX F: 1.9033E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 17 Helmholtz VELX F: 1.3494E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 18 Helmholtz VELX F: 1.7204E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 19 Helmholtz VELX F: 1.2034E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 20 Helmholtz VELX F: 1.5221E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 21 Helmholtz VELX F: 1.0667E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 22 Helmholtz VELX F: 1.3784E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 23 Helmholtz VELX F: 9.4593E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 24 Helmholtz VELX F: 1.2150E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 25 Helmholtz VELX F: 8.3723E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 26 Helmholtz VELX F: 1.0814E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 27 Helmholtz VELX F: 7.4172E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 28 Helmholtz VELX F: 9.5709E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 29 Helmholtz VELX F: 6.5580E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 30 Helmholtz VELX F: 8.4010E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 31 Helmholtz VELX F: 5.8123E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 32 Helmholtz VELX F: 7.4063E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 33 Helmholtz VELX F: 5.1583E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 34 Helmholtz VELX F: 6.4478E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 35 Helmholtz VELX F: 4.5736E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 36 Helmholtz VELX F: 5.5897E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 37 Helmholtz VELX F: 4.0446E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 38 Helmholtz VELX F: 4.7939E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 39 Helmholtz VELX F: 3.5297E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 40 Helmholtz VELX F: 4.0633E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 41 Helmholtz VELX F: 3.0458E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 42 Helmholtz VELX F: 3.4152E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 43 Helmholtz VELX F: 2.5992E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 44 Helmholtz VELX F: 2.8566E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 45 Helmholtz VELX F: 2.2129E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 46 Helmholtz VELX F: 2.3980E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 47 Helmholtz VELX F: 1.8827E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 48 Helmholtz VELX F: 2.0141E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 49 Helmholtz VELX F: 1.6008E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 50 Helmholtz VELX F: 1.6882E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 51 Helmholtz VELX F: 1.3526E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 52 Helmholtz VELX F: 1.4089E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 53 Helmholtz VELX F: 1.1432E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 54 Helmholtz VELX F: 1.1718E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 55 Helmholtz VELX F: 9.7230E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 56 Helmholtz VELX F: 9.7481E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 57 Helmholtz VELX F: 8.3091E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 58 Helmholtz VELX F: 8.0341E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 59 Helmholtz VELX F: 7.1296E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 60 Helmholtz VELX F: 6.5755E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 61 Helmholtz VELX F: 6.0322E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 62 Helmholtz VELX F: 5.2350E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 63 Helmholtz VELX F: 4.9399E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 64 Helmholtz VELX F: 3.9440E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 65 Helmholtz VELX F: 3.8235E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 66 Helmholtz VELX F: 3.0369E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 67 Helmholtz VELX F: 2.9578E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 68 Helmholtz VELX F: 2.1610E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 69 Helmholtz VELX F: 2.3270E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 70 Helmholtz VELX F: 1.8663E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 71 Helmholtz VELX F: 1.5680E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 72 Helmholtz VELX F: 1.3638E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 73 Helmholtz VELX F: 1.4654E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 74 Helmholtz VELX F: 1.0235E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 75 Helmholtz VELX F: 9.1389E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 76 Helmholtz VELX F: 8.8702E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 77 Helmholtz VELX F: 8.6898E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 78 Helmholtz VELX F: 6.0201E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 79 Helmholtz VELX F: 5.0776E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 80 Helmholtz VELX F: 5.8083E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 81 Helmholtz VELX F: 4.2300E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 82 Helmholtz VELX F: 3.6550E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 83 Helmholtz VELX F: 2.9649E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 84 Helmholtz VELX F: 3.1284E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 85 Helmholtz VELX F: 2.9479E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 86 Helmholtz VELX F: 2.5445E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 87 Helmholtz VELX F: 2.2850E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 88 Helmholtz VELX F: 2.3946E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 89 Helmholtz VELX F: 2.5465E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 90 Helmholtz VELX F: 2.5046E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 91 Helmholtz VELX F: 2.2049E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 92 Helmholtz VELX F: 2.1031E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 93 Helmholtz VELX F: 2.3922E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 94 Helmholtz VELX F: 2.4638E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 95 Helmholtz VELX F: 2.2244E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 96 Helmholtz VELX F: 2.1260E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 97 Helmholtz VELX F: 2.0263E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 98 Helmholtz VELX F: 1.9858E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 99 Helmholtz VELX F: 1.8289E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 100 Helmholtz VELX F: 1.4050E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 100 **ERROR**: Failed in HMHOLTZ: VELX 1.4050E-03 5.0372E+00 1.0000E-07 9.99999999999999955E-008 p22 1 1 1 1 Helmholtz VELY F: 5.0372E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 2 Helmholtz VELY F: 3.7548E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 3 Helmholtz VELY F: 2.9811E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 4 Helmholtz VELY F: 3.4244E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 5 Helmholtz VELY F: 2.5126E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 6 Helmholtz VELY F: 3.1289E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 7 Helmholtz VELY F: 2.2675E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 8 Helmholtz VELY F: 2.9705E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 9 Helmholtz VELY F: 2.0600E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 10 Helmholtz VELY F: 2.5881E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 11 Helmholtz VELY F: 1.8571E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 12 Helmholtz VELY F: 2.3955E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 13 Helmholtz VELY F: 1.6764E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 14 Helmholtz VELY F: 2.1144E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 15 Helmholtz VELY F: 1.5109E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 16 Helmholtz VELY F: 1.9033E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 17 Helmholtz VELY F: 1.3494E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 18 Helmholtz VELY F: 1.7204E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 19 Helmholtz VELY F: 1.2034E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 20 Helmholtz VELY F: 1.5221E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 21 Helmholtz VELY F: 1.0667E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 22 Helmholtz VELY F: 1.3784E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 23 Helmholtz VELY F: 9.4593E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 24 Helmholtz VELY F: 1.2150E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 25 Helmholtz VELY F: 8.3723E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 26 Helmholtz VELY F: 1.0814E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 27 Helmholtz VELY F: 7.4172E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 28 Helmholtz VELY F: 9.5709E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 29 Helmholtz VELY F: 6.5580E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 30 Helmholtz VELY F: 8.4010E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 31 Helmholtz VELY F: 5.8123E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 32 Helmholtz VELY F: 7.4063E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 33 Helmholtz VELY F: 5.1583E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 34 Helmholtz VELY F: 6.4478E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 35 Helmholtz VELY F: 4.5736E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 36 Helmholtz VELY F: 5.5897E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 37 Helmholtz VELY F: 4.0446E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 38 Helmholtz VELY F: 4.7939E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 39 Helmholtz VELY F: 3.5297E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 40 Helmholtz VELY F: 4.0633E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 41 Helmholtz VELY F: 3.0458E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 42 Helmholtz VELY F: 3.4152E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 43 Helmholtz VELY F: 2.5992E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 44 Helmholtz VELY F: 2.8566E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 45 Helmholtz VELY F: 2.2129E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 46 Helmholtz VELY F: 2.3980E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 47 Helmholtz VELY F: 1.8827E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 48 Helmholtz VELY F: 2.0141E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 49 Helmholtz VELY F: 1.6008E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 50 Helmholtz VELY F: 1.6882E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 51 Helmholtz VELY F: 1.3526E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 52 Helmholtz VELY F: 1.4089E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 53 Helmholtz VELY F: 1.1432E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 54 Helmholtz VELY F: 1.1718E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 55 Helmholtz VELY F: 9.7230E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 56 Helmholtz VELY F: 9.7481E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 57 Helmholtz VELY F: 8.3091E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 58 Helmholtz VELY F: 8.0341E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 59 Helmholtz VELY F: 7.1297E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 60 Helmholtz VELY F: 6.5755E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 61 Helmholtz VELY F: 6.0322E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 62 Helmholtz VELY F: 5.2350E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 63 Helmholtz VELY F: 4.9399E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 64 Helmholtz VELY F: 3.9440E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 65 Helmholtz VELY F: 3.8235E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 66 Helmholtz VELY F: 3.0369E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 67 Helmholtz VELY F: 2.9578E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 68 Helmholtz VELY F: 2.1610E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 69 Helmholtz VELY F: 2.3270E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 70 Helmholtz VELY F: 1.8663E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 71 Helmholtz VELY F: 1.5680E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 72 Helmholtz VELY F: 1.3638E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 73 Helmholtz VELY F: 1.4654E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 74 Helmholtz VELY F: 1.0235E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 75 Helmholtz VELY F: 9.1389E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 76 Helmholtz VELY F: 8.8702E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 77 Helmholtz VELY F: 8.6898E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 78 Helmholtz VELY F: 6.0201E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 79 Helmholtz VELY F: 5.0776E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 80 Helmholtz VELY F: 5.8083E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 81 Helmholtz VELY F: 4.2300E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 82 Helmholtz VELY F: 3.6550E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 83 Helmholtz VELY F: 2.9649E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 84 Helmholtz VELY F: 3.1284E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 85 Helmholtz VELY F: 2.9480E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 86 Helmholtz VELY F: 2.5445E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 87 Helmholtz VELY F: 2.2850E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 88 Helmholtz VELY F: 2.3946E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 89 Helmholtz VELY F: 2.5465E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 90 Helmholtz VELY F: 2.5046E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 91 Helmholtz VELY F: 2.2049E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 92 Helmholtz VELY F: 2.1031E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 93 Helmholtz VELY F: 2.3923E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 94 Helmholtz VELY F: 2.4639E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 95 Helmholtz VELY F: 2.2244E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 96 Helmholtz VELY F: 2.1260E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 97 Helmholtz VELY F: 2.0263E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 98 Helmholtz VELY F: 1.9858E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 99 Helmholtz VELY F: 1.8289E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 100 Helmholtz VELY F: 1.4050E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 100 **ERROR**: Failed in HMHOLTZ: VELY 1.4050E-03 5.0372E+00 1.0000E-07 9.99999999999999955E-008 p22 1 1 1 1 Helmholtz VELZ F: 9.5292E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 2 Helmholtz VELZ F: 6.4166E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 3 Helmholtz VELZ F: 3.9902E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 4 Helmholtz VELZ F: 4.2872E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 5 Helmholtz VELZ F: 2.1071E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 6 Helmholtz VELZ F: 2.4260E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 7 Helmholtz VELZ F: 1.9539E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 8 Helmholtz VELZ F: 1.6743E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 9 Helmholtz VELZ F: 1.3834E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 10 Helmholtz VELZ F: 1.4388E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 11 Helmholtz VELZ F: 1.0349E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 12 Helmholtz VELZ F: 1.0620E+01 1.0000E-07 5.0000E-04 1.0000E+03 1 13 Helmholtz VELZ F: 9.7298E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 14 Helmholtz VELZ F: 8.8052E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 15 Helmholtz VELZ F: 6.9648E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 16 Helmholtz VELZ F: 7.8431E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 17 Helmholtz VELZ F: 6.3039E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 18 Helmholtz VELZ F: 5.6527E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 19 Helmholtz VELZ F: 5.8296E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 20 Helmholtz VELZ F: 5.0590E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 21 Helmholtz VELZ F: 4.4749E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 22 Helmholtz VELZ F: 4.0941E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 23 Helmholtz VELZ F: 4.2169E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 24 Helmholtz VELZ F: 3.2098E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 25 Helmholtz VELZ F: 3.4489E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 26 Helmholtz VELZ F: 3.0033E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 27 Helmholtz VELZ F: 2.6518E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 28 Helmholtz VELZ F: 2.4470E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 29 Helmholtz VELZ F: 2.3838E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 30 Helmholtz VELZ F: 2.0417E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 31 Helmholtz VELZ F: 1.8710E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 32 Helmholtz VELZ F: 1.7928E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 33 Helmholtz VELZ F: 1.5012E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 34 Helmholtz VELZ F: 1.3838E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 35 Helmholtz VELZ F: 1.3661E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 36 Helmholtz VELZ F: 1.1942E+00 1.0000E-07 5.0000E-04 1.0000E+03 1 37 Helmholtz VELZ F: 9.7837E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 38 Helmholtz VELZ F: 9.8364E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 39 Helmholtz VELZ F: 8.6608E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 40 Helmholtz VELZ F: 8.1446E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 41 Helmholtz VELZ F: 6.9231E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 42 Helmholtz VELZ F: 6.4787E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 43 Helmholtz VELZ F: 5.3081E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 44 Helmholtz VELZ F: 5.5599E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 45 Helmholtz VELZ F: 4.8954E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 46 Helmholtz VELZ F: 3.9745E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 47 Helmholtz VELZ F: 3.6509E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 48 Helmholtz VELZ F: 3.4559E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 49 Helmholtz VELZ F: 3.2388E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 50 Helmholtz VELZ F: 2.6830E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 51 Helmholtz VELZ F: 2.4632E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 52 Helmholtz VELZ F: 2.0555E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 53 Helmholtz VELZ F: 1.9987E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 54 Helmholtz VELZ F: 1.8225E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 55 Helmholtz VELZ F: 1.5536E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 56 Helmholtz VELZ F: 1.3305E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 57 Helmholtz VELZ F: 1.2566E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 58 Helmholtz VELZ F: 1.1498E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 59 Helmholtz VELZ F: 1.0225E-01 1.0000E-07 5.0000E-04 1.0000E+03 1 60 Helmholtz VELZ F: 9.5690E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 61 Helmholtz VELZ F: 8.7744E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 62 Helmholtz VELZ F: 8.1697E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 63 Helmholtz VELZ F: 7.2685E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 64 Helmholtz VELZ F: 6.5719E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 65 Helmholtz VELZ F: 5.3845E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 66 Helmholtz VELZ F: 4.6931E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 67 Helmholtz VELZ F: 3.8194E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 68 Helmholtz VELZ F: 3.2661E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 69 Helmholtz VELZ F: 2.7015E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 70 Helmholtz VELZ F: 2.4773E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 71 Helmholtz VELZ F: 1.9782E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 72 Helmholtz VELZ F: 1.6632E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 73 Helmholtz VELZ F: 1.5315E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 74 Helmholtz VELZ F: 1.4528E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 75 Helmholtz VELZ F: 1.2165E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 76 Helmholtz VELZ F: 1.0963E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 77 Helmholtz VELZ F: 1.0932E-02 1.0000E-07 5.0000E-04 1.0000E+03 1 78 Helmholtz VELZ F: 9.5822E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 79 Helmholtz VELZ F: 7.6653E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 80 Helmholtz VELZ F: 6.2040E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 81 Helmholtz VELZ F: 5.0502E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 82 Helmholtz VELZ F: 4.3983E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 83 Helmholtz VELZ F: 3.8328E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 84 Helmholtz VELZ F: 3.2025E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 85 Helmholtz VELZ F: 2.9368E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 86 Helmholtz VELZ F: 2.8773E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 87 Helmholtz VELZ F: 2.8604E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 88 Helmholtz VELZ F: 2.7054E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 89 Helmholtz VELZ F: 2.4420E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 90 Helmholtz VELZ F: 2.2751E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 91 Helmholtz VELZ F: 2.2877E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 92 Helmholtz VELZ F: 2.2196E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 93 Helmholtz VELZ F: 2.1330E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 94 Helmholtz VELZ F: 2.1739E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 95 Helmholtz VELZ F: 2.5034E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 96 Helmholtz VELZ F: 2.8429E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 97 Helmholtz VELZ F: 2.5881E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 98 Helmholtz VELZ F: 2.0766E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 99 Helmholtz VELZ F: 1.8595E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 100 Helmholtz VELZ F: 1.7923E-03 1.0000E-07 5.0000E-04 1.0000E+03 1 100 **ERROR**: Failed in HMHOLTZ: VELZ 1.7923E-03 9.5292E+01 1.0000E-07 L1/L2 DIV(V) : -3.8266E-06 1.6062E+00 L1/L2 QTL : 0.0000E+00 0.0000E+00 L1/L2 DIV(V)-QTL: -3.8266E-06 1.6062E+00 WARNING: DIV(V)-QTL too large! 1 1.0000E-03 7.2582E+00 Fluid done filt amp 0.0000 0.0000 0.0000 0.0000 0.0100 filt trn 1.0000 1.0000 1.0000 1.0000 0.9900 1 1.0000E-03 7.8540E-01 1.0000E+00 3.3862E+00 ubar schfile:/home/2014012/jpeixi01/neksamples/expansion/expansion.sch CFL, Ctarg! 76.427286262789934 0.50000000000000000 Step 2, t= 2.0000000E-03, DT= 1.0000000E-03, C= 76.427 7.2928E+00 7.2928E+00 Solving for fluid T T T 9.99999999999999955E-008 p22 2 1 2 PRES gmres: 510 1.6457E-03 1.0000E-06 3.5115E+03 4.4956E+00 6.5480E+00 F 9.99999999999999955E-008 p22 2 1 2 100 **ERROR**: Failed in HMHOLTZ: VELX 5.1695E-05 2.5294E+00 1.0000E-07 9.99999999999999955E-008 p22 2 1 2 100 **ERROR**: Failed in HMHOLTZ: VELY 5.1695E-05 2.5294E+00 1.0000E-07 9.99999999999999955E-008 p22 2 1 2 100 **ERROR**: Failed in HMHOLTZ: VELZ 5.2018E-05 6.2429E+01 1.0000E-07 L1/L2 DIV(V) : -3.9826E-06 1.7367E+00 L1/L2 QTL : 0.0000E+00 0.0000E+00 L1/L2 DIV(V)-QTL: -3.9826E-06 1.7367E+00 WARNING: DIV(V)-QTL too large! 2 2.0000E-03 7.2092E+00 Fluid done 2 2.0000E-03 7.8540E-01 1.0000E+00 4.0515E+00 ubar CFL, Ctarg! 88.963918084087751 0.50000000000000000 Step 3, t= 3.0000000E-03, DT= 1.0000000E-03, C= 88.964 1.4519E+01 7.2262E+00 Solving for fluid T T T 9.99999999999999955E-008 p22 3 1 3 PRES gmres: 510 6.8139E-04 1.0000E-06 2.9318E+03 4.5070E+00 6.5629E+00 F 9.99999999999999955E-008 p22 3 1 3 100 **ERROR**: Failed in HMHOLTZ: VELX 1.2233E-04 2.4292E+00 1.0000E-07 9.99999999999999955E-008 p22 3 1 3 100 **ERROR**: Failed in HMHOLTZ: VELY 1.2233E-04 2.4292E+00 1.0000E-07 9.99999999999999955E-008 p22 3 1 3 100 **ERROR**: Failed in HMHOLTZ: VELZ 2.9226E-05 5.2030E+01 1.0000E-07 L1/L2 DIV(V) : -4.0140E-06 1.4886E+00 L1/L2 QTL : 0.0000E+00 0.0000E+00 L1/L2 DIV(V)-QTL: -4.0140E-06 1.4886E+00 WARNING: DIV(V)-QTL too large! 3 3.0000E-03 7.2199E+00 Fluid done 3 3.0000E-03 7.8540E-01 1.0000E+00 3.9996E+00 ubar CFL, Ctarg! 85.072860002093805 0.50000000000000000 Step 4, t= 4.0000000E-03, DT= 1.0000000E-03, C= 85.073 2.1756E+01 7.2366E+00 Solving for fluid T T T 4 PRES gmres: 510 9.5133E-04 1.0000E-06 1.1727E+03 4.4994E+00 6.5550E+00 F 9.99999999999999955E-008 p22 4 1 4 100 **ERROR**: Failed in HMHOLTZ: VELX 2.1953E-04 5.8873E+00 1.0000E-07 9.99999999999999955E-008 p22 4 1 4 100 **ERROR**: Failed in HMHOLTZ: VELY 2.1953E-04 5.8873E+00 1.0000E-07 9.99999999999999955E-008 p22 4 1 4 100 **ERROR**: Failed in HMHOLTZ: VELZ 4.4696E-05 2.0886E+01 1.0000E-07 L1/L2 DIV(V) : -4.1535E-06 1.3150E+00 L1/L2 QTL : 0.0000E+00 0.0000E+00 L1/L2 DIV(V)-QTL: -4.1535E-06 1.3150E+00 WARNING: DIV(V)-QTL too large! 4 4.0000E-03 7.2123E+00 Fluid done 4 4.0000E-03 7.8540E-01 1.0000E+00 3.6760E+00 ubar CFL, Ctarg! 77.948750333544453 0.50000000000000000 Step 5, t= 5.0000000E-03, DT= 1.0000000E-03, C= 77.949 2.8985E+01 7.2291E+00 Solving for fluid T T T 5 PRES alph1n 6.2326E+03 6.0667E+03 1.0274E+00 1 5 PRES halpha 1 -5.8781E-01 5 PRES gmres: 510 7.7739E-04 1.0000E-06 1.2938E+00 4.5036E+00 6.5645E+00 F 9.99999999999999955E-008 p22 5 1 5 100 **ERROR**: Failed in HMHOLTZ: VELX 4.2651E-04 1.8335E+01 1.0000E-07 9.99999999999999955E-008 p22 5 1 5 100 **ERROR**: Failed in HMHOLTZ: VELY 4.2651E-04 1.8335E+01 1.0000E-07 9.99999999999999955E-008 p22 5 1 5 100 **ERROR**: Failed in HMHOLTZ: VELZ 1.6888E-04 3.9705E+00 1.0000E-07 L1/L2 DIV(V) : -4.1210E-06 1.5917E+00 L1/L2 QTL : 0.0000E+00 0.0000E+00 L1/L2 DIV(V)-QTL: -4.1210E-06 1.5917E+00 WARNING: DIV(V)-QTL too large! 5 5.0000E-03 7.2267E+00 Fluid done 5 5.0000E-03 7.8540E-01 1.0000E+00 4.1843E+00 ubar CFL, Ctarg! 96.699964741224690 0.50000000000000000 Step 6, t= 6.0000000E-03, DT= 1.0000000E-03, C= 96.700 3.6228E+01 7.2439E+00 Solving for fluid T T T 6 PRES alph1n 3.1298E+04 3.5048E+04 8.9299E-01 2 6 PRES halpha 2 -3.6429E+00 5.9155E+01 6 PRES gmres: 510 1.5538E-02 1.0000E-06 6.3365E+00 4.4854E+00 6.5337E+00 F 9.99999999999999955E-008 p22 6 1 6 100 **ERROR**: Failed in HMHOLTZ: VELX 2.1875E-03 7.4612E+01 1.0000E-07 9.99999999999999955E-008 p22 6 1 6 100 **ERROR**: Failed in HMHOLTZ: VELY 2.1875E-03 7.4612E+01 1.0000E-07 9.99999999999999955E-008 p22 6 1 6 100 **ERROR**: Failed in HMHOLTZ: VELZ 8.5719E-04 1.8110E+01 1.0000E-07 L1/L2 DIV(V) : -3.6185E-06 6.6419E+00 L1/L2 QTL : 0.0000E+00 0.0000E+00 L1/L2 DIV(V)-QTL: -3.6185E-06 6.6419E+00 WARNING: DIV(V)-QTL too large! 6 6.0000E-03 7.2321E+00 Fluid done 6 6.0000E-03 7.8540E-01 1.0000E+00 2.6491E+01 ubar CFL, Ctarg! 708.00485131999335 0.50000000000000000 Step 7, t= 7.0000000E-03, DT= 1.0000000E-03, C=708.005 4.3477E+01 7.2489E+00 Solving for fluid T T T 7 PRES alph1n 6.0102E+05 6.3402E+05 9.4795E-01 3 7 PRES halpha 3 -4.1770E+01 -2.0720E+03 2.4424E+03 7 PRES gmres: 510 6.9422E-02 1.0000E-06 1.0606E+02 4.6009E+00 6.7103E+00 F 9.99999999999999955E-008 p22 7 1 7 100 **ERROR**: Failed in HMHOLTZ: VELX 3.4033E-02 1.0637E+03 1.0000E-07 9.99999999999999955E-008 p22 7 1 7 100 **ERROR**: Failed in HMHOLTZ: VELY 3.4033E-02 1.0637E+03 1.0000E-07 9.99999999999999955E-008 p22 7 1 7 100 **ERROR**: Failed in HMHOLTZ: VELZ 7.2466E-03 2.1147E+02 1.0000E-07 L1/L2 DIV(V) : 1.3630E-07 6.9253E+01 L1/L2 QTL : 0.0000E+00 0.0000E+00 L1/L2 DIV(V)-QTL: 1.3630E-07 6.9253E+01 WARNING: DIV(V)-QTL too large! 7 7.0000E-03 7.3974E+00 Fluid done 7 7.0000E-03 7.8540E-01 1.0000E+00 2.3804E+02 ubar CFL, Ctarg! 6270.1912022422794 0.50000000000000000 Step 8, t= 8.0000000E-03, DT= 1.0000000E-03, C=******* 5.0892E+01 7.4147E+00 Solving for fluid T T T 8 PRES alph1n 1.1317E+08 1.0988E+08 1.0299E+00 4 8 PRES halpha 4 -1.2174E+04 1.6148E+06 5.1069E+05 3.2658E+05 8 PRES gmres: 510 9.3158E+00 1.0000E-06 2.1635E+04 4.5082E+00 6.5851E+00 F 9.99999999999999955E-008 p22 8 1 8 100 **ERROR**: Failed in HMHOLTZ: VELX 7.7090E+00 2.0641E+05 1.0000E-07 9.99999999999999955E-008 p22 8 1 8 100 **ERROR**: Failed in HMHOLTZ: VELY 7.7090E+00 2.0641E+05 1.0000E-07 9.99999999999999955E-008 p22 8 1 8 100 **ERROR**: Failed in HMHOLTZ: VELZ 2.4990E+00 7.2553E+04 1.0000E-07 L1/L2 DIV(V) : 4.0704E-03 2.4798E+04 L1/L2 QTL : 0.0000E+00 0.0000E+00 L1/L2 DIV(V)-QTL: 4.0704E-03 2.4798E+04 WARNING: DIV(V)-QTL too large! 8 8.0000E-03 7.2441E+00 Fluid done 8 8.0000E-03 7.8540E-01 1.0000E+00 4.9065E+04 ubar CFL, Ctarg! 1360732.6106204169 0.50000000000000000 call outfld: ifpsco: F 9 8.0000E-03 Write checkpoint: 0 9 OPEN: expansion0.f00001 1 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 2 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 3 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 4 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 5 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 6 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 7 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 9 8.0000E-03 done :: Write checkpoint file size = 88. MB avg data-throughput = 873.0MB/s io-nodes = 1 16 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 17 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 18 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 20 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 21 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 22 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 23 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 104 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 8 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 56 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 40 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 105 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 9 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 57 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 41 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 106 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 10 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 58 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 42 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 107 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 11 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 59 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 43 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 108 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 12 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 60 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 44 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 109 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 24 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 13 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 61 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 45 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 110 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 32 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 80 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 25 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 14 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 62 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 46 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 111 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 33 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 81 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 26 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 15 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 63 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 47 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 48 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 34 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 82 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 96 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 27 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 88 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 49 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 35 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 83 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 97 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 28 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 89 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 50 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 36 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 72 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 84 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 120 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 98 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 29 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 90 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 51 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 37 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 73 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 85 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 121 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 99 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 30 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 91 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 52 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 38 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 112 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 74 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 86 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 122 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 100 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 31 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 92 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 53 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 39 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 113 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 75 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 87 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 123 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 101 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 93 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 128 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 54 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 114 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 76 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 124 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 102 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 94 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 129 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 55 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 115 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 77 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 125 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 103 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 19 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 95 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 130 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 116 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 78 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 126 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 131 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 117 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 79 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 127 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 132 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 118 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 133 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 119 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 134 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 0 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 64 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 65 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 66 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 67 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 68 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 69 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 70 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** 71 Emergency exit: 9 time = 8.00000000000000017E-003 Latest solution and data are dumped for post-processing. *** STOP *** runtime statistics: total time 50.892067055566407 inv3 time 16 1.82151794433593750E-003 3.57917854338106657E-005 invc time 9144 7.39049911499023438E-003 1.45219079172416622E-004 mltd time 24 2.52811908721923828E-002 4.96760936131542138E-004 cdtp time 24 1.45916938781738281E-002 2.86718436141371812E-004 eslv time 0 0.0000000000000000 0.0000000000000000 pres time 8 52.796670913696289 1.0374243761026711 crsl time 4080 7.9781398773193359 0.15676588393645749 crsl min 6.8175837993621826 crsl max 8.9708325862884521 crsl avg 7.9962120532989500 hmhz time 27 57.585806846618652 1.1315281571044009 spro time 9 2.29263305664062500E-003 4.50489278444402358E-005 usbc time 8 1.89974308013916016E-002 3.73288645962234001E-004 usbc min 1.78277492523193359E-002 usbc max 2.98836231231689453E-002 usb avg 2.07468492013436791E-002 axhm time 6763 7.5539877414703369 0.14843153714354990 advc time 24 0.11107635498046875 2.18258682358470979E-003 vdss time 25 1.73659324645996094E-002 3.41230637097892872E-004 vdss min 1.43883228302001953E-002 vdss max 2.66075134277343750E-002 vdss avg 1.70246848353633184E-002 dsum time 10797 2.9744241237640381 5.84457322300628648E-002 dsum min 2.2821061611175537 dsum max 4.7710735797882080 dsum avg 3.1347654890131067 dadd time 0 7.8218324184417725 0.15369453180004497 ddsl time 0 0.0000000000000000 0.0000000000000000 solv time 0 0.0000000000000000 0.0000000000000000 prep time 9 0.11356639862060547 2.23151475644736957E-003 # nid tusbc tdadd tcrsl tvdss tdsum tgop qqq F 0 1.8997E-02 7.8218E+00 7.9781E+00 1.7366E-02 2.9744E+00 1.3191E-02 qqq 1 2.0400E-02 7.8423E+00 8.4285E+00 1.5377E-02 2.5655E+00 1.4935E-02 qqq 2 2.0486E-02 8.0273E+00 8.1594E+00 1.5913E-02 3.1643E+00 1.4970E-02 qqq 3 2.0749E-02 8.0644E+00 7.6260E+00 1.7381E-02 3.4140E+00 1.4896E-02 qqq 4 2.0874E-02 7.9990E+00 8.0812E+00 1.5448E-02 3.0351E+00 1.4895E-02 qqq 5 2.0221E-02 7.8341E+00 8.6697E+00 1.5337E-02 2.6427E+00 1.4889E-02 qqq 6 1.9518E-02 7.7666E+00 8.7718E+00 1.5164E-02 2.5814E+00 1.5003E-02 qqq 7 2.0141E-02 7.9038E+00 8.7661E+00 1.4774E-02 2.7152E+00 1.4891E-02 qqq 8 2.0012E-02 7.8714E+00 8.1369E+00 1.7557E-02 3.2791E+00 1.4847E-02 qqq 9 1.9076E-02 7.7617E+00 7.5790E+00 1.7205E-02 3.8093E+00 1.4864E-02 qqq 10 2.1220E-02 7.9835E+00 8.4108E+00 1.7497E-02 3.3731E+00 1.4930E-02 qqq 11 1.8966E-02 7.8573E+00 8.5894E+00 1.7756E-02 2.5395E+00 1.4891E-02 qqq 12 2.1374E-02 7.8222E+00 8.6283E+00 1.6928E-02 2.5987E+00 1.4871E-02 qqq 13 2.1943E-02 7.8546E+00 8.1899E+00 1.7641E-02 3.1074E+00 1.4981E-02 qqq 14 1.8917E-02 7.8214E+00 8.1624E+00 1.6000E-02 2.8441E+00 1.5014E-02 qqq 15 1.9269E-02 7.8290E+00 8.3168E+00 1.6467E-02 3.0379E+00 1.4860E-02 qqq 16 1.9789E-02 7.8505E+00 7.4406E+00 1.5963E-02 3.9149E+00 1.4874E-02 qqq 17 2.3311E-02 7.9581E+00 7.8078E+00 1.6060E-02 3.6824E+00 1.4981E-02 qqq 18 2.0361E-02 7.8511E+00 7.7393E+00 1.6280E-02 3.7710E+00 1.4913E-02 qqq 19 2.0780E-02 7.8792E+00 7.9376E+00 1.7926E-02 3.6091E+00 1.4794E-02 qqq 20 2.0782E-02 7.9088E+00 8.2752E+00 1.6706E-02 3.0404E+00 1.5031E-02 qqq 21 2.0214E-02 7.9419E+00 8.9708E+00 1.5912E-02 2.5444E+00 1.4861E-02 qqq 22 2.0280E-02 8.2116E+00 8.4133E+00 1.6546E-02 3.2354E+00 1.4845E-02 qqq 23 2.0966E-02 8.0479E+00 8.5005E+00 1.5932E-02 3.1699E+00 1.4866E-02 qqq 24 2.0541E-02 7.9525E+00 8.9421E+00 1.5842E-02 2.6931E+00 1.4996E-02 qqq 25 1.9933E-02 7.7591E+00 8.5776E+00 1.5635E-02 2.7921E+00 1.4881E-02 qqq 26 1.9004E-02 7.7590E+00 8.2606E+00 1.5479E-02 3.1258E+00 1.5049E-02 qqq 27 1.8812E-02 7.7192E+00 8.2581E+00 1.7215E-02 3.1503E+00 1.4877E-02 qqq 28 2.0898E-02 7.6397E+00 8.6141E+00 1.4835E-02 2.6966E+00 1.4866E-02 qqq 29 2.0755E-02 7.7618E+00 8.8365E+00 1.5570E-02 2.5246E+00 1.4892E-02 qqq 30 2.0029E-02 7.6449E+00 8.7606E+00 1.5057E-02 2.5247E+00 1.4941E-02 qqq 31 2.0037E-02 7.7788E+00 8.4921E+00 1.6268E-02 2.9307E+00 1.4912E-02 qqq 32 1.8462E-02 8.2316E+00 7.2717E+00 1.8251E-02 4.4490E+00 1.4899E-02 qqq 33 1.9141E-02 8.0813E+00 6.9938E+00 1.7984E-02 4.7711E+00 1.4889E-02 qqq 34 2.0012E-02 7.7092E+00 7.8662E+00 1.6642E-02 3.4381E+00 1.4867E-02 qqq 35 1.9677E-02 7.7678E+00 7.7212E+00 1.7287E-02 3.4074E+00 1.4859E-02 qqq 36 1.9125E-02 7.8100E+00 7.7308E+00 1.6609E-02 3.6359E+00 1.5044E-02 qqq 37 1.8250E-02 7.5938E+00 7.4622E+00 1.5498E-02 3.6643E+00 1.4992E-02 qqq 38 2.1388E-02 7.5545E+00 7.3959E+00 1.8540E-02 3.8832E+00 1.4945E-02 qqq 39 1.8054E-02 7.6125E+00 7.6501E+00 1.6223E-02 3.5646E+00 1.4865E-02 qqq 40 1.7828E-02 7.5351E+00 7.8967E+00 2.0605E-02 3.3991E+00 1.4906E-02 qqq 41 2.1564E-02 7.5247E+00 7.8982E+00 2.0040E-02 3.4334E+00 1.4852E-02 qqq 42 1.9394E-02 7.3486E+00 8.2573E+00 1.7259E-02 2.8028E+00 1.5020E-02 qqq 43 1.9411E-02 7.9668E+00 8.4359E+00 2.0545E-02 2.8271E+00 1.4862E-02 qqq 44 2.0115E-02 8.0115E+00 8.3086E+00 2.0946E-02 3.2512E+00 1.4877E-02 qqq 45 2.0313E-02 7.6972E+00 8.9501E+00 1.7477E-02 2.5241E+00 1.4859E-02 qqq 46 2.2603E-02 7.6659E+00 8.7327E+00 1.7000E-02 2.6030E+00 1.4980E-02 qqq 47 1.9747E-02 7.6953E+00 8.2351E+00 1.8872E-02 3.0575E+00 1.4940E-02 qqq 48 1.9066E-02 7.7223E+00 8.7886E+00 1.5275E-02 2.6701E+00 1.4944E-02 qqq 49 1.8854E-02 7.6590E+00 8.0037E+00 1.6690E-02 2.9829E+00 1.4903E-02 qqq 50 2.1531E-02 7.7643E+00 7.7398E+00 1.5923E-02 3.8266E+00 1.4902E-02 qqq 51 2.0745E-02 7.7767E+00 8.1410E+00 1.5507E-02 3.4132E+00 1.5073E-02 qqq 52 2.1969E-02 7.7002E+00 8.4550E+00 1.6098E-02 3.0825E+00 1.5024E-02 qqq 53 1.9438E-02 7.7663E+00 8.3370E+00 1.5162E-02 3.0130E+00 1.4897E-02 qqq 54 2.0162E-02 7.7306E+00 8.2410E+00 1.6105E-02 3.0812E+00 1.4979E-02 qqq 55 1.8429E-02 7.6938E+00 7.5033E+00 1.5954E-02 3.8848E+00 1.4934E-02 qqq 56 1.8039E-02 7.7522E+00 7.3345E+00 1.8445E-02 4.1581E+00 1.4913E-02 qqq 57 1.9627E-02 7.5588E+00 8.0168E+00 1.6092E-02 3.2113E+00 1.4856E-02 qqq 58 2.1672E-02 7.6851E+00 8.1503E+00 1.5781E-02 3.4153E+00 1.4944E-02 qqq 59 2.1654E-02 7.6977E+00 8.6098E+00 1.5740E-02 2.8923E+00 1.4865E-02 qqq 60 1.8691E-02 7.5917E+00 8.6822E+00 1.6062E-02 2.6721E+00 1.4993E-02 qqq 61 2.1041E-02 7.5110E+00 8.8209E+00 1.5756E-02 2.6777E+00 1.4864E-02 qqq 62 1.9603E-02 7.5547E+00 8.1930E+00 1.6981E-02 3.2808E+00 1.5053E-02 qqq 63 2.0205E-02 7.7503E+00 8.4462E+00 1.6958E-02 3.0425E+00 1.4885E-02 qqq 64 2.0476E-02 7.6624E+00 8.5345E+00 1.6692E-02 2.9758E+00 1.2822E-03 qqq 65 2.0051E-02 7.5855E+00 8.9003E+00 1.5428E-02 2.4402E+00 1.2832E-03 qqq 66 2.6172E-02 7.6515E+00 8.1743E+00 1.5428E-02 3.3445E+00 1.2672E-03 qqq 67 2.0842E-02 7.8336E+00 6.8920E+00 1.6771E-02 4.6180E+00 1.3671E-03 qqq 68 2.1120E-02 7.6297E+00 7.9715E+00 1.6677E-02 3.2010E+00 1.2579E-03 qqq 69 2.0030E-02 7.5890E+00 8.0874E+00 1.4388E-02 2.9543E+00 1.3192E-03 qqq 70 2.2208E-02 7.6305E+00 7.5651E+00 1.5531E-02 3.4561E+00 1.3931E-03 qqq 71 2.1794E-02 7.5165E+00 8.1758E+00 1.4644E-02 2.8983E+00 1.3182E-03 qqq 72 2.0783E-02 7.5779E+00 8.0675E+00 1.7799E-02 3.0843E+00 1.5056E-02 qqq 73 1.9461E-02 7.6131E+00 7.8338E+00 1.6891E-02 3.3133E+00 1.4898E-02 qqq 74 1.9565E-02 7.7259E+00 7.9647E+00 1.5613E-02 2.9911E+00 1.4886E-02 qqq 75 1.9466E-02 7.6783E+00 7.7528E+00 1.5789E-02 3.3698E+00 1.4937E-02 qqq 76 2.3161E-02 7.7640E+00 7.4294E+00 1.7002E-02 3.3511E+00 1.4886E-02 qqq 77 1.9887E-02 7.7196E+00 7.7493E+00 1.6977E-02 3.0271E+00 1.4904E-02 qqq 78 2.1290E-02 7.6826E+00 7.7420E+00 1.6229E-02 3.0658E+00 1.5017E-02 qqq 79 2.0850E-02 7.7486E+00 7.8388E+00 1.7046E-02 3.0475E+00 1.4970E-02 qqq 80 2.0608E-02 7.7420E+00 7.9109E+00 1.7637E-02 2.8251E+00 1.4989E-02 qqq 81 2.0664E-02 7.6271E+00 8.1421E+00 1.6241E-02 2.2821E+00 1.4903E-02 qqq 82 2.0033E-02 7.6962E+00 8.0105E+00 1.6470E-02 2.8799E+00 1.4858E-02 qqq 83 2.0547E-02 7.8031E+00 7.9124E+00 1.7149E-02 2.9909E+00 1.4873E-02 qqq 84 2.0152E-02 7.6862E+00 7.6792E+00 1.6963E-02 3.3774E+00 1.5022E-02 qqq 85 2.1266E-02 7.7688E+00 7.8481E+00 1.8033E-02 3.1863E+00 1.4942E-02 qqq 86 2.0369E-02 7.7784E+00 7.9312E+00 1.7189E-02 3.1171E+00 1.4839E-02 qqq 87 2.2780E-02 7.7185E+00 7.5206E+00 1.7629E-02 3.6009E+00 1.4862E-02 qqq 88 2.1044E-02 7.5326E+00 8.3698E+00 1.5473E-02 2.3604E+00 1.5050E-02 qqq 89 2.1616E-02 7.6924E+00 8.3328E+00 1.5203E-02 2.4429E+00 1.4963E-02 qqq 90 2.2563E-02 7.7314E+00 8.1804E+00 1.4789E-02 2.7396E+00 1.4999E-02 qqq 91 2.5641E-02 7.6883E+00 7.8726E+00 1.5150E-02 3.0380E+00 1.4883E-02 qqq 92 2.1547E-02 7.7893E+00 8.0129E+00 1.6886E-02 2.6964E+00 1.4895E-02 qqq 93 2.0789E-02 7.7934E+00 7.9006E+00 1.6682E-02 3.1198E+00 1.4881E-02 qqq 94 2.2114E-02 7.8954E+00 7.6806E+00 1.6489E-02 3.3018E+00 1.4925E-02 qqq 95 2.0794E-02 7.8142E+00 7.6028E+00 1.7606E-02 3.3341E+00 1.4892E-02 qqq 96 2.0269E-02 7.6726E+00 7.6262E+00 1.8658E-02 3.1795E+00 1.4849E-02 qqq 97 2.1118E-02 7.8171E+00 7.5935E+00 1.8120E-02 3.2789E+00 1.4895E-02 qqq 98 2.0828E-02 7.8392E+00 7.5101E+00 1.6401E-02 3.2291E+00 1.4845E-02 qqq 99 2.1659E-02 7.9864E+00 7.7380E+00 1.7657E-02 3.0656E+00 1.4928E-02 qqq 100 2.2091E-02 7.7389E+00 7.5233E+00 1.7974E-02 3.2232E+00 1.5021E-02 qqq 101 2.3491E-02 7.9507E+00 7.4743E+00 1.7877E-02 3.4506E+00 1.4873E-02 qqq 102 2.3749E-02 7.9925E+00 7.8873E+00 1.6546E-02 3.0521E+00 1.4968E-02 qqq 103 2.1884E-02 8.1659E+00 8.1194E+00 1.6705E-02 2.9187E+00 1.4843E-02 qqq 104 2.0611E-02 8.0316E+00 7.7686E+00 1.6595E-02 3.1612E+00 1.4870E-02 qqq 105 1.9465E-02 7.7357E+00 8.0229E+00 1.5415E-02 2.6442E+00 1.4919E-02 qqq 106 1.9511E-02 7.7926E+00 7.9963E+00 1.5216E-02 2.7670E+00 1.4984E-02 qqq 107 1.9622E-02 7.8786E+00 8.2851E+00 1.4623E-02 2.5351E+00 1.4867E-02 qqq 108 2.0074E-02 7.8709E+00 8.0219E+00 1.6129E-02 2.6274E+00 1.4949E-02 qqq 109 2.1814E-02 7.8046E+00 8.3771E+00 1.5132E-02 2.4123E+00 1.4867E-02 qqq 110 2.0707E-02 7.9058E+00 8.2860E+00 1.5233E-02 2.5040E+00 1.5041E-02 qqq 111 1.9787E-02 7.9039E+00 8.1343E+00 1.8379E-02 2.5311E+00 1.4904E-02 qqq 112 2.3240E-02 7.8169E+00 7.5040E+00 1.7142E-02 3.2910E+00 1.4847E-02 qqq 113 2.1423E-02 7.8949E+00 7.5822E+00 1.6338E-02 3.1599E+00 1.4842E-02 qqq 114 2.0939E-02 7.9426E+00 7.4697E+00 1.9600E-02 3.3665E+00 1.4833E-02 qqq 115 2.0994E-02 7.9038E+00 7.3860E+00 1.5774E-02 3.3223E+00 1.5022E-02 qqq 116 2.2229E-02 7.8166E+00 7.6017E+00 1.5407E-02 2.9932E+00 1.4965E-02 qqq 117 2.0817E-02 8.0007E+00 7.5879E+00 2.0319E-02 3.2945E+00 1.4897E-02 qqq 118 1.8285E-02 7.8040E+00 7.4030E+00 1.8713E-02 3.2069E+00 1.4925E-02 qqq 119 1.9691E-02 7.6913E+00 7.5841E+00 1.9394E-02 3.2251E+00 1.4873E-02 qqq 120 1.9101E-02 7.7747E+00 7.9222E+00 1.8945E-02 2.7652E+00 1.4875E-02 qqq 121 2.1934E-02 7.8507E+00 7.6348E+00 1.8246E-02 3.1548E+00 1.4961E-02 qqq 122 1.9693E-02 7.8769E+00 8.0103E+00 1.9356E-02 2.8332E+00 1.5056E-02 qqq 123 1.9724E-02 7.7084E+00 8.0795E+00 1.8274E-02 2.6286E+00 1.4919E-02 qqq 124 2.0824E-02 7.8590E+00 8.1825E+00 1.8743E-02 2.6317E+00 1.5001E-02 qqq 125 2.1933E-02 7.7790E+00 7.8159E+00 1.8163E-02 2.9538E+00 1.4880E-02 qqq 126 1.9516E-02 7.9766E+00 7.8348E+00 1.8711E-02 2.8852E+00 1.4885E-02 qqq 127 1.9078E-02 7.7655E+00 7.1991E+00 1.9511E-02 3.6562E+00 1.4900E-02 qqq 128 1.9639E-02 8.0022E+00 7.6943E+00 1.6644E-02 3.0015E+00 1.4925E-02 qqq 129 2.9496E-02 9.3101E+00 7.6622E+00 2.4991E-02 3.6170E+00 1.4986E-02 qqq 130 2.9711E-02 9.3588E+00 7.5876E+00 2.5007E-02 3.7197E+00 1.4921E-02 qqq 131 2.0338E-02 7.8577E+00 7.4443E+00 1.7324E-02 3.7143E+00 1.4978E-02 qqq 132 1.9297E-02 7.9682E+00 7.6298E+00 1.7424E-02 3.2416E+00 1.4947E-02 qqq 133 2.2209E-02 7.9547E+00 7.3134E+00 1.7201E-02 3.8384E+00 1.5018E-02 qqq 134 2.9884E-02 9.3604E+00 6.8176E+00 2.6608E-02 4.5788E+00 1.5059E-02 qqq call exitt: dying ... backtrace(): obtained 12 stack frames. /home/2014012/jpeixi01/neksamples/expansion/nek5000(print_stack_+0x1a) [0x5084aa] /home/2014012/jpeixi01/neksamples/expansion/nek5000(exitt_+0x286) [0x5d80d6] /home/2014012/jpeixi01/neksamples/expansion/nek5000(emerxit_+0x1e9) [0x4b98f9] /home/2014012/jpeixi01/neksamples/expansion/nek5000(setdt_+0x3a0) [0x4b17e0] /home/2014012/jpeixi01/neksamples/expansion/nek5000(settime_+0x49) [0x423489] /home/2014012/jpeixi01/neksamples/expansion/nek5000(nek_advance_+0x2a7) [0x41ab77] /home/2014012/jpeixi01/neksamples/expansion/nek5000(nek__multi_advance_+0x20) [0x41a880] /home/2014012/jpeixi01/neksamples/expansion/nek5000(nek_solve_+0xda) [0x41acfa] /home/2014012/jpeixi01/neksamples/expansion/nek5000(MAIN__+0x26) [0x41a846] /home/2014012/jpeixi01/neksamples/expansion/nek5000(main+0x2a) [0x63b9ea] /lib64/libc.so.6(__libc_start_main+0xfd) [0x330ee1ecdd] /home/2014012/jpeixi01/neksamples/expansion/nek5000() [0x419d09] total elapsed time : 6.34487E+01 sec total solver time incl. I/O : 5.08921E+01 sec time/timestep : 5.65467E+00 sec CPU seconds/timestep/gridpt : 4.57344E-04 sec From nek5000-users at lists.mcs.anl.gov Fri Mar 28 10:53:41 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 28 Mar 2014 16:53:41 +0100 Subject: [Nek5000-users] Emergency Exit error In-Reply-To: References: Message-ID: Hi, The information is here: > CFL, Ctarg! 214.75997258291918 0.50000000000000000 Or at each timestep, the C=*** > Step 3, t= 3.0000000E-03, DT= 1.0000000E-03, C= 88.964 1.4519E+01 7.2262E+00 That means that your CFL is too big (214.75 and 88.964), it should be below 0.5 You should decrease the timestep (parameter 12 in the REA file) -- Ismael -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 198 bytes Desc: Digital signature URL: From nek5000-users at lists.mcs.anl.gov Fri Mar 28 11:03:10 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 28 Mar 2014 17:03:10 +0100 Subject: [Nek5000-users] Emergency Exit error In-Reply-To: References: Message-ID: Thanks a lot. I figured it out and fixed it. The problem was with the meshing, with too big elements. On 28/03/2014 16:53, nek5000-users at lists.mcs.anl.gov wrote: > Hi, > The information is here: >> CFL, Ctarg! 214.75997258291918 0.50000000000000000 > Or at each timestep, the C=*** >> Step 3, t= 3.0000000E-03, DT= 1.0000000E-03, C= 88.964 1.4519E+01 7.2262E+00 > That means that your CFL is too big (214.75 and 88.964), it should be > below 0.5 > > You should decrease the timestep (parameter 12 in the REA file) > > > > _______________________________________________ > 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 Mar 28 16:22:34 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 28 Mar 2014 22:22:34 +0100 Subject: [Nek5000-users] Magnetic boundary conditions (especially (quasi) vacuum conditions) Message-ID: Hello, are there any new developments concerning the implementation? I have been running tests with the pure hydrodynamic case, for which Nek5000 seems to be well suited, and would like to go on and add an magnetic field. I think I could get quasi vacuum bc?s for the special case of an ellipsoid working, but a general implementation would be great. Cheers, Jan From nek5000-users at lists.mcs.anl.gov Sun Mar 30 23:07:18 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 31 Mar 2014 04:07:18 +0000 Subject: [Nek5000-users] Unstructured Mesh In-Reply-To: References: Message-ID: Hi JP, I've just committed two changes to the repo - one is in prenek, the other in the nek source. With these changes your import should work ok. The issue was that the VTK file had one or more elements that were numbered clockwise, rather than counterclockwise. Nek will now reverse this ordering when encountered. Please let me know if this doesn't work for you. Best regards, 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: Wednesday, March 26, 2014 1:31 AM To: nek5000-users Subject: Re: [Nek5000-users] Unstructured Mesh Dear Paul/All, Any possible updates on the error? I have tried several other geometries but unfortunately I get the same error. Thank you, JP -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Sun Mar 30 23:16:30 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 31 Mar 2014 04:16:30 +0000 Subject: [Nek5000-users] Problem with passive scalar fields - subroutine outfld In-Reply-To: References: Message-ID: Dear Donato, I've not looked closely at this yet. However, I suspect it might be related to the projection. Can you run your example with param 94 = 0 in the .rea file ? This should turn off the velocity/passive scalar projection. Another question --- in your .rea file, where it says "IFADVC" there should be two lines like: T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields) F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh) For your cases, you probably want them: T T T T T T T T T T T T IFNAV & IFADVC (convection in P.S. fields) F F T F F F F F F F F F F IFTMSH (IF mesh for this field is T mesh) This would do the following: 1) It turns on advection for all fields (fluid, temperature, and passive scalars) 2) It sets the mesh field to be Fluid It sets the velocity field to be Fluid It sets the temperature to be Temperature (i.e., fluid + solid) It sets the remaining passive scalars to be Fluid For conjugate heat transfer, normally only the temperature acts over both the fluid and solid element set, which is why I set those flags as above. Let me know if this resolves the issue for you... Regards, 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: Friday, March 28, 2014 4:34 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] Problem with passive scalar fields - subroutine outfld Dear Neks, Sorry for my previous mail. I hope that the layout will be good now. _________________________________________________________________ Dear users and developers, I?m simulating a 3D conjugate heat transfer case, considering temperature as a passive scalar. When I simulate only one temperature field everything is ok. If I add more passive scalar fields (each one identical to the first one), after the subroutine prepost.f/outfld is called, then some of these fields are not good. I ran up to 8 additional fields, and results are the followings (?--? means that the field is not good): ldimt T PS1 PS2 PS3 PS4 PS5 PS6 PS7 PS8 1 ok 2 ok -- 3 ok -- ok 4 ok -- ok ok 5 ok -- -- -- ok 6 ok -- -- -- ok ok 7 ok -- -- -- ok ok ok 8 ok -- -- -- ok ok ok ok 9 ok -- -- -- ok ok ok ok ok I?m using PN-PN. Nothing changes with the number of mesh points or the number of steps or the spatial polynomial order. I saw differences with Visit, with mean temperature profiles, and with average temperature in fluid and solid domain. Differences are noticeable also in logfile. Consider the following example, in which the first field file is dumped at step 250; at step 251 the changes are noticeable in the number of iterations per field only for PS1, PS2 and PS3. :: Step 250, t= 3.19E+00, DT= 1.09E-02, C= 0.571 7.85E+01 2.6463E-01 Solving for heat 250 Hmholtz TEMP: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 1: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 2: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 3: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 4: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 5: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 6: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 7: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 Hmholtz PS 8: 17 9.2949E-04 4.5939E-03 1.0000E-03 250 3.1937E+00 2.0457E-01 Heat done Solving for fluid :: call outfld: ifpsco: T 250 3.1937E+00 Write checkpoint: 0 250 OPEN: cjh_tc0.f00001 250 3.1937E+00 done :: Write checkpoint file size = 5.6 MB avg data-throughput = 367.7MB/s io-nodes = 1 Step 251, t= 3.20E+00, DT= 1.09E-02, C= 0.566 7.88E+01 2.93E-01 Solving for heat 251 Hmholtz TEMP: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 1: 74 9.1166E-04 1.1888E+02 1.0000E-03 251 Hmholtz PS 2: 85 9.4452E-04 9.1977E+02 1.0000E-03 251 Hmholtz PS 3: 91 8.5561E-04 2.3562E+03 1.0000E-03 251 Hmholtz PS 4: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 5: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 6: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 7: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 Hmholtz PS 8: 17 8.3175E-04 4.6145E-03 1.0000E-03 251 3.2046E+00 3.6832E-01 Heat done :: Looking to the outfld subroutine I wasn?t able to figure out the reason for that. Does anyone know how to solve this problem? Thanks in advance. Relevant parts of logfile, SIZE, .usr, .box and .rea files are attached at the end of the mail. Best Regards, Donato ______________________________________________________________________ SIZE file parameter (ldimt= 9) ! passive scalars + T ______________________________________________________________________ *.box file base.rea 3 spatial dimension ( < 0 --> generate .rea/.re2 pair) 10 number of fields #====== Box -10 14 -5 nelx,nely,nelz for Box 0. 6. 1. x0,x1,gain -1.00 -0.98 -0.94 -0.86 -0.70 -0.50 -0.25 -0.00 0.25 0.50 0.70 0.86 0.94 0.98 1.00 y0,y1,...,yn 0. 3. 1. z0,z1,gain P ,P ,W ,W ,P ,P vel BCs P ,P ,E ,E ,P ,P temp BCs P ,P ,E ,E ,P ,P PS1 BCs P ,P ,E ,E ,P ,P PS2 BCs P ,P ,E ,E ,P ,P PS3 BCs P ,P ,E ,E ,P ,P PS4 BCs P ,P ,E ,E ,P ,P PS5 BCs P ,P ,E ,E ,P ,P PS6 BCs P ,P ,E ,E ,P ,P PS7 BCs P ,P ,E ,E ,P ,P PS8 BCs Box_up -10 7 -5 nelx,nely,nelz for Box_up 0. 6. 1. x0,x1,gain 1.00 1.02 1.06 1.14 1.30 1.50 1.75 2.00 y0,y1,...,yn 0. 3. 1. z0,z1,gain , , , , , vel BCs P ,P ,I ,E ,P ,P temp BCs P ,P ,I ,E ,P ,P PS1 BCs P ,P ,I ,E ,P ,P PS2 BCs P ,P ,I ,E ,P ,P PS3 BCs P ,P ,I ,E ,P ,P PS4 BCs P ,P ,I ,E ,P ,P PS5 BCs P ,P ,I ,E ,P ,P PS6 BCs P ,P ,I ,E ,P ,P PS7 BCs P ,P ,I ,E ,P ,P PS8 BCs Box_down -10 7 -5 nelx,nely,nelz for Box_down 0. 6. 1. x0,x1,gain -2.00 -1.75 -1.50 -1.30 -1.14 -1.06 -1.02 -1.00 y0,y1,...,yn 0. 3. 1. z0,z1,gain , , , , , vel BCs P ,P ,E ,I ,P ,P temp BCs P ,P ,E ,I ,P ,P PS1 BCs P ,P ,E ,I ,P ,P PS2 BCs P ,P ,E ,I ,P ,P PS3 BCs P ,P ,E ,I ,P ,P PS4 BCs P ,P ,E ,I ,P ,P PS5 BCs P ,P ,E ,I ,P ,P PS6 BCs P ,P ,E ,I ,P ,P PS7 BCs P ,P ,E ,I ,P ,P PS8 BCs ______________________________________________________________________ logfile 103 Parameters 1 1.00000 p1 DENSITY 2 -2812.94 p2 Re_b 7 1.00000 p7 RHOCP 8 1.00000 p8 CONDUCT 11 500. p11 NSTEPS 12 0.05 p12 DT 15 250. p15 IOSTEP 21 1.000000E-06 p21 DIVERGENCE 22 1.000000E-03 p22 HELMHOLTZ 23 8 p23 NPSCAL 24 0.100000E-01 p24 TOLREL 25 0.100000E-01 p25 TOLABS 26 0.50000 p26 TARGET COURANT NUMBER 27 3.00000 p27 TORDER 28 0.000000E+00 p28 30 1.000000E+00 p30 1: user-def props via USERVP subr. 52 1. p52 if >1, history point dump freq.(steps) 54 1.000000E+00 p54 direction of flow rate -- 1-X, 2-Y, 3-Z 55 3.000E+00 p55 volumetric flow rate 65 1.00000 p65 66 0.00000 p66 write format 67 0.00000 p67 read format 68 1. p68 iostep for avg_all 93 20.000 p93 if>0, # res-proj 94 5.00000 p94 if>0, start proj vel and p.s. 95 5.00000 p95 if>0, start proj pressure 99 4.000000E+00 p99 dealiasing 102 1.00000 p102 dump out divergence at each time step IFTRAN = T IFFLOW = T IFHEAT = T IFSPLIT = T IFLOMACH = F IFUSERVP = T IFUSERMV = F IFSTRS = F IFCHAR = F IFCYCLIC = F IFAXIS = F IFMVBD = F IFMELT = F IFMODEL = F IFKEPS = F IFMOAB = F IFNEKNEK = F IFSYNC = T IFVCOR = T IFINTQ = F IFCWUZ = F IFSWALL = F IFGEOM = F IFSURT = F IFWCNO = F IFTMSH for field 1 = F IFADVC for field 1 = T IFNONL for field 1 = F IFTMSH for field 2 = T IFADVC for field 2 = T IFNONL for field 2 = F IFTMSH for field 3 = T IFADVC for field 3 = T IFNONL for field 3 = F IFTMSH for field 4 = T IFADVC for field 4 = T IFNONL for field 4 = F IFTMSH for field 5 = T IFADVC for field 5 = T IFNONL for field 5 = F IFTMSH for field 6 = T IFADVC for field 6 = T IFNONL for field 6 = F IFTMSH for field 7 = T IFADVC for field 7 = T IFNONL for field 7 = F IFTMSH for field 8 = T IFADVC for field 8 = T IFNONL for field 8 = F IFTMSH for field 9 = T IFADVC for field 9 = T IFNONL for field 9 = F IFTMSH for field 10 = T IFADVC for field 10 = T IFNONL for field 10 = F Dealiasing enabled, lxd= 6 ______________________________________________________________________ *.rea file ***** OUTPUT FIELD SPECIFICATION ***** 5 SPECIFICATIONS FOLLOW T COORDINATES T VELOCITY T PRESSURE T TEMPERATURE F TEMPERATURE GRADIENT 8 PASSIVE SCALARS T PS 1 T PS 2 T PS 3 T PS 4 T PS 5 T PS 6 T PS 7 T PS 8 ______________________________________________________________________ *.usr file subroutine uservp (ix,iy,iz,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' if (ifield.eq.1) then udiff = param(2) ! fluid kinematic viscosity utrans = param(1) ! fluid density if (ieg .gt. nelgv) then ! properties in the solid utrans = param(1) ! solid density endif else utrans = param(7) ! fluid rhocp udiff = param(8) ! fluid conductivity if (ieg .gt. nelgv) then ! properties in the solid udiff = param(8) ! solid conductivity utrans = param(7) ! solid rhocp endif endif return end c------------------------------------------ subroutine userq (ix,iy,iz,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' qvol = 1. if (ieg.gt.nelgv) qvol = -1. return end c------------------------------------------ subroutine userbc (ix,iy,iz,iside,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' ux=0.0 uy=0.0 uz=0.0 temp=0.0 return end c------------------------------------------ subroutine useric (ix,iy,iz,ieg) include 'SIZE' include 'TOTAL' include 'NEKUSE' one = 1. pi = 4.*atan(one) ux = (9.*(1.-y**8))/8.-0.6283*Cos(x)*Sin(Pi*y)*Sin(z) uy = 0.1*(1. + Cos(Pi*y))*Sin(x)*Sin(z) uz = 0.3141*Cos(z)*Sin(x)*Sin(Pi*y) temp = 0.d0 return end -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Sun Mar 30 23:31:49 2014 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sun, 30 Mar 2014 23:31:49 -0500 (CDT) Subject: [Nek5000-users] Users Meeting, Thessaloniki, Greece, Aug. 21-22 2014 Message-ID: We are pleased to invite you to attend the third Nek5000 User and Development Meeting, to be hosted by the University of Western Macedonia in Thessaloniki, Greece, August 21-22, 2014. Description: The objective of the meeting is to bring developers and users of Nek5000 together to promote collaborative activities, exchange information, and share experiences in using Nek5000 in areas of common interest. The schedule is under development, but we are planning to have technical sessions on Thursday, August. 21 and Friday, August 22, including hands-on sessions in the use of Nek5000 and VisIt. One objective of the meeting is to assist users in defining their simulation problems and to address technical issues that they may face. Attendees are invited to contribute short (15 min.) presentations on their applications of Nek5000. Contributions should consist of an abstract (1-2 pages) and presentation slides. For more information and registration please visit http://Nek5000.mcs.anl.gov/index.php/Usermeeting2014 (Note: Webpage will be going up shortly.)