[Nek5000-users] to Michael--> Nek5000 & Gambit, blocked file

[nek5000-users at lists.mcs.anl.gov]

Hi Michael,

I have been running your code and have hit the error below.
Unfortunately, I don't ever use Matlab.

Cheers,
Frank


ELEM & PTS MATRIX RESTURCTURE:    STARTED ****
??? Index exceeds matrix dimensions.

Error in ==> readingNEU2 at 99
        ElmF(f+1:f+n1(i),1)=BC1m(f+1:f+n1(i),1);

Error in ==> rea_write at 17
[Npts,Nelm,Elem,ElemMidPts1,Pts,Elem1,ElmF,...
 
>> 



On Fri, 2010-06-04 at 18:42 -0500, nek5000-users at lists.mcs.anl.gov
wrote:
> Hi Frank,
> 
>  
> 
> This is achieved by the step 2 in the list I gave:  You must set the
> edge nodes from 2 to 3. Alternatively you could set the face from 4 to
> 8 or the volume from 8 to 20 nodes. The extra nodes that gambit throws
> because of this are the midpoints.  For Nekton, the recent midside
> node support feature which fits a parabola through the two corner
> points and midpoint of the element side is where this comes into
> play.  
> 
>  
> 
> So all you need to do is make sure you have set the mesh nodes as
> above, and the matlab code does the rest. It separates the element's
> corner nodes from the midpoint nodes. The midpoint nodes are written
> into the curved side section of the the REA, whereas the corner points
> are written into the Mesh data section. So...The matlab code only uses
> the flag 'm' for curved side data. 
> 
>  
> 
> Also, in the inputs section you can specify a maximum radius for
> consideration. The code will check each edge's corner and midpoints
> and determine the radius curvature. If it is too large it considers it
> a straight edge and moves on.
> 
>  
> 
> Hope this helps,
> 
>  
> 
> Michael
> 
> 
> ----- Original Message -----
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Sent: Friday, June 4, 2010 6:24:06 PM GMT -06:00 US/Canada Central
> Subject: Re: [Nek5000-users] to Michael--> Nek5000 & Gambit, blocked
> file
> 
> Hi Michael,
> 
> Thanks very much for the code.  One question concerning the curved
> sides: if one had a 3D cylindrical pipe grid from Gambit, how would
> the
> curvature of the outside (boundary) faces of the grid be set in your
> code?
> 
> Cheers,
> Frank
> 
> 
> On Thu, 2010-06-03 at 21:13 -0500, nek5000-users at lists.mcs.anl.gov
> wrote:
> > Hi Frank,
> > 
> >  
> > 
> > Try the attached.
> > 
> >  
> > 
> > - Michael
> > 
> > 
> > ----- Original Message -----
> > From: nek5000-users at lists.mcs.anl.gov
> > To: nek5000-users at lists.mcs.anl.gov
> > Sent: Thursday, June 3, 2010 9:01:14 PM GMT -06:00 US/Canada Central
> > Subject: [Nek5000-users] to Michael--> Nek5000 & Gambit, blocked
> file
> > 
> > Hello Michael,
> > 
> > It seems that you sent an email with a *.zip attachment, which due
> to
> > a
> > spam filter, got blocked.  If it is renamed to say *.zipp, then it
> > should get through.
> > 
> > Cheers,
> > Frank
> > 
> > 
> > -- 
> > 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
> -- 
> 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
-- 
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