How to get vertex(node) number in CFX user Fortran..
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)
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 : 188.8.131.52 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??
What are you trying to do with the external data? Is it during the solver or post processing?
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.
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.
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.
|All times are GMT -4. The time now is 20:51.|