[Nek5000-users] lambda2

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

Hi JC,

L2 does not look to be properly declared?  Perhaps

      real l2(lx1,ly1,lz1,lelt)

fixes this?  (Let me know if otherwise...)

Paul

On Fri, 28 Jan 2011, nek5000-users at lists.mcs.anl.gov wrote:

> Hi Nek's,
>
> I have a few troubles using the lambda-2 subroutine. Here is my usrchk:
>
> subroutine userchk
> include 'SIZE'
> include 'TOTAL'
> include 'RESTART'
>
> character*80 filename(9999)
> parameter(lt = lx1*ly1*lz1*lelv)
> real L2(lx1,ly1,lz1,1)
> integer nfiles
>
> ntot = nx1*ny1*nz1*nelv
>
> ifreguo = .true.
>
> 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))
> call lambda2(L2)
> ifto = .true.
> call outpost(vx,vy,vz,pr,L2,'his')
> enddo
>
> call exitt
>
> return
> end
>
> The thing is that it actually outposts nothing (or at least it looks like to
> me that it doesn't outpost anything). I have tried switching from my l2
> variable to the temperature t variable but it did not changed anything. Any
> clues?
>
> Regards,
> JC
>