fortran log的表示问题 运行错误 意思是说 没有alog这个数组 fortran里不是自定义了这个函数?real::tw(0:n) 湿球温度估算值序列real::etireal::etw 饱和水汽压real,parameter::a=8.15*10**(-4)real::ej(0:n) 水汽压的

fortran log的表示问题 运行错误 意思是说 没有alog这个数组 fortran里不是自定义了这个函数?real::tw(0:n) 湿球温度估算值序列real::etireal::etw 饱和水汽压real,parameter::a=8.15*10**(-4)real::ej(0:n) 水汽压的
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)

fortran log的表示问题 运行错误 意思是说 没有alog这个数组 fortran里不是自定义了这个函数?real::tw(0:n) 湿球温度估算值序列real::etireal::etw 饱和水汽压real,parameter::a=8.15*10**(-4)real::ej(0:n) 水汽压的
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
这句你想做什么?
等号左边是alog10(eti) 是一个函数返回值.
等号右边是一个表达式.
你不能让一个函数返回值等于一个表达式.这从逻辑上就是说不通的.
你是想解方程么?