[Nek5000-users] Solution for Re: genbox error

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Mon May 3 13:36:59 CDT 2010


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




More information about the Nek5000-users mailing list