[Nek5000-users] r,s,t index for given point
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Mon Jul 12 00:46:32 CDT 2010
Hi,
that was not it, I'm afraid (though I thought gllel transforms from
local to global).
But I also realized that the hpts routine does not deliver the correct
values for me. To test this, I set up a simple 3d box with inflow,
outflow and symmetry all around, Re very low, initial condition=boundary
condition and 4 points.
Please find all files attached, the output for the first point is OK,
but then the values seem to shift positions and don't make sense.
Also, when running it for 1 step or more, I get the error message:
"
dump history points
done :: dump history points
hisfile:
/Users/m0s1978/nektonruns/rectest/rectest.his
ERROR: .sch file already exists. ierr= 1
"
And lastly, the newest repo won't compile, here is the error message:
"
Error: /Users/m0s1978/nek5trial/trunk/nek/drive2.f, line 1551: This name
has already been assigned a data type. [NVTOT]
integer*8 ntot,nvtot,nptot
---------------------^
"
Thanks for any help,
Markus
nek5000-users at lists.mcs.anl.gov wrote:
>
> Hi Markus,
>
> It's possible that it's returning the local (to the processor)
> element number. You would get the global (to the computational
> domain, defined by the .rea file) element number via:
>
> integer e,eg
>
> eg = lglel(e)
>
>
> Paul
>
>
> On Sun, 11 Jul 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
>> Hi again,
>>
>> Regarding finding element numbers (such as what gllel would return)
>> for arbitrary history points, I was experimenting with the history
>> point routine hpts a little and it seems that this function call:
>>
>> call findpts(inth_hpts,rcode,1,
>> & proc,1,
>> & elid,1,
>> & rst,ndim,
>> & dist,1,
>> & pts(1,1),ndim,
>> & pts(2,1),ndim,
>> & pts(3,1),ndim,npoints)
>>
>> does not return the same element number (elid) as it can be found in
>> the .rea file, for instance.
>> For a handful of manually set points, I checked that the returned
>> element number does not provide the right coordinates to enclose the
>> desired point.
>>
>> Is there a way to transform the function return or am I missing
>> something else?
>>
>> Thanks,
>> Markus
>>
>> _______________________________________________
>> 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
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rectest.tar.gz
Type: application/x-gzip
Size: 24067 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20100712/2f47e9e7/attachment.bin>
More information about the Nek5000-users
mailing list