[Nek5000-users] Post processing mode

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Wed Mar 2 06:46:52 CST 2016


NSTEPS=FINTIME=0 in .rea is enough to trigger post-processing mode -- i.e. call of userchk and exit after initialization.

Aleks
________________________________
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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20160302/f6e9a280/attachment.html>


More information about the Nek5000-users mailing list