[Nek5000-users] Solution for Re: genbox error

Mon May 3 14:59:33 CDT 2010

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