fortran log的表示问题 运行错误 意思是说 没有alog这个数组 fortran里不是自定义了这个函数?

fortran log的表示问题 运行错误 意思是说 没有alog这个数组 fortran里不是自定义了这个函数?
real::tw(0:n) 湿球温度估算值序列
real::eti
real::etw 饱和水汽压
real,parameter::a=8.15*10**(-4)
real::ej(0:n) 水汽压的计算值
real::err(0:n) 计算误差
real::ermin 计算误差的最小值
real::twmin 合理的湿球温度
do while(.true.)
输入变量
write(*,*)"相对湿度为"
read(*,*)u
write(*,*)"干球温度为"
ead(*,*)t
write(*,*)"本站气压为"
read(*,*)p
计算水汽压 t0=273.15+t
if(t>0)then
alog10(eti)=10.79574*(1-t1/t0)-5.02800*alog10(t0/t1)+1.50475*10**(-4)*(1-10**((-8.2969*(t0/t1-1))))+0.42873*10**(-3)*(10**(4.76955*(1-t1/t0))-1)+0.78614
else
alog10(eti)=-9.09685)*(t1/t0-1)-3.56654*alog10(t1/t0)+0.87682*(1-t0/t1)+0.78614
end if
e=eti*u
算出误差序列
tw(0)=t-15
do i=0,n
if(tw(i)>0)then
alog10(etw)=10.79574*(1-t1/t0)-5.02800*alog10(t0/t1)+1.50475*10**(-4)*(1-10**((-8.2969*(t0/t1-1))))+0.42873*10**(-3)*(10**(4.76955*(1-t1/t0))-1)+0.78614
else
alog10(etw)=-9.09685)*(t1/t0-1)-3.56654*alog10(t1/t0)+0.87682*(1-t0/t1)+0.78614
end if
ej(i)=etw-a*p*(t-tw(i)) 算出水汽压的计算值
err(i)=abs(e-ej(i))
if(i==n) exit
tw(i+1)=tw(i)+0.01
end do
算出最合适的湿球温度
errmin=err(0)
do i=1,n
if(err(i)