[Nek5000-users] "byte swap" + grid to grid subroutine

Tue Sep 7 12:06:46 CDT 2010

Hi Stefan,

Thanks for your help. I tested the grid to grid interpolation routine with
the suggestions that you gave and it worked!!

I guess I was confused about the last 2 statements from the Wiki page,

# run nek using the new geometry (e.g. new_geom.f0000)
# run nek using the old geometry and add this code snipplet to userchk()

Now I see what these 2 statements mean.

Regards,

Harish.

On Tue, 7 Sep 2010 nek5000-users at lists.mcs.anl.gov wrote:

> Harish:
>
> I think you mixed up the old field file (oldfld) and the fld file containing the new geometry (newgfld).
>
> call g2gi(newfld, oldfld, newgfld) ! grid2grid interpolation
>
> Note: You need to run your old case (oldfld.rea, oldfld.map) to interpolate oldfld onto the newgfld mesh.
> The output of this operation results in a new field file (newfld).
>
> hth,
> Stefan
>
> ----- Original Message -----
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Sent: Monday, September 6, 2010 6:22:26 PM GMT -06:00 US/Canada Central
> Subject: Re: [Nek5000-users] "byte swap" + grid to grid subroutine
>
> Hi Stefan,
>
> Thanks a lot for looking into this issue.
>
> Here is the link with bfs.tar file containing *.rea, *.usr, SIZE, old
> field file (bfs10.f00001) and new geometry file (bfs20.f00001).
>
> www2.uic.edu/~hkanch1/bfs.tar
>
> I actually converted a bfs1.fld20 file to bfs10.f00001. I use *.fld format
> because, my Visit 1.12.2 on Windows does not work with *.f???? format
> files. I can only see the mesh and not the velocity contours.
>
> Thanks,
>
> Harish.
>
>
> On Mon, September 6, 2010 4:47 pm, nek5000-users at lists.mcs.anl.gov wrote:
> > Harish:
> >
> > Can you send me your case? I cannot reproduce your problem here.
> >
> > -Stefan
> >
> > On Mon, Sep 6, 2010 at 8:39 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
> >
> >> Thanks Stefan!!
> >>
> >> Harish
> >>
> >> On Mon, 6 Sep 2010 nek5000-users at lists.mcs.anl.gov wrote:
> >>
> >> > Harish,
> >> >
> >> > I will check into that.
> >> >
> >> > Stefan
> >> >
> >> > On Mon, Sep 6, 2010 at 7:52 AM, <nek5000-users at lists.mcs.anl.gov>
> >> wrote:
> >> >
> >> > > Hello,
> >> > >
> >> > > I am testing subroutine g2gi() for a case of backward facing step
> >> (BFS)
> >> > > field file. After starting Nek5k in post-processing mode, the code
> >> stops
> >> > > in the subroutine g2gi() while reading the old field file.
> >> > >
> >> > > Nek5k is compiled with GNU compiler and OpenMPI. Here is the code
> >> snippet
> >> > > from the userchk() subroutine.
> >> > > ***********************************************************************
> >> > >      subroutine userchk
> >> > >      include 'SIZE'
> >> > >      include 'TOTAL'
> >> > >
> >> > >      character*132 newfld, oldfld, newgfld
> >> > >      data newfld, oldfld, newgfld /'new0.f00001','bfs10.f00001',
> >> > >     $                              'bfs20.f00001'/
> >> > >
> >> > >      call g2gi(newfld,oldfld,newgfld) ! grid2grid interpolation
> >> > >      call exitt()
> >> > > ***********************************************************************
> >> > >
> >> > > Here is the part of the logfile which shows the point at which the
> >> > > simulation hangs and I have to kill the job after this point.
> >> > >
> >> > > ***********************************************************************
> >> > >  nsteps=0 -> skip time loop
> >> > >  running solver in post processing mode
> >> > >
> >> > >  call userchk
> >> > >  grid-to-grid interpolation
> >> > >  set initial conditions
> >> > >  Checking restart options: bfs10.f00001
> >> > >
> >> > >  Reading checkpoint data
> >> > >       1 OPEN: bfs10.f00001
> >> > >  byte swap: F   6.5432100     -2.93127722E+35
> >> > > ***********************************************************************
> >> > >
> >> > > I am attaching the logfile for you reference. What is the solution
> >> to
> >> this
> >> > > problem?
> >> > >
> >> > > Thanks for all the help,
> >> > >
> >> > > Harish.
> >> > > _______________________________________________
> >> > > 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
> >
>
>
> _______________________________________________
> 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
>