|
[Sponsors] |
January 11, 2006, 15:27 |
cell distiction
|
#1 |
Guest
Posts: n/a
|
Hi, im working with open chanel flow, my problem is,
i created a variable called ALT, is defined by multiplying rho*1000 in non porous zones and rho*900 in porous media, how can i make that distiction? ive tried with this ground code IF((IX.LT.15).AND.(IX.GT.14))THEN CALL FN0(LBNAME('ALT'), LBNAME('DEN1')) CALL FN25(LBNAME('ALT'), 0.001) ELSE CALL FN0(LBNAME('ALT'), LBNAME('DEN1')) CALL FN25(LBNAME('ALT'), 0.0009) ENDIF The problem is that the condition if is not working, all the domain goes directly to else, why? ix is the cell number in X direction right? help please |
|
January 12, 2006, 02:52 |
Re: cell distiction
|
#2 |
Guest
Posts: n/a
|
Simply because the condition 14<ix<15 is never true for integer ix...
|
|
January 12, 2006, 02:54 |
Re: cell distiction
|
#3 |
Guest
Posts: n/a
|
Simply because the condition 14 < ix < 15 is never true for integer ix...
|
|
January 12, 2006, 12:20 |
Re: cell distiction
|
#4 |
Guest
Posts: n/a
|
Rami, thatks for your answer, buti misstyped, instead of 14 is 34, what i want to know is that if is correct to use ix for that calculation, or what integer should i use
IF((IX.LT.15).AND.(IX.GT.34))THEN CALL FN0(LBNAME('ALT'), LBNAME('DEN1')) CALL FN25(LBNAME('ALT'), 0.001) ELSE CALL FN0(LBNAME('ALT'), LBNAME('DEN1')) CALL FN25(LBNAME('ALT'), 0.0009) ENDIF |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
How to identify boundary cell and internal cell | dbxmcf | OpenFOAM Running, Solving & CFD | 2 | November 6, 2016 09:11 |
cell velocity U and its cell face area | Chinayu | Siemens | 0 | May 6, 2006 22:44 |
unsturcture cell & cylinder cell Difference~ | lenson | Siemens | 1 | January 9, 2006 05:32 |
Cell vertex FV vs Cell centered FV | peter.zhao | Main CFD Forum | 1 | September 15, 2004 09:26 |
Cell-vertex and cell-centered finite volume method | Praveen | Main CFD Forum | 4 | September 8, 2003 23:58 |