[Nek5000-users] Post processing mode
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Tue Mar 1 22:11:20 CST 2016
We should indeed add this more explicitly in our documentation.
The main aspect is to have an output file which is written at which iostep value is written in the .rea file, see parameter 15
http://nek5000.github.io/NekDoc/Nek_usersch7.html#x9-690007.1
This file can be later visualized in Visit, for which you will need to run the visnek script. (visnek is in the tools/scripts folder)
Oana
________________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Tuesday, March 01, 2016 4:50 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Post processing mode
Hi neks,
can someone tell me step by step how can i set Nek5000 in post-processing mode? I read something about it into the mailing list archive, but nothing that explain step by step what i have to do.
I have done the following step:
-create a file, named file.list where in the first row i have to write the total number analyzed. In the following rows i have to write the file's output name, since the first to the last,
- now, i read i have to set in .rea file the parameter related to Nstep equal 0 (param p011), i have to fix P010,P14 and P15 equal 0 too?
I find out this post in your mailing list archive and i added it in my usercheck subroutine, but everytime i have errors in compilation!
SUBROUTINE USERCHK
INCLUDE 'SIZE'
INCLUDE 'TOTAL'
INCLUDE 'RESTART'
character*80 filename(9999)
ntot = nx1*ny1*nz1*nelv
ifreguo = .true. ! dump on regular (uniform) grid instead of GLL
nrg = 16 ! dimension of regular grid (nrg**ndim)
! read file-list
if (nid.eq.0) then
open(unit=199,file='file.list',form='formatted',status='old')
read(199,*) nfiles
read(199,'(A80)') (filename(i),i=1,nfiles)
close(199)
endif
call bcast(nfiles,isize)
call bcast(filename,nfiles*80)
do i = 1,nfiles
call load_fld(filename(i))
! do something
! note: make sure you save the result into arrays which are
! dumped by prepost() e.g. T(nx1,ny1,nz1,nelt,ldimt)
...
! dump results into file
call prepost(.true.,'his')
enddo
! we're done
call exitt
Ifreguo is a logical switch to have a grid without elemnts?Only knots?
Thanks to everyone.
More information about the Nek5000-users
mailing list