[DSMC123]

Hi,

I have a working vectorCodedSource:

momentumSource
{
    type            vectorCodedSource;
    active          on;
	selectionMode   cellSet;
	cellSet        ionicWindRegion;

    vectorCodedSourceCoeffs
    {
		fieldNames (U);
		redirectType    velocitySource;

		codeCorrect
               #{
                      Pout<< "**codeCorrect**" << endl;
                #};

               codeAddSup
               #{
                      Pout << "***codeAddSup***" << endl;
                #};

               codeSetValue
              #{
                        Pout<< "**codeSetValue**" << endl;
			cellSet selectedCells(mesh_, cellSetName_);
			labelList cells = selectedCells.toc();
			vectorField values(cells.size());

			forAll(cells, i)
			{
				scalar ccX = mesh_.cellCentres()[cells[i]].x()*1000.0;
				scalar ccY = mesh_.cellCentres()[cells[i]].y()*1000.0;

				scalar Umax = someEquation;

				scalar Unorm = someOtherEquation;

				scalar Ux = Umax*Unorm;

				values[i] = Ux;
				values[i].y() = 0.0;
				values[i].z() = 0.0;
			}

			eqn.setValues(cells, values);
        #};

        // Dummy entry. Make dependent on above to trigger recompilation
        code
        #{
            $codeInclude
            $codeCorrect
            $codeAddSup
            $codeSetValue
        #};
    }
    
//     sourceTimeCoeffs
    velocitySourceCoeffs
    {
    // Dummy entry
    }
  
}

This works as expected, but ideally I would only constrain the x-component of velocity, I get a nonphysical looking velocity field because it keeps setting the y-component to zero inside the cellSet. Is it possible to only constrain one component of the velocity? If I don't set all 3 components, I of course get runtime errors.

Thanks!