How to get vertex(node) number in CFX user Fortran..

flyply · February 1, 2014, 05:30

Hello.
I try to input the external data into the CFX.
External data was generated by 3rd party software.

Data is exactly same with mesh data, but only contain the X, Y, Z coordinate and one Scalar data. so external data array size is (total node number x 4).

anyway, At present, I match the data like this(very simple &stupid

)
=====================================
PRINT *, 'NLOC = ', NLOC
DO ILOC = 1, NLOC
FR(ILOC) = 0.0
DO LN = 1, NLIN
IF (ABS(X(ILOC) - TABLE(1,LN)) .LE. 0.000001 .AND.
& ABS(Y(ILOC) - TABLE(2,LN)) .LE. 0.000001 .AND.
& ABS(Z(ILOC) - TABLE(3,LN)) .LE. 0.000001) THEN
FR(ILOC) = TABLE(4, LN)
END IF
END DO
END DO
=====================================

This method working well. but very slow, especially in Large Problem.
In my case, using 9 million node, It took 5 hours.

So, I try to match the node(vertex) number.
but there are no information, how to access the vertex number, in User Guide.
but, i find some clue in CFX Modeling Guide : 18.5.5.1 Message

CALL MESAGE( 'WRITE', 'Max Temp is '//CFROMR(TMAX)//’ at vertex’// & CFROMI(IVX)

'IVX' seems to Vertex number, but when i try to like this
=======================================
CALL PEEKI ('/USER/'//CDIR_GLOB//'/IVX', 1, IVX,
& 'STOP', CRESLT, IZ)
=======================================
It show the error message like bellow.
=======================================
Details of erros:--
Error detected by routine PEEKI
CDNAM = /USER/MESH_INFO/IVX
CRESLT = NONE
Current Directory : /FLOW/SOLUTION/TSTEP10/CLOOP1/ZN1/VERTICES
.
.
=======================================
Anyway, Is there anyone who knows about how to get Vertex Number in CFX user Fortran??

ghorrocks · February 2, 2014, 16:18

What are you trying to do with the external data? Is it during the solver or post processing?

flyply · February 4, 2014, 10:44

Hello, ghorrocks
It is for Chemical Reaction and External Data is UV intensity.
for the present, I improved the algorithm and it is worth, but not sufficient.
I take fluent into account.

ghorrocks · February 4, 2014, 17:16

CFX already has a 3D interpolation function available. But if you have millions of nodes it will probably take a while to run as well. But it might be faster than your implementation.

flyply · February 4, 2014, 19:51

Yes. I already know, but 3D interpolation is very slow in my case. so I trying to do with fortran. but i will one more check the user guide about 3D interpolation.
thanks.

February 1, 2014, 05:30	How to get vertex(node) number in CFX user Fortran..	#1
flyply New Member Join Date: Feb 2012 Posts: 6 Rep Power: 14	Hello. I try to input the external data into the CFX. External data was generated by 3rd party software. Data is exactly same with mesh data, but only contain the X, Y, Z coordinate and one Scalar data. so external data array size is (total node number x 4). anyway, At present, I match the data like this(very simple &stupid ) ===================================== PRINT *, 'NLOC = ', NLOC DO ILOC = 1, NLOC FR(ILOC) = 0.0 DO LN = 1, NLIN IF (ABS(X(ILOC) - TABLE(1,LN)) .LE. 0.000001 .AND. & ABS(Y(ILOC) - TABLE(2,LN)) .LE. 0.000001 .AND. & ABS(Z(ILOC) - TABLE(3,LN)) .LE. 0.000001) THEN FR(ILOC) = TABLE(4, LN) END IF END DO END DO ===================================== This method working well. but very slow, especially in Large Problem. In my case, using 9 million node, It took 5 hours. So, I try to match the node(vertex) number. but there are no information, how to access the vertex number, in User Guide. but, i find some clue in CFX Modeling Guide : 18.5.5.1 Message CALL MESAGE( 'WRITE', 'Max Temp is '//CFROMR(TMAX)//’ at vertex’// & CFROMI(IVX) 'IVX' seems to Vertex number, but when i try to like this ======================================= CALL PEEKI ('/USER/'//CDIR_GLOB//'/IVX', 1, IVX, & 'STOP', CRESLT, IZ) ======================================= It show the error message like bellow. ======================================= Details of erros:-- Error detected by routine PEEKI CDNAM = /USER/MESH_INFO/IVX CRESLT = NONE Current Directory : /FLOW/SOLUTION/TSTEP10/CLOOP1/ZN1/VERTICES . . ======================================= Anyway, Is there anyone who knows about how to get Vertex Number in CFX user Fortran??

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
CFX User Fortran: Particle User Sources	hustxinxin	CFX	0	March 8, 2012 08:31
CFX Solver Reynolds Number	haider760	CFX	2	March 4, 2012 22:05
Check particle impaction with User Fortran	Julian K.	CFX	3	January 12, 2012 09:46
Unaligned accesses on IA64	andre	OpenFOAM	5	June 23, 2008 10:37
how to use USER FORTRAN with CFX	cfd_99	Main CFD Forum	1	June 4, 1999 05:42

February 2, 2014, 16:18		#2
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,703 Rep Power: 143	What are you trying to do with the external data? Is it during the solver or post processing?

February 4, 2014, 10:44		#3
flyply New Member Join Date: Feb 2012 Posts: 6 Rep Power: 14	Hello, ghorrocks It is for Chemical Reaction and External Data is UV intensity. for the present, I improved the algorithm and it is worth, but not sufficient. I take fluent into account.

February 4, 2014, 17:16		#4
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,703 Rep Power: 143	CFX already has a 3D interpolation function available. But if you have millions of nodes it will probably take a while to run as well. But it might be faster than your implementation.

February 4, 2014, 19:51		#5
flyply New Member Join Date: Feb 2012 Posts: 6 Rep Power: 14	Yes. I already know, but 3D interpolation is very slow in my case. so I trying to do with fortran. but i will one more check the user guide about 3D interpolation. thanks.