[Nek5000-users] Solution for Re: genbox error

Mon May 3 15:33:53 CDT 2010

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