fortranlog的表示问题运行错误意思是说没有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合理的湿球温度

dowhile(.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

endif

e=eti*u

算出误差序列

tw(0)=t-15

doi=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

endif

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

enddo

算出最合适的湿球温度

errmin=err(0)

doi=1,n

if(err(i)