[Nek5000-users] Changed mesh coordinates

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

Dear Neks,


I'm doing some data processing for my simulations in nek5000. I found for some elements, however, the coordinates are not exactly as specified. For example, I have an element whose vertices are as follows.


     ELEMENT         195 [    1m]    GROUP    0
  0.000000E+00 -0.130526E+00 -0.129129E+00  0.000000E+00
  0.100000E+01  0.991445E+00  0.979179E+00  0.987876E+00
  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
  0.000000E+00 -0.130526E+00 -0.129129E+00  0.000000E+00
  0.100000E+01  0.991445E+00  0.979179E+00  0.987876E+00
  0.263158E+00  0.263158E+00  0.263158E+00  0.263158E+00


When I used the following code to print the coordinates for this element, I got the values below.


      do e=1,nelv
        eg = lglel(e)
        if(eg .eq. 195) then
          write(*,*) xm1(1,1,1,e),ym1(1,1,1,e),zm1(1,1,1,e)
          write(*,*) cbc(1:6,e,1)
        end if
      end do


  7.500381456730263E-017   1.00000000000000       0.000000000000000E+000
  0.000000000000000E+000
 W  E  E  E  P  E


Is there any reason why the x coordinate for the first element vertex is not exactly 0.0e0? Hope someone can help. Thank you very much.


Best regards,


Tony
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20150201/6ec991ab/attachment.html>