[Nek5000-users] Solution for Re: genbox error

Mon May 3 15:39:01 CDT 2010

Thanks for bringing this up. That's a typical situation. We know Nek and the Tools are not straightforward to use and there are many undocumented tricks and options.
It simply takes some time to become familiar with it. But I guess that's true for all software.

Sorry for all the trouble.

Stefan

On May 3, 2010, at 10:33 PM, nek5000-users at lists.mcs.anl.gov wrote:

> Hello Aleks,
> 
> Yes, I see this would work also.  I was (stubbornly) following the
> examples; everything on one line for each coordinate direction.
> 
> Thanks,
> Frank
> 
> 
> On Mon, 2010-05-03 at 14:59 -0500, nek5000-users at lists.mcs.anl.gov
> wrote:
>> Hello,
>> 
>> I disagree that genbox reuires a fix for meshes with large number of elements in one direction -- below is a part of .box file with nelx=27 written with each element coordinate per new line:
>> 
>> Best,
>> Aleks
>> 
>> #
>> Box
>> 27 -106 -54                                 (nelx,nely,nelz  for Box)
>> 0.5378790
>> 0.5463419
>> 0.5556915
>> 0.5663933
>> 0.5789130
>> 0.5936488
>> 0.6109992
>> 0.6329137
>> 0.6599956
>> 0.6928485
>> 0.7320758
>> 0.7782808
>> 0.8320671
>> 0.8933098
>> 0.9533322
>> 1.0101231
>> 1.0637385
>> 1.1142344
>> 1.1616669
>> 1.2060918
>> 1.2475652
>> 1.2890386
>> 1.3305120
>> 1.3719854
>> 1.4134588
>> 1.4549322
>> 1.4964056
>> 1.5378790                                    (x0,x1,gain or x
>> 
>> 
>> 
>> 
>> ----- Original Message -----
>> From: nek5000-users at lists.mcs.anl.gov
>> To: nek5000-users at lists.mcs.anl.gov
>> Sent: Monday, May 3, 2010 1:47:05 PM GMT -06:00 US/Canada Central
>> Subject: Re: [Nek5000-users] Solution for Re:  genbox error
>> 
>> Please don't start distributing some code because it will end up in a mess! That's why we have a central and public available SVN server.
>> 
>> If you think you know how to improve something please contact one of the developers. We will review the code and may add it to the repo.
>> 
>> Stefan
>> 
>> On May 3, 2010, at 8:36 PM, nek5000-users at lists.mcs.anl.gov wrote:
>> 
>>> Hello all,
>>> 
>>> Genbox has an issue in that the number of characters in a line for its
>>> input file was <= 80.  This makes it impossible to have a large number
>>> of elements in one direction if one is specifying the element
>>> locations.  
>>> 
>>> A solution is to extend the line length to 9999999.  A version with this
>>> change is at.
>>> http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/temp/genbox.f
>>> 
>>> The diff between old & new is:
>>> 
>>> diff */genbox.f ../*1/g*/genbox.f
>>> 1157,1158c1157,1158
>>> <       character*9999999 string
>>> <       character*1  string1(9999999)
>>> ---
>>>>     character*80 string
>>>>     character*1  string1(80)
>>> 1167c1167
>>> <          call blank(string,9999999)
>>> ---
>>>>        call blank(string,80)
>>> 1174c1174
>>> <               len = ltrunc(string,9999999)
>>> ---
>>>>             len = ltrunc(string,80)
>>> 1184,1185c1184,1185
>>> <    80 format(a9999999)
>>> <    81 format(9999999a1)
>>> ---
>>>>  80 format(a80)
>>>>  81 format(80a1)
>>> 
>>> 
>>> 
>>> 
>>> 
>>> 
>>> On Mon, 2010-05-03 at 14:13 +0200, nek5000-users at lists.mcs.anl.gov
>>> wrote:
>>>> Hello all,
>>>> 
>>>> I have come across what may be an error in genbox.  It concerns reading
>>>> in data.  If instead of 8 elements, 7 are used or if the trailing zeros
>>>> are removed, all is fine.  It almost seems as if genbox has a limitation
>>>> on the length of the lines defining the element edges.
>>>> 
>>>> Cheers,
>>>> Frank
>>>> 
>>>> Reading      32 =       8       4       4 elements for box   1.
>>>> Reading       2 =
>>>> error:           8           9
>>>> error:           0           5
>>>> 
>>>> 
>>>> The following is the input file:
>>>> 
>>>> openboat2D.rea
>>>> -2                      spatial dimension
>>>> 2                     number of fields
>>>> #========================================================
>>>> #
>>>> #    Open boat with the Marangoni effect at the top surface.
>>>> #
>>>> #========================================================
>>>> #
>>>> Box
>>>> 8 4
>>>> -4.000000 -3.946236 -3.806717 -3.497438 -3.000000 -2.500000 -2.000000
>>>> -1.500000 -1.000000
>>>> -5.000000 -3.347128  0.000000  3.754690  5.000000
>>>> W  ,W  ,W  ,shl,    V bc's  ! NB:  3 characters each !
>>>> t  ,t  ,   ,   ,    T bc's  !      You must have 2 spaces!!
>>>> 
>>>> 
>>>> 
>>>> 
>>> -- 
>>> Frank Herbert Muldoon, Ph.D. Mechanical Engineering
>>> Technische Universität Wien (Technical University of Vienna)
>>> Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid
>>> Mechanics and Heat Transfer)
>>> Resselgasse 3
>>> 1040 Wien
>>> Tel: +4315880132232
>>> Fax: +4315880132299 
>>> Cell:+436765203470
>>> fmuldoo (skype)
>>> http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.html
>>> 
>>> _______________________________________________
>>> 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
> -- 
> Frank Herbert Muldoon, Ph.D. Mechanical Engineering
> Technische Universität Wien (Technical University of Vienna)
> Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid
> Mechanics and Heat Transfer)
> Resselgasse 3
> 1040 Wien
> Tel: +4315880132232
> Fax: +4315880132299 
> Cell:+436765203470
> fmuldoo (skype)
> http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.html
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users