STAR source UDF syntax

Emmanuel Resch

Hi,

Another question about syntax for UDFs...

In SORSCA.F, I'm having a hard time getting the scalar sources to perform the way I expect. I'd like to have cells ICTID 4 & 5 consume/produce scalars according to the following (the main part of my SORSCA.F subroutine is below)...

IF(IS.EQ.1.AND.ICTID.EQ.4)THEN !scalar 1 is hydrogen

S1P=(-1)*Janode/Z/F*MH2*GREAT

S2P=GREAT

END IF

IF(IS.EQ.2.AND.ICTID.EQ.4)THEN !scalar 2 is water

S1P=Janode/Z/F*MH2O*GREAT

S2P=GREAT

ENDIF

IF(IS.EQ.3.AND.ICTID.EQ.5)THEN !scalar 3 is Nitrogen (inert)

S1P=Jcathode/Z/F*MN2*.7/0.3*GREAT !For a given amount of Oxygen in a space, there will be 0.7/0.3 times as much N2

S2P=GREAT

ENDIF

IF(IS.EQ.4.AND.ICTID.EQ.5)THEN !scalar 4 is oxygen

S1P=(-Jcathode)/Z/F*MO2*GREAT

S2P=GREAT

ENDIF

Can anyone notice any obvious errors in the above coding?

Thanks in advance,

-Emmanuel