Hi,

I want to use my own groud to write results according a specific format in a file that I created in ground. The problem is that I don't know the names used for variables and if for example U1, W1, Temp1 are used as vectors or matrix, if yes what are names of arguments,are they IX, IY and IZ.

Thanks.

Mohamed Ouzzane

All the variables are stored as 1-dimension array, F array. You should indicate the location of the variable inside somewhere 3-D grid in the F array. You should use the L0f function to find the starting location of the variable for this Z-slab inside F array. (PHOENICS use slab-by-slab solver)

L0FU1=L0F(U1)

L0FTEM1=L0F(INAME('TEM1 ')) (if you define the variable name in group 7) And then

DO II=1,NX

JI=(II-1)*NY

DO JJ=1,NY

LFU1=L0FU1+JJ+JI

LFTEM1=L0FTEM1+JJ+JI

F(LFU1)=function[F(LFTEM1)] (this means some function expression.)

ENDDO

ENDDO

Richard

Remember that PHOENICS uses slab-by-slab solver as default, not always.

You should add that in PHOENICS you can find "slab wise variables", "whole field variables" ... you have to be carefull to obtain their L0F

I will recommned to Mohamed to read POLIS issue about F array. It is available from CHAM web page. What do you think?

Regards

Kike

Dear Mohamed

You can find a lot of information about the ways by which PHOENICS store data in memory. Go to CHAM web page and enter into POLIS service.

There you will find an index of issues. Search for the "F array" and read the page content. As Richard say the F array is a vector in which PHOENICS store all the data sequentially, ie. a1, a2..., aN, b1, b2..., bN, ...

So you should know where the data for U1 beguin into the array and it can be achieve by certain standard ways in PHOENICS

Ask here for help if you have any doubt after reading.

Regards

Kike

Thanks for all,

I used with succes in a simple example the suggested below procedure to write results in a new file according to a specific format.

C * ------------------- SECTION 6 ---- Finish of iz slab.

IF(ISWEEP.EQ.LSWEEP) THEN

L0T=L0F(LBNAME('TEM1'))

DO 211 I=1, NX*NY C KES=L0T+I 211 WRITE(9,*)I,F(L0T+I)

Now I like to know where I can find the variable names of x, y and z cordinates for cells (like LBNAME('TEMP1') above .

In Additionnal to the precedant procedure , I want to write temperature with space coordinates.

Thank you

Mohamed Ouzzane

Dear Mohamed

You can see under include files (ie. grdloc and satgrd) combining your search with POLIS infos. I insist, read POLIS, at least once, before to come here. You will learn more.

In this case (ie. for temperature) you can use XG2D, YG2D and ZGNZ. You should note that you need to insert lines as: call make(xg2d) in Group 1 into GROUND sub.

Regards

Kike