OpenFOAM后处理操作编程问题

[sendy]

电场E对于时间的偏导数，在OpenFOAM中编写这个后处理的操作，该怎么编写呢？
Foam::fvc::ddt(psi.boundaryField()[pI].snGrad() + phiE.boundaryField()[pI].snGrad())
我尝试这样写，但是编译问题是这样的：

请教一下大家，这个编程怎么处理呢，谢谢大家

[Veen]

先输出里边需要做时间导数的物理量看看呢？之前遇到过OF有的函数在里边计算会报错

[sendy]

感谢 我尝试先创建了一个物理量随时间变化的场，然后再对这个场的求梯度值 可以啦