[gezibash]

Hi,

I am writing my own solver for a problem and I've been able to implement most features I need, however one is proving to be a bit difficult.

The idea is that I have an h field which I can most easily describe as a function

h = f(x,y)

This field is then used in the solver to generate a pressure distribution. Here is the solve block as is atm

        while (simple.correctNonOrthogonal())
        {

          if(runTime.value() < 3e-5){
            h = planarOffset;
          }
          else {
            h = 6.667e-3 * mesh.C().component(vector::Y);
          } 

          solve
          (
            (1/(12*Mu))*(fvm::laplacian(pow(h,3)*p,p)) == fvm::ddt(h,p)
          );
        }

This will compile and it works. My issue is can I somehow pass this portion

          if(runTime.value() < 3e-5){
            h = planarOffset;
          }
          else {
            h = 6.667e-3 * mesh.C().component(vector::Y);
          }

from the problem statement through some sort of boundary condition or something. I can do it like this, but I would like to test out a few different h fields, and I have to wmake all the time. I would prefer if the solver was a bit more general.