[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:47:05 CDT 2010


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




More information about the Nek5000-users mailing list