关于子程序问题，求高手帮忙，谢谢

marc编译子程序时出错，求高手帮忙，谢谢
Log文件：
E:\luosen\circle\hcontrol.f(3) : Error: The same letter must not appear as a single letter or be included in a range of letters more than once in all IMPLICIT statements in a scoping unit.
implicit real*8 (a-h,o-z)
--------^
compile failed for E:\luosen\circle\hcontrol.f
MSC.Marc Exit number 3
子程序为：
subroutine flux(f,temflu,mibody,time)
include '../common/implicit'
implicit real*8 (a-h,o-z)
dimension f(2),mibody(6),temflu(7)
c* * * * * *
c
c user subroutine for non-uniform flux input.
c
c f(1) flux value (to be defined)
c f(2) derivative of flux with respect to temperature
c (to be defined; optional, this might improve the
c convergence behavior)
c
c temflu(1) estimated temperature
c temflu(2) previous volumetric flux
c temflu(3) temperature at beginning of increment
c temflu(4,5,6)integration point coordinates
c mibody(1) element number
c mibody(2) flux type
c mibody(3) integration point number
c mibody(4) flux index
c mibody(5) not used
c mibody(6) =1 : heat transfer
c =2 : joule
c =3 : bearing
c =4 : electrostatic
c =5 : magnetostatic
c =6 : acoustic
c mibody(8) layer number for heat transfer shells elements
c and volume flux
c time time
c
c* * * * * *
pai=3.1415926
f(1)=0
f0=6.5e4
y0=temflu(5)
z0=temflu(6)
ss=sqrt(y0*y0+z0*z0)
th0=acos(y0/ss)
deltat=2*pai/36
th=deltat*time
if(z0.le.0) th0=2*pai-th0
deltath=(th-th0)*180/pai
if(deltath.ge.-5.0.and.deltath.le.5.0) then
if(time.le.0.1) then
f(1)=1.5*f0*time
else
f(1)=1.2*f0
endif
else
f(1)=0
endif
return
end