[himanshu28]

Hi,

I am trying to achieve following using PyFoam. I know it can be done in a very trivial way, but not sure how to go about it.

I have a data(scalarField(Numpy array)), I want to write it as a initial condition of a Passive scalar at a given time-step. Here is how I implemented it using PyFoam.

fname= loc+"T";
f = ParsedParameterFile(fname)
    for j in range(0,len(data)):
        f["internalField"][j] = data[j];
    f.writeFile()

But the bottleneck of this implementation is that the T should have a nonuniform<list> in it so that each entry in that list can be replaced by my data.
Following is I want, I want to use the a general T file of following format

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.4.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0; 
    format      ascii;
    class       volScalarField;
    location    "0"; 
    object      T;   
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //  

dimensions      [0 0 0 1 0 0 0]; 

internalField   uniform 0;


boundaryField
{
    inlet
    {    
        type            zeroGradient;
    }    
    outlet
    {    
        type            zeroGradient;
    }    
    wall.ceiling
    {    
        type            zeroGradient;
    }    
 wall.1
    {
        type            zeroGradient;
    }
    wall.3
    {
        type            zeroGradient;
    }
    emptyZ
    {
        type            empty;
    }
}

and replace the uniform to the non-uniform list of my data. I don't know how to go about it. I looked into various *.py files like PyFoam.RunDictionary.SolutionFile import SolutionFile
and some other too but couldn't accomplish, I am not able to get any detail documentations of using these functions.Any guidance will be really helpful.

Edit:
In simple words I would like to write the data in the openFoam format using the case file as a template.

@gschaider any suggestions...


Regards