User sub routine for Stokes number
 

Hi,

I'm am modeling an airblast atomizer and would like to look at the Stokes number of my particles amongst other things in the subsequent solved flow. According to the CFX help file I have to wrtie a user sub routine in Fortran to do so.

Can anyone point me as to where I could find similar examples or a nice user guide on how to implement this.

Thanks in advacne for any time spent on this.

Cheers

Mat

Ok after doing some research over the web on this issue I found that I have to use the USER_GETVAR function in Fortran.

Can anyone point me to somewhere where an example of a subroutine with this could be found or help me with the issue?

Cheers

Edit-- this is what I have so far:

#include "cfx5ext.h"

dllexport(stokes)

SUBROUTINE STOKES (NLOC, NRET, NARG, RET, ARGS, CRESLT,

& CZ,DZ,IZ,LZ,RZ)

IMPLICIT NONE

C Calculate magnitude of Stokes number

#include "stack_point.h"

C Result flag

CHARACTER CRESLT*(*)

INTEGER NLOC

INTEGER NARG

INTEGER NRET

REAL ARGS(NLOC,NARG)

REAL RET(NLOC,NRET)

INTEGER IZ(*)

CHARACTER CZ(*)*(1)

DOUBLE PRECISION DZ(*)

LOGICAL LZ(*)

REAL RZ(*)

__stack_point__ pST

C Obtain Stokes number

CALL USER_GETVAR('ptstt',CRESLT,
& pST,CZ,DZ,IZ,LZ,RZ)

IF (CRESLT.NE. 'GOOD') STOP


END