Re: [Nek5000-users] "byte swap" + grid to grid subroutine
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: [email protected] To: [email protected] 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, [email protected] wrote:
Harish:
Can you send me your case? I cannot reproduce your problem here.
-Stefan
On Mon, Sep 6, 2010 at 8:39 AM, <[email protected]> wrote:
Thanks Stefan!!
Harish
On Mon, 6 Sep 2010 [email protected] wrote:
Harish,
I will check into that.
Stefan
On Mon, Sep 6, 2010 at 7:52 AM, <[email protected]> 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 [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
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 [email protected] 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: [email protected] To: [email protected] 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, [email protected] wrote:
Harish:
Can you send me your case? I cannot reproduce your problem here.
-Stefan
On Mon, Sep 6, 2010 at 8:39 AM, <[email protected]> wrote:
Thanks Stefan!!
Harish
On Mon, 6 Sep 2010 [email protected] wrote:
Harish,
I will check into that.
Stefan
On Mon, Sep 6, 2010 at 7:52 AM, <[email protected]> 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 [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users _______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
participants (1)
-
nek5000-users＠lists.mcs.anl.gov