如何获得网格坐标？

[danao]

定义了一个矢量，类型为dimensionedVector，=网格坐标*vector(1,0,0)。发现mesh.c()是一个scalerField，所以编译就报错“error: conversion from ‘Foam::tmp<Foam::GeometricField....". mesh.c()这里应该如何修改？

[史浩]

fvMesh::C()返回的是一个矢量场，而不是一个矢量（dimensionedVector）

 const Foam::volVectorField& Foam::fvMesh::C() const
 {
     if (!CPtr_)
     {
         makeC();
     }
 
     return *CPtr_;
 }

[danao]

@史浩 那怎么才能得到网格xyz坐标的标量值？

[史浩]

@danao 这个矢量场中，每一个矢量的三个分量就是网格中心的x、y、z值

[danao]

@史浩 mesh.C().component(0)，是这个吗？我试了好像还是volVectorField 类型，而不是一个scalar

[史浩]

@danao 你说的方法也可以，不过我没有用过，你试一下可以写一下这个函数的用法，我也学习一下。
我个人比较喜欢用是[]运算符，类似于这样

forAll(mesh.C(), celli)
{
    Scalar x = mesh.C()[celli].x();
}

[东岳]

@danao 在 如何获得网格坐标？ 中说：

而不是一个scalar

网格坐标不就是个vector么？

[Dingcy]

@danao 不知道你的最终目的是什么，如果就是为了获得网格坐标，paraview里就可以吧

[danao]

@Dingcy 我要编写源项，公式里需要网格坐标xyz的值

[danao]

@史浩 我试了，提示[celli]没声明，后来我用的mesh.C().component(0)，这个是个volScalarField 类型，有单位m

[danao]

@东岳 是的，我要用到坐标三个方向的具体数值，标量