I am trying to apply different heat fluxes at different pre-defined surfaces but the results does not show any change in temperature. I copied below the subroutine. Can someone please check what is wrong with my subroutine?

SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,

1 TEMP,PRESS,SNAME)

C

INCLUDE 'ABA_PARAM.INC'

C

DIMENSION COORDS(3),FLUX(2),TIME(2)

CHARACTER*80 SNAME

IF (SNAME .EQ. 'ExposedA') THEN

FLUX(1)=-250.*SOL

FLUX(2)=-250.

ELSE IF (SNAME .EQ. 'ExposedB') THEN

FLUX(1)=350.*SOL

FLUX(2)=350.

END IF

RETURN

END