Settler
Structural
- May 22, 2010
- 88
Hello all,
I have this subroutine and I'm getting problems during compilation. Here's 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
C
PARAMETER(PI_OVER_180=3.141592653589793/180)
C SUBMERGE_SPEED: in/sec
SUBMERGE_SPEED = 0.525
C THERMAL ENVIRONMENT
BATH_TEMPERATURE = 828.0
AIR_TEMPERATURE = 65.0
C THERMAL PROPERTY
BATH_CONVECTION_H = 0.000509578
AIR_CONVECTION_H = 0.0000339719
C GEOMETRY
GX0 = 0.0
GY0 = -23.5
GZ0 = 0.0
GX1 = 0.0
GY1 = 16.6875
GZ1 = 21.0478
C PLANE
THETA = 86.0*PI_OVER_180
A = 0.0
B = SIN(THETA)
C = COS(THETA)
GD0 = A*GX0+B*GY0+C*GZ0
GD1 = A*GX1+B*GY1+C*GZ1
D = GD0 + SUBMERGE_SPEED*TIME(2)
C FLUX
F = A*COORDS(1)+B*COORDS(2)+C*COORDS(3)-D
IF (F .LT. 0.0) THEN
FLUX(1) = BATH_CONVECTION_H*(BATH_TEMPERATURE-SOL)
FLUX(2) = -1.0*BATH_CONVECTION_H
ELSE
FLUX(1) = AIR_CONVECTION_H*(AIR_TEMPERATURE-SOL)
FLUX(2) = -1.0*AIR_CONVECTION_H
ENDIF
RETURN
END
Abaqus is linked to a compiler as you can see in the following link:
I have this subroutine and I'm getting problems during compilation. Here's 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
C
PARAMETER(PI_OVER_180=3.141592653589793/180)
C SUBMERGE_SPEED: in/sec
SUBMERGE_SPEED = 0.525
C THERMAL ENVIRONMENT
BATH_TEMPERATURE = 828.0
AIR_TEMPERATURE = 65.0
C THERMAL PROPERTY
BATH_CONVECTION_H = 0.000509578
AIR_CONVECTION_H = 0.0000339719
C GEOMETRY
GX0 = 0.0
GY0 = -23.5
GZ0 = 0.0
GX1 = 0.0
GY1 = 16.6875
GZ1 = 21.0478
C PLANE
THETA = 86.0*PI_OVER_180
A = 0.0
B = SIN(THETA)
C = COS(THETA)
GD0 = A*GX0+B*GY0+C*GZ0
GD1 = A*GX1+B*GY1+C*GZ1
D = GD0 + SUBMERGE_SPEED*TIME(2)
C FLUX
F = A*COORDS(1)+B*COORDS(2)+C*COORDS(3)-D
IF (F .LT. 0.0) THEN
FLUX(1) = BATH_CONVECTION_H*(BATH_TEMPERATURE-SOL)
FLUX(2) = -1.0*BATH_CONVECTION_H
ELSE
FLUX(1) = AIR_CONVECTION_H*(AIR_TEMPERATURE-SOL)
FLUX(2) = -1.0*AIR_CONVECTION_H
ENDIF
RETURN
END
Abaqus is linked to a compiler as you can see in the following link: